;$Id: pcomp_Lu.pro,v 1.1 2012/02/20 05:42: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=8
  thick2=3
endif else begin
  thick1=4
  thick2=1
endelse
;
;  mv idl.ps ~/tex/tina/ObsSigna/pcomp_kf.eps
;  mv idl.ps ~/tex/tina/ObsSigna/pcomp_kf_hel.eps
;
;  mv idl.ps ~/tex/tina/Decay/fig/pcomp_kf_hel.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pcomp_kf2.ps
;
ihel=1 ;(helical)
ihel=0 ;(non-helical)
;
!p.charsize=1.7
!x.margin=[8.8,.5]
!y.margin=[3.2,.5]
;
;  sample to compare time series form different directories
;
;dir1='128c5_noforce3'
;dir2='128c5_noforce3_pm10'
;dir3='128c5_noforce3_pm100'
;
;dir1='128h2_noforce1_abisko'
;dir2='128h2_noforce1_pm10'
;dir3='128h2_noforce1_pm100'
;dir3b='256c5_pm100_noforce1'
;
dir1='short256pm1b_noforce'
dir2='short256pm10b_noforce'
dir3='short256pm100b_noforce'
dir4='short256pm1000b_noforce'
;
dir1='short256pm1a4_noforce'
dir1='short512pm1a3_noforce'
dir2='short256pm10a_noforce'
dir3='short256pm100a_noforce'
dir4='short256pm1000a_noforce'
;
dir1='short256pm1a4_noforce'
dir2='short512pm1a3_noforce'
dir3='short512pm100a3_noforce'
dir3='short512pm1b3_noforce'
;
dir1='short512pm1a3_noforce'
dir2='short512pm1b3_noforce'
;
if ihel eq 1 then begin
  dir1='short512pm1a3_noforce'
  dir2='short512pm100a3_noforce'
  yr=[1e3,1e4]
  print & print,'--> helical case'
  print,"mv idl.ps ~/tex/tina/ObsSigna/pcomp_kf_hel.eps" & print
endif else begin
  dir1='short512pm1b3_noforce'
  dir2='short512pm100b3_noforce'
  yr=[1e3,1e4]
  print & print,'--> non-helical case'
  print,"mv idl.ps ~/tex/tina/ObsSigna/pcomp_kf.eps" & print
endelse
;
!x.range=[1.,180.]
!x.title='!7g!6'
!y.title='!6'
eta=1e-5
siz2=2.0
siz=1.6
;
restore,'../'+dir1+'/pkf.sav'
vA=sqrt(2.*EM)
plot_oo,dt,vA/(eta*kI),yr=yr,/nodata
oplot,dt,vA/(eta*kI),col=122,thick=thick2
;
restore,'../'+dir2+'/pkf.sav'
vA=sqrt(2.*EM)
oplot,dt,vA/(eta*kI),col=55,thick=thick1
;
;restore,'../'+dir3+'/pkf.sav'
;oplot,dt,kI,col=188
;oplot,dt,kd,col=188,li=2
;
;restore,'../'+dir4+'/pkf.sav'
;oplot,dt,kI,col=55
;oplot,dt,kd,col=55,li=2
;
;restore,'../'+dir3+'/pkf.sav'
;oplot,dt,kI1,col=55
;
;restore,'../'+dir3b+'/pkf.sav'
;oplot,dt,kI1,col=55,li=2
;
;  put labels
;
if ihel eq 1 then begin
  ; helical case
  xx=[16.,140.] & oplot,xx,1.3e-2*xx^.667,li=1
  xx=[30.,140.] & oplot,xx,2.4e-3*xx^.4,li=1
  xyouts,siz=siz,80,1.5e-1,'!7g!6!u-2/3!n'
  xyouts,siz=siz,80,2e-2,'!7g!6!u-0.4!n'
  xyouts,siz=siz2,1.8,9e-2,'!7n!d!8M!6!n'
  xyouts,siz=siz2,1.8,6e-3,'!7k!d!8d!6!n'
  xyouts,siz=siz,8.,4e-2,'!6Pr!d!8M!6!n=100',col=55
  xyouts,siz=siz,20.,4e-1,'!6Pr!d!8M!6!n=1',col=122
  xyouts,siz=siz,70.,2e-3,'!6Pr!d!8M!6!n=1',col=122
endif else begin
  ; non-helical case
  xx=[30.,140.] & oplot,xx,8e-4*xx^.5,li=1
  xx=[30.,140.] & oplot,xx,6e-3*xx^.5,li=1
  xyouts,siz=siz,80,2e-2,'!7g!6!u-1/2!n'
  xyouts,siz=siz2,1.8,9e-2,'!7n!d!8M!6!n'
  xyouts,siz=siz2,1.8,5e-3,'!7k!d!8d!6!n'
  xyouts,siz=siz,50.,9e-2,'!6Pr!d!8M!6!n=100',col=55
  xyouts,siz=siz,30.,2e-3,'!6Pr!d!8M!6!n=1',col=122
  xyouts,siz=siz,20.,2e-1,'!6Pr!d!8M!6!n=1',col=122
endelse
;
END