;Id$
;
;  The eigenfunctions are plotted afterwards using ".r pevec.pro"
;  But before that, run on the command line ".r rayleighII.pro".
;  Report problems to Axel Brandenburg <Axel.Brandenburg@Colorado.edu>
;
if !d.name eq 'PS' then begin
  device,xsize=9,ysize=23,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ../fig/pevec_vectors1.ps
;
!p.charsize=1.0
!x.margin=[7.8,.5]
!y.margin=[3.2,.5]
;
!p.multi=[0,2,1]
!p.multi=0
!x.title='!8x!6'
!y.title='!8y!6'
;
nx=n
ny=80 & y=40*(findgen(ny)-ny/2)/(ny-1)
nlev=30 & lev=2.*(findgen(nlev)-nlev/2)/(nlev-1)
;
ik=10
print,k[ik]
psi=float(reform(eveck[*,ik])#exp(ii*k[ik]*y))
ux=float(reform(eveck[*,ik])#(ii*k[ik]*exp(ii*k[ik]*y)))
uy=-xder_6th_bc(psi)
help,psi
!x.range=[-1,1]*4
!y.range=[-1,1]*12
contour,psi,x,y,lev=lev,/fil,/iso
;velovect,ux,uy,x,y,/over
vel_a,ux,uy,x,y,/over,len=.003,nvec=1400,col=255
;
;isamp=4
;x_samp=sample_1d(x,isamp)
;ux_samp=sample_2d(ux,isamp,1)
;uy_samp=sample_2d(uy,isamp,1)
;velovect,ux_samp,uy_samp,x_samp,y,/over
end