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