This zipped file contains a MatLab script and functions to model an isothermal ice cap.
The ice-cap code is adapted from the transient 2-D advective-diffusion code.
Basal sliding constitutes advection, and internal deformation through Glen's law using the shallow Ice Approximation (SIA) produces the nonlinear diffusion term.
The code still has a bug - as the dome approaches a steady state, the surface develops un-physical undulations on the scale of 2 finite volumes.