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