def setup_orbit(): b = Bunch( equatorial_radius = 6378*km, gravitational_parameter = 398600*km**3, ) o = Orbit(body=b) o.e = 0.3 o.h = 60000*km**2 o.θ0 = 120*deg return o