;$Id: pkpm.pro,v 1.1 2026/05/10 08:24:56 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
;
!p.charsize=1.7
!x.margin=[7.8,.5]
!y.margin=[3.2,2.5]
;
xr=[1.,100]
xr=[.6,200]
yr=[1e-5,1e-1]
!x.title='!7n!6!dM!n(!8t!6)'
!y.title='!8v!6!dA!n(!8t!6)'
;
file='/kpm.txt'
;
;-----------------------------------------------------------------------------
dir1='1024_k1em2_k1em0_vA1em1_nu1em4'  ;black, t=2.45E+05
dir2='1024_k1em2_k1em1_vA1em1_nu1em4'  ;red, t=4.64e+04 (did reach it)
dir3='1024_k1em1_k1e1_vA3em3_nu1em5'  ;black, t=2.45E+05
dir3='1024_k1em3_k1em1_vA3em3_nu1em3'  ;black, t=2.45E+05
dir3='1024_k1em3_k1em1_vA7em3_nu1em3'  ;black, t=2.45E+05
dir32='1024_k1em3_k1em1_vA7em3_nu1em3b'  ;black, t=2.45E+05
dir33='1024_k1em3_k1em1_vA7em3_nu1em4b'  ;orange, running
dir34='1024_k1em3_k1em1_vA7em3_nu1em4'  ;should be blue, running
dir4='1024_k1em2_k1em0_vA1em1_nu1em4b'  ;orange, t=6.81e+03
dir42='1024_k1em3_k1em1_vA1em1_nu1em3b'  ;green, running, CB=2.9-2.95
;pc_diffrun 1024_k1em2_k1em1_vA1em1_nu1em4 1024_k1em3_k1em1_vA1em1_nu1em3b
;
;  i   t      EMag       xiMag        CM          Lu
;  & CM=a[4,*] & Lu=a[5,*]
;
a=rtable('../'+dir1+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
plot_oo,xiM,vA,ps=1,xr=xr,yr=yr
;
a=rtable('../'+dir2+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,xiM,vA,ps=1,col=122
;
a=rtable('../'+dir3+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,xiM,vA,ps=1,col=55,th=4
;
a=rtable('../'+dir32+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,xiM,vA,ps=1,col=55,th=4
;
a=rtable('../'+dir33+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,xiM,vA,ps=1,col=188,th=4
;
a=rtable('../'+dir34+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,xiM,vA,ps=1,col=188,th=4
;
a=rtable('../'+dir4+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,1.06*xiM,1.09*vA,ps=1,col=188
;
a=rtable('../'+dir42+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
loadct,6
oplot,xiM,vA,ps=1,col=144
loadct,5
;
;-----------------------------------------------------------------------------
xx=[2.,80.] & oplot,xx,xx/1e4
xx=[2.,80.] & oplot,xx,xx/2e4,li=3
xx=[2.,80.] & oplot,xx,xx/4e4,li=1
xx=[.4,80.] & oplot,xx,1e-2/xx,li=1
;
print,"$mv idl.ps ~/tex/teach/COSMOMAG26/3_decaying_turbulence/fig/pkpm.eps"
print,"$mv idl.ps ~/Overleaf/Pallavi/Decaying-MHD-Turbulence---Need-to-Converge/fig/pkpm.eps"
END