Global Index (short | long) | Local contents | Local Index (short | long)
| This script calls | |
|---|---|
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load RAW_detrend_L1_EOF_yr101-1000.mat; rspc = 1*pcs; rslam = lam;
load RAW_detrend_L1-7_EOF_yr101-1000.mat; rpc = -1*pcs; rlam = lam;
load LP9_detrend_L1-7_EOF_yr101-1000.mat; lppc = -1*pcs; llam = lam;
load HP10_detrend_L1-7_EOF_yr101-1000.mat; hppc = 1*pcs; hlam = lam;
ct = getnc('temp', [165 205 17.5 37.5], 1, 101:1000);
ct = mean(mean(ct, 2), 3);
ct = detrend(ct);
ct2 = getnc('temp', [180 270 -6 6], 1, 101:1000);
ct2 = mean(mean(ct2, 2), 3);
ct2 = detrend(ct2);
ctcorr = [corr(rpc(:,1), ct) corr(rspc(:,1), ct) corr(ct, ct2)]
ctcorr =
-0.5108 -0.4907 -0.4259
[b1, a1] = butter(9, 2/10);
[b2, a2] = butter(9, 2/9);
ctlp = [corr(lppc(:,1), filtfilt(b2, a2, ct)) ...
corr(filtfilt(b2, a2, ct), filtfilt(b2, a2, ct2))]
ctlp =
-0.5735 -0.4173
ctlp = [corr(lppc(:,1), filtfilt(b2, a2, ct), 1) ...
corr(filtfilt(b2, a2, ct), filtfilt(b2, a2, ct2), -1)]
ctlp =
-0.5969 -0.4177
==> ct leads lppc by one year
ct leads ct by one year -- note small difference in correlation
cthp = [corr(hppc(:,1), ct - filtfilt(b1, a1, ct)) ...
corr(ct - filtfilt(b1, a1, ct), ct2 - filtfilt(b1, a1, ct2))]
cthp =
-0.4806 -0.4373