def prepare_earth_position_vel(time):
 from ..solar_system import get_body_barycentric, get_body_barycentric_posvel
 earth_pv = get_body_barycentric_posvel(u'earth', time)
 sun = get_body_barycentric(u'sun', time)
 earth_heliocentric = (earth_pv[0] - sun).get_xyz(xyz_axis=(-1)).to(u.au).value
 earth_pv = np.concatenate((earth_pv[0].get_xyz(xyz_axis=(-1)).to(u.au)[..., np.newaxis, :].value, earth_pv[1].get_xyz(xyz_axis=(-1)).to((u.au / u.d))[..., np.newaxis, :].value), axis=(-2))
 return (earth_pv, earth_heliocentric)

def prepare_earth_position_vel ( time ) : from . . solar_system import get_body_barycentric , get_body_barycentric_posvel earth_pv = get_body_barycentric_posvel ( u'earth' , time ) sun = get_body_barycentric ( u'sun' , time ) earth_heliocentric = ( earth_pv [ 0 ] - sun ) . get_xyz ( xyz_axis = ( - 1 ) ) . to ( u . au ) . value earth_pv = np . concatenate ( ( earth_pv [ 0 ] . get_xyz ( xyz_axis = ( - 1 ) ) . to ( u . au ) [ ... , np . newaxis , : ] . value , earth_pv [ 1 ] . get_xyz ( xyz_axis = ( - 1 ) ) . to ( ( u . au / u . d ) ) [ ... , np . newaxis , : ] . value ) , axis = ( - 2 ) ) return ( earth_pv , earth_heliocentric )
get barycentric position and velocity .