Global Index (short | long) | Local contents | Local Index (short | long)
Plot data
| This script calls | |
|---|---|
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/New_calcs/data
load heat_pcs.mat
lims = [80 310 -80 80];
tim = 101:1000;
[sst, taux] = getnc('temp', 'taux', lims, 1, tim);
[slp, rnd] = getflx('psl', 'rnd', lims, tim);
sst = detrend(sst);
corrtaux = detrend(taux);
slp = detrend(slp);
rnd = detrend(rnd);
tim = hpcs(:,1);
[b, a] = butter(9, 2/10);
sst2 = sst - filtfilt(b, a, sst);
taux2 = taux - filtfilt(b, a, taux);
slp2 = slp - filtfilt(b, a, slp);
rnd2 = rnd - filtfilt(b, a, rnd);
[pat3, c3] = regress_eof(sst2, tim);
[pat2, c2] = regress_eof(taux2, tim);
[pat1, c1] = regress_eof(slp2, tim);
[pat4, c4] = regress_eof(rnd2, tim);
[lat1, lon1, depth, lm] = getll('temp', lims); % sst, rnd, slp, lm
[lat2, lon2] = getll('taux', lims); % tau
tim = -lpcs(:,1);
[b, a] = butter(9, 2/9);
sst2 = filtfilt(b, a, sst);
taux2 = filtfilt(b, a, taux);
slp2 = filtfilt(b, a, slp);
rnd2 = filtfilt(b, a, rnd);
[pat3b, c3b] = regress_eof(sst2, tim);
[pat2b, c2b] = regress_eof(taux2, tim);
[pat1b, c1b] = regress_eof(slp2, tim);
[pat4b, c4b] = regress_eof(rnd2, tim);
cd ~/Papers/mlcsiro/matlab/Data2
save HPC1_reg_DYN_HC.mat pat1 pat2 pat3 pat4 c1 c2 c3 c4 ...
lat1 lon1 lat2 lon2 lm lims
save LPC1_reg_DYN_HC.mat pat1b pat2b pat3b pat4b c1b c2b c3b c4b ...
lat1 lon1 lat2 lon2 lm lims
clean
cd ~/Papers/mlcsiro/matlab/Data2
load HPC1_reg_DYN_HC.mat
load LPC1_reg_DYN_HC.mat
figure(1); fo(1); clf;
sz = get(gcf, 'PaperSize');
hw = 1/sz(1); vw = 1/sz(2); hsz = hw*(2.5); vsz = vw*1.8;
bmarg = (sz(2) - 8.5)*0.5; tmarg = (sz(2) - 8.5)*0.5;
mid = 0.45;
[tem, tem, tem, lm] = getll('temp', lims);
global GRDX_SPACING GRDY_SPACING FRAME
GRDX_SPACING = 45; GRDY_SPACING = 20;
cint = [0.15 0.015 0.05 0.075];
FRAME = [95 300 -75 75];
% MIX run
dg2(lat1, lon1);
subplot('position', [0.5-hsz-0.04 ...
vw*(sz(2)-tmarg)-vsz hsz vsz]);
gcont(pat1, cint(1));
dc2(lm, 0.3, -1000);
color_shade(squeeze(c1.^2), 0.25, 0.7*[1 1 1])
set(gca, 'fontsize', 8);
t1 = title(['< HPC1, SLP >']);
set(t1, 'units', 'normalized', 'HorizontalAlignment', 'left', ...
'VerticalAlignment', 'bottom', 'fontsize', 9);
pos2 = get(t1, 'position');
set(t1, 'position', [0 pos2(2)-0.01 0]);
tx1 = text(mean(FRAME(1:2)), 110, 'HPC1 of 0-270m Heat Content');
set(tx1, 'fontsize', 10, 'HorizontalAlignment', 'center');
dg2(lat2, lon2);
subplot('position', [0.5-hsz-0.04 ...
vw*(sz(2)-tmarg-mid)-2*vsz hsz vsz]);
gcont(pat2, cint(2));
dg2(lat1, lon1);
dc2(lm, 0.3, 1000);
dg2(lat2, lon2);
color_shade(squeeze(c2.^2), 0.25, 0.7*[1 1 1])
set(gca, 'fontsize', 8);
t1 = title(['< HPC1, TAUX >']);
set(t1, 'units', 'normalized', 'HorizontalAlignment', 'left', ...
'VerticalAlignment', 'bottom', 'fontsize', 9);
pos2 = get(t1, 'position');
set(t1, 'position', [0 pos2(2)-0.01 0]);
dg2(lat1, lon1);
subplot('position', [0.5-hsz-0.04 ...
vw*(sz(2)-tmarg-2*mid)-3*vsz hsz vsz]);
gcont(pat3, cint(3));
dc2(lm, 0.3, 1000);
color_shade(squeeze(c3.^2), 0.25, 0.7*[1 1 1])
set(gca, 'fontsize', 8);
t1 = title(['< HPC1, SST >']);
set(t1, 'units', 'normalized', 'HorizontalAlignment', 'left', ...
'VerticalAlignment', 'bottom', 'fontsize', 9);
pos2 = get(t1, 'position');
set(t1, 'position', [0 pos2(2)-0.01 0]);
subplot('position', [0.5-hsz-0.04 ...
vw*(sz(2)-tmarg-3*mid)-4*vsz hsz vsz]);
gcont(pat4, cint(4));
dc2(lm, 0.3, -1000);
color_shade(squeeze(c4.^2), 0.25, 0.7*[1 1 1])
set(gca, 'fontsize', 8);
t1 = title(['< HPC1, PREC >']);
set(t1, 'units', 'normalized', 'HorizontalAlignment', 'left', ...
'VerticalAlignment', 'bottom', 'fontsize', 9);
pos2 = get(t1, 'position');
set(t1, 'position', [0 pos2(2)-0.01 0]);
% DYN run
dg2(lat1, lon1);
subplot('position', [0.5+0.04 ...
vw*(sz(2)-tmarg)-vsz hsz vsz]);
gcont(pat1b, cint(1));
dc2(lm, 0.3, -1000);
color_shade(squeeze(c1b.^2), 0.25, 0.7*[1 1 1])
set(gca, 'fontsize', 8);
t1 = title(['< LPC1, SLP >']);
set(t1, 'units', 'normalized', 'HorizontalAlignment', 'left', ...
'VerticalAlignment', 'bottom', 'fontsize', 9);
pos2 = get(t1, 'position');
set(t1, 'position', [0 pos2(2)-0.01 0]);
tx1 = text(mean(FRAME(1:2)), 110, 'LPC1 of 0-270m Heat Content');
set(tx1, 'fontsize', 10, 'HorizontalAlignment', 'center');
subplot('position', [0.5+0.04 ...
vw*(sz(2)-tmarg-mid)-2*vsz hsz vsz]);
dg2(lat2, lon2);
gcont(pat2b, cint(2));
dg2(lat1, lon1);
dc2(lm, 0.3, 1000);
dg2(lat2, lon2);
color_shade(squeeze(c2b.^2), 0.25, 0.7*[1 1 1])
set(gca, 'fontsize', 8);
t1 = title(['< LPC1, TAUX >']);
set(t1, 'units', 'normalized', 'HorizontalAlignment', 'left', ...
'VerticalAlignment', 'bottom', 'fontsize', 9);
pos2 = get(t1, 'position');
set(t1, 'position', [0 pos2(2)-0.01 0]);
dg2(lat1, lon1);
subplot('position', [0.5+0.04 ...
vw*(sz(2)-tmarg-2*mid)-3*vsz hsz vsz]);
gcont(pat3b, cint(3));
dc2(lm, 0.3, 1000);
color_shade(squeeze(c3b.^2), 0.25, 0.7*[1 1 1])
set(gca, 'fontsize', 8);
t1 = title(['< LPC1, SST >']);
set(t1, 'units', 'normalized', 'HorizontalAlignment', 'left', ...
'VerticalAlignment', 'bottom', 'fontsize', 9);
pos2 = get(t1, 'position');
set(t1, 'position', [0 pos2(2)-0.01 0]);
subplot('position', [0.5+0.04 ...
vw*(sz(2)-tmarg-3*mid)-4*vsz hsz vsz]);
gcont(pat4b, cint(4));
dc2(lm, 0.3, -1000);
color_shade(squeeze(c4b.^2), 0.25, 0.7*[1 1 1])
set(gca, 'fontsize', 8);
t1 = title(['< LPC1, PREC >']);
set(t1, 'units', 'normalized', 'HorizontalAlignment', 'left', ...
'VerticalAlignment', 'bottom', 'fontsize', 9);
pos2 = get(t1, 'position');
set(t1, 'position', [0 pos2(2)-0.01 0]);
cd /home/disk/tao/dvimont/matlab/CSIRO/Tony/Figs_28June2000
print -dps2 HP_LP_PC1_regressions.ps