@classmethod
def true_phase(cls, phase=None, phase_shift=None):
    """
        returns shifted phase of the orbit by the amount phase_shift

        :param phase: numpy.array
        :param phase_shift: numpy.float
        :return:
        """
    return phase + phase_shift
