Global Index (short | long) | Local contents | Local Index (short | long)
num2str(round(per(num))) ...
'% Variance Explained; (Lamda = ' ...
num2str(round(1e-16*lam(num))/100) ')']);
| This script calls | |
|---|---|
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load RAW_detrend_L1-7_EOF_yr101-1000.mat; pcs3 = pcs;
load LP9_detrend_L1-7_EOF_yr101-1000.mat; pcs3 = pcs;
load HP10_detrend_L1-7_EOF_yr101-1000.mat; pcs3 = -1*pcs;
lims = [-0.1 360 -90 90];
tim = 101:1000;
lev = 1;
heat = getnc('taux', lims, 1, tim);
[lath, lonh] = getll('taux', lims);
temp = getnc('z250', lims, 1, tim);
[latt, lont] = getll('z250', lims);
lm = getnc('temp', lims, 1, 1);
[latlm, lonlm] = getll('temp', lims);
[ntim, nlat, nlon] = size(heat);
heat = reshape(heat, ntim, nlat*nlon);
temp = reshape(temp, ntim, nlat*nlon);
temp = detrend(temp);
heat = detrend(heat);
temp = reshape(temp, ntim, nlat, nlon);
heat = reshape(heat, ntim, nlat, nlon);
num = 1;
[hreg, hcoef] = regress_eof(heat, -1*pcs3(:,num), 0);
[treg, tcoef] = regress_eof(temp, -1*pcs3(:,num), 0);
default_global; FRAME = [0 360 -90 90]; XAX = lonh; YAX = lath;
figure(1); figure_orient;
cint = 1; clev = [-10:cint:-cint cint:cint:10];
subplot(2,1,2);
gcont(1e2*hreg, clev);
XAX = lonlm; YAX = latlm;
dc2(lm)
XAX = lonh; YAX = lath;
color_shade(squeeze(hcoef.^2), 0.25, [.8 .8 .8]);
title(['Zonal Wind Sress anomaly regressed on RAW PC1']);% -- '...
xlabel(['Contour Interval: ' num2str(cint) 'x10^-^3 N m^-^2 std^-^1']);
% ylabel(['EOF over [110E 60W 60S 60N]'])
default_global; FRAME = [0 360 -90 90]; XAX = lont; YAX = latt;
cint = 2; clev = [-20:cint:-cint cint:cint:20];
subplot(2,2,1); FRAME = [0 360 15 90];
mcont(treg, clev, 'stereo');
greyshd(squeeze(tcoef.^2), 0.25, [.8 .8 .8])
axis(1.66*[-1 1 -1 1])
title(['NH 250mb Height Regressed on RAW PC' num2str(num)]);
xlabel(['Contour Interval: ' num2str(cint) ' m std^-^1']);
subplot(2,2,2); FRAME = [0 360 -90 -15];
mcont(treg, clev, [-90 270]);
greyshd(squeeze(tcoef.^2), 0.25, [.8 .8 .8])
setm(gca, 'FLatLimit', [15 90], 'MapLatLimit', [-15 -90]);
axis(1.66*[-1 1 -1 1])
title(['SH 250mb Height Regressed on RAW PC' num2str(num)]);
xlabel(['Contour Interval: ' num2str(cint) ' m std^-^1']);
cd /home/disk/tao/dvimont/Thesis/Chap2
figure(2); figure_orient;
cint = 0.1; clev = [-1:cint:-cint cint:cint:10];
sp(1);
gcont(hcoef, clev);
dc2(hcoef)
title(['Correlation Coefficient between Unfiltered PC' num2str(num) ...
' and Unfiltered 0-270m Heat Content'])
xlabel(['Contour Interval: ' num2str(cint) 'x10^8 J m^-^2 std^-^1']);
ylabel(['EOF over [110E 60W 60S 60N]'])
cint = 0.1; clev = [-1:cint:-cint cint:cint:10];
sp(2);
gcont(tcoef, clev);
dc2(tcoef)
title(['Correlation Coefficient between Unfiltered PC' num2str(num) ...
' and Unfiltered 12.5m Temperature'])
xlabel(['Contour Interval: ' num2str(cint) ' K std^-^1']);
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap1/Plots
figure(1);
figure(2);