;$Id: pres.pro,v 1.11 2026/10/11 20:23:44 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=16,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,.5]
;print,"$mv idl.ps ~/Overleaf/Axel/Transfer/fig/p1.eps"
;
file='res.txt' & nhead=15
file='res_eOT.txt' & nhead=1
file='res_x2.txt' & nhead=1
;
a=rtable(file,head=nhead,8)
nu=reform(a[0,*])
epsK=reform(a[2,*])
tepsK=reform(a[3,*])
epsM=reform(a[4,*])
tepsM=reform(a[5,*])
nx=fix(reform(a[6,*]))
;nu       eta    epsK   tepsK  epsM   tepsM   nx    tmax  run
;
!p.multi=[0,1,2]
!x.title='!7m!6'
!y.title='!7e!6!dK!n  and  !7e!6!dM!n'
xr=[3e-5,7e-3]
xr=[1e-6,2e-3]
xr=[3e-3,1e-1]
;
;-----------------------------------------------------------------------------
yr=[1e-4,1e-0]
yr=[3e-3,2e-2]
yr=[3e-2,2e-0]
circ_sym,1.3,1
plot_oo,xr,yr,ps=8,/nodata
;
good=where(nx eq 1024)
oplot,nu(good),epsK(good),ps=8;,col=55
;
;good=where(nx eq 2048)
;oplot,nu(good),epsK(good),ps=8,col=122
;
;loadct,6
;good=where(nx eq 4096)
;oplot,nu(good),epsK(good),ps=8,col=144
;loadct,5
;
;good=where(nx eq 8192)
;oplot,nu(good),epsK(good),ps=8,col=155
;
;good=where(nx eq 16384)
;oplot,nu(good),epsK(good),ps=8,col=55
;
circ_sym,1.3,0
good=where(nx eq 1024)
oplot,nu(good),epsM(good),ps=8;,col=55
;
;good=where(nx eq 2048)
;oplot,nu(good),epsM(good),ps=8,col=122
;
;loadct,6
;good=where(nx eq 4096)
;oplot,nu(good),epsM(good),ps=8,col=144
;loadct,5
;
;good=where(nx eq 8192)
;oplot,nu(good),epsM(good),ps=8,col=155
;
;good=where(nx eq 16384)
;oplot,nu(good),epsM(good),ps=8,col=55
;
;-----------------------------------------------------------------------------
!y.title='!8t!6!dK!n  and  !8t!6!dM!n'
yr=[1e-2,1e-0]
yr=[.17,.33]
yr=[0e-2,2e0]
circ_sym,1.3,1
plot_oi,xr,yr,ps=8,/nodata
;
good=where(nx eq 1024)
oplot,nu(good),tepsK(good),ps=8;,col=55
;
;good=where(nx eq 2048)
;oplot,nu(good),tepsK(good),ps=8,col=122
;
;loadct,6
;good=where(nx eq 4096)
;oplot,nu(good),tepsK(good),ps=8,col=144
;loadct,5
;
;good=where(nx eq 8192)
;oplot,nu(good),tepsK(good),ps=8,col=155
;
;good=where(nx eq 16384)
;oplot,nu(good),tepsK(good),ps=8,col=55
;
;
circ_sym,1.3,0
good=where(nx eq 1024)
oplot,nu(good),tepsM(good),ps=8;,col=55
;
;good=where(nx eq 2048)
;oplot,nu(good),tepsM(good),ps=8,col=122
;
;loadct,6
;good=where(nx eq 4096)
;oplot,nu(good),tepsM(good),ps=8,col=144
;loadct,5
;
;good=where(nx eq 8192)
;oplot,nu(good),tepsM(good),ps=8,col=155
;
;good=where(nx eq 16384)
;oplot,nu(good),tepsM(good),ps=8,col=55
;
!p.multi=0
print,"$mv idl.ps ../tex/fig/pres.eps"
END