Ice sheet weathering crust evolution code (1D, saturated)

This repository contains the code used to produce the results in the article "Modelling the evolution of an ice sheet's weathering crust" by Tilly Woods and Ian J. Hewitt, 2024, IMA Journal of Applied Mathematics, https://doi.org/10.1093/imamat/hxae031. The code finds time-dependent solutions for the porosity and temperature profiles in and below an ice sheet's weathering crust. The model is based on mass conservation, energy conservation, internal shortwave radiation and a surface energy balance. It is solved using an ethalpy method with finite volumes and semi-implicit timestepping. The code was written in MATLAB R2024a.