def smooth(topg, topg_smoothed, usurf, theta): "Smooth the bed topography." grid = topg.grid() smoother = PISM.BedSmoother(grid, 1) smoother.preprocess_bed(topg) smoother.theta(usurf, theta) topg_smoothed.copy_from(smoother.smoothed_bed())