Documentation of std_ct

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

Cross-Reference Information

Listing of script std_ct

clear
cd ~/matlab/CSIRO/Thesis/Data
load seas_ct2.mat
DJF CT:
std(detrend(ct2(1:4:3600)))
    0.2989
MAM CT:
std(detrend(ct2(2:4:3600)))
    0.2903
JJA CT:
std(detrend(ct2(3:4:3600)))
    0.3375
SON CT:
std(detrend(ct2(4:4:3600)))
    0.2891
for i = 1:12;
  std_ct(i) = std(detrend(ct(i:12:10800)));
end
mon = ['jan'; 'feb'; 'mar'; 'apr'; 'may'; 'jun'; ...
       'jul'; 'aug'; 'sep'; 'oct'; 'nov'; 'dec'];
figure(1); clf; fo
subplot(2,1,1)
  h = plot(1:12, std_ct, 'x-k')
  axis([0 13 0 1]);
  grid on
  set(gca, 'XTick', 1:12)
  set(gca, 'XTickLabel', mon);
  title('Annual Cycle of the Standard Deviation of the CT index -- CSIRO model');
  ylabel('K')
  xlabel('Month')
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);
  coads = nc{'data'}((104*12+1):1704, yk, xk);
  mv = nc{'data'}.missing_value(:);
  ao = nc{'data'}.add_offset(:);
  sf = nc{'data'}.scale_factor(:);
nc = close(nc);
coads(coads == mv) = NaN;
coads = squeeze(mean2(mean2(shiftdim(coads, 1))));
coads = coads*sf + ao;
for i = 1:12;
  tem = coads(i:12:length(coads));
  std_ct2(i) = std(detrend(tem(~isnan(tem))));
end
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);
  nmc = nc{'air'}(:, yk, xk);
  mv = nc{'air'}.missing_value(:);
  ao = nc{'air'}.add_offset(:);
  sf = nc{'air'}.scale_factor(:);
nc = close(nc);
nmc(nmc == mv) = NaN;
nmc = squeeze(mean2(mean2(shiftdim(nmc, 1))));
nmc = nmc*sf + ao;
for i = 1:12;
  tem = nmc(i:12:length(nmc));
  std_ct3(i) = std(detrend(tem(~isnan(tem))));
end
figure(1); clf; fo
subplot(2,1,1)
  h = plot(1:12, std_ct, 'x-k', 1:12, std_ct2, 'o-k', 1:12, std_ct3, 's-k')
  axis([0 13 0 1]);
  grid on
  set(gca, 'XTick', 1:12)
  set(gca, 'XTickLabel', mon);
  title('Annual Cycle of the Standard Deviation of the CT index');
  ylabel('K')
  xlabel('Month')
  legend(h, 'CSIRO (101 - 1000)', 'COADS (1900 - 1995)', 'NMC (1958 - 1997)')
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Seasonal
print -dps2 CSIRO_COADS_CT_STD.ps