eps3a = rhoShelfIce*SHELFICEheatCapacity_Cp
eps3 = eps3a/pLoc(I,J)
& + eps3a*( b
_RL eps1, eps2, eps3, eps3a, eps4, eps5, eps6, eps7, eps8