Global Index (short | long) | Local contents | Local Index (short | long)
tser = lpcs(:,1); lags = -8:2:6; biff = 2;
This script calls | |
---|---|
clear loadpcs; lims = [100 300 -35 35]; tim = 101:1000; ut = getheat(lims, 1:3, tim, 'u'); vt = getheat(lims, 1:3, tim, 'v'); ub = getheat(lims, 4:7, tim, 'u'); vb = getheat(lims, 4:7, tim, 'v'); wt = getheat(lims, 1:3, tim, 'wl'); wb = getheat(lims, 4:7, tim, 'wl'); [lat, lon] = getll('u', lims); [lat1, lon1, tem, lm] = getll('temp', lims); [lat2, lon2] = getll('wl', lims); tser = hpcs(:,1); lags = -1:1; biff = 1; nlag = length(lags); utr = regress_eof(ut, tser, lags); vtr = regress_eof(vt, tser, lags); ubr = regress_eof(ub, tser, lags); vbr = regress_eof(vb, tser, lags); wtr = regress_eof(wt, tser, lags); wbr = regress_eof(wb, tser, lags); default_global; FRAME = [105 240 -30 30]; figure(biff); fo(1); clf; for i = 1:nlag XAX = lon; YAX = lat; subplot(3,2,2*i-1); gquiv(utr(i,:,:), vtr(i,:,:), -1, 2, ''); hold on; XAX = lon2; YAX = lat2; gcont(wtr(i,:,:), 0.05); hold off; XAX = lon1; YAX = lat1; dc2(lm); xlabel(''); ylabel(['Lag = ' num2str(lags(i))]); XAX = lon; YAX = lat; subplot(3,2,2*i-1); gquiv(ubr(i,:,:), vbr(i,:,:), -.5, 2, ''); hold on; XAX = lon2; YAX = lat2; gcont(wbr(i,:,:), 0.025); hold off; XAX = lon1; YAX = lat1; dc2(lm); xlabel(''); end figure(biff+2); fo(1); clf; for i = 1:nlag XAX = lon2; YAX = lat2; subplot(6,2,2*i-1); hold on; gcont(wtr(i,:,:), 0.05); % XAX = lon1; YAX = lat1; % dc2(lm); % ylabel(['Lag = ' num2str(lags(i))]); XAX = lon2; YAX = lat2; subplot(6,2,2*i); gcont(wbr(i,:,:), 0.015); XAX = lon1; YAX = lat1; dc2(lm); end tser = hpcs(:,1); lags = -1:1; biff = 1; nlag = length(lags); ubr = regress_eof(ub, tser, lags); vbr = regress_eof(vb, tser, lags); wbr = regress_eof(wb, tser, lags); tser = lpcs(:,1); lags = -4:3:2; biff = 1; nlag = length(lags); ubr2 = regress_eof(ub, tser, lags); vbr2 = regress_eof(vb, tser, lags); wbr2 = regress_eof(wb, tser, lags); figure(1); fl; clf; default_global; FRAME = [105 270 -30 30]; for i = 1:nlag; XAX = lon; YAX = lat; subplot(3,2,2*i-1); % gquiv(ubr(i,:,:), vbr(i,:,:), -0.5, 2, ''); % hold on; XAX = lon2; YAX = lat2; gcont(wbr(i,:,:), 0.025); hold off; XAX = lon1; YAX = lat1; dc2(lm); xlabel(''); ylabel(['Lag = ' num2str(lags(i))]); end for i = 1:nlag; XAX = lon; YAX = lat; subplot(3,2,2*i); % gquiv(ubr2(i,:,:), vbr2(i,:,:), -0.5, 2, ''); % hold on; XAX = lon2; YAX = lat2; gcont(wbr2(i,:,:), 0.015); hold off; XAX = lon1; YAX = lat1; dc2(lm); xlabel(''); ylabel(['Lag = ' num2str(lags(i))]); end subplot(3,2,1); title('< 80-270m W, HPPC1 >'); subplot(3,2,2); title('< 80-270m W, LPPC1 >'); cd ~/Thesis/Talk print -dps2 HP10_LP9_W_eq_lags.ps % Look at wind stress tx = getnc('taux', lims, 1, tim); ty = getnc('tauy', lims, 1, tim); tser = hpcs(:,1); lags = -1:1; biff = 1; nlag = length(lags); ubr = regress_eof(tx, tser, lags); vbr = regress_eof(ty, tser, lags); wbr = regress_eof(wb, tser, lags); tser = lpcs(:,1); lags = -4:3:2; biff = 1; nlag = length(lags); ubr2 = regress_eof(tx, tser, lags); vbr2 = regress_eof(ty, tser, lags); wbr2 = regress_eof(wb, tser, lags); figure(1); fl; clf; default_global; FRAME = [105 270 -30 30]; for i = 1:nlag; XAX = lon; YAX = lat; subplot(3,2,2*i-1); gquiv(ubr(i,:,:), vbr(i,:,:), -0.05, 1, ''); hold on; XAX = lon2; YAX = lat2; % gcont(wbr(i,:,:), 0.025); hold off; XAX = lon1; YAX = lat1; dc2(lm); xlabel(''); ylabel(['Lag = ' num2str(lags(i))]); end for i = 1:nlag; XAX = lon; YAX = lat; subplot(3,2,2*i); gquiv(ubr2(i,:,:), vbr2(i,:,:), -0.05, 1, ''); hold on; XAX = lon2; YAX = lat2; % gcont(wbr2(i,:,:), 0.015); hold off; XAX = lon1; YAX = lat1; dc2(lm); xlabel(''); ylabel(['Lag = ' num2str(lags(i))]); end