Documentation of walker_eof

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

Cross-Reference Information

Listing of script walker_eof

cd /home/disk/tao/data/nmc.reanalysis/monthly
filin = 'omega.mon.mean.nc';
nc = netcdf(filin, 'nowrite');
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  lev = nc{'level'}(:);
  [xk, yk] = keep_var([0 360 -7 7], lon, lat);
  omega = nc{'omega'}(:,:,yk,xk);
  add_offset = nc{'omega'}.add_offset(:);
  scale_factor = nc{'omega'}.scale_factor(:);
nc = close(nc);
omega = omega * scale_factor;
omega = omega + add_offset;
omega = shiftdim(squeeze(mean(shiftdim(omega, 2))), 1);
[ntim, nlev, nlon] = size(omega);
omega = reshape(omega, ntim, nlev*nlon);
[omega, clim] = annave(omega);
clim = reshape(clim, 12, nlev, nlon);
c = omega*omega';
[lam, pcs, per] = eof(c);
pc10 = pcs(:,1:10) ./ (ones(ntim, 1) * std(pcs(:,1:10)));
lds = pc10' * omega / ntim;
num = 1;
tem = reshape(lds(num,:), nlev, nlon);
figure(1)
sd(1)
     pncont(lon, -1*lev, tem, [-.02:.004:.03]);
     set(gca, 'YTickLabel', [1000:-100:100])
sd(2)
     plot(pc10(:,num));
     axis([0 ntim+1 -3 3])
     grid
     set(gca, 'XTick', 25:60:500, 'XTickLabel', 60:5:100);
tem = squeeze(mean(clim([1 2 12], :, :)));
figure(2)
pncont(lon, -1*lev, tem, [-.1:.02:.1])