Documentation of HP10_LP9_wreg

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

Cross-Reference Information

Listing of script HP10_LP9_wreg

clear
loadpcs;
lims = [100 300 -35 35];
tim = 101:1000;
wt = getheat(lims, 1:3, tim, 'wl');
wb = getheat(lims, 4:7, tim, 'wl');
[lat1, lon1, tem, lm] = getll('temp', lims);
[lat2, lon2] = getll('wl', lims);
tser = hpcs(:,1); lags1 = -1:1; biff = 1;
nlag = length(lags);
wbr = regress_eof(wb, tser, lags1);
tser = lpcs(:,1); lags2 = -4:3:2; biff = 1;
nlag = length(lags);
wbr2 = regress_eof(wb, tser, lags2);
figure(1); fo; clf;
default_global; FRAME = [105 295 -30 30];
for i = 1:nlag;
  sptalk(5,2,2*i-1);
    XAX = lon2; YAX = lat2;
    gcont(wbr(i,:,:), 0.025);
    XAX = lon1; YAX = lat1;
    dc2(lm);
    yl(i) = ylabel(['Lag = ' num2str(lags1(i))]);
    set(gca, 'fontsize', 9);
end
for i = 1:nlag;
  sptalk(5,2,2*i);
    XAX = lon2; YAX = lat2;
    gcont(wbr2(i,:,:), 0.015);
    XAX = lon1; YAX = lat1;
    dc2(lm);
    yl(i+nlag) = ylabel(['Lag = ' num2str(lags2(i))]);
    set(gca, 'fontsize', 9);
end
sptalk(5,2,1);
t(1) = title('< 80-270m W, HPPC1 >');
sptalk(5,2,2);
t(2) = title('< 80-270m W, LPPC1 >');
set(yl, 'fontsize', 11)
set(t, 'fontsize', 10);
cd ~/Thesis/Talk
print -dps2 HP10_LP9_W_80-270_eq_lags.ps