;$Id: pjz_hist_comp.pro,v 1.3 2026/06/06 18:48:08 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=11,yoffset=3
  !p.charthick=2 & !p.thick=2 & !x.thick=2 & !y.thick=2
end
;
siz=1.7
si2=1.0
!p.charsize=1.7
!x.margin=[7.9,.6]
!y.margin=[3.2,.3]
!x.title='!8x!6'
!y.title='!8y!6'
!p.position=0
;
file1='../8192_k1em2_k1em0_vA1em1_nu1em5/jz_ivar1000.sav'   ;panel (c), Run A at 1e5 (-> 1.7e5) + 1100
file2='../8192_k1em3_k1em1_vA9em3_nu1em5u2/jz_ivar200.sav'  ;panel (b), Run C at 
file3='../8192_k1em3_k1em1_vA9em3_nu1em5u2/jz_ivar10.sav'   ;panel (a), Run C
;
;  full domain, but downsampled with aver option
;
default,lread,1
if lread then begin
  restore,file1
  print,file1,',  t=',t
  s=size(jz)
  nx=s[1]
  ny=s[2]
  jz1=jz
  ;
  restore,file2
  print,file2,',  t=',t
  s=size(jz)
  nx=s[1]
  ny=s[2]
  jz2=jz
  ;
  restore,file3
  print,file3,',  t=',t
  s=size(jz)
  nx=s[1]
  ny=s[2]
  jz3=jz
  lread=0
  ;
  n=500
  pdf,jz1,xx1,yy1,n=n
  pdf,jz2,xx2,yy2,n=n
  pdf,jz3,xx3,yy3,n=n
endif
;
!p.multi=[0,1,3]
!p.multi=0
!x.range=[-1,1]*27.
!y.range=[3e-7,2]
!x.title='!8J!dz!n!6/!8J!s!dz!n!6!r!urms!n'
!y.title='!8p!6(!8J!dz!n!6/!8J!s!dz!n!6!r!urms!n)'
siz=1.7
;
;s1=sqrt(total(xx1^2*yy1)/total(yy1))
dx1=xx1[1]-xx1[0] & s1=sqrt(total(xx1^2*yy1)*dx1)
plot_io,xx1/s1,yy1*s1
xyouts,-20.,4e-5,siz=siz,'!6(c)'
;
dx2=xx2[1]-xx2[0] & s2=sqrt(total(xx2^2*yy2)*dx2)
oplot,xx2/s2,yy2*s2,col=122
xyouts,-17.,4e-6,siz=siz,'!6(b)',col=122
;
dx3=xx3[1]-xx3[0] & s3=sqrt(total(xx3^2*yy3)*dx3)
oplot,xx3/s3,yy3*s3,col=55
xyouts,-6.9,4e-6,siz=siz,'!6(a)',col=55
;
print,'minmax(jz)=',minmax(jz)
is=16 & is0=is & is1=16 & is2=2
;
print,"$mv idl.ps ~/Overleaf/Mattia/Isochrones/fig/pjz_hist_comp.eps"
END