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

pn = 1;
x = (0:Lx/200:Lx)*(pi*RADUS/180)/ae;
Lxn = Lx*pi*RADUS/(180*ae);
fx = exp(ij*k*x);
%fx = exp(ij*[0:1/200:1]*2*pi);

u = -0.5*pn*vc(1:4:end, ind(ii))*fx*c;
v = -0.5*pn*vc(2:4:end, ind(ii))*fx*c;
h = -0.5*pn*vc(3:4:end, ind(ii))*fx*He;
T =      pn*vc(4:4:end, ind(ii))*fx*To;
%T = T .* ((Kq'/max(Kq)) * ones(1,length(x)));

mT = max(max(real(T)));
mh = max(max(real(h)));
lev = [-0.9:.3:-0.3 0.3:.3:0.9];

%figure_tall(3); clf;
figure(1); orient tall; wysiwyg;
scl = 1.25;
%global_axes(scl*3, scl*1.5, 1, 1, 1.5);

subplot('Position', [0.2794    0.6932    0.4412    0.1705]);
%subplot2(1,1); cla;
  [c2a, h2a] = contour(x, y, real(T), lev(lev<0)*mT, '--r');
  hold on;
  [c2b, h2b] = contour(x, y, real(T), lev(lev>0)*mT, '-r');
  [c1a, h1a] = contour(x, y, real(h), lev(lev<0)*mh, '--k');
  [c1b, h1b] = contour(x, y, real(h), lev(lev>0)*mh, '-k');
  hold off;
  set([h2a; h2b], 'linewidth', 2, 'color', 'r');
  set([h1a; h1b], 'linewidth', 2);
  
  set(gca, 'XTick', 0:Lxn/4:Lxn, 'XTickLabel', []);
  set(gca, 'YTick', 3*[-1:.5:1], 'YTickLabel', []);
  
  
  nthin = [ceil((1/dy)/5), 10];
  temu = thin(real(u), nthin);
  temv = thin(real(v), nthin);
  global_latlon(y, x); global XAX2 YAX2
  temx = thin(XAX2, nthin);
  temy = thin(YAX2, nthin);
  hold on;
  q = quiver(temx, temy, temu, temv, 1.5);
  hold off;
  set(q, 'color', 'k');
  
  axis([0 Lxn -3 3]);
  
set(gca, 'XTickLabel', [0:.25:1]);
set(gca, 'YTickLabel', 3*[-1:.5:1]);
xl = xlabel('X / Lx', 'fontsize', 10);
yl = ylabel('Y / ae', 'fontsize', 10);
lam1 = real(lamc(ind(ii))); lam2 = -1*imag(lamc(ind(ii)));
lam = lam1 + sqrt(-1)*lam2;
lamtxt = num2str(round(1000*3600*24*100*c/ae*lam)/1000);
t1 = title(['{\bfSWM Mode ' num2str(ii) ':}  \lambda = ' lamtxt ...
	     ' (100 days)^-^1']);
set(t1, 'fontsize', 12);

if 1;
%global_axes(scl*.7, scl*1.5, scl*3+.1, 1, 1.5);
%subplot2(2,2); cla;
subplot('Position', [0.7265    0.6932    0.1029    0.1705]); cla;
pp2 = plot(wes*rho*co*Ho*To/ae, y, 'b', ...
     Kq*(max(wes*rho*co*Ho*To/ae)/max(Kq)), y, 'r');
axis([[-1.15 1.15]*max(abs(wes))*rho*co*Ho*To/ae [-1 1]*3]);
set(pp2, 'linewidth', 2);
set(gca, 'YTick', -3:(3/2):3, 'YTickLabel', []);
grid on;
t2 = title2('{\bfWES}, Kq', 'fontsize', 12);
xl = xlabel('(W s m^-^3)', 'fontsize', 10);
set(pp2(1), 'color', 0*[1 1 1], 'linewidth', 2);
set(pp2(2), 'color', [0 0 0], 'linewidth', 1);
%ll = legend(pp2, 'WES', 'Kq');
%set(ll, 'fontsize', 10);
%delete(pp2(2));
end