specfun = interpolate.interp1d(w, S1, kind='cubic')