Solving the equation for thickness evolution

From Interactive System for Ice sheet Simulation
Revision as of 10:18, 5 August 2009 by Sprice (Talk | contribs)

Jump to: navigation, search

A generic conservation equation

Evolution of the ice sheet geometry is simply the result of mass conservation. Consider a generic conservation equation for some scalar value φ


\frac{\partial \varphi }{\partial t}=-\nabla \cdot Q+S_{\varphi }


On the left-hand side we have the change in φ per unit time, which is given by the negative of the divergence (the \nabla \cdot operation) of the flux, Q, plus any contribution from sources, S_{\varphi }.