eps3 = rUnit2mass*ICEFRONTheatCapacity_Cp
aqe = a0 *(-eps1+eps3)
bqe = - eps1*eps6 -sLoc*a0*eps3 + eps3*eps7 + eps2
cqe = -(eps2+eps3*eps7)*sLoc
_RL eps1, eps2, eps3, eps4, eps5, eps6, eps7
eps3 = eps3a/pLoc(I,J)
aqe = a0 *(eps1+eps3-eps8)
print *, 'ml-shelfice: eps3 = ', eps3
& ( eps3
_RL eps1, eps2, eps3, eps3a, eps4, eps5, eps6, eps7, eps8