;$Id: pppvar_9panels.pro,v 1.1 2026/06/08 11:40:41 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=17.8,yoffset=3
  !p.charthick=2 & !p.thick=2 & !x.thick=2 & !y.thick=2
end
;
!p.charsize=1.7
!x.margin=[5.9,.4]
!y.margin=[3.2,.2]
!p.multi=[0,3,3]
nlev=17 & color,/sgi,steps=nlev
tilde='!9!s!aA!n!r!6'
;
restore,'b0perp.sav'
;
lev=grange(0,1,nlev)*1
lev2=grange(-1,1,nlev)*!pi
xlab=5.2
ylab=5.0
;
for j=0,2 do begin
  contour,clip(abs(b0perp[*,*,j]),minmax(lev)),x,y,lev=lev,/fil,/iso
  xyouts,xlab,ylab,'!8r!6!d'+str(j+1)+'!n',col=255
endfor
;
for j=0,2 do begin
  contour,arg(b0perp1[*,*,j]),x,y,lev=lev2,/fil,/iso
  xyouts,xlab,ylab,'!7h!6!d'+str(j+1)+'!n',col=255
endfor
;
for j=0,2 do begin
  contour,arg(b0perp2[*,*,j]),x,y,lev=lev2,/fil,/iso
  xyouts,xlab,ylab,tilde+'!7h!6!d'+str(j+1)+'!n',col=255
endfor
;
;
!x.title='!6'
ytickn=['!6-!7p!6','!6-!7p!6/2','!60','!7p!6/2','!7p!6']
yy0=.08 & yy1=yy0+2./3. & dxx=.01 & dyy=.22 & xx0=.085
colorbar,pos=[xx0,yy1,xx0+dxx,yy1+dyy], range=minmax(lev),/bot,/ver,$
form='(f4.1)',charsize=2.0,div=5,ytit='!6',col=255,/right
;
yy1=yy0+1./3.
colorbar,pos=[xx0,yy1,xx0+dxx,yy1+dyy], range=minmax(lev2),/bot,/ver,$
form='(f4.1)',charsize=2.0,div=4,ytit='!6',col=255,/right,ytickn=ytickn
;
yy1=yy0
colorbar,pos=[xx0,yy1,xx0+dxx,yy1+dyy], range=minmax(lev2),/bot,/ver,$
form='(f4.1)',charsize=2.0,div=4,ytit='!6',col=255,/right,ytickn=ytickn
;
cwd,run
!p.multi=0
print,"$mv idl.ps ~/Overleaf/Axel/Roberts-Flow-Test/fig/pppvar_"+run+".eps"
END