dRho = rhoNil-rhoConst
rhoLoc=rhoNil* (
rhoLoc(i,j)=rhoNil*(
rholoc = rhoNil*(
& ) + rhoNil
IF ( rhoConst .EQ. UNSET_RL ) rhoConst=rhoNil
IF ( rhoNil .LE. 0. ) THEN
recip_rhoNil = 1. _d 0 / rhoNil
rhoNil = 999.8 _d 0
& *rhoNil*recip_rhoConst*gravity
$ *HeatCapacity_Cp*rhoNil
& * rhoNil/1000.