changed initial_zenith convention to be ra/dec in radians. test passes. trying to deal with the spherical harmonics problem in scipy