FUNCTION EShift3, ext_file, obj eout = FltArr(101,195) ef = ext_file eft = ef[22,720:1690] wv0 = FindGen(195) * 0.005 + 1.52 wv2 = FindGen(1711) * 0.001 + 0.80 obf = obj[1,140:334] e1 = Interpol(eft, 9701, /SPLINE) o1 = Interpol(obf, 9701, /SPLINE) lag = IndGen(201) - 100 c = C_Correlate(o1, e1, lag) m = Max(c) offs = lag[!C] * 0.0001 print, 'Region 3 wavelength offset =', offs w_offs = wv0 + offs FOR j=1, 101 DO BEGIN aa = ef[j-1,*] ab = Interpol(aa, wv2, w_offs, /SPLINE) eout[j-1,*] = ab ENDFOR RETURN, eout END