Documentation of ct_hgt_nmc

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

Function Synopsis

pat = ct_hgt_nmc(lev);

Help text

ctstar = ctstar(1:396);

Cross-Reference Information

Listing of function ct_hgt_nmc

function pat = ct_hgt_nmc(lev);
cur_dir = cd;
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
load ct_gr.mat

cd /home/disk/tao/data/nmc.reanalysis/monthly

filin = 'hgt.mon.mean.nc'
nc = netcdf(filin, 'nowrite');
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  levnmc = nc{'level'}(:);
  lind = find(levnmc == lev);
  if isempty(lind)
    error(['Value for lev = ' num2str(lev) ' is not found'])
  end
  hgt500 = squeeze(nc{'hgt'}(37:444,lind,:,:));
%  hgt500 = squeeze(nc{'hgt'}(37:432,lind,:,:));
nc = close(nc);

kp = sort([1:12:408 2:12:408 12:12:408]);
%kp = sort([1:12:396 2:12:396 12:12:396]);
hgt500 = hgt500(kp,:,:);
ctstar = myrunning_ave(ctstar,1);
ctstar = ctstar(kp);
ctstar = (ctstar-mean(ctstar))/std(ctstar);

[ntim, nlat, nlon] = size(hgt500);

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

pat = reshape((ctstar' * reshape(hgt500, ntim, nlat*nlon)/ntim),nlat,nlon);