Documentation of sst_ceof

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

Help text

  data; save ML_NH_ceof.mat lam lds pcs per lims; back
  data; save COUP_NH_ceof.mat lam lds pcs per lims; back

Cross-Reference Information

Listing of script sst_ceof

clear
csirod; load ML_ANN_sst.mat; back;
lims = [110 270 15 67.5];
[xk, yk] = keep_var(lims, lon, lat);
sst = sst2(:, yk, xk);
sst = getnc('temp', lims, 1, 101:1000);
[lat2, lon2, depth, lm] = getll('temp', lims);
lm(~isnan(lm)) = 1;
[ntim, nlat, nlon] = size(sst);
for i = 1:ntim;
  sst(i,:,:) = squeeze(sst(i,:,:)) .* lm;
end
clim = squeeze(mean(sst));
sst = detrend(sst);
sst = cosweight(sst, lat2);
kp = find(~isnan(clim));
sst = sst(:, kp);
[lam, lds, pcs, per] = complex_eof(sst);

%  Look at patterns

csirod; load ML_ANN_sst.mat; back;
data; load ML_NH_ceof.mat; back;

sst2 = getnc('temp', [-0.1 360 -90 90], 1, 101:1000);
data; load COUP_NH_ceof.mat; back;

num = 2;
[reg1, temtim] = regress_ceof(sst2, pcs(:,num), 6);
[lat, lon, depth, lm] = getll('temp', [-0.1 360 -90 90]);

default_global; FRAME = [0 360 -90 90];

figure(1); fo(1);
sptalk(2,1,1);
  gcont(reg1(1,:,:), 0.1);
  dc2(lm, 0.55, -1);
sptalk(2,1,2);
  gcont(reg1(2,:,:), 0.1);
  dc2(lm, 0.55, -1);

figure(2); fo(1);
sptalk(2,1,1);
  gcont(reg1(3,:,:), 0.1);
  dc2(lm, 0.55, -1);
sptalk(2,1,2);
  gcont(reg1(4,:,:), 0.1);
  dc2(lm, 0.55, -1);

figure(3); fo(1);
sptalk(2,1,1);
  gcont(reg1(5,:,:), 0.1);
  dc2(lm, 0.55, -1);
sptalk(2,1,2);
  gcont(reg1(6,:,:), 0.1);
  dc2(lm, 0.55, -1);