;$Id: ppdf2d_comp_BE_3times.pro,v 1.2 2026/08/03 09:17:55 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=7,yoffset=3
  !p.charthick=2 & !p.thick=2 & !x.thick=2 & !y.thick=2
  col0=0
endif else begin
  col0=255
end
;
;  Used in paper to compare 2 times.
;
;  note that the idl save files are produced with pvar.pro, which, in turn,
;  calls pdf2d,lr,lb,xx,yy,ww,nx,ny,xr=[-1,1]*1.8,yr=[-16.,2.],/unnorm
;
siz=0.9
!p.charsize=1.5
!x.margin=[5.9,0.2]
!y.margin=[3.2,1.8]
!p.multi=[0,3,1]
;
dir1='Hy512h_noshock2'
dir2=dir1
dir3=dir1
file1='Larger_pdf2d_lr_lo_26'
file1='Larger_pdf2d_lr_lo_0'
file1='Larger_pdf2d_lr_lo_1'
file2='Larger_pdf2d_lr_lo_28'
file2='Larger_pdf2d_lr_lo_26'
file3='Larger_pdf2d_lr_lo_30'
offset1=+2.7
offset2=5.1
;
xr1=[-1.5,+4.5] & yr1=[-7.0,+2.]
xr2=[-1.5,+4.5] & yr2=[-7.0,+2.]
lev2=grange(-12.,+0.,11)
lev1=lev2
;
;  Panel 1
;
!x.title='!6ln!7q!6'
!y.title='!6ln!7x!6'
!p.title='!8t!6=0.1'
restore,'../'+dir1+'/'+file1+'.sav'
contour,-alog(ww>1.),xx,yy,lev=lev1,/fil,xr=xr1,yr=yr1
;xr=[-1.5,4.0] & oplot,xr,xr*2./3.-offset1,col=0
;xr=[-1.5,1.5] & oplot,xr,xr*2.   -offset2,col=0,li=2
;xr=[-1.5,1.5] & oplot,xr,xr*3.   -offset2,col=0,li=2
;
;  Panel 2
;
!y.title='!6'
!p.title='!8t!6=2.6'
restore,'../'+dir2+'/'+file2+'.sav'
contour,-alog(ww>1.),xx,yy,lev=lev2,/fil,xr=xr2,yr=yr2
xr=[-1.5,4.0] & oplot,xr,xr*2./3.-offset1,col=0
;xr=[-1.5,1.5] & oplot,xr,xr*2.   -offset2,col=0,li=2
;xyouts,3.1,0.7,siz=siz,'!7x!9A!7q!6!u2/3!n'
;xyouts,0.1,-6.3,siz=siz,'!7x!9A!7q!6!u2!n'
xyouts,2.8,0.4,siz=siz,'!7x!9A!7q!6!u2/3!n'
;xyouts,0.1,-6.3,siz=siz,'!7x!9A!7q!6!u2!n'
;
;  Panel 3
;
!y.title='!6'
!p.title='!8t!6=3.0'
restore,'../'+dir3+'/'+file3+'.sav'
contour,-alog(ww>1.),xx,yy,lev=lev2,/fil,xr=xr2,yr=yr2
xr=[-1.5,4.0] & oplot,xr,xr*2./3.-offset1,col=0
;xr=[-1.5,1.5] & oplot,xr,xr*2.   -offset2,col=0,li=2
;xyouts,3.1,0.7,siz=siz,'!7x!9A!7q!6!u2/3!n'
;xyouts,0.1,-6.3,siz=siz,'!7x!9A!7q!6!u2!n'
;
col=0
!p.title='!6ln!8P!6'
!y.title='!6'
!x.title='!6'
xx0=.28     & yy1=.19 & dxx=.03 & dyy=.20
colorbar,pos=[xx0,yy1,xx0+dxx,yy1+dyy], range=minmax(lev1),/bot,/ver,$
form='(f5.1)',charsize=1.2,div=2,ytit='!6',col=col,/rev
;
xx0=.28+.33 & yy1=.19 & dxx=.03 & dyy=.20
colorbar,pos=[xx0,yy1,xx0+dxx,yy1+dyy], range=minmax(lev2),/top,/ver,$
form='(f5.1)',charsize=1.2,div=2,ytit='!6',col=col,/rev
;
xx0=.28+.67 & yy1=.19 & dxx=.03 & dyy=.20
colorbar,pos=[xx0,yy1,xx0+dxx,yy1+dyy], range=minmax(lev2),/top,/ver,$
form='(f5.1)',charsize=1.2,div=2,ytit='!6',col=col,/rev
;
print,"$mv idl.ps ~/Overleaf/Eva/BE_collapse/fig/ppdf2d_comp_BE_3times.eps"
END