def fref(rad,cosb,height,t_setx,r_setx,s_setx,sang): ref=reflectance(rad,cosb,t_setx,height,r_setx,s_setx,sang) return interpolate.RectBivariateSpline(t_setx,r_setx,ref)