fenics_ice

fenics_ice is a Python code intended for Bayesian quantification of uncertainty of large-scale parameter sets within an ice-sheet flow model. the code contains a dynamic solver for ice-sheet flow which implements the Shallow-Shelf Approximation (MacAyeal et al, 1989, ) and approximates parameter covariance using a low-rank approximation to the inverse of the cost-function Hessian. The code uses Algorithmic Differentiation to find sensitivity of time-evolving Quantities of Interest to parameter sets, allowing projections of parameter uncertainty forward in time.

MacAyeal, D. R. (1989), Large‐scale ice flow over a viscous basal sediment: Theory and application to ice stream B, Antarctica, J. Geophys. Res., 94( B4), 4071– 4087.