Documentation of precip_ct_nmc

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

Help text

filin = 'prate.mon.mean.nc';
var_nam = 'prate';

Cross-Reference Information

Listing of script precip_ct_nmc

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
load ct_gr.mat
cd /home/disk/tao/data/nmc.reanalysis/monthly
filin = 'vwnd.mon.mean.nc';
var_nam = 'vwnd';
lev = 1000;

nc = netcdf(filin, 'nowrite');
  levs = nc{'level'}(:);

  lind = find(levs == lev);

  var = squeeze(nc{var_nam}(37:444,lind,:,:));
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  add_offset = nc{var_nam}.add_offset(:);
  scale_factor = nc{var_nam}.scale_factor(:);
nc = close(nc);

var = var * scale_factor;
var = var + add_offset;
get_global;

kp = sort([1:12:408 2:12:408 12:12:408]);
var = var(kp,:,:);
[ntimw, nlat, nlon] = size(var);

i = 0;
clear clim
for m = 1:3;
  i = i+1;
  clim(i,:,:) = mean(var(m:3:ntimw,:,:));
end
x = zeros(ntimw,nlat,nlon);
for m = 1:ntimw;
  l=rem(m-1,3)+1;
  var(m,:,:) = var(m,:,:) - clim(l,:,:);
end;

%temct = ctstar(kp);
temct = -1*gr(kp);
temct = (temct-mean(temct))/std(temct);
pat = reshape((temct' * reshape(var, ntimw, nlat*nlon)/ntimw),nlat,nlon);

lims = [0 360 -90 90];

%top = squeeze(mean(clim)) * (3600 * 24 * 1000 / 1e4);
%bot = pat * (3600 * 24 * 1000 / 1e4);
top = squeeze(mean(clim));
bot = pat;

figure(2)
subplot(2,1,1)
     gcont(top, [-20:2:20]);
     dc
subplot(2,1,2)
     gcont(bot, [-5:.5:5]);
     dc