ConvertFac = HeatCapacity_Cp*rhoConst
ConvertFac = rhoConst
plumetend(i,j,k) = ConvertFac *
IF(k.eq.1) plumetend(i,j,k) = ConvertFac *
& SurfVal, SEAICE_Tfrz, ConvertFac, recip_ConvertFac