Serguei pointed out that the Alvi term is really only meant to be valid for equatorial orbits. Adjust as necessary.