Global Index (short | long) | Local contents | Local Index (short | long)
This script calls | |
---|---|
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