Documentation of ct_statistics

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

Help text

  Monthly

Cross-Reference Information

Listing of script ct_statistics

clear
[b1, a1] = butter(9, 2/9);
[b2, a2] = butter(9, 2/10);
ct = getct;
ct = detrend(ct);
[std(ct) std(filtfilt(b1, a1, ct)) std(ct - filtfilt(b2, a2, ct))]
ans =
    0.2799    0.2041    0.1918
cd /home/disk/tao/data/coads/sstanom4by6
nc = netcdf('sstcoadsanom4by6.18541995.nc', 'nowrite');
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  [xk, yk] = keep_var([180 270 -6 6], lon, lat);
  ct2 = nc{'data'}(553:1704, yk, xk);
  sf = nc{'data'}.scale_factor(:);
  ao = nc{'data'}.add_offset(:);
  mv = nc{'data'}.missing_value(:);
nc = close(nc);
ct2(ct2 == mv) = NaN;
ct2 = ct2*sf;
ct2 = squeeze(mean2(mean2(shiftdim(ct2, 1))));
for i = 1:(length(ct2)/12);
  ind = 12*(i-1) + [1:12];
  ct3(i) = mean2(ct2(ind));
end
[std(ct3) std(filtfilt(b1, a1, ct3)) std(ct3 - filtfilt(b2, a2, ct3))]
ans =
    0.5580    0.3300    0.4533
ct3 = detrend(ct3);
[std(ct3) std(filtfilt(b1, a1, ct3)) std(ct3 - filtfilt(b2, a2, ct3))]
ans =
    0.5137    0.2582    0.4533
cd /home/disk/tao/data/nmc.reanalysis/monthly
nc = netcdf('sst.mon.mean.nc', 'nowrite');
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  [xk, yk] = keep_var([180 270 -6 6], lon, lat);
  ct2 = nc{'air'}(:, yk, xk);
  sf = nc{'air'}.scale_factor(:);
  ao = nc{'air'}.add_offset(:);
  mv = nc{'air'}.missing_value(:);
nc = close(nc);
ct2(ct2 == mv) = NaN;
ct2 = ct2*sf;
ct2 = squeeze(mean2(mean2(shiftdim(ct2, 1))));
for i = 1:(length(ct2)/12);
  ind = 12*(i-1) + [1:12];
  ct4(i) = mean2(ct2(ind));
end
[std(ct4) std(filtfilt(b1, a1, ct4)) std(ct4 - filtfilt(b2, a2, ct4))]
ans =
    0.6061    0.3465    0.5766
ct4 = detrend(ct4);
[std(ct4) std(filtfilt(b1, a1, ct4)) std(ct4 - filtfilt(b2, a2, ct4))]
ans =
    0.5965    0.3156    0.5766

clear

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load seas_ct2.mat
std(ct)

ans =

    0.3268

>> std(detrend(ct))

ans =

    0.3123


cd /home/disk/tao/data/coads/sstanom4by6
nc = netcdf('sstcoadsanom4by6.18541995.nc', 'nowrite');
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  [xk, yk] = keep_var([180 270 -6 6], lon, lat);
 
  ct2 = nc{'data'}(553:1704, yk, xk);
  sf = nc{'data'}.scale_factor(:);
  ao = nc{'data'}.add_offset(:);
  mv = nc{'data'}.missing_value(:);
nc = close(nc);

ct2(ct2 == mv) = NaN;
ct2 = ct2*sf;
ct3 = squeeze(mean2(mean2(shiftdim(ct2, 1))));
[ct3, clim] = annave(ct3);

std(ct3(find(~isnan(ct3)))) 

ans =

     0.7373


tem = rave(ct3, 5);
ct3(find(isnan(ct3))) = tem(find(isnan(ct3)));
std(detrend(ct3))

ans =

    0.7051



cd /home/disk/tao/data/nmc.reanalysis/monthly
nc = netcdf('sst.mon.mean.nc', 'nowrite');
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  [xk, yk] = keep_var([180 270 -6 6], lon, lat);
 
  ct2 = nc{'air'}(:, yk, xk);
  sf = nc{'air'}.scale_factor(:);
  ao = nc{'air'}.add_offset(:);
  mv = nc{'air'}.missing_value(:);
nc = close(nc);

ct2(ct2 == mv) = NaN;
ct2 = ct2*sf;
ct4 = squeeze(mean2(mean2(shiftdim(ct2, 1))));
[ct4, clim] = annave(ct4);

std(ct4)

ans =

    0.7929

std(detrend(ct4))

ans =

    0.7858