;$Id: pxy_slice.pro,v 1.1 2022/01/24 05:51:06 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=6.1,yoffset=3
  !p.charthick=1.6 & !p.thick=1.6 & !x.thick=1.6 & !y.thick=1.6
  col0=0
endif else begin
  col0=255
end
;
;  this is to be run in: 2048_1e2_1e4_4e8_1em6aD_k10_bes
;
!p.charsize=1.4
!x.margin=[7.0,0.3]
!y.margin=[3.2,0.6]
tilde='!9!s!aA!n!r!6'
siz=.9
;
!p.multi=[0,3,1]
!x.title='!8k!dx!n!6'
!y.title='!8k!dy!n!6'
dk=5.
lev=grange(-1.,1.,18)*!pi
lev=grange(0.,1.,18)*1e-5
lev=grange(0.,1.,18)*1e-5
fo='(f3.1)'
;
file='xy_slice_0.sav'
restore,file
s=size(xy_slice1)
k=dk*findgen(s[1])
print,minmax(xy_slice1)
!p.title='!8t!6 = '+str(t,fo=fo)
contour,sqrt(xy_slice1^2+xy_slice2^2),k,k,lev=lev,/fil,/iso
;
file='xy_slice_1.sav'
restore,file
!p.title='!8t!6 = '+str(t,fo=fo)
contour,sqrt(xy_slice1^2+xy_slice2^2),k,k,lev=lev,/fil,/iso
;
file='xy_slice_2.sav'
restore,file
!p.title='!8t!6 = '+str(t,fo=fo)
contour,sqrt(xy_slice1^2+xy_slice2^2),k,k,lev=lev,/fil,/iso
;
!x.title='!9!!!6 '+tilde+'!8T !9!!!6'
!y.title='!6'
!p.title='!6'
xx0=.97 & yy1=.35 & dxx=.01 & dyy=.5
colorbar,pos=[xx0,yy1,xx0+dxx,yy1+dyy], range=minmax(lev),/top,/ver,$
form='(e8.1)',charsize=1.2,div=4,ytit='!6',col=255
;
;contour,atan(xy_slice2,xy_slice1),k,k,lev=lev,/fil,/iso
print,'$mv idl.ps ~/GitHub/Axel/GWspectra_from_CME/fig/pxy_slice.eps'
END