tsubshelf_freezeon.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
---
tsubshelf_freezeon.mac (706B)
---
1 /* This corresponds to the three-equation model for the shelf-base freeze-on,
2 case (equations 26 and 31) 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 /* Temperature gradient at the base of the ice; zero in the freeze-on
18 case. */
19 grad_T : 0 $
20
21 /* solve salinity and heat flux balance equations and find
22 coefficients of the quadratic equation for the shelf base salinity */
23 load("flux_balance.mac")$