Global Index (short | long) | Local contents | Local Index (short | long)
Plot negative lags
This script calls | |
---|---|
clear cd ~/matlab/CSIRO/Thesis/Data load LP9_detrend_L1-7_EOF_yr101-1000.mat; pcs = 1*pcs; lims = [150 205 -44 44]; tim = 101:1000; lags = -15:15; cd ~/matlab/CSIRO/Thesis/Chap5 [vdtdy, wdtdz] = ... vert_vdtdy_wdtdz(pcs, lims, lags); tit1 = 'Vbar T''_Y'; tit2 = 'Wbar T''_Z'; [vdtdy, wdtdz] = ... vert_vptb_wptb(pcs, lims, lags); tit1 = 'V'' Tbar_Y'; tit2 = 'W'' Tbar_Z'; [vdtdy, wdtdz] = vert_ubtp_uptb(pcs, lims, lags); tit1 = 'Ubar T''_X'; tit2 = 'U'' Tbar_X'; [latv, lonv, depthv] = getll('u', lims); [latw, lonw, depthw] = getll('u', lims); [latv, lonv, depthv] = getll('v', lims); [latw, lonw, depthw] = getll('wl', lims); vdtdy2 = 1e9*squeeze(mean2(shiftdim(vdtdy, 3))); wdtdz2 = 1e9*squeeze(mean2(shiftdim(wdtdz, 3))); cint = 0.5; clev1 = sort([-cint:-cint:min(min(min(vdtdy2))) cint:cint:max(max(max(vdtdy2)))]); clev2 = sort([-cint:-cint:min(min(min(wdtdz2))) cint:cint:max(max(max(wdtdz2)))]); clev = [-10 -2:.1:2 20]; figure(5); fo(1); lagplot = -14:2:-2; lp = length(lagplot); for i = 1:lp; subplot(lp, 2, 2*i-1); tind = find(lags == lagplot(i)); contourf(latv, -1*depthv, squeeze(1*vdtdy2(tind,:,:)), clev); caxis([-2 2]); % pncont(latv, -1*depthv, squeeze(1*vdtdy2(tind,:,:)), clev1, 0, 'k'); set(gca, 'YTick', -500:125:0, 'YTickLabel', 500:-125:0, ... 'XTick', -40:10:40); axis([-40 40 -545 0]); shading flat ylabel(['Lag = ' num2str(lagplot(i))]); hold on; vline(-40:10:40, ':k'); hline(-500:125:-125, ':k'); hold off; end for i = 1:lp; subplot(lp, 2, 2*i); tind = find(lags == lagplot(i)); contourf(latw, -1*depthw, squeeze(1*wdtdz2(tind,:,:)), clev); caxis([-2 2]); % pncont(latw, -1*depthw, squeeze(1*wdtdz2(tind,:,:)), clev2, 0, 'k'); set(gca, 'YTick', -500:125:0, 'YTickLabel', 500:-125:0, ... 'XTick', -40:10:40); axis([-40 40 -545 0]); shading flat hold on; vline(-40:10:40, ':k'); hline(-500:125:-125, ':k'); hold off; end subplot(lp, 2, 1); title(tit1); subplot(lp, 2, 2); title(tit2); subplot(lp, 2, 2*lp-1); % xlabel(['Contour Interval: ' num2str(cint) 'x10^9 K s^-^1 std^-^1']); subplot(lp, 2, 2*lp); % xlabel(['Contour Interval: ' num2str(cint) 'x10^9 K s^-^1 std^-^1']); %%%%%%%%% Positive lags figure(6); fo(1); lagplot = 0:2:12; lp = length(lagplot); for i = 1:lp; subplot(lp, 2, 2*i-1); tind = find(lags == lagplot(i)); contourf(latv, -1*depthv, squeeze(1*vdtdy2(tind,:,:)), clev); caxis([-2 2]); % pncont(latv, -1*depthv, squeeze(1*vdtdy2(tind,:,:)), clev1, 0, 'k'); set(gca, 'YTick', -500:125:0, 'YTickLabel', 500:-125:0, ... 'XTick', -40:10:40); axis([-40 40 -545 0]); shading flat ylabel(['Lag = ' num2str(lagplot(i))]); hold on; vline(-40:10:40, ':k'); hline(-500:125:-125, ':k'); hold off; end for i = 1:lp; subplot(lp, 2, 2*i); tind = find(lags == lagplot(i)); contourf(latw, -1*depthw, squeeze(1*wdtdz2(tind,:,:)), clev); caxis([-2 2]); % pncont(latw, -1*depthw, squeeze(1*wdtdz2(tind,:,:)), clev2, 0, 'k'); set(gca, 'YTick', -500:125:0, 'YTickLabel', 500:-125:0, ... 'XTick', -40:10:40); axis([-40 40 -545 0]); shading flat hold on; vline(-40:10:40, ':k'); hline(-500:125:-125, ':k'); hold off; end subplot(lp, 2, 1); title(tit1); subplot(lp, 2, 2); title(tit2); subplot(lp, 2, 2*lp-1); % xlabel(['Contour Interval: ' num2str(cint) 'x10^9 K s^-^1 std^-^1']); subplot(lp, 2, 2*lp); % xlabel(['Contour Interval: ' num2str(cint) 'x10^9 K s^-^1 std^-^1']); cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap5/Plots % Plot Vbar and Wbar vbar = getnc('v', lims, 1:10, tim); wbar = getnc('wl', lims, 1:9, tim); vbar = mean(vbar); wbar = mean(wbar); vbar = squeeze(mean2(shiftdim(vbar, 2))); wbar = squeeze(mean2(shiftdim(wbar, 2))); figure(3); clf; fo(1); subplot(3, 1, 1); pncont(latv, -1*depthv, vbar, [-20:0.25:20], 0, 'k'); axis([-40 40 -545 0]); set(gca, 'YTick', -500:125:0, 'YTickLabel', 500:-125:0, ... 'XTick', -40:10:40); title('Vbar'); xlabel('Contour Interval: 0.25 cm s^-^1') subplot(3, 1, 2); pncont(latw, -1*depthw, wbar, [-20:0.25:20], 0, 'k'); axis([-40 40 -545 0]); set(gca, 'YTick', -500:125:0, 'YTickLabel', 500:-125:0, ... 'XTick', -40:10:40); title('Wbar'); xlabel('Contour Interval: 2.5 x 10^-^5 cm s^-^1') temp = getnc('temp', lims, 1:10, tim); temp = detrend(temp); treg = regress_eof(temp, pcs, 0); treg = squeeze(mean2(shiftdim(treg, 3))); [latt, lont, deptht] = getll('temp', lims); subplot(3, 1, 3); pncont(latt, -1*deptht, -treg, [-.3:0.025:.3], 0, 'k'); axis([-40 40 -545 0]); set(gca, 'YTick', -500:125:0, 'YTickLabel', 500:-125:0, ... 'XTick', -40:10:40); title('Wbar'); xlabel('Contour Interval: 0.025 K std^-^1')