tdiffusion_only.mac - pism - [fork] customized build of PISM, the parallel ice sheet model (tillflux branch)
HTML git clone git://src.adamsgaard.dk/pism
DIR Log
DIR Files
DIR Refs
DIR LICENSE
---
tdiffusion_only.mac (1165B)
---
1 /* This corresponds to the three-equation model for the diffusion-only
2 case (equation 21) in
3 @article{HollandJenkins1999,
4 title={Modeling thermodynamic ice-ocean interactions at the base of an ice shelf},
5 author={Holland, David M and Jenkins, Adrian},
6 journal={Journal of Physical Oceanography},
7 volume={29},
8 number={8},
9 pages={1787--1800},
10 year={1999}
11 }
12 */
13
14 /* define all necessary variables */
15 load("common_definitions.mac")$
16
17 /* add some shortcuts */
18 h : ice_thickness $
19 c_pW : ocean_heat_capacity $
20 c_pI : ice_heat_capacity $
21 rho_I : ice_density $
22 rho_W : ocean_density $
23 kappa : ice_thermal_diffusivity $
24
25 /* Temperature gradient at the base of the ice; assume linear
26 temperature profile through the thickness of the shelf */
27 grad_T : delta_T / h $
28
29 B_simplified : gamma_s*L + gamma_t*c_pW*(Theta_W - b[2]*h - b[1]) + rho_I*c_pI*kappa*(T_S - a[2]*h - a[1]) / (h * rho_W) $
30
31 /* solve salinity and heat flux balance equations and find
32 coefficients of the quadratic equation for the shelf base salinity */
33 load("flux_balance.mac")$
34
35 print("Checking the simplified form of B; this should print 0:", expand(B - B_simplified))$
36
37