def METHOD_NAME(self):
    # Find a conjunction which we have confirmed by GBT data and Shapiro delay
    mJ0737 = m.get_model("0737A_latest.par")
    x = mJ0737.conjunction(55586.25)
    assert np.isclose(x, 55586.29643451057), "J0737 conjunction time is bad"
    # And now make sure we calculate the true anomaly for it correctly
    nu = mJ0737.orbital_phase(x, anom="true").value
    omega = mJ0737.components["BinaryDD"].binary_instance.omega().value
    # Conjunction occurs when nu + OM == 90 deg
    assert np.isclose(
        np.degrees(np.remainder(nu + omega, 2 * np.pi)), 90.0
    ), "J0737 conjunction time gives bad true anomaly"
    # Now verify we can get 2 results from .conjunction
    x = mJ0737.conjunction([55586.0, 55586.2])
    assert len(x) == 2, "conjunction is not returning an array"
    # make sure true anomaly before T0 is positive
    assert mJ0737.orbital_phase(52000.0, anom="true").value > 0.0