Modification of this model, with heat flow added.
It sovles equations describing firn compaction, grain growth, heat flow, and age on an Eulerian grid using the model of lines.
The addition of heat flow is a work in progress.
The python directory contains a pythin version of the isothermal model (fcm.py).