_RL cpIce
& cpIce, cpWater,
qbot = -cpIce *Tf + Lfresh
& + cpIce *(Tmlt1-Tf) + Lfresh*(1. _d 0-Tmlt1/Tf)
qicen(2)= -cpIce *Tf + Lfresh
Qice2(i,j,bi,bj) = -cpIce *Tf + Lfresh
WRITE(iUnit,*) 'ThSI: cpIce =',cpIce
cpIce = 2106. _d 0
a1 = cpIce
b1 = qicen(1) + (cpWater-cpIce )*Tmlt1 - Lfresh
Tice(2) = (Lfresh-qicen(2)) / cpIce
a10 = rhoi*cpIce *hi/(2. _d 0*dt) +
& k32 * (4. _d 0*dt*k32 + rhoi*cpIce *hi)
& / (6. _d 0*dt*k32 + rhoi*cpIce *hi)
& (rhoi*cpIce*Tice(1)+rhoi*Lfresh*Tmlt1/Tice(1))
& - k32 * (4. _d 0*dt*k32*Tf+rhoi*cpIce *hi*Tice(2))
& / (6. _d 0*dt*k32 + rhoi*cpIce *hi) - fswint
& + rhoi*cpIce *hi*Tice(2))
& /(6. _d 0*dt*k32 + rhoi*cpIce *hi)
qicen(1) = -cpWater*Tmlt1 + cpIce *(Tmlt1-Tice(1))
qicen(2) = -cpIce *Tice(2) + Lfresh