;$Id: pkpm_comp.pro,v 1.11 2026/10/24 18:59:19 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.multi=[0,1,3]
!p.charsize=1.7
!x.margin=[6.2,.4]
!y.margin=[3.2,.4]
;
xxx0=8e-4
xr=[xxx0,3e2]
yrB=[1e-0,3e1]
!x.title='!8t!6'
!y.title='!8B!6(!8t!6) !8t!6/!7n!6!s!dM!n!r!u2!n'
;
file='/kpm.txt'
;
;-----------------------------------------------------------------------------
dir1='2048_e1em5_k180a'  ;black, t=2.45E+05
dir2='2048_e1em5_k180b'  ;red, t=4.64e+04 (did reach it)
dir3='2048_e1em5_k180c' ;blue, t=6.81e+04 (didn't reach 1e4)
dir4='2048_e1em5_k360a'  ;orange, t=6.81e+03
dir4='2048_k1_e1em5_k60a'  ;orange, t=6.81e+03
dir3='2048_k1_e1em5_k360b'  ;orange, t=6.81e+03
dir1='2048_k1_e1em5_k360a'  ;orange, t=6.81e+03
dir2='2048_k1_e1em5_k180a'  ;orange, t=6.81e+03
;
dir1='2048_k1_e1em6_k60c'  ;black
dir2='2048_k1_e2em6_k60c'  ;red
dir3='2048_k1_e5em6_k60c'  ;blue
dir4='2048_k1_e1em5_k60a'  ;orange
dir5='2048_k1_e1em5_k60b'  ;yellow
dir6='2048_k1_e1em5_k60c'  ;green
;
dir1='2048_k1_e1em4_k60c'  ;black
dir2='2048_k1_e2em5_k60c'  ;red
dir3='2048_k1_e5em6_k60c'  ;blue
dir4='2048_k1_e1em7_k60c'  ;orange
dir5='2048_k1_e1em8_k60c'  ;yellow
dir5='4096_k1_e1em9_k60c'  ;yellow
dir6='2048_k1_e1em9_k60c'  ;green
;
dir3='2048_k1_e1em6_k60a_h2em14'  ;blue
dir4='2048_k1_e1em7_k60a_h1em14'  ;orange
dir5='2048_k1_e1em5_k60a_h1em14'  ;yellow
dir6='2048_k1_e1em6_k60a_h1em14'  ;green
dir6='2048_k1_e1em6_k60c_h1em14'  ;green
dir2='2048_k1_e1em6_k180a_h1em14'  ;green
dir6='2048_k1_e1em3_k60a_h2em13'  ;green  (dt=1e-5)
dir6='2048_k1_e1em3_k60a_h2em13'  ;green  (dt=1e-5)
dir6='2048_k1_e1em4_k60a_h1em14'  ;green  (dt=1e-5)
dir6='2048_k1_e1em6_k60a_h2em13_hel'  ;green  (dt=1e-5)
;
dir1='2048_k1_e1em3_k60a_h2em13'  ;green  (dt=1e-5)
dir2='2048_k1_e1em4_k60a_h1em14'  ;green  (dt=1e-5)
dir3='2048_k1_e1em5_k60a_h1em14'  ;yellow
dir4='2048_k1_e1em6_k60a_h1em14'  ;green
dir5='2048_k1_e1em7_k60a_h1em14'  ;orange  (dt=3e-5)
dir6='2048_k1_e1em8_k60a_h2em13'  ;orange  (dt=3e-5)
;
eta1=1e-3
eta2=1e-4
eta3=1e-5
eta4=1e-6
eta5=1e-7
eta6=1e-8
;
; 2048_k1_e1em5_k360a 2048_e1em5_k360a
; 2048_k1_e1em5_k180a 2048_e1em5_k180c
;  i   t      EMag       xiMag        CM          Lu
;  & CM=a[4,*] & Lu=a[5,*]
;
yr=yrB
tcrit=30.
a=rtable('../'+dir1+file,head=1,6)
t=reform(a[1,*]) & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
plot_oo,t,vA*t/xiM^2,ps=1,xr=xr,yr=yr
;
;circ_sym,1.3,1 & it=findex(tcrit,t) & oplot,xiM[it]*[1,1],vA[it]*[1,1],ps=8
;
a=rtable('../'+dir2+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,t,vA*t/xiM^2,ps=1,col=122
;circ_sym,1.3,1 & it=findex(tcrit,t) & oplot,xiM[it]*[1,1],vA[it]*[1,1],ps=8,col=122
;
a=rtable('../'+dir3+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,t,vA*t/xiM^2,ps=1,col=55
;circ_sym,1.3,1 & it=findex(tcrit,t) & oplot,xiM[it]*[1,1],vA[it]*[1,1],ps=8,col=55
;
a=rtable('../'+dir4+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,t,vA*t/xiM^2,ps=1,col=155
;circ_sym,1.3,1 & it=findex(tcrit,t) & oplot,xiM[it]*[1,1],vA[it]*[1,1],ps=8,col=155
;
a=rtable('../'+dir5+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,t,vA*t/xiM^2,ps=1,col=188
;circ_sym,1.3,1 & it=findex(tcrit,t) & oplot,xiM[it]*[1,1],vA[it]*[1,1],ps=8,col=188
;
a=rtable('../'+dir6+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
loadct,6
oplot,t,vA*t/xiM^2,ps=1,col=144
;circ_sym,1.3,1 & it=findex(tcrit,t) & oplot,xiM[it]*[1,1],vA[it]*[1,1],ps=8,col=144
loadct,5
;
xx=[xxx0,1.] & oplot,xx,0.*xx+8.5,li=1
xx=[xxx0,10.] & oplot,xx,0.*xx+15.,li=1
xx=[xxx0,100.] & oplot,xx,0.*xx+21.,li=1
;
;-----------------------------------------------------------------------------
yr=[1e-0,5e7]
a=rtable('../'+dir1+file,head=1,6)
t=reform(a[1,*]) & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
plot_oo,t,vA/eta1,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,t,vA/eta2,ps=1,col=122
;
a=rtable('../'+dir3+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,t,vA/eta3,ps=1,col=55
;
a=rtable('../'+dir4+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,t,vA/eta4,ps=1,col=155
;
a=rtable('../'+dir5+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,t,vA/eta5,ps=1,col=188
;
a=rtable('../'+dir6+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
loadct,6
oplot,t,vA/eta6,ps=1,col=144
loadct,5
;
;-----------------------------------------------------------------------------
xr=[1e-0,5e7]
yr=yrB
a=rtable('../'+dir1+file,head=1,6)
t=reform(a[1,*]) & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
plot_oo,vA/eta1,vA*t/xiM^2,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,vA/eta2,vA*t/xiM^2,ps=1,col=122
;
a=rtable('../'+dir3+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,vA/eta3,vA*t/xiM^2,ps=1,col=55
;
a=rtable('../'+dir4+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,vA/eta4,vA*t/xiM^2,ps=1,col=155
;
a=rtable('../'+dir5+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,vA/eta5,vA*t/xiM^2,ps=1,col=188
;
a=rtable('../'+dir6+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
loadct,6
oplot,vA/eta6,vA*t/xiM^2,ps=1,col=144
loadct,5
;
;-----------------------------------------------------------------------------
print,"$mv idl.ps ~/Overleaf/Vinay/Helicity-Gauge/figs/pkpm_comp.eps"
print,"$mv idl.ps ~/Overleaf/Vinay/Helicity-Gauge/figs/pkpm_comp_eta.eps"
!p.multi=0
END