Documentation of plot_wbtp_tcline_reg

Global Index (short | long) | Local contents | Local Index (short | long)

Help text

  Plot wbT'_z

Cross-Reference Information

Listing of script plot_wbtp_tcline_reg

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP9_detrend_L1-7_EOF_yr101-1000.mat
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP9_hctend_ALL.mat
newlag = [-6:2:6];
newlag = [-4 -1 2];
wbtp = wbtp13(find(ismember(lags, newlag)), :, :);
lims2 = [lims(1:2) -3 3];
temp = getnc('temp', lims2, 1:10, tim);
temp2 = squeeze(shiftdim(mean2(shiftdim(temp, 2)), 2));
tc = regress_eof(temp2, pcs, newlag);

figure(1); clf;

[lat, lon] = getll('wl', lims);
[lat2, lon2, depth2] = getll('temp', lims);
lm = getnc('temp', lims, 1, 1);

default_global;  FRAME = [110 300 -30 30];
cint = 0.25; clev = [-10:cint:-cint cint:cint:10];

for i = 1:7;
subplot(7,2,2*i-1);
  XAX = lon; YAX = lat;
  gcont(1*wbtp(i,:,:), cint);
  ylabel(['Lag = ' num2str(newlag(i))]);
  XAX = lon2; YAX = lat2;
    dc2(lm);
end

subplot(7,2,1)
xtic = get(gca, 'XTick');
xtlab = get(gca, 'XTickLabel');

cint2 = 0.025; clev2 = [-5:cint2:-cint2 0 cint2:cint2:5];
for i = 1:3;
subplot(6,2,2*i-1);
  pncont(lon2, -1*depth2, -1*squeeze(tc(i,:,:)), clev2, 'k');
  dcv(lon2, -1*depth2, squeeze(tc(i,:,:)));
  axis([110 300 -500 0]);
  set(gca, 'YTick', [-500:100:0], 'YTickLabel', [500:-100:0]);
  set(gca, 'XTick', xtic, 'XTickLabel', xtlab);
  ylabel(['Lag = ' num2str(newlag(i))]);
end
subplot(6,2,1);
title('Equatorial Temperature');

subplot(7,2,1);
  title('HP10:  0-80m Wbar T''_z');
subplot(7,2,2);
  title('HP10:  Equatorial T''');

subplot(7,2,13);
  xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2 std^-^1']);
subplot(7,2,14);
  xlabel(['Contour Interval:  ' num2str(cint2) ' K std^-^1']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap3/Plots
% print -dps2 wbtp_tcline_reg.ps