;$Id: pkpm.pro,v 1.5 2026/10/25 20:59:33 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
;
siz=1.5
si2=1.0
!p.charsize=1.8
!x.margin=[8.0,2.1]
!y.margin=[3.2,0.3]
;
;  good examples:
;  2048_k1_e1em8_k60a_h2em13 & nless=4
;  2048_k1_e1em7_k60a_h1em14
;
nless=4
file='kpm.txt'
a=rtable(file,6,head=1)
t=reform(a[1,*])
EMag=reform(a[2,*])
xiMag=reform(a[3,*])
;
lnt=alog(t)
p=-deriv(lnt,alog(EMag))
q=+deriv(lnt,alog(xiMag))
;plot_io,t,xiMag,yr=[1e-7,1e-1]
;
!x.title='!8q!6'
!y.title='!8p!6'
xr=[0,1] & yr=[0,2.2]
xr=[0,.5] & yr=[0,1.0]
plot,q,p,xr=xr,yr=yr,ps=-1,li=1,/nodata
;
;oplot,[0,1],[2,0]
oplot,[0,.5],[2,0],li=1
oplot,[0,1],[0,1]
oplot,[0,1],[0,1]*2
;oplot,[0,1],[0,1]*2.5,li=2
oplot,[0,1],[0,1]*3.
;oplot,[0,1],[0,1]*4.
oplot,[0,1],[0,1]*5.
;
oplot,[1,1]/3.,yr,li=3
oplot,xr,[1,1]*2./3.,li=3
xyouts,.02,.70,siz=siz,'!8p!6=2/3'
xyouts,.25,.04,siz=siz,'!8q!6=1/3'
;
xyouts,.44,.50,siz=siz,'!7b!6=0'
xyouts,.44,.83,siz=siz,'!7b!6=1'
xyouts,.21,.78,siz=siz,'!7b!6=2'
xyouts,.12,.86,siz=siz,'!7b!6=4'
;
;oplot,[1,1]*2./5.,yr,li=3
;oplot,xr,[1,1]*2./5.,li=3
;
for j=1,n_elements(q)-nless do begin
  s=1.4*t[j]^.2
  circ_sym,s,1
  oplot,[1,1]*q[j],[1,1]*p[j],ps=8,li=0,col=122
  wait,.01
endfor
;
cwd,run
print,"$mv idl.ps ~/Overleaf/Vinay/Helicity-Gauge/figs/pkpm_"+run+".eps"
END