def METHOD_NAME(self):
    """
    Test orbit method to ensure proper output
    """
    t_ref = Time(2000.0, format="jyear")
    for mod in self.allmods:
        with RedirectStreams(stdout=self.dev_null):
            if "SotoStarshade" in mod.__name__:
                obj = mod(f_nStars=4, **copy.deepcopy(self.spec))
            else:
                obj = mod(**copy.deepcopy(self.spec))
        r_sc = obj.orbit(t_ref)
        # the r_sc attribute is set and is a 3-tuple of astropy Quantity's
        self.assertEqual(type(r_sc), type(1.0 * u.km))
        self.assertEqual(r_sc.shape, (1, 3))