Global Index (short | long) | Local contents | Local Index (short | long)
| This script calls | |
|---|---|
clear
cd /home/disk/hayes2/dvimont/data
load ct_gr.mat
load nmc_logistics.mat
load nmc_noice_sst.mat
ctemp = ctstar;
ctemp = (ctemp - mean(ctemp))/std(ctemp);
a1 = (1/(ctemp'*ctemp))*ctemp'*ts;
resid = ts - ctemp*a1;
ctemp = grhp;
ctemp = (ctemp - mean(ctemp))/std(ctemp);
a1 = (1/(ctemp'*ctemp))*ctemp'*ts;
resid = resid - ctemp*a1;
[ntim, nkeep] = size(ts);
tsnew = zeros(ntim, nlat*nlon);
tsnew(:,kpt) = resid;
tsnew = reshape(tsnew, ntim, nlat, nlon);
tsnew = cosweight(tsnew, lat);
tsnew = reshape(tsnew, ntim, nlat*nlon);
resid = tsnew(:,kpt);
clear tsnew
c = (resid*resid')/nkeep;
[lam, pc, per] = eof(c);
pc10 = ones(ntim,1) * (1./(std(pc(:,1:10)))) .* pc(:,1:10);
lds = pc10'*resid/ntim;
lds = pc10'*ts/ntim;
get_global
default_global
num = 1;
tem = NaN * ones(1, nlat*nlon);
tem(kpt) = lds(num,:);
tem = reshape(tem, nlat, nlon);
sd(1)
gcont(1.*tem, [-4:.1:4]);
dc
title('GR -- Covariance Matrix (Latitudes between 60S and 60N)')
xlabel('Contour Interval -- 0.1 K/std')
sd(2)
plot(1.*pc10(:,num))
title('GR')
xlabel('Year')
ylabel('STD')
set(gca, 'XTick',25:60:ntim,'XTickLabel',[1960:5:2000])
axis([1 408 -2.5 2.5])
grid
grlp2 = myrunning_ave(myrunning_ave(pc10(:,1),25),37);
grlp2 = (grlp2 - mean(grlp2))/std(grlp2);
grlpat2 = grlp2'*resid/ntim;
grlpat2 = grlp2'*ts/ntim;
tem = NaN * ones(1, nlat*nlon);
tem(kpt) = grlpat2;
tem = reshape(tem, nlat, nlon);
sd(1)
gcont(1.*tem, [-4:.1:4]);
dc
title('GR -- Covariance Matrix (Latitudes between 60S and 60N)')
xlabel('Contour Interval -- 0.1 K/std')
sd(2)
plot(1.*grlp2)
title('GR')
xlabel('Year')
ylabel('STD')
set(gca, 'XTick',25:60:ntim,'XTickLabel',[1960:5:2000])
axis([1 408 -2.5 2.5])
grid