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;
soipac = getflx('psl', [179 287 -20 20], 101:1000);
[tem, tem, tem, lmpac] = getll('temp', [179 287 -20 20]);
soiind = getflx('psl', [287 360 -20 20], 101:1000);
[tem, tem, tem, lmind] = getll('temp', [287 360 -20 20]);
soiind2 = getflx('psl', [-0.1 179 -20 20], 101:1000);
[tem, tem, tem, lmind2] = getll('temp', [-0.1 179 -20 20]);
for i = 1:900;
soipac(i,:,:) = squeeze(soipac(i,:,:)).*lmpac;
soiind(i,:,:) = squeeze(soiind(i,:,:)).*lmind;
soiind2(i,:,:) = squeeze(soiind2(i,:,:)).*lmind2;
end
soiind = shiftdim([shiftdim(soiind, 2); shiftdim(soiind2, 2)], 1);
soipac = mean2(mean2(shiftdim(soipac, 1)));
soiind = mean2(mean2(shiftdim(soiind, 1)));
soi = squeeze(soipac - soiind);
soi = detrend(soi);
ct = getnc('temp', [180 270 -6 6], 1, 101:1000);
ct = mean(mean(ct, 2), 3);
ct = detrend(ct);
taux = getnc('taux', [150 180 -6 6], 1, 101:1000);
taux = mean(mean(taux, 2), 3);
taux = detrend(taux);
ctcorr = [corr(ct, lppc(:,1)) corr(ct, hppc(:,1)) ...
corr(ct, rpc(:,1)) corr(ct, rspc(:,1)) ...
corr(ct, soi) corr(ct, taux)]
ctcorr =
0.6785 0.6306 0.9014 0.9587 -0.8926 0.7776
rawcor = [corr(lppc(:,1), rpc(:,1)) corr(hppc(:,1), rpc(:,1)) ...
corr(rpc(:,1), rspc(:,1)) corr(rpc(:,1), soi) ...
corr(rpc(:,1), taux) ...
corr(rspc(:,1), soi)]
rawcor =
0.8240 0.5772 0.9453 -0.8096 0.7201 -0.8637
[b1, a1] = butter(9, 2/10);
[b2, a2] = butter(9, 2/9);
lpcorr = [corr(filtfilt(b2, a2, ct), lppc(:,1)), ...
corr(filtfilt(b2, a2, rpc(:,1)), lppc(:,1)), ...
corr(filtfilt(b2, a2, taux), lppc(:,1)) ...
corr(filtfilt(b2, a2, soi), lppc(:,1))]
lpcorr =
0.8992 0.9850 0.7378 -0.8447
lpcorrp1 = [corr(filtfilt(b2, a2, ct), lppc(:,1), 1), ...
corr(filtfilt(b2, a2, rpc(:,1)), lppc(:,1), 1), ...
corr(filtfilt(b2, a2, taux), lppc(:,1), 1) ...
corr(filtfilt(b2, a2, soi), lppc(:,1), 1)]
lpcorrp1 =
0.7915 0.9142 0.6385 -0.7322
==> var2 leads var1
lpcorrm1 = [corr(filtfilt(b2, a2, ct), lppc(:,1), -1), ...
corr(filtfilt(b2, a2, rpc(:,1)), lppc(:,1), -1), ...
corr(filtfilt(b2, a2, taux), lppc(:,1), -1) ...
corr(filtfilt(b2, a2, soi), lppc(:,1), -1)]
lpcorrm1 =
0.9022 0.9504 0.7451 -0.8577
==> var1 leads var2
hpcorr = [corr(ct-filtfilt(b1, a1, ct), hppc(:,1)), ...
corr(rpc(:,1)-filtfilt(b1, a1, rpc(:,1)), hppc(:,1)), ...
corr(taux - filtfilt(b1, a1, taux), hppc(:,1)), ...
corr(soi - filtfilt(b1, a1, soi), hppc(:,1))]
hpcorr =
0.8865 0.9496 0.8153 -0.7984
ctreg = [ct' * lppc(:,1), ct' * hppc(:,1) ct'*rpc(:,1)]./900
ctreg =
0.1897 0.1763 0.2521
tauxcorr = [corr(taux, lppc(:,1)) corr(taux, hppc(:,1)) corr(taux, rpc(:,1)) ...
corr(taux, soi)]
tauxcorr =
0.4265 0.6982 0.7201 -0.8585
tauxcorr2 = [corr(filtfilt(b2, a2, taux), lppc(:,1)), ...
corr(taux - filtfilt(b1, a1, taux), hppc(:,1)), ...
corr(taux, rpc(:,1)) ...
corr(filtfilt(b2, a2, taux), filtfilt(b2, a2, soi)) ...
corr(taux - filtfilt(b1, a1, taux), soi - filtfilt(b1, a1, soi))]
tauxcorr2 =
0.7378 0.8153 0.7201 -0.9051 -0.8436
cttaux = [corr(ct, taux), ...
corr(filtfilt(b2, a2, ct), filtfilt(b2, a2, taux)), ...
corr(ct - filtfilt(b1, a1, ct), taux - filtfilt(b1, a1, taux))]
cttaux =
0.7776 0.8659 0.7503
ctsoi = [corr(ct, soi), ...
corr(filtfilt(b2, a2, ct), filtfilt(b2, a2, soi)), ...
corr(ct - filtfilt(b1, a1, ct), soi - filtfilt(b1, a1, soi))]
ctsoi =
-0.8926 -0.9496 -0.8466
pc_per = [rper(1:3)', lper(1:3)', hper(1:3)']
pc_per =
24.8994 28.6640 26.5249
13.5284 15.2824 9.7593
6.2328 6.3732 7.4138
pc_lam = [rlam(1:3)', llam(1:3)', hlam(1:3)']
pc_lam =
1.0e+19 *
1.2378 0.8621 0.5214
0.6725 0.4596 0.1919
0.3098 0.1917 0.1457
pc_lam_sum = [sum(rlam) sum(llam) sum(hlam)]
pc_lam_sum =
1.0e+19 *
4.9711 3.0075 1.9658