Documentation of define_1950-93_clim

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

Cross-Reference Information

Listing of script define_1950-93_clim

clear
cd /home/disk/hayes/dvimont/ccm3.6/data/OI_RUN
sst = getnc('sst_jan1978-sep1993.nc', 'SST');
[lat, lon] = getll('sst_jan1978-sep1993.nc');
sstclim = sst(1:12,:,:);
sstcur = sst(25:189,:,:);
[sstcur, climcur] = annave(sstcur);
cd /home/disk/hayes/dvimont/ccm3.6/data/GR
sst5079 = getnc('sst5079t31.nc', 'SST');
[lat31, lon31] = getll('sst5079t31.nc');
for i = 1:12;
  sst8093(i,:,:) = interp2(lon, lat, squeeze(climcur(i,:,:)), lon31, lat31');
end
default_global
tem = squeeze(mean(climcur - sstclim));
w = (79 - 50 + 1) / (93 - 50 + 1);
clim = w * sst5079 + (1 - w) * sst8093;
tem = squeeze(mean(clim - sstclim));
cd /home/disk/hayes/dvimont/ccm3.6/data/GR
nc = netcdf('sst5093t31.nc', 'write');
  sstold = nc{'SST'}(:,:,:);
  icpt = find(clim <= -1.799);
  clim(icpt) = -1.8 * ones(size(icpt));
  nc{'SST'}(:,:,:) = clim;
nc = close(nc);