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