_RL hThinIce
& hThinIce, hThickIce, hNewIceMax
aMlt = (c1Mlt+c2Mlt)/(hThickIce-hThinIce)
hMlt = hThinIce+c2Mlt/aMlt
aFrz = (c1Frz+c2Frz)/(hThickIce-hThinIce)
hFrz = hThinIce+c2Frz/aFrz
aMlt = -1. _d 0 /(hThickIce-hThinIce)
aFrz = -1. _d 0 /(hThickIce-hThinIce)
ELSEIF (hi.GE.hThinIce) THEN
ELSEIF (hi.LT.hThinIce) THEN
iceThick = MIN(hThinIce,newIce/iceMaskMin)
hThinIce = 0.2 _d 0
WRITE(iUnit,*) 'ThSI: hThinIce =',hThinIce
& fracEnMelt, fracEnFreez, hThinIce, hThickIce, hNewIceMax