b1 = qicen(1) + (cpwater-cpice )*Tmlt1 - Lfresh
qicen(1) = -cpwater*Tmlt1 + cpice *(Tmlt1-Tice(1)) +
& cpice, cpwater,
_RL cpwater
cpchr =cpwater*rhosw*transcoef
WRITE(iUnit,*) 'ThSI: cpwater =',cpwater
cpwater = HeatCapacity_Cp
qicen(1)= -cpwater*Tmlt1
cphm = cpwater*rhosw*hOceMxL(i,j,bi,bj)
heatFac = ocean_deltaT/(cpwater*rhosw)