;$Id: pkt_hel_difftimes.pro,v 1.10 2012/11/17 15:26:44 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=24,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_sig003_noforce2.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
pc_read_param,/param2,obj=param
nu=param.nu
;
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=[7.6,.8]
!y.margin=[3.3,.3]
!p.charsize=3.2
!p.multi=[0,1,5]
tilde='!9!s!aA!n!r!6'
@parameters
;
;  compute k1=1/k
;
k1=k
k1(0)=1.
k1=1./k1
k1(0)=1.
;
;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 3 then !x.title=xtit else !x.title='!6'
  if i eq 2 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=4
  oplot,k/k0,s*spec1(*,it),li=1,thick=2
  kM=total(spec2(*,it))/total(k1*spec2(*,it))
  urms=sqrt(total(spec1(*,it))*2.)
  ratio=sqrt(total(spec2(*,it))*2.)/sqrt(total(spec1(*,it))*2.)
  Rey=urms/(nu*kM)
  print,it,t(it)-t(0),t(it),urms,ratio,Rey
  xx1=[1,4] 
  xx2=[14,50]
  if i eq 1 then begin
    oplot,xx1/k0,1.0e-6*xx1^4,col=122
    xyouts,.10,1.4e-4,'!7g!6!u4!n',siz=siz,col=122
  endif
  if i eq 2 then begin
    oplot,xx2/k0,2e-3/xx2^1.6667,col=122
    xyouts,2.6,2.4e-5,'!7g!6!u-5/3!n',siz=siz,col=122
  endif
  xyouts,7.,4e-4,'!7g!6='+str(fix(.5+t(it)-t(0))),siz=siz
  wait,w
endfor
;
print,'mv idl.ps ~/tex/tina/QCD/fig/short512pm1a3_sig001b_noforce.eps'
END