//  [python]
def test_orbit_method_inputro_quantity():
    from astropy import units

    from galpy.orbit import Orbit
    from galpy.potential import MWPotential2014

    o = Orbit([1.1, 0.1, 1.1, 0.2, 0.3, 0.3])
    ro = 11.0
    assert (
        numpy.fabs(
            o.E(pot=MWPotential2014, ro=ro * units.kpc)
            - o.E(pot=MWPotential2014, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method E does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.ER(pot=MWPotential2014, ro=ro * units.kpc)
            - o.ER(pot=MWPotential2014, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method ER does not return the correct value as Quantity"
    assert (
        numpy.fabs(
            o.Ez(pot=MWPotential2014, ro=ro * units.kpc)
            - o.Ez(pot=MWPotential2014, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method Ez does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.Jacobi(pot=MWPotential2014, ro=ro * units.kpc)
            - o.Jacobi(pot=MWPotential2014, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method Jacobi does not return the correct value when input ro is Quantity"
    assert numpy.all(
        numpy.fabs(
            o.L(pot=MWPotential2014, ro=ro * units.kpc)
            - o.L(pot=MWPotential2014, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method L does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.rap(pot=MWPotential2014, analytic=True, ro=ro * units.kpc)
            - o.rap(pot=MWPotential2014, analytic=True, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method rap does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.rperi(pot=MWPotential2014, analytic=True, ro=ro * units.kpc)
            - o.rperi(pot=MWPotential2014, analytic=True, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method rperi does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.zmax(pot=MWPotential2014, analytic=True, ro=ro * units.kpc)
            - o.zmax(pot=MWPotential2014, analytic=True, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method zmax does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.jr(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro * units.kpc)
            - o.jr(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method jr does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.jp(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro * units.kpc)
            - o.jp(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method jp does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.jz(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro * units.kpc)
            - o.jz(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method jz does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.wr(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro * units.kpc)
            - o.wr(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method wr does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.wp(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro * units.kpc)
            - o.wp(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method wp does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.wz(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro * units.kpc)
            - o.wz(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method wz does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.Tr(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro * units.kpc)
            - o.Tr(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method Tr does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.Tp(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro * units.kpc)
            - o.Tp(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method Tp does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.Tz(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro * units.kpc)
            - o.Tz(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method Tz does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.Or(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro * units.kpc)
            - o.Or(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro)
        )
        < 10.0**-8.0
    ), "Orbit method Or does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.Op(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro * units.kpc)
            - o.Op(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro)
        )
        < 10.0**-8.0
    ), "Opbit method Or does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(
            o.Oz(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro * units.kpc)
            - o.Oz(pot=MWPotential2014, type="staeckel", delta=0.5, ro=ro)
        )
        < 10.0**-8.0
    ), "Ozbit method Or does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.time(ro=ro * units.kpc) - o.time(ro=ro)) < 10.0**-8.0
    ), "Orbit method time does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.R(ro=ro * units.kpc) - o.R(ro=ro)) < 10.0**-8.0
    ), "Orbit method R does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.vR(ro=ro * units.kpc) - o.vR(ro=ro)) < 10.0**-8.0
    ), "Orbit method vR does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.vT(ro=ro * units.kpc) - o.vT(ro=ro)) < 10.0**-8.0
    ), "Orbit method vT does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.z(ro=ro * units.kpc) - o.z(ro=ro)) < 10.0**-8.0
    ), "Orbit method z does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.vz(ro=ro * units.kpc) - o.vz(ro=ro)) < 10.0**-8.0
    ), "Orbit method vz does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.phi(ro=ro * units.kpc) - o.phi(ro=ro)) < 10.0**-8.0
    ), "Orbit method phi does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.vphi(ro=ro * units.kpc) - o.vphi(ro=ro)) < 10.0**-8.0
    ), "Orbit method vphi does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.x(ro=ro * units.kpc) - o.x(ro=ro)) < 10.0**-8.0
    ), "Orbit method x does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.y(ro=ro * units.kpc) - o.y(ro=ro)) < 10.0**-8.0
    ), "Orbit method y does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.vx(ro=ro * units.kpc) - o.vx(ro=ro)) < 10.0**-8.0
    ), "Orbit method vx does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.vy(ro=ro * units.kpc) - o.vy(ro=ro)) < 10.0**-8.0
    ), "Orbit method vy does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.theta(ro=ro * units.kpc) - o.theta(ro=ro)) < 10.0**-8.0
    ), "Orbit method theta does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.vtheta(ro=ro * units.kpc) - o.vtheta(ro=ro)) < 10.0**-8.0
    ), "Orbit method vtheta does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.vr(ro=ro * units.kpc) - o.vr(ro=ro)) < 10.0**-8.0
    ), "Orbit method vr does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.ra(ro=ro * units.kpc) - o.ra(ro=ro)) < 10.0**-8.0
    ), "Orbit method ra does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.dec(ro=ro * units.kpc) - o.dec(ro=ro)) < 10.0**-8.0
    ), "Orbit method dec does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.ll(ro=ro * units.kpc) - o.ll(ro=ro)) < 10.0**-8.0
    ), "Orbit method ll does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.bb(ro=ro * units.kpc) - o.bb(ro=ro)) < 10.0**-8.0
    ), "Orbit method bb does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.dist(ro=ro * units.kpc) - o.dist(ro=ro)) < 10.0**-8.0
    ), "Orbit method dist does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.pmra(ro=ro * units.kpc) - o.pmra(ro=ro)) < 10.0**-8.0
    ), "Orbit method pmra does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.pmdec(ro=ro * units.kpc) - o.pmdec(ro=ro)) < 10.0**-8.0
    ), "Orbit method pmdec does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.pmll(ro=ro * units.kpc) - o.pmll(ro=ro)) < 10.0**-8.0
    ), "Orbit method pmll does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.pmbb(ro=ro * units.kpc) - o.pmbb(ro=ro)) < 10.0**-8.0
    ), "Orbit method pmbb does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.vlos(ro=ro * units.kpc) - o.vlos(ro=ro)) < 10.0**-8.0
    ), "Orbit method vlos does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.vra(ro=ro * units.kpc) - o.vra(ro=ro)) < 10.0**-8.0
    ), "Orbit method vra does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.vdec(ro=ro * units.kpc) - o.vdec(ro=ro)) < 10.0**-8.0
    ), "Orbit method vdec does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.vll(ro=ro * units.kpc) - o.vll(ro=ro)) < 10.0**-8.0
    ), "Orbit method vll does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.vbb(ro=ro * units.kpc) - o.vbb(ro=ro)) < 10.0**-8.0
    ), "Orbit method vbb does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.helioX(ro=ro * units.kpc) - o.helioX(ro=ro)) < 10.0**-8.0
    ), "Orbit method helioX does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.helioY(ro=ro * units.kpc) - o.helioY(ro=ro)) < 10.0**-8.0
    ), "Orbit method helioY does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.helioZ(ro=ro * units.kpc) - o.helioZ(ro=ro)) < 10.0**-8.0
    ), "Orbit method helioZ does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.U(ro=ro * units.kpc) - o.U(ro=ro)) < 10.0**-8.0
    ), "Orbit method U does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.V(ro=ro * units.kpc) - o.V(ro=ro)) < 10.0**-8.0
    ), "Orbit method V does not return the correct value when input ro is Quantity"
    assert (
        numpy.fabs(o.W(ro=ro * units.kpc) - o.W(ro=ro)) < 10.0**-8.0
    ), "Orbit method W does not return the correct value when input ro is Quantity"
    return None