Published December 12, 2021 | Version v1.0.0
Software Open

agstub/subglacial-inversion: First release of subglacial parameter inversion code

Description

This python code inverts ice-sheet altimetry data for two types of subglacial parameters: the basal vertical velocity or basal drag coefficient field. The forward model is based on a small-perturbation approximation of the Stokes equations, and the least-squares inverse problem is solved with the conjugate gradient method. The primary operations involve map-plane Fourier transforms and convolution over time, which are implemented with SciPy's FFT methods. The complete details are included in a forthcoming manuscript. The Jupyter notebooks in the "notebooks" directory reproduce the computational examples from the manuscript.

Files

agstub/subglacial-inversion-v1.0.0.zip

Files (375.3 kB)

Name Size Download all
md5:a092656abed89f0834b37f127771c155
375.3 kB Preview Download

Additional details