& gravity, recip_gravity, gBaro,
_RL recip_gravity
scalingFactor=rhoConst*gravitySign*recip_gravity
scalingFactor = recip_rhoConst*recip_gravity
scalingFactor = tRef(k)*recip_gravity
& shelficeLoadAnomaly(i,j,bi,bj)*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
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
& )*recip_gravity
zHalf3d(:,:,k) = phiRef(2*ki-1)*recip_gravity
& *recip_gravity
dm=100.*(p0-pt)*recip_gravity
dz_bottom = dz_bottom*recip_gravity*recip_rhoConst