;$Id: pkt_hel_difftimes2.pro,v 1.1 2013/07/20 03:32:43 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=18,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ~/tex/tina/QCD/fig/short512pm1a3_sig01b_noforce.eps
;  mv idl.ps ~/tex/tina/QCD/fig/short512pm1a3_sig001b_noforce.eps
;
default,iread,0
if iread eq 0 then begin
  power,'_kin','_mag',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot & print,n
  power,'hel_kin','hel_mag',k=k,spec1=spec1h,spec2=spec2h,i=n,tt=t,/noplot & print,n
  iread=1
endif
;
siz=1.5
!p.multi=0
xtit='!8k!6/!8k!6!d0!n'
ytit2='!8E!6!dM!n(!8k!6,!7g!6),  !8E!6!dK!n(!8k!6,!7g!6),  and  !8kH!6!dM!n(!8k!6,!7g!6)/2'
!x.margin=[9.4,.8]
!y.margin=[3.3,.3]
!p.charsize=2.7
!p.multi=[0,1,3]
tilde='!9!s!aA!n!r!6'
@parameters
;
;w=0.
default,w,.1
s=1./k^2.
s=k^2.5
s=1. & default,yr,[2e-9,2e-2]
k0=15.
E=t
EK=t
D=t
II=t
IIK=t
default,it3,80
nit=n_elements(iit)
for i=0,nit-1 do begin
  if i eq 2 then !x.title=xtit else !x.title='!6'
  if i eq 1 then !y.title=ytit2 else !y.title='!6'
  it=iit(i)
  plot_oo,k/k0,s*spec2(*,it),xr=[1,256]/k0,yr=yr
  oplot,k/k0,s*spec2h(*,it)*k/2.,li=2,thick=.5
  oplot,k/k0,s*spec1(*,it),li=1,thick=.5
  print,it,t(it)-t(0)
  if i eq 1 then begin
    xx=[1,4] & oplot,xx/k0,1.0e-6*xx^4,col=122
    xx=[20,80] & oplot,xx/k0,1.6e-2/xx^1.6667,col=122
  endif
  xyouts,8.,2e-3,'!7g!6='+str(fix(t(it)-t(0))),siz=siz
  wait,w
endfor
;
print,'mv idl.ps ~/tex/tina/QCD/fig/short512pm1a3_sig001b_noforce.eps'
END