_RL cpWater
& cpIce, cpWater,
cpchr =cpWater*rhosw*bMeltCoef
qicen(1)= -cpWater*Tmlt1
Qice1(i,j,bi,bj) = -cpWater*Tmlt1
WRITE(iUnit,*) 'ThSI: cpWater =',cpWater
cpWater = HeatCapacity_Cp
b1 = qIc1(i,j) + (cpWater-cpIce )*Tmlt1 - Lfresh
qIc1(i,j) = -cpWater*Tmlt1 + cpIce *(Tmlt1-tIc1(i,j))