SNOWPACK is a multi-purpose snow and land-surface model, which focuses on a detailed description of the mass and energy exchange between the snow, the atmosphere and optionally with the vegetation cover and the soil.
snowpackBuoyantPimpleFoam is a two-phase solver implemented to model convection of water vapor with phase change in snowpacks. This new solver is based on the standard solver of buoyantPimpleFoam in OpenFOAM 5.
StreamFlow is an extension to the spatially-distributed snow model Alpine3D which allows the user to perform hydrological simulations. On top of discharge and water height, StreamFlow can also compute stream temperature at any point along the stream.
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.