[work in progress] this is an implementation of a 1-d coupled radiative-thermodynamic heat transfer model applicable to glacier ice written by Glen Liston, UCAR with updates by Matt Cooper, PNNL. Stay tuned for a completely overhauled version written in Matlab.
Liston, G. E., Bruland, O., Elvehøy, H., & Sand, K. (1999). Below-surface ice melt on the coastal Antarctic ice sheet. Journal of Glaciology, 45(150), 273–285. https://doi.org/10.3189/002214399793377130