_RL newIce
newIce = esurp*deltaTice/qicAv
IF ( newIce.GT.hIceMin*iceMaskMin ) THEN
iceThick = MIN(hThinIce,newIce/iceMaskMin)
iceThick = MAX(iceThick,newIce/iceMaskMax)
iceFrac = newIce/iceThick
iceFormed = newIce
iceFrac = MIN(icFrac(i,j)+newIce/hNewIce,iceMaskMax)
iceThick = MIN(hiMax,(iceVol+newIce)/iceFrac)
& iceThick, newIce, iceFrac-icFrac(i,j)