Global Index (short | long) | Local contents | Local Index (short | long)
Heat content and ocean currents
| This script calls | |
|---|---|
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat; pcs = -1*pcs;
lims = [100 300 -65 65];
tim = 101:1000;
lev = 1;
lags = -10:-2;
sst = getnc('temp', lims, lev, tim);
[lat, lon] = getll('temp', lims);
taux = getnc('taux', lims, 1, tim);
[latx, lonx] = getll('taux', lims);
[treg, tcoef] = regress_eof(sst, pcs, lags);
[xreg, xcoef] = regress_eof(taux, pcs, lags);
default_global; FRAME = [105 295 -40 60];
dx = mean(diff(XAX));
dy = mean(diff(YAX));
XAX2 = [XAX(1)-dx; XAX; XAX(nx)+dx];
YAX2 = [YAX(1)-dy; YAX; YAX(ny)+dy];
ind = 4;
cint = 0.05; clev = [-2:cint:-cint cint:cint:2];
for i = 1:4;
subplot(4,2,2*i-1);
treg2 = squeeze(treg(i+ind,:,:));
[ny, nx] = size(treg2);
treg2 = [ones(1, nx+2); [ones(ny, 1) treg2 ones(ny, 1)]; ones(1, nx+2)];
contourf2(XAX2, YAX2, treg2, [-0.25:.025:.275]);
hold on;
pncont(XAX, YAX, squeeze(treg(i+ind,:,:)), [-0.25:.05:.3], 0, 'k');
hold off;
ylabel(['Lag = ' num2str(lags(i+ind))]);
axis([105 295 -60 60]);
colorbar2
caxis([-0.25 .275])
set(gca, 'XTick', 120:30:270, 'YTick', -60:30:60);
dc3(0);
end
xlabel(['Contour Interval: ' num2str(cint) ' K std^-^1']);
XAX = lonx; YAX = latx;
dx = mean(diff(XAX));
dy = mean(diff(YAX));
XAX2 = [XAX(1)-dx; XAX; XAX(nx)+dx];
YAX2 = [YAX(1)-dy; YAX; YAX(ny)+dy];
cint = 2; clev = [-2:cint:-cint cint:cint:2];
for i = 1:4;
subplot(4,2,2*i);
xreg2 = squeeze(xreg(i+ind,:,:));
[ny, nx] = size(xreg2);
xreg2 = [ones(1, nx+2); [ones(ny, 1) xreg2 ones(ny, 1)]; ones(1, nx+2)];
contourf2(XAX2, YAX2, 100*xreg2, [-8.5:.5:5]);
hold on;
pncont(XAX, YAX, squeeze(100*xreg(i+ind,:,:)), [-8:5], 0, 'k');
hold off;
ylabel(['Lag = ' num2str(lags(i+ind))]);
axis([105 295 -60 60]);
colorbar2
caxis([-8.5 5])
dc3(0);
end
xlabel(['Contour Interval: 1x10^-^3 N m^-^2 std^-^1']);
subplot(4,2,1);
title('LP10: 12.5m TEMP');
colorbar2
subplot(4,2,2);
title('LP10: ZONAL WIND STRESS');
colorbar2
cd /home/disk/tao/dvimont/matlab/CSIRO/IUGG_talk/Plots
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat; pcs = -1*pcs;
lims = [100 300 -65 65];
tim = 101:1000;
lev = 1;
lags = -6:-3;
hc = getheat(lims, 4:7, tim);
[lat, lon] = getll('temp', lims);
u = getheat(lims, 4:7, tim, 'u');
v = getheat(lims, 4:7, tim, 'v');
[latv, lonv] = getll('v', lims);
[hreg, hcoef] = regress_eof(hc, pcs, lags);
[ureg, ucoef] = regress_eof(u, pcs, lags);
[vreg, vcoef] = regress_eof(v, pcs, lags);
default_global; FRAME = [105 295 -40 60];
ind = 0;
cint = 0.15; clev = [-3:cint:-cint cint:cint:3];
for i = 1:4;
subplot(4,2,2*i-1);
gcont(1e-8*hreg(i+ind,:,:), clev);
dc2(hreg(i,:,:));
% color_shade(squeeze(tcoef(i+ind,:,:).^2), 0.1, [.9 .9 .9]);
ylabel(['Lag = ' num2str(lags(i+ind))]);
end
xlabel(['Contour Interval: ' num2str(cint) ' J m¯^2 std^-^1']);
yl = get(gca, 'YLim');
FRAME = [105 295 -30 60];
ind = 0;
for i = 1:4;
subplot(4,2,2*i);
XAX = lonv; YAX = latv;
gquiv(ureg(i+ind,:,:), vreg(i+ind,:,:), -.15, 2, 'cm s^-^1 std^-^1');
XAX = lon; YAX = lat;
dc2(hreg(i,:,:));
set(gca, 'YLim', yl);
% color_shade(squeeze(tcoef(i+ind,:,:).^2), 0.1, [.9 .9 .9]);
ylabel(['Lag = ' num2str(lags(i+ind))]);
end
subplot(4,2,1);
title('LP10: 80-270m HEAT CONTENT');
subplot(4,2,2);
title('LP10: 80-270m OCEAN VELOCITY');
cd /home/disk/tao/dvimont/matlab/CSIRO/IUGG_talk/Plots
% Heat Flux and ocean currents
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat; pcs = -1*pcs;
lims = [100 300 -65 65];
tim = 101:1000;
lev = 1;
lags = -6:-3;
hc = getnc('heat', lims, 1, tim);
[lat, lon] = getll('temp', lims);
u = getheat(lims, 1:3, tim, 'u');
v = getheat(lims, 1:3, tim, 'v');
[latv, lonv] = getll('v', lims);
[hreg, hcoef] = regress_eof(hc, pcs, lags);
[ureg, ucoef] = regress_eof(u, pcs, lags);
[vreg, vcoef] = regress_eof(v, pcs, lags);
default_global; FRAME = [105 295 -40 60];
ind = 0;
cint = 0.5; clev = [-5:cint:-cint cint:cint:5];
for i = 1:4;
subplot(4,2,2*i-1);
gcont(hreg(i+ind,:,:), clev);
dc2(hreg(i,:,:));
% color_shade(squeeze(tcoef(i+ind,:,:).^2), 0.1, [.9 .9 .9]);
ylabel(['Lag = ' num2str(lags(i+ind))]);
end
xlabel(['Contour Interval: ' num2str(cint) ' W m^-^2 std^-^1']);
yl = get(gca, 'YLim');
FRAME = [105 295 -30 60];
ind = 0;
for i = 1:4;
subplot(4,2,2*i);
XAX = lonv; YAX = latv;
gquiv(ureg(i+ind,:,:), vreg(i+ind,:,:), -.25, 2, 'cm s^-^1 std^-^1');
XAX = lon; YAX = lat;
dc2(hreg(i,:,:));
set(gca, 'YLim', yl);
% color_shade(squeeze(tcoef(i+ind,:,:).^2), 0.1, [.9 .9 .9]);
ylabel(['Lag = ' num2str(lags(i+ind))]);
if i ~= 4; xlabel(''); end;
end
subplot(4,2,1);
title('LP10: NET HEAT FLUX (Pos. Down)');
subplot(4,2,2);
title('LP10: 0-80m OCEAN VELOCITY');
cd /home/disk/tao/dvimont/matlab/CSIRO/IUGG_talk/Plots