oggm.tasks.filter_inversion_output#

oggm.tasks.filter_inversion_output(gdir, n_smoothing=5, min_ice_thick=1.0, max_depression=5.0)[source]#

Filters the last few grid points after the physically-based inversion.

For various reasons (but mostly: the equilibrium assumption), the last few grid points on a glacier flowline are often noisy and create unphysical depressions. Here we try to correct for that. It is not volume conserving, but area conserving. If a parabolic shape factor is getting smaller than the minimum defined one (cfg.PARAMS[‘mixed_min_shape’]) the grid point is changed to a trapezoid, similar to what is done during the actual physically-based inversion.

Parameters
gdiroggm.GlacierDirectory

the glacier directory to process

n_smoothingint

number of grid points which should be smoothed. Default is 5

min_ice_thickfloat

the minimum ice thickness after the smoothing. Default is 1 m

max_depressionfloat

the limit allowed bed depression without smoothing. Default is 5 m

Notes

Files written to the glacier directory:
inversion_output.pkl

List of dicts containing the output data from the inversion.