fig, axes = plt.subplots(1, 2, figsize=(10, 5),
                         sharex=True, sharey=True,
                         constrained_layout=True)

i1 = np.abs(w0.v_z.value - 120).argmin()
i2 = np.abs(w0.v_z.value - 100).argmin()
orbits[:, i1].cylindrical.plot(['rho', 'z'], alpha=0.5, marker=',', axes=[axes[0]]);
orbits[:, i2].cylindrical.plot(['rho', 'z'], alpha=0.5, marker=',', axes=[axes[1]]);