icepack/icepack-paper

icepack

icepack is a library for modeling the flow of ice sheets and glaciers using the finite element method. For more information and installation instructions, see the project webpage. Once you've installed icepack, the directory notebooks/tutorials/ contains several Jupyter notebooks with demonstrations of the main features of icepack to get you started. The directory notebooks/how-to contains more advanced material for when you've got a good feel for how things work. Some of these demos use real data, for which you'll need to have an account with NASA EarthData.

See also firedrake, a finite element modeling package used to implement the solvers in icepack.

icepack-paper's Issues

MISMIP+ demo is slow

The MISMIP+ demo requires very high resolution near the grounding line and probably also the shear margins. At coarser resolutions the simulation exhibits weird oscillatory artifacts and obvious mesh imprinting. Right now, the demo spins up to steady state on a coarse mesh, then uniformly refines and spins up again several times to get to the desired resolution. This gets a good simulation result, but it takes about a day and half to run to completion. We could probably get a way with a mesh that's much coarser upstream than near the grounding line.

Just to save the trouble of running that long simulation again, I've attached a tar archive containing the checkpoint files with the results. You can fool Make into thinking that these are more up to date than, say, the simulation script by doing touch *.h5.

mismip+.tar.gz

update the grants / funding acknowledgments section

We're missing a CryoSAT grant and possibly others.