;$Id: ptab_k0.pro,v 1.11 2026/08/03 09:53:35 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=11,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
@parameters
!p.charsize=1.7
!x.margin=[7.3,.5]
!y.margin=[3.2,.5]
!x.title='!8k!6!d0!n/!8k!6!d1!n'
!x.title='!8k!6!d0!n'
!y.title='!6'
k1=.5
;
;-----------------------------------------------------------------------------
file='tab_k0.txt'
a=rtable(file,head=1,4)
k0=reform(a[0,*])  ;(already in units of k1, need to multiply by k1 to get actual k0)
divu2m=reform(a[1,*])
orms2=reform(a[2,*])
epsK=reform(a[3,*])
xr=[1.,20.] & yr=[4e-3,3.]
xr=[.7,34.] & yr=[4e-5,9.]
circ_sym,1.0,1
plot_oo,k0*k1,sqrt(orms2),ps=8,xr=xr,yr=yr,/nodata
oplot,k0*k1,sqrt(orms2),ps=8,col=122
oplot,k0*k1,sqrt(divu2m),ps=8,col=55
xx=[1.8,14.]
xx=[.9,14.]
oplot,xx,2.15/xx^.5,col=55
oplot,xx,.38/xx,col=122
;
loadct,6
oplot,k0*k1,epsK,ps=8,col=144
xx2=[1.8,7.]
oplot,xx2,.085/xx2,col=144
loadct,5
;
xxx=.9 & dxxx=.4 & yyy=2e-3 & dyyy=2.7 & si2=1.4
legend,/log,siz=si2,xxx,dxxx,yyy,0,ps=8,'!8 G!6!dN!n!9=!60, Runs N, O, P'
;
;-----------------------------------------------------------------------------
file='tab_k0_d.txt'
a=rtable(file,head=1,4)
k0=reform(a[0,*])  ;(already in units of k1, need to multiply by k1 to get actual k0)
divu2m=reform(a[1,*])
orms2=reform(a[2,*])
epsK=reform(a[3,*])
xr=[.01,.6]
xr=[1.,20.]
yr=[1e-2,3.]
circ_sym,1.0,0
oplot,k0*k1,sqrt(orms2),ps=8,col=122
oplot,k0*k1,sqrt(divu2m),ps=8,col=55
oplot,xx,1.4/xx^.5,li=2,col=55
oplot,xx,.19/xx,li=2,col=122
yyy=yyy/dyyy & legend,/log,siz=si2,xxx,dxxx,yyy,0,ps=8,'!8 G!6!dN!n!9=!60, Runs F, L, M'
;
;-----------------------------------------------------------------------------
file='tab_k0_G0.txt'
a=rtable(file,head=1,4)
k0=reform(a[0,*])  ;(already in units of k1, need to multiply by k1 to get actual k0)
divu2m=reform(a[1,*])
orms2=reform(a[2,*])
epsK=reform(a[3,*])
xr=[.01,.6]
xr=[1.,20.]
yr=[1e-2,3.]
ps=4
symsiz=2
oplot,k0*k1,sqrt(orms2),ps=ps,symsiz=symsiz,col=122
oplot,k0*k1,sqrt(divu2m),ps=ps,symsiz=symsiz,col=55
;oplot,xx,1.4/xx^.5,li=2,col=55
;oplot,xx,.19/xx,li=2,col=122
yyy=yyy/dyyy & legend,/log,siz=si2,xxx,dxxx,yyy,0,ps=ps,'!8 G!6!dN!n!6=0, !8k!6!dcut!n=25'
;
;-----------------------------------------------------------------------------
file='tab_k0_G0_kcut100.txt'
a=rtable(file,head=1,4)
k0=reform(a[0,*])  ;(already in units of k1, need to multiply by k1 to get actual k0)
divu2m=reform(a[1,*])
orms2=reform(a[2,*])
epsK=reform(a[3,*])
xr=[.01,.6]
xr=[1.,20.]
yr=[1e-2,3.]
ps=1
symsiz=3
oplot,k0*k1,sqrt(orms2),ps=ps,symsiz=symsiz,col=122
oplot,k0*k1,sqrt(divu2m),ps=ps,symsiz=symsiz,col=55
;oplot,xx,1.4/xx^.5,li=2,col=55
;oplot,xx,.19/xx,li=2,col=122
yyy=yyy/dyyy & legend,/log,siz=si2,xxx,dxxx,yyy,0,ps=ps,'!8 G!6!dN!n!6=0, !8k!6!dcut!n=50'
;
;=============================================================================
;kcut=.7*.5*50.
;xxx=grange(xx[0],xx[1],50)
;oplot,xxx,1e-2*(kcut^1.333-xxx^1.333),li=2,col=55,th=4
;print,1e-2*(kcut^1.333-xxx^1.333)
;
;=============================================================================
;
siz=1.6
xyouts,5.52,1.30,siz=siz,'!6(!9G.!8u!6)!drms!n!9A !8k!6!s!d0!n!r!u-1/2!n',col=55
xyouts,5.46,.09,siz=siz,'!7x!6!drms!n!9A !8k!6!s!d0!n!r!u-1!n',col=122
;
loadct,6
xyouts,1.70,.01,siz=siz,'!7e!6!dK!n!9A !8k!6!s!d0!n!r!u-1!n',col=144
oplot,k0*k1,epsK,ps=8,col=144
oplot,xx2,.073/xx2,li=2,col=144
loadct,5
;
print,"$mv idl.ps ~/Overleaf/Eva/BE_collapse/fig/ptab_k0.eps"
END