;$Id: pres_eOT.pro,v 1.2 2026/10/12 02:10:29 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"
;
;  E_B = B^2/2, so B=sqrt(2*EB)=sqrt(2*.5)=1
;                           or =sqrt(2*2)=2
;
file='res.txt' & nhead=15
file='res_eOT.txt' & nhead=1 & B0=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='!6Lu = !8B!6!d0!n/!7g!8k!6!d0!n'
!y.title='!7e!6!dK!n  and  !7e!6!dM!n'
xr=1./[3e-5,7e-3]
xr=1./[1e-6,2e-3]
;
;-----------------------------------------------------------------------------
yr=[1e-4,1e-0]
yr=[1e-3,3e-2]
circ_sym,1.3,1
plot_oo,xr,yr,ps=8,/nodata
;
good=where(nx eq 1024)
oplot,B0/nu(good),epsK(good),ps=8,col=55
;
good=where(nx eq 2048)
oplot,B0/nu(good),epsK(good),ps=8,col=122
;
loadct,6
good=where(nx eq 4096)
oplot,B0/nu(good),epsK(good),ps=8,col=144
loadct,5
;
;good=where(nx eq 8192)
;oplot,B0/nu(good),epsK(good),ps=8,col=155
;
;good=where(nx eq 16384)
;oplot,B0/nu(good),epsK(good),ps=8,col=55
;
;
circ_sym,1.3,0
good=where(nx eq 1024)
oplot,B0/nu(good),epsM(good),ps=8,col=55
;
good=where(nx eq 2048)
oplot,B0/nu(good),epsM(good),ps=8,col=122
;
loadct,6
good=where(nx eq 4096)
oplot,B0/nu(good),epsM(good),ps=8,col=144
loadct,5
;
;good=where(nx eq 8192)
;oplot,B0/nu(good),epsM(good),ps=8,col=155
;
;good=where(nx eq 16384)
;oplot,B0/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=[6e0,9e0]
circ_sym,1.3,1
plot_oi,xr,yr,ps=8,/nodata
;
good=where(nx eq 1024)
oplot,B0/nu(good),tepsK(good),ps=8,col=55
;
good=where(nx eq 2048)
oplot,B0/nu(good),tepsK(good),ps=8,col=122
;
loadct,6
good=where(nx eq 4096)
oplot,B0/nu(good),tepsK(good),ps=8,col=144
loadct,5
;
;good=where(nx eq 8192)
;oplot,B0/nu(good),tepsK(good),ps=8,col=155
;
;good=where(nx eq 16384)
;oplot,B0/nu(good),tepsK(good),ps=8,col=55
;
;
circ_sym,1.3,0
good=where(nx eq 1024)
oplot,B0/nu(good),tepsM(good),ps=8,col=55
;
good=where(nx eq 2048)
oplot,B0/nu(good),tepsM(good),ps=8,col=122
;
loadct,6
good=where(nx eq 4096)
oplot,B0/nu(good),tepsM(good),ps=8,col=144
loadct,5
;
;good=where(nx eq 8192)
;oplot,B0/nu(good),tepsM(good),ps=8,col=155
;
;good=where(nx eq 16384)
;oplot,B0/nu(good),tepsM(good),ps=8,col=55
;
!p.multi=0
print,"$mv idl.ps ~/Overleaf/Liubov/Dissipation-code-comparison/fig/pres_eOT.eps"
END