if np.abs(omega) > M_PI:
        raise ValueError(
            """longitude of sun ({} deg) not between -180 and """
            """+180 degrees""".format(omega*180/M_PI))