Documentation of plot_wbtp_tcline_reg


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


Help text

temp2 = squeeze(shiftdim(mean2(shiftdim(temp, 2)), 2));

Cross-Reference Information

This script calls

Listing of script plot_wbtp_tcline_reg


clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load HP10_detrend_L1-7_EOF_yr101-1000.mat
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load HP10_hctend_ALL.mat
uptb47=uptb17-uptb13;
varn = 'wl';
lims2 = [lims(1:2) -3 3];
temp = getnc(varn, lims2, 1:9, tim);
temp2 = temp;

tc = regress_eof(temp2, pcs, nfrm);

tc = squeeze(mean(tc, 3));

%  Plot wbT'_z

[lat, lon] = getll('wl', lims);
[lat2, lon2, depth2] = getll(varn, 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:5;
subplot(7,2,2*i-1);
  XAX = lon2; YAX = lat2;
  gcont(-1*wptb13(i+8,:,:), clev);
  ylabel(['Lag = ' num2str(lags(i+8))]);
  XAX = lon; YAX = lat;
    dc2(lm);
end

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

cint2 = 0.01; clev2 = [-.5:cint2:-cint2 0 cint2:cint2:.5];
for i = 1:5;
subplot(7,2,2*i);
  pncont(lon2, -1*depth2, squeeze(tc(i+1,:,:)), clev2, 'k');
  dcv(lon2, -1*depth2, squeeze(tc(i+1,:,:)));
  axis([110 300 -270 0]);
  set(gca, 'YTick', [-500:100:0], 'YTickLabel', [500:-100:0]);
  set(gca, 'XTick', xtic, 'XTickLabel', xtlab);
end

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

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

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