& EnergyInNewTotalIceVolume(i,j)
_RL EnergyInNewTotalIceVolume (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
EnergyInNewTotalIceVolume (I,J) = 0.0
EnergyInNewTotalIceVolume(I,J) =
& EnergyInNewTotalIceVolume(I,J)
& EnergyInNewTotalIceVolume(I,J)*ENERGY_TO_DELTA_TEMP