%  This program is executed AFTER running lin_gill
%  See ~/Research/alex/matlab3/lin_gill for variable definitions

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  Calculate Modes
%
%  Note:  real(lam) gives the growth rate (> 0 ==> growth, < 0 ==> damped)
%         imag(lam) gives phase speed (> 0 ==> westward, < 0 ==>
%         eastward)
%  Note also:  I've switched the polarity of the phase below, so < 0 on
%  the PLOT is westward.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%  Now, plot the growth rates
%[tem, ind] = sort(real(lamc));  %  Order from least stable to most stable
%[tem, indu] = sort(real(lamu));

%figure_tall(2); clf;
figure(2); orient tall; wysiwyg;
clf;
scl = 24*3600*sqrt(c*beta);
%global_axes(scl*3, scl*1.5, 1, 1, 1.5);

subplot('Position', [0.2794    0.6932    0.4412    0.1705]);
pp2 = plot(-scl*imag(lamu), scl*real(lamu), '.b');
hold on;
  pp1 = plot(-scl*imag(lamc), scl*real(lamc), 'or');
hold off;
%set(pp2, 'color', 0.7*[1 1 1]);
%set(pp1, 'color', 'k')

axis(scl*[-8 8 -.7 0.05]);
xl = xlabel('Phase (day^-^1)', 'fontsize', 10)
yl = ylabel('Damping (day^-^1)', 'fontsize', 10);
grid on;
t2 = title('\bfEigenspectrum:  Coupled (circles), Uncoupled (dots)', 'fontsize', 12);

if 0;
hold on;
ax = axes('position', [0.33 0.75 0.07 0.08]);
pp1 = plot(3600*24*c/ae*real(om), 3600*24*c/ae*imag(om), 'or');
set(pp1, 'color', 'k')
hold on;
  pp2 = plot(3600*24*c/ae*real(omu), 3600*24*c/ae*imag(omu), '.b');
hold off;
set(pp2, 'color', 0.7*[1 1 1]);
axis([-.25 .25 -.125 0.075])
grid on;
set(gca, 'YAxisLocation', 'right');
end;