oggm.tasks.compute_inversion_velocities

oggm.tasks.compute_inversion_velocities(gdir, settings_filesuffix='', input_filesuffix='', output_filesuffix='', glen_a=None, fs=None, with_sliding=None)

Surface velocities along the flowlines from inverted ice thickness.

Computed following the methods described in Cuffey and Paterson (2010) Eq. 8.35, pp 310:

u_s = u_basal + (2A/n+1)* tau^n * H

In the case of no sliding (or if with_sliding=False, which is a justifiable simplification given uncertainties on basal sliding):

u_z/u_s = [n+1]/[n+2] (= 0.8 if n = 3).

The output is written in ‘inversion_output.pkl’ in m yr-1

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

with_slidingbool

if set to True, we will compute velocities the sliding component. the default is to check if sliding was used for the inversion (fs != 0)

glen_afloat

Glen’s deformation parameter A. Defaults to PARAMS[‘inversion_glen_a’]

fsfloat

sliding paramter, defaults to PARAMS[‘inversion_fs’]

Notes

Files written to the glacier directory:

inversion_output.pkl

List of dicts containing the output data from the inversion.