& gravity, recip_gravity, gBaro,
_RL gravity, recip_gravity
& *recip_gravity*recip_gravFacC(k)
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
gravFacC(k) = gravFacC(k)*recip_gravity
deepFacC(k) = ( rSphere + phiRef(2*k)*recip_gravity
deepFacF(k) = ( rSphere + phiRef(2*k-1)*recip_gravity
& phiRef(k)*recip_gravity
& )*recip_gravity
& *recip_gravity
zHalf3d(:,:,k) = phiRef(2*ki-1)*recip_gravity
dm=100.*(p0-pt)*recip_gravity
dz_bottom = dz_bottom*recip_gravity*recip_rhoConst
if(gencost_name(k).EQ.'gmbp') fac = 1. _d 2 * recip_gravity