fig = orbit.cylindrical.plot(['rho', 'z'])