ICEFRONTthetaSurface = SHELFICEthetaSurface
& SHELFICEthetaSurface,
_RL SHELFICEthetaSurface
CALL WRITE_0D_RL( SHELFICEthetaSurface, INDEX_NONE,
SHELFICEthetaSurface = - 20.0 _d 0
eps7 = eps4 - SHELFICEthetaSurface
& + ( c0 - SHELFICEthetaSurface )/pLoc(I,J) )
& SHELFICEthetaSurface
& * ( thetaFreeze - SHELFICEthetaSurface )