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