Documentation of covar_gr

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

Cross-Reference Information

Listing of script covar_gr

clear
cd /home/disk/tao/data/coads/sstanom4by6
filin = 'sstcoadsanom4by6.18541993.nc';
nc = netcdf(filin, 'nowrite');
  ts = nc{'data'}(1153:1680,:,:);
  datets = nc{'time'}(1153:1680);
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  add_offset = nc{'data'}.add_offset(:);
  scale_factor = nc{'data'}.scale_factor(:);
  mv = nc{'data'}.missing_value(:);
nc = close(nc)
ts = ts * scale_factor;
ts = ts + add_offset;
mv = mv * scale_factor;
mv = mv + add_offset;
[ntim, nlat, nlon] = size(ts);
ts = reshape(ts, ntim, nlat*nlon);
tsave = mean(ts);
kp = find(tsave ~= mv);
clear ts tsave
cd /home/disk/hayes2/dvimont/data
load coads_sst.mat
ntim = size(ts, 1);
tsnew = NaN * ones(ntim, nlat*nlon);
tsnew(:, kp) = ts;
tsnew = reshape(tsnew, ntim, nlat, nlon);
tsnew = cosweight(tsnew, lat);
tsnew = reshape(tsnew, ntim, nlat*nlon);
tsnew = tsnew(:, kp);
c = covar_nan(tsnew', tsnew');
cd /home/disk/tao/dvimont/CCM/GR/DATA
save covar.mat c
end