oggm.tasks.compute_velocities¶
-
oggm.tasks.
compute_velocities
(gdir, glen_a=None, fs=None, filesuffix='')[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 * HIn the case of no 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
You’ll need to call prepare_for_inversion with the add_debug_var=True kwarg for this to work!
Parameters: - gdir : Glacier directory
- with_sliding : bool
default is True, if set to False will not add the sliding component.
- filesuffix : str
add a suffix to the output file
Notes
- Files writen to the glacier directory:
- inversion_output.pkl
- List of dicts containing the output data from the inversion.