;$Id: pcomp_EM_EK.pro,v 1.7 2012/02/18 14:09:17 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
end
;
;  mv idl.ps ~/tex/tina/Decay/fig/pcomp_EM_EK_hel.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pcomp_EM_EK.ps
;
!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='short256pm1a4_noforce'
dir2='short256pm10a_noforce'
dir3='short256pm100a_noforce'
dir4='short256pm1000a_noforce'
;
dir1='short256pm1b_noforce'
dir2='short256pm10b_noforce'
dir3='short256pm100b_noforce'
dir4='short256pm1000b_noforce'
;
dir1='short256pm1a4_noforce'
;
dir1='short512pm1a3_noforce'
dir2='short512pm1b3_noforce'
dir3='short512pm100a3_noforce'
dir4='short512pm100b3_noforce'
;
!x.range=[1.,600]
!x.title='!8t!6'
!y.title='!8k!6!df!n/!8k!6!d1!n'
;
restore,'../'+dir1+'/pkf.sav'
plot_oo,dt,EM,yr=[1e-4,2e-1]
oplot,dt,EK,li=2
;
restore,'../'+dir2+'/pkf.sav'
oplot,dt,EM,col=122
oplot,dt,EK,col=122,li=2
;
restore,'../'+dir3+'/pkf.sav'
oplot,dt,EM,col=55
oplot,dt,EK,col=55,li=2
;
xx=[10.,100.] & oplot,xx,5e-1/xx^.667,li=3,col=188
xx=[2.,200.] & oplot,xx,8e-2/xx^1,li=3
xx=[2.,200.] & oplot,xx,5e-2/xx^1.2,li=3
xyouts,10,3,'!8t!6!u-1/2!n'
;
END