Documentation of som_setup

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

Help text

tsnew = tsnew(:,keep);

Cross-Reference Information

Listing of script som_setup

clear
cd /home/disk/hayes2/dvimont/ccm/som/data
[TS, ORO] = ncread;
whos
nya = (-87.2:(2*87.2/47):87.2);
nxa = [0:7.5:352.5];
lims = [120 290 -60 60];
oronew = subset(ORO, nxa, nya, lims);
tsnew = subset(TS, nxa, nya, lims);
[ntstp,m,n]=size(oronew);
oronew = reshape(oronew,768,m*n);
tsnew = reshape(tsnew,768,m*n);
keep = find(mean(oronew)==0);
keepnya = find(nya>=lims(3) & nya <= lims(4));
keepnxa = find(nxa>=lims(1) & nxa <=lims(2));
[tsa, clim] = annave(tsnew);
tsa = cosweight(tsa, nya(keepnya));
tsa = tsa(:,keep);
c = tsa(180:768,:)' * tsa(180:768,:);
%c = tsa' * tsa;
[lam, loads, per] = eof(c);
pc = loads' * tsa(180:768,:)';
whos

temp = zeros(10,m*n);
temp(:,keep) = loads(:,1:10)';
gcont(reshape(temp(1,:),m,n),lims)



temp = reshape(tsnew,768,m,n);
tempa = reshape(tsa,768,m,n);
tempc = reshape(clim,12,m,n);
gcont(squeeze(TS(1,:,:)))
gcont(squeeze(temp(400,:,:)), lims)
gcont(squeeze(tempa(400,:,:)), lims)
gcont(squeeze(clim(1,:,:), lims)