;$Id: ptab_ampl.pro,v 1.4 2026/07/29 16:53:57 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,4)
ampl=reform(a[0,*])
divu2m=reform(a[1,*])
orms2=reform(a[2,*])
epsK=reform(a[3,*])
xr=[.01,.6]
yr=[2e-3,3.]
circ_sym,1.3,1
plot_oo,ampl,sqrt(orms2),ps=8,xr=xr,yr=yr,/nodata
oplot,ampl,sqrt(orms2),ps=8,col=122
oplot,ampl,sqrt(divu2m),ps=8,col=55
oplot,ampl,1.3333*divu2m*.02*1.15,ps=8,col=155
xx=[.04,.3]
oplot,xx,7.6*xx,col=55
oplot,xx,2.1*xx^2,col=122
;
siz=1.8
xyouts,.02,.12,siz=siz,'!6(!9G.!8u!6)!drms!n!9A !8u!6!dini!n',col=55
xyouts,.03,.02,siz=siz,'!7x!6!drms!n!9A !8u!6!s!dini!n!r!u2!n',col=122
;
loadct,6
oplot,ampl,epsK,ps=8,col=144
oplot,xx,1.72*xx^2,col=144
xyouts,.16,.02,siz=siz,'!7e!6!dK!n!9A !8u!6!s!dini!n!r!u2!n',col=144
loadct,5
;
print,"$mv idl.ps ~/Overleaf/Eva/BE_collapse/fig/ptab_ampl.eps"
END