snowmodel_fenics

This python code uses the Finite Element library FENICS (via docker) to solve the one dimensional partial differential equations for heat and mass transfer in snow. The results are written in vtk format from which the paper figure is reproduced.