Documentation of eof_sst

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

Help text

lims = [-0.1 360 -90 90];

Cross-Reference Information

Listing of script eof_sst

clear
csirod
load ML_ANN_sst.mat; sst = sst2; clear sst2;
back
lims = [105 295 -60 60];
[xk, yk] = keep_var(lims, lon, lat);
sst = sst(:,yk,xk);
[lat, lon, 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

sst2 = detrend(sst);

if 0;
  sst2 = 0.25*(sst2(:,1:2:nlat,1:2:nlon) + sst2(:,2:2:nlat,1:2:nlon) + ...
               sst2(:,1:2:nlat,2:2:nlon) + sst2(:,2:2:nlat,2:2:nlon));
  lat2 = 0.5*(lat(1:2:nlat)+lat(2:2:nlat));
else
  lat2 = lat;
end

clim = squeeze(mean(sst2));
kp = find(~isnan(clim));
sst2 = sst2(:, kp);

stdsst = std(sst2);
sst2 = sst2 ./ (ones(ntim,1)*stdsst);

tem = NaN*ones(ntim, nlat*nlon);
tem(:,kp) = sst2;
tem = cosweight(tem, lat2);
sst2 = tem(:,kp);

[lam, lds, pcs, per] = eof_dan(sst2);

%lamc = lam; ldsc = lds; pcsc = pcs; perc = per;
%lamr = lam; ldsr = lds; pcsr = pcs; perr = per;

data;
save ML_SST_pac_corr.mat lamc lamr ldsc ldsr pcsc pcsr perc perr ...
                         lat lon lims kp 
back

%  Look at results

%clear

%data; load ML_SST_pac_corr.mat; back;
csirod
load ML_ANN_sst.mat; sst = sst2; clear sst2;
back

lims2 = [-0.1 360 -90 90];
[tem, tem, tem, lm] = getll('temp', lims2);

sst2 = detrend(sst);
[reg1, c1] = regress_eof(sst2, -pcsc(:,1), 0);
[reg2, c2] = regress_eof(sst2, pcsc(:,2), 0);
[reg3, c3] = regress_eof(sst2, pcsc(:,3), 0);
[reg4, c4] = regress_eof(sst2, pcsc(:,4), 0);

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

figure(2); fo(1); clf;
sptalk(2,1,1);
gcont(c1, .1);
dc2(lm);
sptalk(2,1,2);
gcont(c2, .1);
dc2(lm);

[regr1, cr1] = regress_eof(sst2, -pcsr(:,1), 0);
[regr2, cr2] = regress_eof(sst2, -pcsr(:,2), 0);
[regr3, cr3] = regress_eof(sst2, pcsr(:,3), 0);
[regr4, cr4] = regress_eof(sst2, pcsr(:,4), 0);

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

figure(2); fo(1); clf;
sptalk(2,1,1);
gcont(regr1, .05);
dc2(lm);
sptalk(2,1,2);
gcont(regr2, .05);
dc2(lm);

figure(3); fo(1); clf;
sptalk(2,1,1);
gcont(regr3, .05);
dc2(lm);
sptalk(2,1,2);
gcont(regr4, .05);
dc2(lm);

%  Determine whether the NH and SH components of pcsc1 are related