Global Index (short | long) | Local contents | Local Index (short | long)
tim = 551:1000; tim = 101:1000;
| This script calls | |
|---|---|
clear
lev = 1;
tim = 101:550;
vortlims = [165 190 37.5 50];
%vortlims = [150 180 37.5 52.5];
[u, v] = getnc('taux', 'tauy', vortlims, lev, tim);
[latu, lonu] = getll('taux', vortlims);
ctlim = [180 270 -6 6];
[temp] = getnc('temp', ctlim, lev, tim);
wstplim = [190 205 30 37.5];
heat = getnc('temp', wstplim, 1, tim);
[ntim, nlat, nlon] = size(u);
clear vor;
for i = 1:ntim;
vor(i,:,:) = sph_curl1(u(i,:,:), v(i,:,:), latu, lonu, 0);
end;
vor = squeeze(mean(mean(shiftdim(vor, 1))));
ct = squeeze(mean(mean(shiftdim(temp, 1))));
wstp = squeeze(mean(mean(shiftdim(heat, 1))));
vor = detrend(vor);
ct = detrend(ct);
wstp = detrend(wstp);
% Lowpass filter
[b, a] = butter(6, 2/10);
vorlp = filtfilt(b, a, vor);
ctlp = filtfilt(b, a, ct);
wstplp = filtfilt(b, a, wstp);
lags = -25:25;
for i = 1:length(lags);
cv(i) = corr(ctlp, vorlp, lags(i));
cw(i) = corr(ctlp, wstplp, lags(i));
cc(i) = corr(ctlp, ctlp, lags(i));
end
figure(3); figure_orient
subplot(3,1,1);
bar(lags, cc);
axis([-26 26 -.6 1.1])
set(gca, 'YTick', -1:.25:1);
grid on
title('Lagged Autocorrelation for CT');
xlabel('Lag')
subplot(3,1,2);
bar(lags, cv);
axis([-26 26 -.6 .6])
set(gca, 'YTick', -1:.25:1);
grid on
title('Lagged Correlation between CT and CURL(TAU) [165W to 165E, 37.5N to 50N]');
xlabel('Lag: Lag < 0 ==> CT leads; Lag > 0 ==> CURL(TAU) leads');
subplot(3,1,3);
bar(lags, cw);
axis([-26 26 -.6 .6])
set(gca, 'YTick', -1:.25:1);
grid on
title('Lagged Correlation between CT and NP [190W to 205E, 30N to 37.5N]');
xlabel('Lag: Lag < 0 ==> CT leads; Lag > 0 ==> NP leads');
% Look at western boundary current
lev = 1:3;
wbclims = [105 155 15 35];
[wbc, lat, lon] = getheat(3:5, tim, wbclims, 'v');
[ntim, nlat, nlon] = size(wbc);
for tind = 1:ntim;
for yind = 1:nlat;
wbc2(tind,yind) = wbc(tind,yind,min(find(~isnan(wbc(tind,yind,:)))));
end
end
wbc = squeeze(mean(shiftdim(wbc2, 1)));
wbc = detrend(wbc);
wbc = filtfilt(b, a, wbc);