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