Global Index (short | long) | Local contents | Local Index (short | long)
This script calls | |
---|---|
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])