;$Id: pexponents.pro,v 1.2 2012/02/04 22:50:01 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/pexponents.ps
;
!p.charsize=1.7
!x.margin=[6.8,.8]
!y.margin=[3.2,.8]
!x.title='!8R!6'
siz=2.0
;
file='exponents.dat'
a=rtable(file,head=1,4)
R=reform(a(0,*))
m=reform(a(1,*))
n=reform(a(2,*))
s=reform(a(3,*))
;
circ_sym,1.,1
plot_oo,R,m,xr=[.5,800],yr=[.05,3],ps=-2
;
circ_sym,1.3,1
oplot,R,n,ps=8
;
ds=1.4
circ_sym,1.8,0
oplot,R,s,li=1,ps=8
errplot,R,s/ds,s*ds
;
;  overplot lines
;
xx=grange(6.,800,20)
oplot,xx,8./xx,li=1
;
xx=grange(.5,100,20)
oplot,xx,2.3/xx^.25,li=2
;
xyouts, 2.,.38,'!8m!6',siz=siz
xyouts,20.,.20,'!8s!6',siz=siz
xyouts,20.,1.4,'!8n!6',siz=siz
;
xx=grange(.5,800,80)
nn=(2.3/xx^.25 > .667)
ss=(8.0/xx < nn)
mm=(nn-ss) >  .5
;oplot,xx,mm
;
END