Documentation of ct_coads_hgt_nmc

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

Function Synopsis

pat = ct_coads_hgt_nmc(lev);

Cross-Reference Information

Listing of function ct_coads_hgt_nmc

function pat = ct_coads_hgt_nmc(lev);
cur_dir = cd;
cd /home/disk/tao/dvimont/pub
filin = 'CT_58.bin';
fid = fopen(filin, 'r','b');
ctstar = fread(fid, 'float');
ctstar = ctstar(37:432);
cd /home/disk/tao/data/nmc.reanalysis/monthly
filin = 'hgt.mon.mean.nc'
nc = netcdf(filin, 'nowrite');
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  hgt500 = squeeze(nc{'hgt'}(37:432,lev,:,:));
nc = close(nc);
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);