Documentation of plot_ceof_heat_sft

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

Help text

[temp, lat, lon, depth, middepth] = getheat(lev, tim, lims);

Cross-Reference Information

Listing of script plot_ceof_heat_sft

clear
varn = 'z250';
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load HP10_detrend_L1-7_CEOF_yr101-1000.mat
lev = 1:7;
tim = 101:1000;
lims = [105 300 -65 65];
nfrm = 6;  lags = 180/nfrm * [0:nfrm-1];
temp = getnc(varn, lims, 1, tim);
[lat, lon, depth] = getll(varn, lims);

treg = regress_ceof(temp, pcs, nfrm);

%  Plot the various phases

default_global; FRAME = [105 295 -60 60];

cint = 1; clev = [-10:cint:-cint cint:cint:10];
for i = 1:3;
  subplot(3,2,2*i-1);
    gcont(1e2*treg(2*i-1,:,:), clev);
    dc2(treg(2*i-1,:,:));
    title(['Lag = ' num2str(lags(2*i-1))]);
  subplot(3,2,2*i);
    gcont(1e2*treg(2*i,:,:), clev);
    dc2(treg(2*i,:,:));
    title(['Lag = ' num2str(lags(2*i))]);
end
for i = 5:6;
  subplot(3,2,i);
%  xlabel(['Contour Interval = ' num2str(cint) ' K std^-^1']);
%  xlabel(['Contour Interval = ' num2str(10*cint) 'x10^7 J m^-^2 std^-^1']);
  xlabel(['Contour Interval = ' num2str(cint) 'x10^-^3 N m^-^2 std^-^1']);
end

subplot(3,2,3);
  ylabel(['HP10 Complex EOF of 0-270m Heat Content;  '  num2str(round(per(1))) ...
          '% Variance Explained']);
  ylabel('12.5m Temperature Regressed on CPC1 of 0-270m Heat Content');
  ylabel('Zonal Wind Stress Regressed on CPC1 of 0-270m Heat Content');


cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap2/Plots
%print -dps2 HC0-270m_CEOF.ps