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