Documentation of lag_regress_on_ct

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

Help text

cint = [0.1 0.1 0.1];

Cross-Reference Information

Listing of script lag_regress_on_ct

clean
ct = getct;
ct = detrend(ct);
[b1, a1] = butter(9, 2/4);
[b2, a2] = butter(9, 2/12);
ct2 = filtfilt(b1, a1, ct)-filtfilt(b2, a2, ct);
[b1, a1] = butter(9, 2/7);
ct2 = filtfilt(b1, a1, ct);
lims = [90 305 -80 80];
tim = 101:1000;
varn = 'temp';
dat1 = getnc(varn, lims, 1, tim);
[lat1, lon1, depth, lm] = getll(varn, lims);
varn = 'psl';
dat2 = getflx(varn, lims, tim);
[lat2, lon2] = getll('temp', lims);
varn = 'taux';
dat3 = getnc(varn, lims, 1, tim);
[lat3, lon3] = getll(varn, lims);
varn = 'tauy';
dat4 = getnc(varn, lims, 1, tim);
[lat4, lon4] = getll(varn, lims);
[dat3, lat3, lon3] = getheat(lims, 1:3, tim, 'u');
[dat4, lat4, lon4] = getheat(lims, 1:3, tim, 'v');
dat1 = detrend(dat1);
dat2 = detrend(dat2);
dat3 = detrend(dat3);
dat4 = detrend(dat4);
dat1 = filtfilt(b1, a1, dat1)-filtfilt(b2, a2, dat1);
dat2 = filtfilt(b1, a1, dat2)-filtfilt(b2, a2, dat2);
dat3 = filtfilt(b1, a1, dat3)-filtfilt(b2, a2, dat3);
dat1 = dat1 - filtfilt(b1, a1, dat1);
dat2 = dat2 - filtfilt(b1, a1, dat2);
dat3 = dat3 - filtfilt(b1, a1, dat3);
dat1 = filtfilt(b1, a1, dat1);
dat2 = filtfilt(b1, a1, dat2);
dat3 = filtfilt(b1, a1, dat3);
dat4 = filtfilt(b1, a1, dat4);
lags = -5:5;
[reg1, c1] = regress_eof(dat1, ct2, lags);
[reg2, c2] = regress_eof(dat2, ct2, lags);
[reg3, c3] = regress_eof(dat3, ct2, lags);
[reg4, c4] = regress_eof(dat4, ct2, lags);
cint = [0.05 0.1 0.1 0.1];
pn = [1000 -1000 1000 1000];

lind = find(ismember(lags, -2:3));
for i = 3:4;
  eval(['tem = reg' num2str(i) ';']);
  figure(i+1); fo(1); clf;
  for j = 1:length(lind);
    subplot(3,2,j);
      eval(['dg2(lat' num2str(i) ', lon' num2str(i) ');']);
      FRAME = [95 300 -60 60];
      gcont(tem(lind(j),:,:), cint(i));
      dg2(lat1, lon1);
      dc2(lm, 0.7, pn(i));
      title(['Lag ' num2str(j)]);
  end
end