Documentation of warm_WPAC.define

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

Cross-Reference Information

Listing of script warm_WPAC.define

cd /home/disk/hayes/dvimont/ccm3.6/data/GR
[sst] = getnc('warm_WPAC.t31x15.nc', 'SST');
[lat, lon] = getll('warm_WPAC.t31x15.nc');
sstw = getnc('wgr_TPAC_t31.nc', 'SST');
sstc = getnc('cgr_TPAC_t31.nc', 'SST');
ssta = (sstw-sstc);
[ntim, nlat, nlon] = size(sst);
ssta2 = zeros(nlat, nlon);
[xk, yk] = keep_var([150 210 -10 10], lon, lat);
ssta2(yk, xk) = ones(length(yk), length(xk));
ssta2 = myrunning_ave(myrunning_ave(ssta2, 3), 3);
ssta2 = myrunning_ave(myrunning_ave(ssta2', 3), 3);
for i = 1:12;
  ssta(i,:,:) = squeeze(ssta(i,:,:)) .* ssta2;
end;
nc = netcdf('warm_WPAC.t31x15.nc', 'write');
  sst = nc{'SST'}(:,:,:);
  sst = sst + ssta;
  nc{'SST'}(:,:,:) = sst;
nc = close(nc);
nc = netcdf('cold_WPAC.t31x15.nc', 'write');
  sst = nc{'SST'}(:,:,:);
  sst = sst - ssta;
  nc{'SST'}(:,:,:) = sst;
nc = close(nc);