b1 = qicen(1) + (cpwater-cpice )*Tmlt1 - Lfresh
c1 = Lfresh * Tmlt1
& (rhoi*cpice*Tice(1)+rhoi*Lfresh*Tmlt1/Tice(1))
c1 = rhoi*Lfresh*hi*Tmlt1 / (2. _d 0*dt)
qicen(1) = -cpwater*Tmlt1 + cpice *(Tmlt1-Tice(1)) +
& Lfresh*(1. _d 0-Tmlt1/Tice(1))
if (Tice(1) .ge. Tmlt1) then
& 'BBerr - Bug: IceT(1) > Tmlt',i,j,bi,bj,Tice(1),Tmlt1
_RL Tmlt1
& Tf0kel, Tmlt1,
Tmlt1=-mu_Tf*S_winton
WRITE(iUnit,*) 'ThSI: Tmlt1 =',Tmlt1
qicen(1)= -cpwater*Tmlt1
& + cpice *(Tmlt1-Tf) + Lfresh*(1. _d 0-Tmlt1/Tf)