Global Index (short | long) | Local contents | Local Index (short | long)
[reg1, c1] = regress_eof(sst, tim, 0);
| This script calls | |
|---|---|
clear
csirod;
load ML_ANN_sst.mat; sst = sst2;
load ML_ANN_slp.mat;
load prec_ML_annave.mat;
back
lims = [100 305 -65 65];
[xk, yk] = keep_var(lims, lon, lat);
sst = sst(:, yk, xk); slp = slp(:, yk, xk); prec = prec(:, yk, xk);
[lat, lon, depth, lm] = getll('temp', lims);
cd ~/matlab/CSIRO/Thesis/Data; load ML_SLP_eof_npac.mat; back;
cd ~/matlab/CSIRO/Thesis/Data; load ML_SLP_eof_spac.mat; back;
default_global; FRAME = [105 300 -60 60];
num = 1;
figure(1); fo(1); clf;
for i = 1:3;
if i == 1;
tim = rpcs(:,num);
var1 = slp;
var2 = prec;
elseif i == 2;
tim = hpcs(:,num);
[b, a] = butter(9, 2/10);
var1 = slp - filtfilt(b, a, slp);
var2 = prec - filtfilt(b, a, prec);
elseif i == 3;
tim = lpcs(:,num);
[b, a] = butter(9, 2/9);
var1 = filtfilt(b, a, slp);
var2 = filtfilt(b, a, prec);
end
[reg1, c1] = regress_eof(var1, tim, 0);
[reg2, c2] = regress_eof(var2, tim, 0);
% pos = [.2 (9.8/11)-i*7/33 .3 .18];
% subplot('position', pos);
sptalk(3,2,2*i-1);
gcont(reg1, .1);
dc2(lm, .5, 1);
color_shade(squeeze(c1.^2), .1, .87*[1 1 1]);
if i == 1; title('SST'); end;
% pos = [.58 (9.8/11)-i*7/33 .3 .18];
% subplot('position', pos);
sptalk(3,2,2*i);
gcont(reg2, .05);
dc2(lm, .5, -1);
color_shade(squeeze(c2.^2), .1, .87*[1 1 1]);
if i == 1; title('SLP'); end;
end
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap5/Plots
rtem = rpcs;
ltem = lpcs;
htem = hpcs;
sptalk(4,1,i)