;$Id: pjz_comp.pro,v 1.5 2026/06/07 03:03:42 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=17.2,yoffset=3
  !p.charthick=2 & !p.thick=2 & !x.thick=2 & !y.thick=2
end
;
siz=1.7
si2=1.0
!p.charsize=1.0
!x.margin=[8.9,.6]
!y.margin=[3.2,.3]
!x.title='!8x!6'
!y.title='!8y!6'
!p.position=0
loadct,0
;
file1='../8192_k1em2_k1em0_vA1em1_nu1em5/jz_ivar1000.sav'
file1='../8192_k1em2_k1em0_vA1em1_nu1em5/jz_ivar1125.sav'
file2='../8192_k1em3_k1em1_vA9em3_nu1em5u2/jz_ivar200.sav'
file3='../8192_k1em3_k1em1_vA9em3_nu1em5u2/jz_ivar10.sav'
xout=105.
yout=130.
;
;  full domain, but downsampled with aver option
;
restore,file2
print,file2,',  t=',t & tb=t & tb01=69000. & tb02=76000.
s=size(jz)
nx=s[1]
ny=s[2]
print,'minmax(jz)=',minmax(jz)
is=16 & is0=is & is1=16 & is2=2
jmax=.002
contour,sample_2d(jz,is0,is0,/aver),lev=grange(-1,1,30)*jmax*.5,/fill,/iso,sample_1d(x,is0),sample_1d(y,is0)
;
;-----------------------------------------------------------------------------
;  same file, but subdomain
;
loadct,5
nx2=nx/is
ny2=ny/is
dp=.3
xp1=.23 & yp1=.2 & yp1b=.20
xp3=.65 & yp3=.64
xp2=xp1+dp & yp2=yp1+dp
xp4=xp3+dp & yp4=yp3+dp
;
x2=x[nx-(is/2+.5)*nx2:nx-(is/2+.5)*nx2+nx2-1]
y2=x[nx-(is/2+.5)*nx2:nx-(is/2+.5)*nx2+nx2-1]
jz2=jz[nx-(is/2+.5)*nx2:nx-(is/2+.5)*nx2+nx2-1,nx-(is/2+.5)*nx2:nx-(is/2+.5)*nx2+nx2-1]
xr=minmax(x2)
yr=minmax(y2)
;
xx=[xr[0],xr[1],xr[1],xr[0],xr[0]]
yy=[yr[0],yr[0],yr[1],yr[1],yr[0]]
oplot,xx,yy,col=255
oplot,[xr[0], 650.],[yr[1],2750.],li=1,col=255
oplot,[xr[1],2700.],[yr[0], 700.],li=1,col=255
;
!x.title='!8J!dz!n!6'
!x.title='!6'
!y.title='!6'
xx0=.26 & yy1=.20 & dxx=.26 & dyy=.03
xx0=.46 & yy1=.10 & dxx=.26 & dyy=.03
colorbar,pos=[xx0,yy1,xx0+dxx,yy1+dyy], range=[-1,1],/top,/hor,$
form='(f4.1)',charsize=.009,div=2,ytit='!6',col=255
;
!x.title='!8x!6'
!y.title='!8y!6'
!p.position=[xp3,yp3,xp4,yp4]
;
help,jz2
print,minmax(jz2)
contour,sample_2d(jz2,is2,is2),lev=grange(-1,1,30)*jmax,/fill,/iso,sample_1d(x2,is2),sample_1d(y2,is2),col=255,/noerase
xyouts,xout,yout,siz=siz,col=255,'!6(b)'
;
;-----------------------------------------------------------------------------
!p.position=[xp1,yp3,xp2,yp4]
restore,file3
print,file3,',  t=',t & ta=t & ta01=69000. & ta02=76000.
jz2=jz[nx-(is/2+.5)*nx2:nx-(is/2+.5)*nx2+nx2-1,nx-(is/2+.5)*nx2:nx-(is/2+.5)*nx2+nx2-1]
help,jz2
is2=2
print,minmax(jz2)
contour,sample_2d(jz2,is2,is2),lev=grange(-1,1,30)*jmax*2,/fill,/iso,sample_1d(x2,is2),sample_1d(y2,is2),col=255,/noerase
xyouts,xout,yout,siz=siz,col=255,'!6(a)'
;
;-----------------------------------------------------------------------------
!p.position=[xp3,yp1,xp4,yp2]
restore,file1
print,file1,',  t=',t & tc=t & tc01=1100. & tc02=740.
contour,sample_2d(jz,is1,is1,/aver),lev=grange(-1,1,30)*jmax*1.,/fill,/iso,sample_1d(x,is1),sample_1d(y,is1),/noerase,xr=xr,yr=yr,col=255
xyouts,xout,yout,siz=siz,col=255,'!6(c)'
;
;-----------------------------------------------------------------------------
!p.position=[xp1+.02,yp1b,xp2,yp2]
!x.title='!7n!6!dM!n(!8t!6)'
!y.title='!8v!6!dA!n(!8t!6)'
;
dir4='8192_k1em2_k1em0_vA1em1_nu1em5' & istart4=2+6 ;green, running, CB=2.9-2.95
dir5='8192_k1em3_k1em1_vA9em3_nu1em5u2' & istart5=0+6*0 & i1_dir2=21 & i2_dir2=29  ;Orange
;
xr=[.4,40]
yr=[1e-4,2.2e-2]
!x.title='!7n!6!dM!n(!8t!6)'
!y.title='!8v!6!dA!n(!8t!6)'
file='/kpm.txt'
;
th=5.
siz=.8
circ_sym,1.1,1
dir=dir5 & i1=istart5
a=rtable('../'+dir+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
plot_oo,xiM,vA,xr=xr,yr=yr,xst=9,yst=9,/noerase,col=255,/nodata
oplot,xiM,vA,col=155,th=th
it=findex(ta+ta02,t) & oplot,[1,1]*xiM[it],[1,1]*vA[it],col=155,ps=8
it=findex(tb+tb02,t) & oplot,[1,1]*xiM[it],[1,1]*vA[it],col=155,ps=8
xyouts,4.,.0012,siz=siz,'(a)',col=155
xyouts,7.,.0006,siz=siz,'(b)',col=155
;
loadct,6
circ_sym,.9,1
dir=dir4 & i1=istart4
a=rtable('../'+dir+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,xiM,vA,col=144,th=th
it=findex(tc+tc01,t) & oplot,[1,1]*xiM[it],[1,1]*vA[it],col=144,ps=8
xyouts,12.,.0014,siz=siz,'(c)',col=144
loadct,5
;
print
print,ta,ta01,ta+ta01
print,tb,tb01,tb+tb01
print,tc,tc01,tc+tc01
print
;!p.position=0
print,"$mv idl.ps ~/Overleaf/Mattia/Isochrones/fig/pjz_comp.eps"
END