;$Id: reb_spec.pro,v 1.2 2018/05/08 10:40:22 brandenb Exp $ if !d.name eq 'PS' then begin device,xsize=18,ysize=9,yoffset=3 !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3 end ; ; mv idl.ps ~/tex/tina/prop/dustEB/fig/peb.ps ; loadct,5 !p.position=0 !p.charsize=1.2 !p.multi=[0,1,2] !x.margin=[2.0,0.2] !y.margin=[1.8,0.2] ; ; background ; !p.position=0 !p.multi=[0,2,1] restore,'eb.sav' s=size(eb) nx=s[1] ny=s[2] x=2*!pi*findgen(nx)/nx y=2*!pi*findgen(ny)/ny ; ; compute E and B vectors ; E=float(eb) B=imaginary(eb) ; ; plot in full size ; power2d_hel0,E,E,wav=k,spe=EEspec,/noplo,/quiet power2d_hel0,B,B,wav=k,spe=BBspec,/noplo,/quiet power2d_hel0,E,B,wav=k,spe=EBspec,/noplo,/quiet ; save,file='EBspec.sav',k,EEspec,BBspec,EBspec ; END