Global Index (short | long) | Local contents | Local Index (short | long)
load LP10_detrend_L1-7_yr101-1000.mat; pcs3 = pcs; load HP9_detrend_L1-7_EOF_yr101-1000.mat; pcs3 = -1*pcs;
This script calls | |
---|---|
clear cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data load RAW_detrend_L1-7_EOF_yr101-1000.mat; pcs3 = pcs; lims = [-0.1 360 -90 90]; tim = 101:1000; lev = 1:7; [heat, lath, lonh, depth, middepth] = getheat(lims, lev, tim); temp = getnc('temp', lims, 1, tim); [latt, lont] = 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); colr = 0.5; pn = 1; default_global; FRAME = [0 360 -90 90]; XAX = lonh; YAX = lath; figure(1); figure_orient; cint = 0.4; clev = [-10:cint:-cint cint:cint:10]; %subplot(2,1,2); spthes(2); gcont(1e-8*hreg, clev); dc2(hreg, colr, pn); color_shade(squeeze(hcoef.^2), 0.25, [0.87 0.87 0.87]); title(['< PC1 of HC270, HC270 >']); % title(['HC270 anomaly regressed on HP10 PC1']);% -- '... % num2str(round(per(num))) ... % '% Variance Explained; (Lamda = ' ... % num2str(round(1e-16*lam(num))/100) ')']); % xlabel(['Contour Interval: ' num2str(10*cint) 'x10^7 J m^-^2 std^-^1']); % ylabel(['EOF over [110E 60W 60S 60N]']) h = line([110 300 300 110 110], [60 60 -60 -60 60], [100 100 100 100 100]); set(h, 'linewidth', 2, 'linestyle', '-', 'Color', [0 0 0]); default_global; FRAME = [0 360 -90 90]; XAX = lont; YAX = latt; cint = 0.05; clev = [-1:cint:-cint cint:cint:1]; %subplot(2,1,1); spthes(1); gcont(treg, clev); dc2(treg, colr, pn); title(['< PC1 of HC270, SST >']) % title(['12.5m Temperature Regressed on HP10 PC' num2str(num)]); % xlabel(['Contour Interval: ' num2str(cint) ' K std^-^1']); % dc2(treg) color_shade(squeeze(tcoef.^2), 0.25, [.87 .87 .87]); % xlabel(['Contour Interval: ' num2str(cint) ' K 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);