eps4 = b*pLoc + c0
eps6 = eps4 - tLoc
eps7 = eps4 - thetaIce
thetaFreeze = a0*saltFreeze + eps4
_RL eps1, eps2, eps3, eps4, eps5, eps6, eps7
eps4 = b*pLoc(I,J) + c0
eps6 = eps4 - tLoc(I,J)
eps7 = eps4 - SHELFICEthetaSurface
print *, 'ml-shelfice: eps4 = ', eps4
_RL eps1, eps2, eps3, eps3a, eps4, eps5, eps6, eps7, eps8