pro stb_summ_plot_spec, filename

get_timespan,tspan

load_stb_l1_mag,/polar,res=60.
load_stb_l1_let
load_stb_l1_sweadist
load_stb_l1_ste

;ylim,'stb_let_H_UNSEC_FLUX',1.,1e2
ylim,'stb_ste_D0',1.,1e2
ylim,'stb_B',-20.,20.
ylim,'stb_B_mag',0.,20.
zlim,'stb_let_H_UNSEC_FLUX',2.82e-3,1e4

;options,'stb_let_H_UNSEC_FLUX','spec',0
;options,'stb_let_H_UNSEC_FLUX','ytitle','LET H+ Flux'
;get_data,'stb_let_H_UNSEC_FLUX',ptr=plet
;options,'stb_let_H_UNSEC_FLUX','labels',string(*plet.v,format='(f4.1)')+'MeV'
;options,'stb_let_H_UNSEC_FLUX','ztitle','H+ Flux'
options,'stb_let_H_UNSEC_FLUX','zlog',1
options,'stb_let_H_UNSEC_FLUX','labflag',-1
;letbins = replicate(1b,12)
;letbins[11] = 0b
;options,'stb_let_H_UNSEC_FLUX','bins',letbins
options,'stb_let_H_UNSEC_FLUX','ystyle',1

;options,'stb_ste_D0','spec',0
;options,'stb_ste_D0','ytitle','STE e- Cnts'
;get_data,'stb_ste_D0',ptr=pste
;options,'stb_ste_D0','labels',reform(string((*pste.v)[0,*],format='(f5.1)')+'keV')
options,'stb_ste_D0','ztitle','e- Cnts'
options,'stb_ste_D0','zlog',1
options,'stb_ste_D0','labflag',-1
;options,'stb_ste_D0','labsize',1.
;stebins = bytarr(32)
;for i = 0,10 do stebins[3*i] = 1b
;options,'stb_ste_D0','bins',stebins
options,'stb_ste_D0','ystyle',1

options,'stb_B','ytitle','B (nT)'
options,'stb_B','labflag',-1

get_data,'stb_swea_dist.DISTRIBUTION',ptr=pdist
distsp = total(*pdist.y,3)

get_data,'stb_swea_dist.ENERGY',ptr=penergy
energy = reform((*penergy.y)[0,*])

store_data,'sweaspec',data={x: (*pdist.x), y: distsp, v:energy}
ylim,'sweaspec',1.,2e3,1
;options,'sweaspec','labels',string(energy,format='(f6.1)')+'eV'
options,'sweaspec','ytitle','SWEA'
options,'sweaspec','ztitle','e- Cnts'
;options,'sweaspec','labflag',-1
;options,'sweaspec','labsize',1.
;sweabins = bytarr(16)
;for i=0,7 do sweabins[i*2] = 1b
;options,'sweaspec','bins',sweabins
options,'sweaspec','spec',1
options,'sweaspec','zlog',1

tplot_options,'title','STEREO-B IMPACT Summary Plot'

set_plot,'z'
device,set_res=[780,780]
loadct2,34
!p.background = 255
!p.color = 0

tplot,['stb_let_H_UNSEC_FLUX','stb_ste_D0','sweaspec','stb_B','stb_B_mag']

tevent1 = time_double('2006-12-05/10:18')
tevent2 = time_double('2006-12-13/02:14')
tevent3 = time_double('2006-12-06/18:29')

timebar,tevent1,linestyle=2
timebar,tevent2,linestyle=2
timebar,tevent3,linestyle=2

mevent1 = time_double('2006-12-12/01:48')
mevent2 = time_double('2006-12-12/22:00')

timebar,mevent1,linestyle=1
timebar,mevent2,linestyle=1

makegif,filename
device,/close

return
end