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. The dataset contains the code and the output data to reproduce the key Figure 5 from the related publication: _Schürholt, K., Kowalski, J., Löwe, H.; Elements of future snowpack modeling - Part 1: A physical instability arising from the non-linear coupling of transport and phase changes, The Cryosphere, 2022_ The code and potential updates can be accessed directly through git via: https://gitlabext.wsl.ch/snow-physics/snowmodel_fenics