oggm.tasks.filter_inversion_output

oggm.tasks.filter_inversion_output(gdir, settings_filesuffix='', input_filesuffix='', output_filesuffix='', n_smoothing=5, min_ice_thick=1.0, max_depression=5.0)

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 (gdir.settings[‘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

settings_filesuffix: str

You can use a different set of settings by providing a filesuffix. This is useful for sensitivity experiments. Code-wise the settings_filesuffix is set in the @entity-task decorater.

input_filesuffix: str

the filesuffix of the ‘inversion_output’ to use

output_filesuffix:str

the filesuffix of the resulting ‘inversion_output’ to save back into the gdir

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.