Documentation of get_taux_std

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

Help text

  Start with FSU

Cross-Reference Information

Listing of script get_taux_std

clear
txlim = [150 230 -6 6];
ctlim = [180 270 -6 6];
cd /home/disk/tao/dvimont/Data/20dISO
tim1 = 1:444;
[taux, lat1, lon1] = getnc2('fsu_taux_data.cdf', 'taux', txlim, ...
					  1, tim1);
%  Get COADS for the same years (1961--1995)
cd /home/disk/tao/data/coads/coads1a
tim2 = 133:576;
[sst, lat2, lon2] = getnc2('sst.mean.5097.nc', 'sst', ctlim, ...
					 1, tim2);


%  Form indices
ct = squeeze(mean2(mean2(shiftdim(sst, 1))));
tx = squeeze(mean2(mean2(shiftdim(taux, 1))));

ct = annave(ct);
tx = annave(tx);

ct = detrend(ct);
tx = detrend(tx);

a1 = ct'*tx./((length(ct)-1)*var(ct));
tx2 = tx - a1*ct;

Tdamp = -30/log(corr(tx2, tx2, 1));

%  results
>> var(tx)
    0.0330
>> var(tx2)
    0.0191

>> std(tx) 
    0.1816
>> std(tx2)
    0.1381

>> sqrt(var(tx)-var(tx2)) 
    0.1179
>> sqrt(var(a1*ct))
    0.1179
>> var(a1*ct)
    0.0139