;$Id: pslice.pro,v 1.5 2014/04/08 10:15:42 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=11,ysize=20,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
;
;  mv idl.ps ~/tex/tina/Decay/fig/pslice.ps
;
!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,1,2]
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,-.14,.67,siz=siz2,'(b)'
;
!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=.165 & x1=.63 & 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