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
& cpice, cpwater,
_RL cpice
qbot = -cpice *Tf + Lfresh
WRITE(iUnit,*) 'ThSI: cpice =',cpice
cpice = 2106. _d 0
& + cpice *(Tmlt1-Tf) + Lfresh*(1. _d 0-Tmlt1/Tf)
qicen(2)= -cpice *Tf + Lfresh