_RL fracEnFreez
& fracEnMelt, fracEnFreez,
c1Frz = fracEnFreez**rec_pLaw
c2Frz = (1. _d 0 -fracEnFreez)**rec_pLaw
xxFrz = fracEnFreez-xxFrz
Fbot(i,j) = (1. _d 0 - fracEnFreez)*fzMlOc(i,j)
fracEnFreez= 0. _d 0
WRITE(iUnit,*) 'ThSI: fracEnFreez=',fracEnFreez
& fracEnMelt, fracEnFreez, hThinIce, hThickIce, hNewIceMax