bqe = - eps1*eps6 -sLoc*a0*eps3 + eps3*eps7 + eps2
cqe = -(eps2+eps3*eps7)*sLoc
_RL eps1, eps2, eps3, eps4, eps5, eps6, eps7
eps2 = rUnit2mass*ICEFRONTlatentHeat*ICEFRONTsaltTransCoeff
eps2 = rUnit2mass*SHELFICElatentHeat*SHELFICEsaltTransCoeff
bqe = eps1*eps6 + eps3*eps7 - eps2
cqe = eps2*sLoc(I,J)
print *, 'ml-shelfice: eps2 = ', eps2