;  $Id: pu_xz.pro,v 1.5 2005/07/24 12:30:41 brandenb Exp $
;
;  plots toroidal velocity contours together with velocity vectors
;  in the meridional plane.
;
;  for growth rate, calculate rms velocity in z
;
uzrms=sqrt(mean(uu(l1:l2,m1:m2,n1:n2,2)^2))
print,t,uzrms
;
openw,11,'uzrms.tmp'
printf,11,t,uzrms
close,11
spawn,'cat uzrms.tmp >> uzrms.dat'
;
;  one should see nice Taylor vortices after t=50
;
m=ny/2+2
print,y(m)
xxx=x(l1:l2)
zzz=z(n1:n2)
ux=means(reform(uu(l1:l2,m:m+1,n1:n2,0)),2)
uy=means(reform(uu(l1:l2,m:m+1,n1:n2,1)),2)
uz=means(reform(uu(l1:l2,m:m+1,n1:n2,2)),2)
contour,uy,xxx,zzz,/fil,nlev=30
vel_a,ux,uz,xxx,zzz,nvec=1000,/over,len=.02
;
oplot,-par2.r_int*[1,1],[par.xyz0(0),par.xyz1(0)],thick=5
oplot,+par2.r_int*[1,1],[par.xyz0(0),par.xyz1(0)],thick=5
oplot,-par2.r_int*[1,1],[par.xyz0(0),par.xyz1(0)],thick=4,col=0,li=2
oplot,+par2.r_int*[1,1],[par.xyz0(0),par.xyz1(0)],thick=4,col=0,li=2
;
oplot,-par2.r_ext*[1,1],[par.xyz0(0),par.xyz1(0)],thick=5
oplot,+par2.r_ext*[1,1],[par.xyz0(0),par.xyz1(0)],thick=5
oplot,-par2.r_ext*[1,1],[par.xyz0(0),par.xyz1(0)],thick=4,col=0,li=2
oplot,+par2.r_ext*[1,1],[par.xyz0(0),par.xyz1(0)],thick=4,col=0,li=2
;
;  save
;
save,file=varfile+'.slice',ux,uy,uz,xxx,zzz
;
END