tgrounded_cell_fraction.mac - pism - [fork] customized build of PISM, the parallel ice sheet model (tillflux branch)
 (HTM) git clone git://src.adamsgaard.dk/pism
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) LICENSE
       ---
       tgrounded_cell_fraction.mac (680B)
       ---
            1 /* This Maxima script derives H[0] and H[1], ice thicknesses at two
            2 adjacent grid points needed to achieve the "LI" grounding line
            3 position parameterization lambda defined below.
            4 
            5 z[0] is the sea level elevation, mu = rho_ice / rho_ocean
            6 
            7 I assume that the bed is flat and that the bed elevation is zero.
            8 */
            9 
           10 /* Assume that H[1] is half of H[0] */
           11 H[1] : H[0] / 2;
           12 
           13 /* definitions of alpha and beta */
           14 alpha : mu * H[0] - z[0];
           15 beta  : mu * H[1] - z[0];
           16 
           17 /* lambda, as derived from the linear approximation of ice thickness. */
           18 lambda : alpha / (alpha - beta);
           19 
           20 /* Redefine lambda_target and re-run to recompute H[0]. */
           21 lambda_target : 0.5;
           22 
           23 solve(lambda = lambda_target, H[0]);