Global Index (short | long) | Local contents | Local Index (short | long)
| This script calls | |
|---|---|
clear
lims = [180 270 -6 6];
lev = 1:10;
tim = 101:1000;
ct = getnc('temp', lims, lev, tim);
ct = detrend(ct);
ct = squeeze(mean2(mean2(shiftdim(ct, 2))));
[b, a] = butter(9, 2/9);
ct = filtfilt(b, a, ct);
[lat, lon, depth] = getll('temp', lims);
taux = getnc('taux', [165 210 -6 6], 1, tim);
taux = detrend(taux);
taux = squeeze(mean2(mean2(shiftdim(taux, 1))));
taux = filtfilt(b, a, taux);
newlev = [12.5 100:100:500]';
ctind = ct(:, 1);
for i = 1:length(tim);
ct2(:, i) = interp1(depth, ct(i,:)', newlev);
end;
ct2 = ct2';
lags = -15:15;
for i = 1:length(lags);
for j = 1:length(newlev);
ct_corr(i, j) = corr(ct2(:, j), ctind, lags(i));
end
ct_corr(i, length(newlev)+1) = corr(taux, ctind, lags(i));
end
figure(1); fl;
subplot(1,1,1);
h1 = plot( ...
lags, ct_corr(:, 1), '-k', ...
lags, ct_corr(:, 2), '^-k', ...
lags, ct_corr(:, 3), 'o-k', ...
lags, ct_corr(:, 4), 's-k', ...
lags, ct_corr(:, 5), '+-k', ...
lags, ct_corr(:, 6), 'v-k', ...
lags, ct_corr(:, 7), 'x--k');
set(h1([1 7]), 'linewidth', 2);
axis([-10 10 -0.5 1.1]);
grid on
set(gca, 'XTick', -20:2:20, 'YTick', -1:.2:1);
title(['Lagged Correlation of LP9 Equatorial Temperature (CT region)' ...
' at Various Depths']);
xlabel('Lag (years)');
ylabel('Correlation Coefficient');
legend(h1, '12.5m Temperature', '100m Temperature', '200m Temperature', ...
'300m Temperature', '400m Temperature', '500m Temperature', ...
'Zonal Wind Stress');
cd /home/disk/tao/dvimont/Thesis/Chap5