Exemple #1
0
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())