FUNCTION EShift1, ext_file, obj eout = FltArr(101,107) ef = ext_file eft = ef[22,20:550] wv0 = FindGen(107) * 0.005 + 0.82 wv2 = FindGen(1711) * 0.001 + 0.80 obf = obj[1,0:106] e1 = Interpol(eft, 5301, /SPLINE) o1 = Interpol(obf, 5301, /SPLINE) lag = IndGen(201) - 100 c = C_Correlate(o1, e1, lag) m = Max(c) offs = lag[!C] * 0.0001 print, 'Region 1 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