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