Effect of snow on soil erosion - snowpack dynamics

The aim of the presented code is to calculate monthly snowpack dynamics and to identify monthly snowmelt to later incorporate its erosive forces.
The results are written as georeferenced raster files (.tif format) and can be accessed by all common GIS programs.
Optionally, the user can plot various statistical values (e.g. snow cover) over time.
If only a subarea is of interest, this can be defined by a shapefile. The previously calculated snowfall (mm) and the identified snow cover serve as input data.
The following figure provides a highly simplified overview. Details can be found in the code diagram and the code description.