& rhoNil, rhoConst, recip_rhoConst, thetaConst,
_RL rhoNil
CALL WRITE_0D_RL( rhoNil, INDEX_NONE, 'rhoNil =',
dRho = rhoNil-rhoConst
rhoLoc(i,j)=rhoNil*(
rholoc = rhoNil*(
& ) + rhoNil
IF ( rhoConst .EQ. UNSET_RL ) rhoConst=rhoNil
IF ( eosType.EQ.'LINEAR' .AND. rhoNil.LE.0. ) THEN
rhoNil = 999.8 _d 0
& *rhoNil*recip_rhoConst*gravity
& * (HeatCapacity_Cp*rhoNil*drF(1))*nearIce
$ *HeatCapacity_Cp*rhoNil
& HeatCapacity_Cp*rhoNil*drF(1)
& - (epsilonTemp*HeatCapacity_Cp*rhoNil*drF(1)))
& HeatCapacity_Cp*rhoNil*drF(1) )
& * rhoNil/1000.