MATLAB code for firn thermodynamic and hydrological modelling

The default configuration of the subsurface model is for a 35-m firn column, with 0.1-m thick layers from 0-0.6 m, 0.2-m thick layers from 0.6-2 m, 0.4-m thick layers from 2-10 m, and 1-m thick layers below that. The model requires forcing from mean daily automatic weather station, climate model, or climate reanalysis data, with the default inputs being: minimum, mean, and maximum daily air temperature, mean daily relative or specific humidity, wind speed, air pressure, and longwave radiation, and mean and maximum daily incoming incoming shortwave.

The model internally calculates snow albedo and snow surface temperature (hence, outgoing longwave radiation), as well as conductive heat flux to the snow surface, based on temperature gradients in the top three layers of the subsurface snow/firn model. See the referenced papers for further details on the model physics and parameterizations.