Documentation of look_prec_ITCZ_run

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

Cross-Reference Information

Listing of script look_prec_ITCZ_run

clear
lims = [0 360 -90 90];
tim = 13:132;
lev = 1;
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ITCZ/prec
nc = netcdf('warm_precip.nc', 'nowrite');
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  prec1 = nc{'PRECL'}(:,:,:)+nc{'PRECC'}(:,:,:);
nc = close(nc);
nc = netcdf('cold_precip.nc', 'nowrite');
  prec2 = nc{'PRECL'}(:,:,:)+nc{'PRECC'}(:,:,:);
nc = close(nc);
[ntim, nlat, nlon] = size(prec1);
tim = sort([1:12:ntim 2:12:ntim 3:12:ntim ...
	    4:12:ntim 11:12:ntim 12:12:ntim]);
tim = sort([5:12:ntim 6:12:ntim 7:12:ntim ...
	    8:12:ntim 9:12:ntim 10:12:ntim]);
var1 = prec1(tim,:,:);
var2 = prec2(tim,:,:);
ntim = length(tim);
tem = 1000*24*3600*mean(var1 - var2);
[tstat, dof] = ttest(reshape(var1, ntim, nlat*nlon), ...
		     reshape(var2, ntim, nlat*nlon));
tstat = reshape(tstat, nlat, nlon);
score = tscore(dof, 5);
figure(1); fl(1); clf;
default_global;
gcont(tem, 0.5);
dc
color_shade(abs(tstat), score, 0.7*[1 1 1]);
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ITCZ/data
nc = netcdf('warmclim.nc', 'nowrite');
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  var1a = nc{'TAUX'}(:,:,:);
  var1b = nc{'TAUY'}(:,:,:);
nc = close(nc);
nc = netcdf('coldclim.nc', 'nowrite');
  var2a = nc{'TAUX'}(:,:,:);
  var2b = nc{'TAUY'}(:,:,:);
nc = close(nc);
tema = mean(var1a-var2a);
temb = mean(var1b-var2b);
default_global;
figure(2); fl(1); clf;
gquiv(tema, temb, -0.025, 1, '');