This is a project on mass-conserving layer models which have a nonnegativity constraint on the conserved thickness variable. This situation is especially important for climate applications such as ice sheets, sea ice, subglacial water, shallow water, etc.
This repository includes a paper, several talks, and a 1D example C code which uses PETSc. The paper has been accepted by SIAM J. Applied Math. It is available in preprint form:
- E. Bueler (2020). Conservation laws for free-boundary fluid layers, arxiv:2007.05625
A nontrivial ice sheet example was moved to a separate repo and published:
- E. Bueler (2016). Stable finite volume element schemes for the shallow-ice approximation, J. Glaciol. 62 (232), 230-242, doi:10.1017/jog.2015.3
To run the C code see petsc/README.md.