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
- gdir
oggm.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
- gdir
Notes
- Files written to the glacier directory:
- inversion_output.pkl
List of dicts containing the output data from the inversion.