Documentation of get_cti

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

Cross-Reference Information

Listing of script get_cti

clear
ctlims = [180 270 -6 6];
tim = 1:651; 
cd /home/disk/tao/data/nmc.reanalysis/monthly
[ct, latct, lonct] = getnc2('skt.mon.mean.nc', 'skt', ctlims, 1, tim);
[ntim, nlat, nlon] = size(ct);
scale = ones(2, nlat, nlon);
ct = cosweight(ct, latct);
ct = cosweight(ct, latct);
scale = cosweight(scale, latct);
scale = cosweight(scale, latct);
ct = squeeze(sum(sum(shiftdim(ct, 1))));
scale = sum(sum(squeeze(scale(1,:,:))));
ct = ct/scale;
ct_anncyc = ct;
ct2 = annave(ct);
desc = ['CTI, from NCEP reanalysis SKT.  Years given in tim.' ...
	'The data are cosine-weighted prior to averaging.   '];
cd /home/disk/tao/dvimont/matlab/NMC/Data
save cti2.mat ct ct_anncyc ct2 latct lonct ctlims tim desc