;Id$
;
;  The eigenfunctions are plotted afterwards using ".r pevec_vectors2.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=5.8,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ../fig/pevec_vectors2.ps
;
!p.charsize=2.0
!x.margin=[4.2,0.5]
!y.margin=[2.2,2.0]
;
!p.multi=0
!p.multi=[0,3,1]
!x.title='!6'
!y.title='!6'
;
nx=n
ny=80 & y=40*(findgen(ny)-ny/2)/(ny-1)
nlev=30 & lev=20.*(findgen(nlev)-nlev/2)/(nlev-1)
;
!x.ticks=2
!x.range=[-1,1]*8
;!y.range=[-1,1]*12
;
ik=10
print,k[ik]
intx=alog(cosh(x))
psi0=rebin(reform(intx,nx,1),nx,ny)
psi1=float(reform(eveck[*,ik])#exp(ii*k[ik]*y))
;
!p.title='!7e!6 = 0.05'
contour,psi0+.05*psi1,x,y,lev=lev,/iso
;
!p.title='!7e!6 = 0.5'
contour,psi0+.5*psi1,x,y,lev=lev,/iso
;
!p.title='!7e!6 = 2'
contour,psi0+2.*psi1,x,y,lev=lev,/iso
;
end