scipy.interpolate.RectSphereBivariateSpline
scipy.interpolate.RectSphereBivariateSpline.ev
Evaluate the spline at the grid ponts defined by the coordinate arrays theta, phi.