rho_ice 910[kg/m^3] "Density of ice"
mu_ice 5e12[Pa*s] "Dynamic viscosity of ice"
LSlip 50[m] "Slip length"
T_init -10[degC] "Initial temperature"
betaCC 9.8e-8[K/Pa] "Clapeyron constant"
n_ice 3 "Rheological exponent for ice"
T_tp 0.01[degC] "Temperature at triple point of water"
p_tp 611.657[Pa] "Pressure at triple point of water"
q_geo 120[mW/m^2] "Geothermal heat flux"