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')