Documentation of plot_reg_HEAT_SFC

Global Index (short | long) | Local contents | Local Index (short | long)

Help text

load LP10_detrend_L1-7_yr101-1000.mat; pcs3 = pcs;
load HP9_detrend_L1-7_EOF_yr101-1000.mat; pcs3 = -1*pcs;

Cross-Reference Information

Listing of script plot_reg_HEAT_SFC

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);