Global Index (short | long) | Local contents | Local Index (short | long)
Get lower layer
| This script calls | |
|---|---|
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat
load HP10_detrend_L1-7_EOF_yr101-1000.mat
lev = 1:3;
lims = [178 272 -3 3]; % EQ reg
lims = [151 180.5 -3 3]; % West EQ reg
lims = [118 152 7.5 24]; % WSTP reg
lims = [179.9 214 30 37]; % NP cold anom
lims = [168 209 27 40]; % NP cold anom (2nd attempt);
lims = [168 203 20 33.5]; % NP cold anom (HP10 attempt);
lims = [213 242 30 37]; % CALIF warm anom
nfrm = -12:12; lags = nfrm;
tim = 101:1000;
cd ~/matlab/CSIRO/Heat/Old_routines
[utem, vtem, lat_out, lon_out, depth_out] = ...
get_ubar_gradtprime(pcs, lims, nfrm, tim, lev);
ubtp13 = squeeze(mean2(mean2(shiftdim(utem, 1))));
vbtp13 = squeeze(mean2(mean2(shiftdim(vtem, 1))));
[utem, vtem, lat_out, lon_out, depth_out] = ...
get_uprime_gradtbar(pcs, lims, nfrm, tim, lev);
uptb13 = squeeze(mean2(mean2(shiftdim(utem, 1))));
vptb13 = squeeze(mean2(mean2(shiftdim(vtem, 1))));
[wtem, lat_out, lon_out] = ...
get_wbar_dtprimedz(pcs, lims, nfrm, tim, lev);
wbtp13 = squeeze(mean2(mean2(shiftdim(wtem, 1))))';
[wtem, lat_out, lon_out] = ...
get_wprime_dtbardz(pcs, lims, nfrm, tim, lev);
wptb13 = squeeze(mean2(mean2(shiftdim(wtem, 1))))';
tim2 = 100:1000; ntim = length(tim2);
hc13 = getheat(lims, lev, tim2);
dhc13dt = (hc13(3:ntim,:,:) - hc13(1:(ntim-2),:,:)) / (2*3600*24*365);
hc13 = hc13(2:ntim,:,:); ntim = length(tim);
hc13 = regress_eof(hc13, pcs, nfrm);
dhc13dt = regress_eof(dhc13dt, pcs(1:(ntim-1), :), nfrm);
hc13 = squeeze(mean2(mean2(shiftdim(hc13, 1))))';
dhc13dt = squeeze(mean2(mean2(shiftdim(dhc13dt, 1))))';
hflx = getnc('heat', lims, 1, tim);
hflx = regress_eof(hflx, pcs, nfrm);
hflx = squeeze(mean2(mean2(shiftdim(hflx, 1))))';
cd ~/matlab/CSIRO/Heat/Old_routines
lev = 1:7;
[utem, vtem, lat_out, lon_out, depth_out] = ...
get_ubar_gradtprime(pcs, lims, nfrm, tim, lev);
ubtp17 = squeeze(mean2(mean2(shiftdim(utem, 1))));
vbtp17 = squeeze(mean2(mean2(shiftdim(vtem, 1))));
[utem, vtem, lat_out, lon_out, depth_out] = ...
get_uprime_gradtbar(pcs, lims, nfrm, tim, lev);
uptb17 = squeeze(mean2(mean2(shiftdim(utem, 1))));
vptb17 = squeeze(mean2(mean2(shiftdim(vtem, 1))));
[wtem, lat_out, lon_out] = ...
get_wbar_dtprimedz(pcs, lims, nfrm, tim, lev);
wbtp17 = squeeze(mean2(mean2(shiftdim(wtem, 1))))';
[wtem, lat_out, lon_out] = ...
get_wprime_dtbardz(pcs, lims, nfrm, tim, lev);
wptb17 = squeeze(mean2(mean2(shiftdim(wtem, 1))))';
tim2 = 100:1000; ntim = length(tim2);
hc17 = getheat(lims, lev, tim2);
dhc17dt = (hc17(3:ntim,:,:) - hc17(1:(ntim-2),:,:)) / (2*3600*24*365);
hc17 = hc17(2:ntim,:,:); ntim = length(tim);
hc17 = regress_eof(hc17, pcs, nfrm);
dhc17dt = regress_eof(dhc17dt, pcs(1:(ntim-1), :), nfrm);
hc17 = squeeze(mean2(mean2(shiftdim(hc17, 1))))';
dhc17dt = squeeze(mean2(mean2(shiftdim(dhc17dt, 1))))';
taux = getnc('taux', [165 210 -3 3], 1, tim);
taux = regress_eof(taux, pcs, nfrm);
taux = squeeze(mean2(taux'));
taux = getnc('taux', lims, 1, tim);
taux = regress_eof(taux, pcs, nfrm);
tauy = getnc('tauy', lims, 1, tim);
tauy = regress_eof(tauy, pcs, nfrm);
[lat, lon] = getll('tauy', lims);
for i = 1:length(nfrm);
tx = squeeze(taux(i,:,:)); ty = squeeze(tauy(i,:,:));
curlt(i,:,:) = sph_curl1(tx, ty, lat, lon, 0);
end
curlt = squeeze(mean2(mean2(shiftdim(curlt, 1))));
% Plot all the terms
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_eq_terms.mat
load LP10_wstp_terms.mat
load LP10_npcold_anom_terms.mat
load LP10_calif_anom_terms.mat
load HP10_west_eq_anom_terms.mat
load HP10_np_anom_terms.mat
ubtp47 = ubtp17 - ubtp13;
vbtp47 = ubtp17 - vbtp13;
wbtp47 = wbtp17 - wbtp13;
uptb47 = uptb17 - uptb13;
vptb47 = vptb17 - vptb13;
wptb47 = wptb17 - wptb13;
hc47 = hc17 - hc13;
dhc47dt = dhc17dt - dhc13dt;
figure(1); figure_orient(1);
subplot(2,1,1);
h2 = plot(...
nfrm, ubtp17, 'x-k', ...
nfrm, vbtp17, 'o-k', ...
nfrm, wbtp17, 'v-k', ...;
nfrm, uptb17, '*--k', ...
nfrm, vptb17, 's--k', ...
nfrm, wptb17, '^--k', ...
nfrm, -hflx, '>-.k', ...
nfrm, -dhc17dt, 'd-.k', ...
nfrm, -1e-8*hc17, '-k');%, ...
% nfrm, -40*taux, '--k');
set(h2(9), 'linewidth', 2);
grid on
xlabel('Lag (in years)');
ylabel('W m^-^2');
% title(['LP10: NP Cold Anomaly [170E 150W 27.5N 40N]'...
title(['HP10: NP Cold Anomaly [170E 155W 20N 32.5N]'...
' 0-270m Integrated Anomalous Heat Content'...
' Tendency Terms']);
axis([-12 12 -3 3]);
set(gca, 'XTick', -20:2:20, 'YTick', -5:.5:5);
legend(h2, 'Ubar dHC''/dx', 'Vbar dHC''/dy', 'Wbar dHC''dz',...
'U'' dHCbar/dx', 'V'' dHCbardy', 'W'' dHCbardz',...
'Net Heat Flux', 'dHC''/dt', 'HC (scaled)');%
% (1x10^8 J m^-^2)');%, 'TAUX (scaled)');
cd /home/disk/tao/dvimont/matlab/CSIRO/IUGG_talk/Plots
figure(1); figure_orient(1);
subplot(2,1,2);
h1 = plot(...
nfrm, ubtp47, 'x-k', ...
nfrm, vbtp47, 'o-k', ...
nfrm, wbtp47, 'v-k', ...;
nfrm, uptb47, '*--k', ...
nfrm, vptb47, 's--k', ...
nfrm, wptb47, '^--k', ...
nfrm, -dhc47dt, 'd-.k', ...
nfrm, -1.5e-8*hc47, '-k', ...;
nfrm, -.5e2*taux, '--k')
set(h1(8:9), 'linewidth', 2);
grid on
xlabel('Lag (in years)');
ylabel('W m^-^2');
title('LP10: EQ 80-270m Integrated Anomalous Heat Content Tendency Terms');
axis([-12 12 -1 2]);
set(gca, 'XTick', -20:2:20);
legend(h1, 'Ubar dHC''/dx', 'Vbar dHC''/dy', 'Wbar dHC''dz',...
'U'' dHCbar/dx', 'V'' dHCbardy', 'W'' dHCbardz',...
'dHC''/dt', 'HC (scaled)', 'CURL(TAU) (scaled)');
% 'dHC''/dt', 'HC (scaled)', 'TAUX (scaled)');
cd /home/disk/tao/dvimont/matlab/CSIRO/IUGG_talk/Plots
% Cleaner plots:
figure(1); figure_orient(1);
subplot(2,1,2);
h1 = plot(...
nfrm, (ubtp17+vbtp17+wbtp17'), 's-k', ...
nfrm, (uptb17+vptb17+wptb17'), '*--k', ...
nfrm, -hflx, '>-.k', ...
nfrm, -dhc17dt, 'd-.k', ...
nfrm, -1e-8*hc17, '-k');%, ...;
% nfrm, -.5e2*taux, '--k')
set(h1(5), 'linewidth', 2);
grid on
xlabel('Lag (in years)');
ylabel('W m^-^2');
% title(['LP10: NP Cold Anomaly [170E 150W 27.5N 40N]'...
% title(['HP10: NP Cold Anomaly [170E 155W 20N 32.5N]'...
title(['LP10: NP Cold Anomaly [170E 150W 27.5N 40N]'...
' 0-270m Integrated Anomalous Heat Content'...
' Tendency Terms']);
axis([-12 12 -2.5 2.5]);
set(gca, 'XTick', -20:2:20);
legend(h1, 'Ubar Grad(T'')', 'U'' Grad(Tbar)', 'Net Heat Flux', ...
'dHC''/dt', 'HC (scaled)');
cd /home/disk/tao/dvimont/matlab/CSIRO/IUGG_talk/Plots
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
save HP10_np_anom_terms.mat uptb13 ubtp13 vptb13 vbtp13 wptb13 wbtp13 ...
uptb17 ubtp17 vptb17 vbtp17 wptb17 wbtp17 ...
hc13 hc17 dhc13dt dhc17dt hflx taux lags lims tim nfrm
% hc13 hc17 dhc13dt dhc17dt hflx curlt lags lims tim nfrm
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_eq_terms.mat