bqe = - eps1*eps6 -sLoc*a0*eps3 + eps3*eps7 + eps2
& eps1 * ( thetaFreeze - tLoc ) /
_RL eps1, eps2, eps3, eps4, eps5, eps6, eps7
eps1 = rUnit2mass*HeatCapacity_Cp*ICEFRONTheatTransCoeff
aqe = a0 *(-eps1+eps3)
eps1 = rUnit2mass*HeatCapacity_Cp*SHELFICEheatTransCoeff
aqe = a0 *(eps1+eps3)
bqe = eps1*eps6 + eps3*eps7 - eps2
print *, 'ml-shelfice: eps1 = ', eps1
eps1 = 1.0
eps1 = gnorm / gnorm0
if (eps1 .lt. epsg) then
epsg = eps1
double precision r1, tmin, tmax, tact, gnorm, gnorm0, eps1