Documentation of nao_ccm

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

Cross-Reference Information

Listing of script nao_ccm

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
filin = 'psl.nc';
nc = netcdf(filin, 'nowrite')
  psl = nc{'PSL'}(:);
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  date = nc{'DATE'}(:);
nc = close(nc);
latk = find(lat >= 20);
psn = psl(:, latk, :);
[ntim, nlat, nlon] = size(psn);
XAX = lon;
YAX = lat(latk);
psn = reshape(psn, ntim, nlat*nlon);
[psn, clim] = annave(psn);
psn = reshape(psn, ntim, nlat, nlon);
psn = cosweight(psn, lat(latk));
psn = reshape(psn, ntim, nlat*nlon);
c = psn*psn'/(nlat*nlon);
[lam, pc, per] = eof(c);
pc = ones(ntim,1) * (1./(std(pc))) .* pc;
lds = pc(:,1:10)' * psn / ntim;
lds = reshape(lds, 10, nlat, nlon);
num = 1;
tem = squeeze(lds(num,:,:));
mcont(tem, [-1000:50:1000],'stereo')