& gravity, recip_gravity, gBaro,
_RL recip_gravity
scalingFactor=rhoConst*gravitySign*recip_gravity
rC(k) = -pLoc(kcheck)*recip_rhoConst*recip_gravity
dzLoc = dPiHvR(k)*thetaHvR(k)*recip_gravity
& * thetaHvR(kLev)*recip_gravity
hRef(K) = phiRef(K)*recip_gravity
rUnit2mass = recip_gravity
recip_gravity = 1. _d 0 / gravity
deepFacC(k) = ( rSphere + phiRef(2*k)*recip_gravity
deepFacF(k) = ( rSphere + phiRef(2*k-1)*recip_gravity
& phiRef(k)*recip_gravity
dz_bottom = dz_bottom*recip_gravity*recip_rhoConst