;$Id: pcomp_kfk_free.pro,v 1.5 2012/05/13 22:32:20 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=8
  thick2=3
  thick3=3
endif else begin
  thick1=4
  thick2=1
  thick3=3
endelse
;
ihel=0 ;(non-helical)
ihel=1 ;(helical)
;
!p.charsize=1.7
!x.margin=[8.8,.5]
!y.margin=[3.2,.5]
;
;  sample to compare time series form different directories
;
k0=15.
;
if ihel eq 1 then begin
  dir1='short512pm1a3_noforce'
  dir2='short512pm1a3_noforcek1'
  dir3='short512pm1a3_noforcek01'
  ;
  yr=k0*[8e-4,8e-1]
  print & print,'--> helical case'
  print,"mv idl.ps ~/tex/tina/ObsSigna/pcomp_kft_hel_free.eps" & print
endif else begin
  dir1='short512pm1b3_noforce'
  dir2='short512pm1b3_noforcek1'
  dir3='short512pm1b3_noforcek01'
  ;
  yr=k0*[1e-3,7e-1]
  print & print,'--> non-helical case'
  print,"mv idl.ps ~/tex/tina/ObsSigna/pcomp_kft_nohel_free.eps" & print
endelse
;
!x.range=[1.,240.]
!x.title='!7g!6'
!y.title='!6'
siz2=2.0
siz=1.6
;
restore,'../'+dir1+'/pkf.sav'
plot_oo,dt,k0/kI,yr=yr,/nodata
oplot,dt,k0/kI,col=122,thick=thick2
oplot,dt,k0/kM,col=122,thick=thick2,li=2
;
restore,'../'+dir2+'/pkf.sav'
oplot,dt,k0/kI,col=55,thick=thick1
oplot,dt,k0/kM,col=55,thick=thick1,li=2
;
restore,'../'+dir3+'/pkf.sav'
oplot,dt,k0/kI,col=188,thick=thick3
oplot,dt,k0/kM,col=188,thick=thick3,li=2
;
;  put labels
;
if ihel eq 1 then begin
  ; helical case
  xx=[16.,140.] & oplot,xx,k0*3.3e-2*xx^.667,li=1,col=122
  xx=[10.,140.] & oplot,xx,k0*2.8e-2*xx^.500,li=1,col=55,thick=5
  xx=[30.,60.] & oplot,xx,k0*0.8e-2*xx^.667,li=1,col=188,thick=5
  ;xx=[10.,140.] & oplot,xx,k0*5.0e-2*xx^.4,li=1
  ;xx=[12.,140.] & oplot,xx,k0*5.0e-4*xx^1.4,li=1,col=55
  ;xx=[8.,140.] & oplot,xx,k0*1.1e-2*xx^.333,li=1,col=188
  ;xx=[15.,140.] & oplot,xx,k0*8.0e-3*xx^.4,li=1,col=188
  ;xx=[25.,140.] & oplot,xx,k0*8.0e-4*xx^.4,li=1,col=122
  ;xx=[8.,140.] & oplot,xx,k0*5.0e-2*xx^.3333,li=1
  ;xyouts,siz=siz,80,k0*1.5e-1,'!7g!6!u2/3!n'
  ;xyouts,siz=siz,80,k0*2e-2,'!7g!6!u0.4!n'
  xyouts,siz=siz2,1.5,k0*1e-1,'!7n!d!8M!6!n/!7k!d!80!6!n'
  ;xyouts,siz=siz2,1.8,k0*4e-3,'!7k!d!8d!6!n/!7k!d!80!6!n'
  xyouts,siz=siz2,1.8,k0*4e-3,'!7k!d!8M!6!n/!7k!d!80!6!n'
  xyouts,siz=siz,50.,k0*2.0e-1,'!7a!6=1',col=55
  xyouts,siz=siz,50.,k0*9e-2,'!7a!6=10',col=188
  xyouts,siz=siz,10.,k0*2.5e-1,'!7a!6=0',col=122
endif else begin
  ; non-helical case
  ;xx=[10.,100.] & oplot,xx,k0*18e-4*xx^1.2,li=1
  xx=[30.,140.] & oplot,xx,k0*3.4e-2*xx^.333,li=1,col=55,thick=6
  xx=[30.,140.] & oplot,xx,k0*2.4e-2*xx^.333,li=1,col=188,thick=6
  xx=[30.,140.] & oplot,xx,k0*8e-4*xx^.5,li=1
  xx=[30.,140.] & oplot,xx,k0*6e-3*xx^.5,li=1
  xyouts,siz=siz,80,k0*2e-2,'!7g!6!u1/2!n'
  xyouts,siz=siz2,1.5,k0*1e-1,'!7n!d!8M!6!n/!7k!d!80!6!n'
  ;xyouts,siz=siz2,1.8,k0*5e-3,'!7k!d!8d!6!n/!7k!d!80!6!n'
  xyouts,siz=siz2,1.8,k0*5e-3,'!7k!d!8M!6!n/!7k!d!80!6!n'
  xyouts,siz=siz,110.,k0*1.9e-1,'!7a!6=1',col=55
  xyouts,siz=siz,30.,k0*2e-3,'!6Pr!d!8M!6!n=1',col=122
  xyouts,siz=siz,20.,k0*2e-1,'!6Pr!d!8M!6!n=1',col=122
endelse
;
END