FEniCS solver for coupled evolution of firn density and velocity in 1D (vertical column):
= 0)
)] - \nabla p = \rho \bf{g})
where ))
is the compressible (porous) rheology as documented in Gagliardini and Meyssonnier (1997), Lüthi and Funk (2000,2001), Boutillier (2004), Zwinger et al. (2007).