_RL Tf
Tf = tFrz(i,j)
Fbot = cpchr*(Tf-oceTs)*ustar
qbot = -cpIce *Tf + Lfresh
& + cpIce *(Tmlt1-Tf) + Lfresh*(1. _d 0-Tmlt1/Tf)
qicen(2)= -cpIce *Tf + Lfresh
Tf = -1.70 _d 0
Qice2(i,j,bi,bj) = -cpIce *Tf + Lfresh
& - k32 * (4. _d 0*dt*k32*Tf+rhoi*cpIce *hi*Tice(2))
Tice(2) = (2. _d 0*dt*k32*(Tice(1)+2. _d 0*Tf)
flxCnB(i,j) = 4. _d 0*kIce *(Tice(2)-Tf)/hi