Global Index (short | long) | Local contents | Local Index (short | long)
temp2 = squeeze(shiftdim(mean2(shiftdim(temp, 2)), 2));
This script calls | |
---|---|
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