Documentation of ml_coup_slp_eof

Global Index (short | long) | Local contents | Local Index (short | long)

Help text

[reg3, c3] = regress_eof(sstc, lpseudopcs(:,num), 0);
[reg3, c3] = regress_eof(sst, mlpcs(:,num), 0);

Cross-Reference Information

Listing of script ml_coup_slp_eof

clear
data
load ML_SLP_eof_npac.mat
back
mrpcs = rpcs; mlpcs = lpcs; mhpcs = hpcs;
mrper = rper; mlper = lper; mhper = hper;
mrlds = rlds; mllds = llds; mhlds = hlds;
mrlam = rlam; mllam = llam; mhlam = hlam;
sstc = getnc('temp', [-0.1 360 -90 90], 1, 101:1000);
slpc = getflx('psl', [-0.1 360 -90 90], 101:1000);
csirod
load ML_ANN_slp.mat
load ML_ANN_sst.mat; sst = sst2;
back
[lat, lon, depth, lm] = getll('temp', [-0.1 360 -90 90]);
[xk, yk] = keep_var(lims, lon, lat);
for i = 1:3;
  temslp = detrend(reshape(slpc(:,yk,xk), 900, length(xk)*length(yk)));
  if i == 1;
    rpseudopcs = temslp*mrlds;
  elseif i == 2;
    [b, a] = butter(9, 2/9);
    temslp = filtfilt(b, a, temslp);
    lpseudopcs = temslp*mllds;
  elseif i == 3;
    [b, a] = butter(9, 2/10);
    temslp = temslp - filtfilt(b, a, temslp);
    hpseudopcs = temslp*mhlds;
  end
end
data
load slp_eof_npac.mat
back
default_global; FRAME = [105 295 -60 60];
figure(1); fo(1); clf;
for num = 1:2; 
[reg1, c1] = regress_eof(sstc, rpseudopcs(:,num), 0);
[reg2, c2] = regress_eof(slpc, rpseudopcs(:,num), 0);
[reg1m, c1m] = regress_eof(sst, mrpcs(:,num), 0);
[reg2m, c2m] = regress_eof(slp, mrpcs(:,num), 0);

cint = [0.05 0.2]; pm = [1 -1];
for i = 1:2;
  eval(['tem1 = pm(num)*reg' num2str(i) ';']);
  eval(['tem2 = pm(num)*c' num2str(i) ';']);
  sptalk(4,2,(2*i-1+4*(num-1)));
  gcont(tem1, cint(i));
  dc2(lm, 0.5, pm(i)*10);
  color_shade(squeeze(tem2.^2), .1, .87*[1 1 1]);
  set(gca, 'fontsize', 9);
end

for i = 1:2;
  eval(['tem1 = pm(num)*reg' num2str(i) 'm;']);
  eval(['tem2 = pm(num)*c' num2str(i) 'm;']);
  sptalk(4,2,2*i+4*(num-1));
  gcont(tem1, cint(i));
  dc2(lm, 0.5, pm(i)*10);
  color_shade(squeeze(tem2.^2), .1, .87*[1 1 1]);
  set(gca, 'fontsize', 9);
end

end

sptalk(4,2,1);
title('COUP:  Regression on PPC1');
sptalk(4,2,2);
title('ML:  Regression on PC1');
sptalk(4,2,5);
title('COUP:  Regression on PPC2');
sptalk(4,2,6);
title('ML:  Regression on PC2');


cd ~/matlab/CSIRO/Thesis/Chap5/Plots