Global Index (short | long) | Local contents | Local Index (short | long)
tsnew = tsnew(:,keep);
This script calls | |
---|---|
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)