@classmethod
def true_phase(cls, phase, phase_shift):
    """
        Returns shifted phase of the orbit by the amount phase_shift.

        :param phase: Union[numpy.array, float];
        :param phase_shift: float;
        :return: Union[numpy.array, float];
        """
    return phase + phase_shift
