;


for i=800,1000 do begin
	
	;correct the DC point and spectral images
	spc=spectra[*,*,i]
	for h=0,134 do begin
		spc[128,h]=(spc[127,h]+spc[129,h])/2.
		spc[*,h]=specImage(spc[*,h],10,imagedown=22)
		noi=specnoise(spc[*,h],10)
		pk=specpeak(spc[*,h],noi)
		nm=n_elements(pk)
		if nm gt 0 then begin
			cpsp=spc[*,h]*0+10

	mx=max(spc[pk,h],mxind)
	offs=pk[mxind]-128
	cpsp[pk-offs]=spc[pk,h]
;			cpsp[128:128+nm-1]=spc[pk,h]
			spc[*,h]=cpsp
		endif		
	endfor	

	;contour,spectra[*,*,i],velbins,height,levels=[0.5e3,1e3,1.5e3,1e4,1.5e4,1e5,1.5e5,1e6,1.5e6,1e7,1.5e7],/fill,yrange=[0,1.5],/ystyle
	contour,alog(spc),velbins,height,/fill,xrange=[-2,2],yrange=[0,1.5],$
		levels=[8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25]
	   title='i = '+strtrim(i,2)+', Time = '+strtrim(time[i],2)
	wait,0.05
endfor

end