_RL recip_Cp
recip_Cp = 1. _d 0 / HeatCapacity_Cp
& * recip_Cp * mass2rUnit
tmpFac = recip_Cp * mass2rUnit
& *recip_Cp*mass2rUnit
& recip_Cp*recip_rhoConst
& -Qsw(i,j,bi,bj)*recip_Cp*recip_rhoConst
& ) *recip_Cp*mass2rUnit
_RL thetaFreeze, saltFreeze, recip_Cp
& *recip_Cp*recip_rhoConst