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