- oggm.tasks.compute_inversion_velocities(gdir, glen_a=None, fs=None, filesuffix='', with_sliding=None)[source]#
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
the glacier directory to process
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’s deformation parameter A. Defaults to PARAMS[‘inversion_glen_a’]
sliding paramter, defaults to PARAMS[‘inversion_fs’]
add a suffix to the output file (optional)
- Files written to the glacier directory:
List of dicts containing the output data from the inversion.