def METHOD_NAME(self):
    """
    start with a pulsar in equatorial coordinates
    convert to Galactic and make sure the positions are consistent
    then apply the space motion to the equatorial object & convert to Galactic
    compare that to the Galactic object w/ space motion
    """
    # make a test SkyCoord object
    # make sure it has obstime and distance supplied
    # to use it for conversions as well
    J0613_icrs = self.modelJ0613.coords_as_ICRS()
    J0613_icrs_now = utils.add_dummy_distance(J0613_icrs)
    J0613_galactic = self.modelJ0613.coords_as_GAL()
    J0613_galactic_now = utils.add_dummy_distance(J0613_galactic)
    newepoch = self.modelJ0613.POSEPOCH.quantity.mjd + 100
    # what I get converting within astropy
    J0613_galactic_comparison = J0613_icrs_now.transform_to(
        astropy.coordinates.Galactic
    )
    sep = J0613_galactic_now.separation(J0613_galactic_comparison)
    msg = (
        "Equatorial to Galactic conversion for now failed with separation %.1e arcsec"
        % sep.arcsec
    )
    assert sep < 1e-9 * u.arcsec, msg
    J0613_icrs = self.modelJ0613.coords_as_ICRS(epoch=newepoch)
    # what I get converting within astropy
    J0613_galactic_comparison = J0613_icrs.transform_to(
        astropy.coordinates.Galactic
    )
    J0613_galactic_then = utils.remove_dummy_distance(
        J0613_galactic_now.apply_space_motion(
            new_obstime=astropy.time.Time(newepoch, scale="tdb", format="mjd")
        )
    )
    sep = J0613_galactic_then.separation(J0613_galactic_comparison)
    msg = (
        "Equatorial to Galactic conversion for +100d failed with separation %.1e arcsec"
        % sep.arcsec
    )
    assert sep < 1e-9 * u.arcsec, msg