;$Id: pslice_Pm10.pro,v 1.2 2015/12/03 05:56:11 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=18,yoffset=3
  !p.charthick=2.2 & !p.thick=2.2 & !x.thick=2.2 & !y.thick=2.2
  col=255
endif else begin
  col=0
end
;
;  convert idl.ps ~/tex/prop/NSF/2015/machine/pslice_Pm10.png
;
!p.charsize=1.2
!x.margin=[6.0,.5]
!y.margin=[3.2,.5]
siz=2.4
;
restore,'data/slice_ubz.sav'
!p.multi=0
print,minmax(bz4)
print,minmax(uz4)
print,minmax(bz128)
print,minmax(uz128)
lev=grange(-1.,1.,20)
lev_bz4=lev*.02
lev_uz4=lev*.01
lev_bz128=lev*.01
lev_uz128=lev*.005
xbox=[0,1,1,0,0]*128./2304.*2.*!pi
ybox=[0,0,1,1,0]*128./2304.*2.*!pi
;
x4=findgen(576)/576.*2*!pi
x128=findgen(128)/2304.*2.*!pi
;
k0=60
kMk0=8.
kM=kMk0/k0
!x.title='!8k!6!dM!n !8x!6'
!y.title='!8k!6!dM!n !8y!6'
!x.ticks=2 & !x.tickv=[0,.4,.8]
!y.ticks=2 & !y.tickv=[0,.4,.8]
;contour,clip(bz4,minmax(lev_bz4)),lev=lev_bz4,/fill,kM*x4,kM*x4
;oplot,kM*xbox,kM*ybox,col=255
;oplot,li=3,[0,.44],[kM*max(xbox),.78],col=col,thick=5
;oplot,li=3,[kM*max(xbox),.78],[.0,.44],col=col,thick=5
;xyouts,-.14,.67,siz=siz2,'(a)'
;
contour,clip(uz4,minmax(lev_uz4)),lev=lev_uz4,/fill,kM*x4,kM*x4
oplot,kM*xbox,kM*ybox,col=255
oplot,li=3,[0,.44],[kM*max(xbox),.78],col=col,thick=5
oplot,li=3,[kM*max(xbox),.78],[.0,.44],col=col,thick=5
xyouts,.15,.75,siz=3.7,'(b)',col=255,charthick=9
;
!x.title='!6'
!y.title='!6'
!p.charsize=0.8
!x.ticks=1 & !x.tickv=[0,.04]
!y.ticks=1 & !y.tickv=[0,.04]
dy=.385 & x1=.58 & x2=.96
y2=.97
y1=y2-dy
!p.position=[x1,y1,x2,y2]
;contour,clip(bz128,minmax(lev_bz128)),lev=lev_bz128,/fill,kM*x128,kM*x128,/noeras,col=col
;
;y2=.47
;y1=y2-dy
;!p.position=[x1,y1,x2,y2]
contour,clip(uz128,minmax(lev_uz128)),lev=lev_uz128,/fill,kM*x128,kM*x128,/noeras,col=col
;
!p.position=0
END