def setup_orbit(): b = Bunch( equatorial_radius = 6378*km, gravitational_parameter = 398600*km**3, ) o = Orbit(body=b) o.pe = 9600*km o.ap = 21000*km o.epoch = 0 o.θ0 = 0 return o