;$Id: pkt_comp.pro,v 1.1 2020/06/06 09:09:54 brandenb Exp $
if !d.name eq 'PS' then begin
  default,ysize,12
  default,thick,4
  device,xsize=18,ysize=ysize,yoffset=3
  !p.charthick=thick & !p.thick=thick & !x.thick=thick & !y.thick=thick
end
;
dir1='N1_tt5em4_k180a'
dir1='N10_tt5em4_k180a' & s1=10.
dir2='N01_t5em4_k60a' & s2=.1
;
default,siz,1.7
!p.multi=0
half='!s!u 1!n!r!s-!r!d 2!n'
!x.title='!8k!6/!8k!6!d*!n'
!y.title='!8E!6!dM!n!6(!8k,t!6)  [!8v!6!s!u2!n!r!dA!n!8k!6!d*!n]'
!x.margin=[7.4,1.3]
!y.margin=[3.2,.4]
!p.charsize=siz
tilde='!9!s!aA!n!r!6'
;
yr1=[2e-14,5e-4]
yr2=[2e-14,5e-2]
;
k0=180.
vA0=.1
s=1./(vA0^2*k0) & default,yr1,[2e-9,2e-2] & default,xr,[.9,576]/k0
tauA=1./(k0*vA0)
print,'tauA=',tauA
;
default,xstyle,1
if xstyle then begin
  plot_oo,xr,yr1,/nodata,xtickf='logticks_exp',ytickf='logticks_exp'
endif else begin
  plot_oo,xr,yr1,/nodata
endelse
;
;  plot 1
;
restore,'../'+dir1+'/spec.sav'
nt=n_elements(t)
for it=0,nt-1 do oplot,k/k0,s*s1*spec1(*,it),thick=thick2,col=55
;
;  plot 2
;
s=s*10
restore,'../'+dir2+'/spec.sav'
nt=n_elements(t)
for it=0,nt-1 do oplot,k/k0,s*s2*spec1(*,it),thick=thick2,col=122
;
oplot,k/k0,.01/k^(11./3.)
siz=1.7
default,xlab,5.
default,ylab,1e-6
;if label ne '' then xyouts,xlab,ylab,label,siz=siz
;
default,texdir,'heldecay'
print,'$mv idl.ps ~/tex/tina/'+texdir+'/fig/pkt_XX.ps'
!p.multi=0
END