Evaluate the spline at the points (theta[i], phi[i]), i=0,...,len(theta)-1
scipy.interpolate.RectSphereBivariateSpline.__call__
scipy.interpolate.RectSphereBivariateSpline.get_coeffs