meanall =fltarr(length)

for i=0, length-1 do meanall[i] = mean([px35[i], px36[i], px37[i], px38[i],$
px39[i], px40[i], px41[i], px42[i], px43[i], px44[i], $
px45[i], px46[i], px47[i], px48[i], px49[i], px50[i],$
px51[i], px52[i], px53[i], px54[i], px55[i], px56[i], $
px57[i], px58[i], px59[i], px60[i], px61[i], px62[i], $
px63[i], px64[i], px65[i], px70[i], px71[i], px72[i], $
px73[i], px74[i], px75[i]]) 


; madscale all of the individual spectra to the average one.

allSpectra = reform([px35, px36, px37, px38,$
px39, px40, px41, px42, px43, px44, px45, px46, px47, $
px48, px49, px50, px51, px52, px53, px54, px55, px56, $
px57, px58, px59, px60, px61, px62, px63, px64, px65, $
px70, px71, px72, px73, px74, px75], [length, 37])


sz = size(allSpectra, /Dim)

a = fltarr(sz[1])
for i=0, sz[1]-1 do a[i]=madscale(meanall, allSpectra[*,i])

print, mean(a)
;     0.988155
 
sigmaall = fltarr(length)     
for i=0, length-1 do sigmaAll[i] = stddev([1/a[0]*px35[i], 1/a[1]*px36[i], 1/a[2]*px37[i], 1/a[3]*px38[i],$
1/a[4]*px39[i], 1/a[5]*px40[i], 1/a[6]*px41[i], 1/a[7]*px42[i], 1/a[8]*px43[i], 1/a[9]*px44[i], $
1/a[10]*px45[i], 1/a[11]*px46[i], 1/a[12]*px47[i], 1/a[13]*px48[i], 1/a[14]*px49[i], 1/a[15]*px50[i],$
1/a[16]*px51[i], 1/a[17]*px52[i], 1/a[18]*px53[i], 1/a[19]*px54[i], 1/a[20]*px55[i], 1/a[21]*px56[i], $
1/a[22]*px57[i], 1/a[23]*px58[i], 1/a[24]*px59[i], 1/a[25]*px60[i], 1/a[26]*px61[i], 1/a[27]*px62[i], $
1/a[28]*px63[i], 1/a[29]*px64[i], 1/a[30]*px65[i], 1/a[31]*px70[i], 1/a[32]*px71[i], 1/a[33]*px72[i], $
1/a[34]*px73[i], 1/a[35]*px74[i], 1/a[36]*px75[i]], /double)/sqrt(37.)

plot, meanall
errplot, meanall-sigmaall, meanall+sigmaall, color=5

dir = '/Users/colkin/Work04/Pluto_3um/PlutoSpectra3/'

writefits, dir+'meanP.fits', meanall


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

bin, w, meanall, sigmaall, 10., newwall, newmeanall, newsigmaall
for i=0, n_elements(newwall)-1 do print, newwall[i], newmeanall[i], newsigmaall[i]


























;;
;stopped here


sigma198 =fltarr(959)
 
for i=0, 958 do sigma198[i] = stddev([px37[i], px38[i],$
px41[i], px42[i], px45[i], px46[i], $
px49[i], px50[i], px53[i], px54[i], $
px57[i], px58[i], px61[i], px62[i], px65[i], $
px70[i], px73[i], px74[i]], /double) 

plot, wall, pxall, yrange=[-1000., 60000.]
errplot, wall, pxall-sigmaAll/sqrt(37.), pxall+sigmaAll/sqrt(37.), color=4


plot, w60, pave60, yrange=[-1000., 60000.]
errplot, w60, pave60-sigma60/sqrt(19.), pave60+sigma60/sqrt(19.), color=4
oplot, w198, pave198, color=5



plot, w60, px35, yrange=[-10000,90000.], xrange=[3.4, 3.6]
errplot, w60, px35-err35, px35+err35

errall = (1/(19.+18.))*sqrt($
err37^2 + err38^2 + err41^2 + err42^2 + err45^2 + $
err46^2 + err49^2 + err50^2 + err53^2 + err54^2 + err57^2 + err58^2 + err61^2 + $
err62^2 + err65^2 + err70^2 + err73^2 + err74^2 + $
err35^2 + err36^2 + err39^2 + err40^2 + err43^2 + err44^2 + $
err47^2 + err48^2 + err51^2 + err52^2 + err55^2 + err56^2 + $
err59^2 + err60^2 + err63^2 + err64^2 + err71^2 + err72^2 + err75^2)


plot, w60, px2, yrange=[-1000,60000.]
errplot, w60, px2-errMean60, px2+errMean60
plot, w60, px35, yrange=[-10000,90000.], xrange=[3.4, 3.5]
errplot, w60, px35-err35, px35+err35, color=5  

