dir = '/Users/colkin/Data/12aug01_CD1/spec/'

s22 = readfits(dir+'12aus0022.fits')
s23 = readfits(dir+'12aus0023.fits')
s32 = readfits(dir+'12aus0032.fits')
s33 = readfits(dir+'12aus0033.fits')
s68 = readfits(dir+'12aus0068.fits')
s69 = readfits(dir+'12aus0069.fits')
s84 = readfits(dir+'12aus0084.fits')
s85 = readfits(dir+'12aus0085.fits')

superf = readfits('/Users/colkin/Work04/Pluto_3um/Reduction/Flatten/superFlat.fits')

.compile /Users/colkin/Work04/Pluto_3um/IDL_Pro/traceGaussians.pro 

tracegaussians, s22, s23, superf, res22, yFit22, a22, gauss22
tracegaussians, s23, s22, superf, res23, yFit23, a23, gauss23
tracegaussians, s32, s33, superf, res32, yFit32, a32, gauss32
tracegaussians, s33, s32, superf, res33, yFit33, a33, gauss33
tracegaussians, s68, s69, superf, res68, yFit68, a68, gauss68
tracegaussians, s69, s68, superf, res69, yFit69, a69, gauss69
tracegaussians, s84, s85, superf, res84, yFit84, a84, gauss84
tracegaussians, s85, s84, superf, res85, yFit85, a85, gauss85

plot, yFit22, yr=[0,1023]
oplot, yFit23
oplot, yFit32
oplot, yFit33
oplot, yFit68
oplot, yFit69
oplot, yFit84
oplot, yFit85

print, [yFit22[511], yFit23[511], yFit32[511], yFit33[511], yFit68[511], yFit69[511], yFit84[511], yFit85[511]]

print, [a22[1,511], a23[1,511], a32[1,511], a33[1,511], a68[1,511], a69[1,511], a84[1,511], a85[1,511]]

pts=fltarr(8192,3)

xArr = findgen(1024)

pts[0:1023, 0] = yfit22[511]
pts[0:1023, 1] = findgen(1024)
pts[0:1023, 2] = yfit22

frame = 1024

pts[frame:1023+frame, 0] = yfit23[511]
pts[frame:1023+frame, 1] = xArr
pts[frame:1023+frame, 2] = yfit23

pts[2*frame:1023+2*frame, 0] = yfit32[511]
pts[2*frame:1023+2*frame, 1] = xArr
pts[2*frame:1023+2*frame, 2] = yfit32

pts[3*frame:1023+3*frame, 0] = yfit33[511]
pts[3*frame:1023+3*frame, 1] = xArr
pts[3*frame:1023+3*frame, 2] = yfit33

pts[4*frame:1023+4*frame, 0] = yfit68[511]
pts[4*frame:1023+4*frame, 1] = xArr
pts[4*frame:1023+4*frame, 2] = yfit68

pts[5*frame:1023+5*frame, 0] = yfit69[511]
pts[5*frame:1023+5*frame, 1] = xArr
pts[5*frame:1023+5*frame, 2] = yfit69

pts[6*frame:1023+6*frame, 0] = yfit84[511]
pts[6*frame:1023+6*frame, 1] = xArr
pts[6*frame:1023+6*frame, 2] = yfit84

pts[7*frame:1023+7*frame, 0] = yfit85[511]
pts[7*frame:1023+7*frame, 1] = xArr
pts[7*frame:1023+7*frame, 2] = yfit85

.compile /Users/colkin/Work04/Pluto_3um/IDL_Pro/waveCalCubic.pro 

wavecalcubic, pts, par, resid, chiSq, pmapcubic

writefits, '/Volumes/Vesta/CathyWork04/Pluto_3um/Reduction/PMap/pMapCubic.fits', pmapcubic
writefits, '/Volumes/Vesta/CathyWork04/Pluto_3um/Reduction/PMap/pMapParamsCubic.fits', par


d=readfits('/Volumes/Vesta/12aug01_CD1/spec/12aus0032.fits', h)

tvscl, rebin(rotate(d, 3), 512, 512)
tmp = getclix(pmap/2)