Global Index (short | long) | Local contents | Local Index (short | long)
Define weight
This script calls | |
---|---|
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);