Documentation of z500_setup

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

Help text

pc = pc * diag(1./sqrt(lam));

Cross-Reference Information

Listing of script z500_setup

cd /home/disk/hayes2/dvimont/csm/300yr/tz/
load yrs0-99.mat
z1 = z;
load yrs100-159.mat
whos
z = [z1(13:1200,:,:); z(1:612,:,:)];
clear z1
whos
[ntmz, nyz, nxz] = size(z);
z = reshape(z,ntmz,nyz*nxz);
[z, climz] = annave(z);
z3yr = myrunning_ave(z, 36);
z3yr = z(1:3:ntmz,:);
regrsz = pc1(:,1:10)' * z3yr / 600;
tempz = reshape(regrsz,10,nyz,nxz);
regrsz = regrsz' - ones(nxz*nyz,1) * mean(regrsz');
var = diag(regrsz' * regrsz / (nxz*nyz));
regrsz = regrsz * diag(1./sqrt(var));
regts = z3yr * regrsz;
regts = myrunning_ave(regts, 12);
regts = regts * diag(1./std(regts));


gcont(squeeze(tempz(1,:,:)),lims,[-500:50:500]);
dc