& + rhoCen(k)*gravity*gravFacF(k)*dRlocM
& + rhoCen(k)*gravity*gravFacF(k+1)*dRlocP
U rhoCen,
dRho = rhoLoc(k)-rhoCen(k)
rhoCen(k) = rhoLoc(k)
_RL rhoCen(Nr)
& + rhoCen(k)*gravity*gravFacC(k)*drF(k)*halfRL
& + rhoCen(k)*gravity*gravFacC(k)*drF(k)