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