Documentation of define_weight_netcdf

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

Help text

  Define weight

Cross-Reference Information

Listing of script define_weight_netcdf

clear
cd /home/disk/hayes/dvimont/ccm3.6/data/SOM/cl_ice
nc = netcdf('sst5093t31.nc', 'nowrite');
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
nc = close(nc);
nc = netcdf('biff2.nc', 'write');
  nc{'lat'}(:) = lat;
  nc{'lon'}(:) = lon;
nc = close(nc);
nlat = length(lat); nlon = length(lon);
get_global
default_global

ssta = zeros(nlat, nlon);
[xk, yk] = keep_var([110 260 0 25], lon, lat);
ssta(yk, xk) = ones(length(yk), length(xk));
[xk, yk] = keep_var([145 295 -25 0], lon, lat);
ssta(yk, xk) = ones(length(yk), length(xk));

[xk, yk] = keep_var([260 295 0 25]);
for i = 1:length(yk);
  xk = find(lon >= 260 & lon <= (295 - (35/25)*lat(yk(i))));
  ssta(yk(i), xk) = ones(1, length(xk));
end

[xk, yk] = keep_var([110 145 -25 0], lon, lat);
for i = length(yk):-1:1;
  xk = find(lon <= 145 & lon >= (110 - (35/25)*lat(yk(i))));
  ssta(yk(i), xk) = ones(1, length(xk));
end

ssta = myrunning_ave(myrunning_ave(ssta, 3), 3);
ssta = (myrunning_ave(myrunning_ave(ssta', 3), 3))';

weight = zeros(12, nlat, nlon);
for i = 1:12
  weight(i,:,:) = ssta;
end

%  Define SST

cd /home/disk/hayes/dvimont/ccm3.6/data/GR
nc = netcdf('cgr_TPAC_t31.nc', 'nowrite');
  sst = nc{'SST'}(:);
nc = close(nc);

%  Define QO and MLD

cd /home/disk/hayes/dvimont/ccm3.6/data/SOM/cl_ice
nc = netcdf('qomld5093.nc', 'nowrite');
  qo = nc{'QO'}(:);
  mld = nc{'MLD'}(:);
nc = close(nc);

sst = reshape(sst, 12, 1, 48, 48);
weight = reshape(weight, 12, 1, 48, 48);

%  Write variables out to biff2.nc

nc = netcdf('biff2.nc', 'write');
  nc{'WEIGHT'}(:) = weight;
  nc{'SST'}(:) = sst;
  nc{'QO'}(:) = qo;
  nc{'MLD'}(:) = mld;
nc = close(nc);