ICEFRONTthetaSurface = SHELFICEthetaSurface
& SHELFICEthetaSurface,
_RL SHELFICEthetaSurface
SHELFICEthetaSurface = - 20.0 _d 0
eps7 = eps4 - SHELFICEthetaSurface
& SHELFICEthetaSurface
& ( eps3*( thetaFreeze - SHELFICEthetaSurface )