Documentation of plot_ceof_heat_sft

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

Help text

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

Cross-Reference Information

Listing of script plot_ceof_heat_sft

clear
varn = 'heat';
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat
lev = 4:7;
tim = 101:1000;
lims = [105 300 -65 65];
lags = -8:5; nfrm = lags;
temp = getnc(varn, lims, 1, tim);
[lat, lon, depth] = getll(varn, lims);

treg = regress_eof(temp, pcs, nfrm);

%  Plot the various phases

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

cint = 1; clev = [-15:cint:-cint cint:cint:15];
for i = 1:7;
  subplot(7,2,2*i-1);
    gcont(-1e-0*treg(i,:,:), clev);
    dc2(treg(i,:,:));
    ylabel(['Lag = ' num2str(lags(i))]);
  subplot(7,2,2*i);
    gcont(-1e-0*treg(i+7,:,:), clev);
    dc2(treg(i+7,:,:));
    ylabel(['Lag = ' num2str(lags(i+7))]);
end
for i = 13:14
  subplot(7,2,i);
  cint2 = 1e0*cint;
%    xlabel(['Contour Interval = ' num2str(cint2) ' K std^-^1']);
%    xlabel(['Contour Interval = ' num2str(cint2) 'x10^7 J m^-^2 std^-^1']);
    xlabel(['Contour Interval = ' num2str(cint2) 'x10^-^4 N m^-^2 std^-^1']);
end

for i = 1:2;
subplot(7,2,i);
  title(['LP10:  < HC80-270, PC1 >  '  num2str(round(per(1))) ...
          '% Var Expl.']);
  title('LP10:  < 12.5m Temp, PC1 >');
  title('LP10:  < TAUX, PC1 >');
end

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap3/Plots
%print -dps2 HC270m_PC1_EQ.ps