Global Index (short | long) | Local contents | Local Index (short | long)
lims2 = [180 270 -3 3]; tit = 'EQ'; % Equatorial index lims2 = [150 210 -6 6]; tit = 'EQ'; % Equatorial index lims2 = [120 180 -10 10]; tit = 'WEQ'; % WSTP3 index
| This script calls | |
|---|---|
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/New_calcs2/data
load HP12_hctend_hpcs_new.mat; pn = -1;
load heat_pcs.mat; pcs = bpcs1540;
load LP7_hctend_CT_new.mat; pn = 1;
load HP6_hctend_CT_new.mat; pn = 1;
load HP12_hctend_CT_new.mat; pn = 1;
load BP1840_hctend_CT_new.mat; pn = 1;
load BP412_hctend_CT_new.mat; pn = 1;
load BP2030_hctend_CT_new.mat; pn = 1;
order = ['ubtp'; 'uptb'; 'vbtp'; 'vptb'; 'wbtp'; 'wptb'];
levind = ['13'; '47'];
lims2 = [180 270 -6 6]; tit = 'EQ'; % Equatorial index
for j = 1:2;
[lat, lon] = getll('u', lims);
[xk, yk] = keep_var(lims2, lon, lat);
for i = 1:4;
eval(['tem = ' order(i,:) levind(j,:) '(:,yk,xk);']);
if ~iselement(size(tem), 1);
tem = squeeze(mean2(mean2(shiftdim(tem, 1))));
else
tem = squeeze(mean2(squeeze(shiftdim(tem, 1))))';
end
eval([order(i,:) '_' levind(j,:) ' = tem;']);
end
[lat, lon] = getll('wl', lims);
[xk, yk] = keep_var(lims2, lon, lat);
for i = 5:6
eval(['tem = ' order(i,:) levind(j,:) '(:,yk,xk);']);
if ~iselement(size(tem), 1);
tem = squeeze(mean2(mean2(shiftdim(tem, 1))));
else
tem = squeeze(mean2(squeeze(shiftdim(tem, 1))));
end
eval([order(i,:) '_' levind(j,:) ' = tem;']);
end
[lat, lon] = getll('temp', lims);
[xk, yk] = keep_var(lims2, lon, lat);
eval(['tem = dhcdt' levind(j,:) '(:,yk,xk);']);
tem = squeeze(mean2(mean2(shiftdim(tem, 1))));
eval(['dhcdt_' levind(j,:) ' = tem;']);
end
[lat, lon] = getll('temp', lims);
[xk, yk] = keep_var(lims2, lon, lat);
hflx = squeeze(mean2(mean2(shiftdim(hflx_reg(:,yk,xk), 1))));
res_13 = -dhcdt_13 - (ubtp_13+vbtp_13+wbtp_13+uptb_13+vptb_13+wptb_13-hflx);
res_47 = -dhcdt_47 - (ubtp_47+vbtp_47+wbtp_47+uptb_47+vptb_47+wptb_47);
hc13 = getheat(lims2, 1:3, tim);
hc47 = getheat(lims2, 4:7, tim);
hc_13 = regress_eof(hc13, pcs, nfrm);
hc_47 = regress_eof(hc47, pcs, nfrm);
hc_13 = squeeze(mean2(mean2(shiftdim(hc_13, 1))));
hc_47 = squeeze(mean2(mean2(shiftdim(hc_47, 1))));
taux = getnc('taux', [150 210 -6 6], 1, tim);
taux = regress_eof(taux, pcs, nfrm);
taux = squeeze(mean2(mean2(shiftdim(taux, 1))));
figure(1); fo(1); clf;
sz = get(gcf, 'PaperSize');
hw = 1/sz(1); vw = 1/sz(2); hsz = hw*(2+7/8); vsz = vw*1.75;
bmarg = (sz(2) - 8.5)*0.5; tmarg = (sz(2) - 8.5)*0.5;
mid = 0.5;
subplot(2,1,1);
h1 = plot(...
nfrm, pn*dhcdt_13, 'd-.k', ...
nfrm, -1*pn*ubtp_13, 'x-k', ...
nfrm, -1*pn*vbtp_13, 'o-k', ...
nfrm, -1*pn*wbtp_13, 'v-k', ...
nfrm, -1*pn*uptb_13, '*--k', ...
nfrm, -1*pn*vptb_13, 'd--k', ...
nfrm, -1*pn*wptb_13, '^--k', ...
nfrm, pn*hflx, '>-.k', ...
nfrm, -1*pn*res_13, '.:k', ...
nfrm, pn*2.5e-8*hc_13, '-k', ...
nfrm, pn*60*taux, 'x--k');
set(h1(9:11), 'linewidth', 2);
set(gca, 'XTick', nfrm, 'YTick', -10:.5:10);
axis([min(nfrm) max(nfrm) get(gca, 'YLim')]);
grid on
subplot(2,1,2);
h1 = plot(...
nfrm, pn*dhcdt_47, 'd-.k', ...
nfrm, -1*pn*ubtp_47, 'x-k', ...
nfrm, -1*pn*vbtp_47, 'o-k', ...
nfrm, -1*pn*wbtp_47, 'v-k', ...
nfrm, -1*pn*uptb_47, '*--k', ...
nfrm, -1*pn*vptb_47, 'd--k', ...
nfrm, -1*pn*wptb_47, '^--k', ...
nfrm, -1*pn*res_47, '.:k', ...
nfrm, pn*2.5e-8*hc_47, '-k')
set(h1(8:9), 'linewidth', 2);
set(gca, 'XTick', nfrm, 'YTick', -10:.5:10);
axis([min(nfrm) max(nfrm) get(gca, 'YLim')]);
grid on
hc13 = getheat(lims2, 1:3, tim);
hc47 = getheat(lims2, 4:7, tim);
ct = getct; ct = detrend(ct);
[b, a] = butter(9, 2/15); [b2, a2] = butter(9, 2/40);
bpct = filtfilt(b, a, ct)-filtfilt(b2,a2,ct);
hc_13 = regress_eof(hc13, bpct, [-20:20]);
hc_47 = regress_eof(hc47, bpct, [-20:20]);
hc_13 = squeeze(mean2(mean2(shiftdim(hc_13, 1))));
hc_47 = squeeze(mean2(mean2(shiftdim(hc_47, 1))));
[latu, lonu] = getll('u', lims);
[latw, lonw] = getll('wl', lims);
[lat, lon, depth, lm] = getll('temp', lims);
llord = ['uuuuww'];
figure(2); fo(1); clf;
figure(3); fo(2); clf;
default_global; FRAME = [100 295 -60 60];
lind = find(nfrm == 0);
lev = 0.5;
for i = 1:6;
figure(2);
eval(['tem = -1*' order(i,:) '13;']);
eval(['XAX = lon' llord(i) '; YAX = lat' llord(i) ';']);
subplot(3,2,i); dg(latu, lonu);
gcont(tem(lind,:,:), lev);
dc2(lm, 0.7, 1000);
ylabel([order(i,:) '13']);
figure(3);
eval(['tem = -1*' order(i,:) '47;']);
eval(['XAX = lon' llord(i) '; YAX = lat' llord(i) ';']);
subplot(3,2,i); dg(latu, lonu);
gcont(tem(lind,:,:), lev);
dc2(lm, 0.7, 1000);
ylabel([order(i,:) '47']);
end
% Another way of looking: 6 figures
[latu, lonu] = getll('u', lims);
[latw, lonw] = getll('wl', lims);
[lat, lon, depth, lm] = getll('temp', lims);
llord = ['uuuuww'];
default_global; FRAME = [100 295 -20 20];
lind = find(ismember(lags, -2:2:6));
lev = 0.25;
for i = 6:6;
figure(1); fo(1); clf;
eval(['tem1 = -1*' order(i,:) '13(lind,:,:);']);
eval(['tem2 = -1*' order(i,:) '47(lind,:,:);']);
eval(['XAX = lon' llord(i) '; YAX = lat' llord(i) ';']);
for j = 1:5;
subplot(5,2,2*j-1);
gcont(pn*tem1(j,:,:), lev);
dc2(lm, 0.7, 1000);
ylabel([order(i,:) '13']);
subplot(5,2,2*j);
gcont(pn*tem2(j,:,:), lev);
dc2(lm, 0.7, 1000);
ylabel([order(i,:) '47']);
end
end