Global Index (short | long) | Local contents | Local Index (short | long)
title(['Lagged correlations between 10yr LP filtered,'... ' detrended equatorial indices, and the CT index']);
This script calls | |
---|---|
clear cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data load LP10_detrend_L1-7_yr101-1000.mat; pcs = -1*pcs; tim = 101:1000; ct = getnc('temp', [180 270 -6 6], 1, tim); ct = squeeze(mean(mean(ct, 2), 3)); [b, a] = butter(6, 2/10); ct = filtfilt(b, a, detrend(ct)); taux = getnc('taux', [140 210 -6 6], 1, tim); taux = squeeze(mean2(mean2(shiftdim(taux, 1)))); [b, a] = butter(6, 2/10); taux = filtfilt(b, a, detrend(taux)); hc = getheat([180 270 -6 6], 4:7, tim); hc = squeeze(mean(mean(hc, 2), 3)); [b, a] = butter(6, 2/10); hc = filtfilt(b, a, detrend(hc)); pcs = ct; lags = -15:15; nlag = length(lags); for i = 1:nlag; ctl(i) = corr(ct, pcs(:,1), lags(i)); txl(i) = corr(taux, pcs(:,1), lags(i)); hcl(i) = corr(hc, pcs(:,1), lags(i)); end figure(2); figure_orient(1); subplot(2,1,1); h1 = plot( ... lags, hcl, 'o-k', ... lags, txl, 'x-k', ... lags, ctl, '-k'); grid on set(h1(3), 'linewidth', 2); set(gca, 'XTick', -20:2:20, 'YTick', -1:.2:1); axis([-15 15 -.5 1.1]); xlabel('Lag in Years (negative lag implies the index leads LP10 PC1)'); ylabel('Correlation Coefficient'); box on legend(h1, '80-270m Heat Content', 'Zonal Wind Stress', 'CT 12.5m Temperature'); cd /home/disk/tao/dvimont/matlab/CSIRO/IUGG_talk