& rhoNil, rhoConst, recip_rhoConst, rho1Ref,
_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*gravFacF(k)
& * (HeatCapacity_Cp*rhoNil*drF(1))*nearIce
const1=HeatCapacity_Cp*rhoNil*drF(1)