;$Id: pVAR.pro,v 1.2 2026/07/03 03:05:37 brandenb Exp $
if !d.name eq 'PS' then begin
  ;device,xsize=18,ysize=11.7,yoffset=3
  device,xsize=18,ysize=17.7,yoffset=3
  !p.charthick=1.6 & !p.thick=1.6 & !x.thick=1.6 & !y.thick=1.6
end
;
siz=1.3
!p.charsize=1.7
!x.margin=[6.0,.3]
!y.margin=[3.2,.1]
!p.multi=[0,3,3]
nlev=17 & color,/sgi,steps=nlev
levr=grange(-1,1,nlev)*4.
levd=grange(-1,1,nlev)*4.
levo=grange(-1,1,nlev)*.4
xout=+2.5
yout=-5.5
xoud=-1.6
youd=+4.9
col=0
;
!x.title='!8x!6'
!y.title='!8y!6'
restore,'VAR26_divu.sav'
restore,'VAR26_lnrho.sav'
contour,clip(lnrho,minmax(levr)),lev=levr,/iso,/fil,x,y 
xyouts,xout,yout,siz=siz,col=col,'!8t!6=2.6'
;
!y.title='!6'
restore,'VAR28_lnrho.sav'
contour,clip(lnrho,minmax(levr)),lev=levr,/iso,/fil,x,y 
xyouts,xout,yout,siz=siz,col=col,'!8t!6=2.8'
xyouts,xoud,youd,siz=siz,col=col,'!6ln!7q!6'
;
restore,'VAR30_lnrho.sav'
contour,clip(lnrho,minmax(levr)),lev=levr,/iso,/fil,x,y 
xyouts,xout,yout,siz=siz,col=col,'!8t!6=3.0'
;
!x.title='!6'
!y.title='!8y!6'
restore,'VAR26_divu.sav'
contour,clip(divu,minmax(levd)),lev=levd,/iso,/fil,x,y 
xyouts,xout,yout,siz=siz,col=col,'!8t!6=2.6'
;
!y.title='!6'
restore,'VAR28_divu.sav'
contour,clip(divu,minmax(levd)),lev=levd,/iso,/fil,x,y 
xyouts,xout,yout,siz=siz,col=col,'!8t!6=2.8'
xyouts,xoud,youd,siz=siz,col=col,'!6div!8u!6'
;
restore,'VAR30_divu.sav'
contour,clip(divu,minmax(levd)),lev=levd,/iso,/fil,x,y 
xyouts,xout,yout,siz=siz,col=col,'!8t!6=3.0'
;
!x.title='!8x!6'
!y.title='!8y!6'
restore,'VAR26_oz.sav'
contour,clip(oz,minmax(levo)),lev=levo,/iso,/fil,x,y 
xyouts,xout,yout,siz=siz,col=col,'!8t!6=2.6'
;
!y.title='!6'
restore,'VAR28_oz.sav'
contour,clip(oz,minmax(levo)),lev=levo,/iso,/fil,x,y 
xyouts,xout,yout,siz=siz,col=col,'!8t!6=2.8'
xyouts,xoud,youd,siz=siz,col=col,'!7x!d!8z!6!n'
;
restore,'VAR30_oz.sav'
contour,clip(oz,minmax(levo)),lev=levo,/iso,/fil,x,y 
xyouts,xout,yout,siz=siz,col=col,'!8t!6=3.0'
;
!x.title='!6'
!y.title='!6'
xx0=.08 & yy1=.1+.667 & dxx=.02 & dyy=.2
colorbar,pos=[xx0,yy1,xx0+dxx,yy1+dyy], range=minmax(levd),/rig,/ver,$
form='(f4.1)',charsize=1.2,div=4,ytit='!6',col=col
;
xx0=.08 & yy1=.1+.333 & dxx=.02 & dyy=.2
colorbar,pos=[xx0,yy1,xx0+dxx,yy1+dyy], range=minmax(levd),/rig,/ver,$
form='(f4.1)',charsize=1.2,div=4,ytit='!6',col=col
;
xx0=.08 & yy1=.1 & dxx=.02 & dyy=.2
colorbar,pos=[xx0,yy1,xx0+dxx,yy1+dyy], range=minmax(levo),/rig,/ver,$
form='(f4.1)',charsize=1.2,div=4,ytit='!6',col=col
;
print,"$mv idl.ps ~/Overleaf/Eva/BE_collapse/fig/pVAR.eps"
END