agstub/subglacial-inversion: First release of subglacial parameter inversion code
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 co...
Main Author: | |
---|---|
Format: | Software |
Language: | unknown |
Published: |
Zenodo
2021
|
Subjects: | |
Online Access: | https://dx.doi.org/10.5281/zenodo.5775177 https://zenodo.org/record/5775177 |
Summary: | 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. |
---|