Documentation of tao_windstress_std

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

Help text

    1.6199 :  -6 to 6;
    1.7751 :  -3 to 3;

Cross-Reference Information

Listing of script tao_windstress_std

clean
cd /home/disk/tao/data/toga_tao/taohome/data/winds
nc = netcdf('uwnd_anom_hf_sm.cdf', 'nowrite');
  lat = nc{'lat', 1}(:);
  lon = nc{'lon', 1}(:);
  time = nc{'time', 1}(:);
  [xk, yk] = keep_var([150 220 -6 6], lon, lat);
  datu = nc{'WU_422', 1}(:,:,yk,xk);
  mv = nc.missing_value(:);
nc = close(nc);
datu(datu == mv) = NaN;
datu = squeeze(datu);
nc = netcdf('vwnd_anom_hf_sm.cdf', 'nowrite');
  datv = nc{'WV_423', 1}(:,:,yk,xk);
nc = close(nc);
datv(datv == mv) = NaN;
datv = squeeze(datv);
dat2 = squeeze(mean2(mean2(shiftdim(datu, 1))));
ind1 = max(find(isnan(dat2)))+1;
ind2 = max(find(~isnan(dat2)));
dat3 = dat2(ind1:ind2);
std(dat3)

std(dat3 - rave(rave(dat3, 13), 7))
%    0.7076 :  -6 to 6;
%    0.7444 :  -3 to 3;