;$Id: ptab_ampl.pro,v 1.2 2026/06/29 18:12:17 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=10,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
@parameters
!p.charsize=1.7
!x.margin=[5.3,.5]
!y.margin=[3.2,.5]
!x.range=[0,3.3]
!x.title='!8u!6!dini!n'
!y.title='!6'
;
file='tab_ampl.txt'
a=rtable(file,head=1,3)
ampl=reform(a[0,*])
divu2m=reform(a[1,*])
orms2=reform(a[2,*])
xr=[.01,.6]
yr=[2e-3,3.]
circ_sym,1.3,1
plot_oo,ampl,sqrt(orms2),ps=8,xr=xr,yr=yr
oplot,ampl,sqrt(divu2m),ps=8
xx=[.04,.3]
oplot,xx,7.6*xx
oplot,xx,2.1*xx^2
;
siz=1.8
xyouts,.02,.12,siz=siz,'!6(!9G.!8u!6)!drms!n!9A !8u!6!dini!n'
xyouts,.16,.02,siz=siz,'!7x!6!drms!n!9A !8u!6!s!dini!n!r!u2!n'
;
print,"$mv idl.ps ~/Overleaf/Eva/BE_collapse/fig/ptab_ampl.eps"
END