get heliocentric positions and velocities of planets
if moon is being integrated simultaneously, the integrated quantity
is the geocentric moon position, and we still have the heliocentric
position of the embary as the "planet" position being integrated

this block of code used when earth and moon
effects on other bodies are considered separately
here, we set the moon position and velocity vectors and replace the
embary position and velocity with the corresponding earth quantities