;$Id: pcomp_kft_QCD2.pro,v 1.1 2013/07/20 03:32:43 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
  thick1=2
  thick2=5
  thick3=2
  thick4=5
  thickt=6
endif else begin
  thick1=2
  thick2=5
  thick3=2
  thick4=5
  thickt=6
endelse
;
ihel=0 ;(non-helical)
ihel=1 ;(helical)
;
!p.charsize=1.7
!x.margin=[8.8,.9]
!y.margin=[3.2,.9]
;
;  sample to compare time series form different directories
;
k0=15.
;
if ihel eq 1 then begin
  dir1='short512pm1a3_noforce'
  dir2='short512pm1a3_sig01_noforce'
  dir3='short512pm1a3_sig01b_noforce'
  dir2='short512pm1a3_sig01b_noforce'
  dir4='short512pm1a3_sig01_k30_noforce'
  dir4='short512pm1a3_sig001b_noforce'
  ;dirt='short512pm1a3_noforce_tdep'
  ;
  yr=[.04,11.]
  print & print,'--> helical case'
  print,"mv idl.ps ~/tex/tina/QCD/fig/pcomp_kft_QCD.eps" & print
endif else begin
  stop,"nok"
  dir1='short512pm1b3_noforce'
  dir2='short512pm1b3_noforcek1'
  dir3='short512pm1b3_noforcek01'
  dir4='short512pm100b3_noforce'
  dirt='short512pm1b3_noforce_tdep'
  ;
  yr=[.6,10.]
  print & print,'--> non-helical case'
  print,"mv idl.ps ~/tex/tina/QCD/pcomp_kft_nohel_all.eps" & print
endelse
;
!x.range=[1.,200.]
!x.title='!7g!6'
!y.title='!7n!d!8M!6!n/!7k!d!80!6!n'
siz2=2.0
sizs=1.3
siz=1.6
;
restore,'../'+dir1+'/pkf.sav'
plot_oo,dt,k0/kI,yr=yr,/nodata
oplot,dt,k0/kI,thick=thick1,li=2
kmax=2.*EM/HM
oplot,dt,k0/kmax,thick=thick1,li=0
;
restore,'../'+dir2+'/pkf.sav'
oplot,dt,k0/kI,col=55,thick=thick2,li=2
kmax=2.*EM/HM
oplot,dt,k0/kmax,col=55,thick=thick1,li=0
;
;restore,'../'+dir3+'/pkf.sav'
;oplot,dt,k0/kI,col=155,thick=thick3,li=2
;kmax=2.*EM/HM
;oplot,dt,k0/kmax,col=155,thick=thick3,li=0
;
restore,'../'+dir4+'/pkf.sav'
oplot,dt,k0/kI,col=122,thick=thick4,li=2
kmax=2.*EM/HM
oplot,dt,k0/kmax,col=122,thick=thick4,li=0
;
;restore,'../'+dirt+'/pkf.sav'
;oplot,dt,k0/kI,col=188,thick=thickt,li=3
;
;  put labels
;
if ihel eq 1 then begin
  ; helical case
  xx=[16.,140.] & oplot,xx,k0*2.0e-2*xx^.667,li=1
  xx=[20.,80.] & oplot,xx,k0*6.4e-2*xx^.500,li=1
  xx=[2.,10.] & oplot,xx,k0*1.0e-1*xx^.250,li=1
  ;xx=[20.,160.] & oplot,xx,0.7*xx^.400,li=1
  ;xyouts,siz=siz,28.,6.3,'!6Pr!d!8M!6!n=1'
  ;xyouts,siz=siz,40.,2.8,'!7a!6=1',col=55
  ;xyouts,siz=siz,50.,k0*9e-2,'!7a!6=10',col=155
  ;xyouts,siz=sizs,85.,4.95,'!6Pr!d!8M!6!n=100',col=122
  ;xyouts,siz=sizs,75.,3.4,'!6Pr!d!8M!6!n=!7g!6!u3/2!n',col=188
endif else begin
  ; non-helical case
  xx=[30.,140.] & oplot,xx,k0*2.4e-2*xx^.520,li=1
  xx=[30.,140.] & oplot,xx,k0*3.4e-2*xx^.333,li=1
  xx=[30.,140.] & oplot,xx,k0*2.4e-2*xx^.333,li=1
  xx=[20.,140.] & oplot,xx,.7*xx^.250,li=1
  xyouts,siz=siz,28.,3.0,'!6Pr!d!8M!6!n=1'
  xyouts,siz=siz,110.,2.3,'!7a!6=1',col=55
  xyouts,siz=siz,100.,1.5,'!7a!6=10',col=155
  xyouts,siz=sizs,75.,2.8,'!6Pr!d!8M!6!n=100',col=122
  xyouts,siz=sizs,75.,1.9,'!6Pr!d!8M!6!n=!7g!6!u3/2!n',col=188
endelse
;
END