_RL Fbot
Fbot = frzmlt
Fbot = enFrcFrz*frzmlt
Fbot = 0. _d 0
Fbot = (1. _d 0 - fracEnFreez)*frzmlt
Fbot = cpchr*(Tf-oceTs)*ustar
Fbot = max(Fbot,frzmlt)
Fbot = min(Fbot,0. _d 0)
& 'ThSI_CALC_TH: evpAtm, frzmlt, Fbot =', frwAtm(i,j),frzmlt,Fbot
ebot = (flxCnB(i,j)-Fbot) * dt
qleft = qleft + (Fbot+(esurp+etop+ebot)/dt)
& Qleft,Fbot,(etope+ebote)/dt
fzMlOc(i,j) = frzmlt - Fbot*iceMask(i,j)