Global Index (short | long) | Local contents | Local Index (short | long)
Load PCS:
| This script calls | |
|---|---|
clear
%cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
%load LP10_detrend_L1-7_yr101-1000.mat
%pcs = -1*pcs;
tim = 101:1000;
lags = -30:30; nlag = length(lags);
lims = [180 270 -6 6];
ct = getnc('temp', lims, 1, tim);
ct = squeeze(mean(mean(ct, 2), 3));
lims1 = [160 180 30 37.5];
lims2 = [170 190 45 52.5];
tx1 = getnc('taux', lims1, 1, tim);
tx2 = getnc('taux', lims2, 1, tim);
txind = squeeze(mean(mean((tx1 - tx2), 2), 3));
tx = detrend(txind); pcs = detrend(ct);
[b1, a1] = butter(9, 2/9);
[b2, a2] = butter(9, 2/15);
tx = filtfilt(b1, a1, tx)-filtfilt(b2, a2, tx);
pcs = filtfilt(b1, a1, pcs)-filtfilt(b2, a2, pcs);
clear ac r;
for i = 1:nlag;
ac(i) = corr(pcs(:,1), pcs(:,1), lags(i));
r(i) = corr(tx, pcs(:,1), lags(i));
end
figure(1); figure_orient;
subplot(3,1,1);
bar(lags, ac);
xlabel('Lag in Years');
title(['Autocorrelation function for the CT index']);
axis([-15 15 -1.1 1.1]);
set(gca, 'YTick', -1:.2:1, 'XTick', -20:2:20);
grid on
ylabel('Correlation Coefficient')
r_exp = r((nlag+1)/2).*ac;
fishz = 0.5 * log((1+r)./(1-r));
fishmu = 0.5 * log((1+r_exp)./(1-r_exp));
dof = dof_corr(tx, pcs(:,1));
dof = (dof/2) - 3;
sigz = 1 ./ sqrt(dof);
subplot(3,1,3);
bar(lags, fishz);
hold on;
plot(lags, fishmu, '-k', ...
lags, fishmu+1.96*sigz, '--k', ...
lags, fishmu-1.96*sigz, '--k');
hold off
axis([-15 15 -.6 .6]);
set(gca, 'YTick', -2:.2:2, 'XTick', -20:2:20);
grid on
title('Fisher''s Z-Statistic for < CT, NPTX >, with 95% Confidence Intervals');
xlabel('Lag in years (negative lag implies NPTX leads CT)');
ylabel('Correlation Coefficient')
% Now do the same for the WSTP
tx = getheat([125 150 7.5 17.5], 4:7, tim);
tx = squeeze(mean(mean(tx, 2), 3));
[b, a] = butter(6, 2/10);
tx = filtfilt(b, a, detrend(tx));
pcs = filtfilt(b, a, detrend(pcs));
clear ac r;
for i = 1:nlag;
ac(i) = corr(pcs(:,1), pcs(:,1), lags(i));
r(i) = corr(tx, pcs(:,1), lags(i));
end
r_exp = r((nlag+1)/2).*ac;
fishz = 0.5 * log((1+r)./(1-r));
fishmu = 0.5 * log((1+r_exp)./(1-r_exp));
dof = dof_corr(tx, pcs(:,1));
dof = (dof/2) - 3;
sigz = 1 ./ sqrt(dof);
subplot(3,1,2);
bar(lags, fishz);
hold on;
plot(lags, fishmu, '-k', ...
lags, fishmu+1.96*sigz, '--k', ...
lags, fishmu-1.96*sigz, '--k');
hold off
axis([-15 15 -1.2 .6]);
set(gca, 'YTick', -2:.2:2, 'XTick', -20:2:20);
grid on
title('Fisher''s Z-Statistic for < CT, WSTP >, with 95% Confidence Intervals');
xlabel('Lag in years (negative lag implies WSTP leads CT)');
ylabel('Correlation Coefficient')
colormap([.65 .65 .65]);
lims = [135 200 35 52.5]
[tx, ty] = getnc('taux', 'tauy', lims, 1, tim);
[lat, lon] = getll('taux', lims);
rtx = regress_eof(tx, pcs, lags);
rty = regress_eof(ty, pcs, lags);
clear curlt;
for i = 1:length(lags);
tx = squeeze(rtx(i,:,:)); ty = squeeze(rty(i,:,:));
[curlt(i,:,:), lat2, lon2] = sph_curl1(tx, ty, lat, lon, 0);
end