%%%%% Begin file driver_eigs.m
Globals;
% physical parameters: (see text for definition)
m=1; n=3; ggamma=1e-2; delta=1e-2; bbeta = 0.5; eta=0.133;
% numerical parameters:
num_elmts = 4;                              % number of DG elements
N = 24; Np = N+1;                           % DG order per cell
z_bottom = 0; z_top =.75*1/delta*(1-1/(1+(ggamma*(1-bbeta))/eta));
[Nv, VZ, K, EToV] = MeshGen1D(z_bottom,z_top,num_elmts); % Generate mesh
Startup_eigs; Z=z(:); % Initialize solver and construct grid and metric
kx = 4.4; % wavenumber of x-perturbation
[eigval,eigvec] = compute_eigs(kx,m,n);
plot(Z/z_top,real(eigvec),'.-'); xlabel('z / H');ylabel('Eigenfunction');
title(['Eigenvalue = ' num2str(eigval)]);
%%%%% End file driver_eigs.m