wall= (w60+w198)/2

plot, wall, pxall, yrange=[-1000,60000.], xrange=[3.4, 3.5]
errplot, wall, pxall-errAll, pxall+errAll, color=4

scaleFactor = 5.

wnew = fltarr((n_elements(wall)-4)/scaleFactor )

for i=0, ((n_elements(wall)-4)/scaleFactor)-1 do begin &$
	wnew[i]=TOTAL(wall[scaleFactor*i:scaleFactor*i+(scaleFactor-1)])/scaleFactor &$
end

pspec = fltarr((n_elements(wall)-4)/scaleFactor )

for i=0, ((n_elements(wall)-4)/scaleFactor)-1 do begin &$
	pspec[i]=TOTAL(pxall[scaleFactor*i:scaleFactor*i+(scaleFactor-1)])/scaleFactor &$
end

var = errall^2

newvar = fltarr((n_elements(wall)-4)/scaleFactor)

for i=0, ((n_elements(wall)-4)/scaleFactor)-1 do begin &$
	newvar[i]=TOTAL(var[scaleFactor*i:scaleFactor*i+(scaleFactor-1)])/scaleFactor &$
end

newErr = sqrt(newvar)

plot, wall, pxall, yrange=[-1000,60000.], xrange=[3.4, 3.5]
errplot, wall, pxall-errAll, pxall+errAll, color=4

plot, wnew, pspec, yrange=[-1000,60000.], xrange=[3.4, 3.5]
errplot, wnew, pspec-newErr, pspec+newErr, color=4

plot, wall, pxall, yrange=[-1000,60000.]
errplot, wall, pxall-errAll, pxall+errAll, color=4

plot, wnew, pspec, yrange=[-1000,60000.]
errplot, wnew, pspec-newErr, pspec+newErr, color=4

for i=0, 958 do print, wall[i], pxall[i], errall[i]

; binning by 25
scaleFactor = 25.

wnew = fltarr((n_elements(wall)-9)/scaleFactor )

for i=0, ((n_elements(wall)-9)/scaleFactor)-1 do begin &$
	wnew[i]=TOTAL(wall[scaleFactor*i:scaleFactor*i+(scaleFactor-1)])/scaleFactor &$
end

pspec = fltarr((n_elements(wall)-9)/scaleFactor )

for i=0, ((n_elements(wall)-9)/scaleFactor)-1 do begin &$
	pspec[i]=TOTAL(pxall[scaleFactor*i:scaleFactor*i+(scaleFactor-1)])/scaleFactor &$
end

var = errall^2

newvar = fltarr((n_elements(wall)-9)/scaleFactor)

for i=0, ((n_elements(wall)-9)/scaleFactor)-1 do begin &$
	newvar[i]=TOTAL(var[scaleFactor*i:scaleFactor*i+(scaleFactor-1)])/scaleFactor &$
end

newErr = sqrt(newvar)

plot, wnew, pspec, yrange=[-1000,60000.], background=255, color=0   
errplot, wnew, pspec-newErr, pspec+newErr, color=0                  
oplot, rebin(w60[0:949], 950/25), rebin(px1[0:949], 950/25), color=3
oplot, rebin(w198[0:949], 950/25), rebin(px2[0:949], 950/25), color=5

; binning by 10
scaleFactor = 10.

wnew = fltarr((n_elements(wall)-9)/scaleFactor )

for i=0, ((n_elements(wall)-9)/scaleFactor)-1 do begin &$
	wnew[i]=TOTAL(wall[scaleFactor*i:scaleFactor*i+(scaleFactor-1)])/scaleFactor &$
end

pspec = fltarr((n_elements(wall)-9)/scaleFactor )

for i=0, ((n_elements(wall)-9)/scaleFactor)-1 do begin &$
	pspec[i]=TOTAL(pxall[scaleFactor*i:scaleFactor*i+(scaleFactor-1)])/scaleFactor &$
end

var = errall^2

newvar = fltarr((n_elements(wall)-9)/scaleFactor)

for i=0, ((n_elements(wall)-9)/scaleFactor)-1 do begin &$
	newvar[i]=TOTAL(var[scaleFactor*i:scaleFactor*i+(scaleFactor-1)])/scaleFactor &$
end

newErr = sqrt(newvar)

plot, wnew, pspec, yrange=[-1000,60000.], background=255, color=0   
errplot, wnew, pspec-newErr, pspec+newErr, color=0                  
oplot, rebin(w60[0:949], 950/scaleFactor), rebin(px1[0:949], 950/scaleFactor), color=3
oplot, rebin(w198[0:949], 950/scaleFactor), rebin(px2[0:949], 950/scaleFactor), color=5