Documentation of plot_lp_pac_regressions

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

Help text

  Plot the data

Cross-Reference Information

Listing of script plot_lp_pac_regressions

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP9_detrend_L1-7_EOF_yr101-1000.mat
cd ~/matlab/CSIRO/Heat/Old_routines
nfrm = [-4:3:2]; lags = nfrm;
tim = 101:1000;
lims = [100 305 -65 65];
top = getheat(lims, 1:3, tim);
bot = getheat(lims, 4:7, tim);
[lat, lon, depth, lm] = getll('temp', lims);
top = detrend(top); bot = detrend(bot);
[b, a] = butter(9, 2/9);
top = filtfilt(b, a, top);
bot = filtfilt(b, a, bot);
[treg, tcoef] = regress_eof(top, pcs, nfrm);
[breg, bcoef] = regress_eof(bot, pcs, nfrm);

figure(1); fo;
default_global;  FRAME = [105 300 -60 60];

cintt = 0.15; clevt = [-2:cintt:-cintt  cintt:cintt:2];
cintb = 0.3; clevb = [-5:cintb:-cintb  cintb:cintb:5];

for i = 1:3;
%  subplot(3,3,3*i-2);
  pos = [.2 (9.8/11)-i*7/(3*11) .3 .18];
  subplot('position', pos);
  gcont(-1e-8*(treg(i,:,:)), cintt);
  dc2(tcoef(i,:,:));
  color_shade(squeeze(tcoef(i,:,:).^2), 0.25, [0.9 0.9 0.9]);
  ylabel(['Lag = ' num2str(lags(i))]);
  if i == 1; 
    title('LP9:  0-80m HC');
  end
end
for i = 1:3;
%  subplot(3,3,3*i-1);
  pos = [.58 (9.8/11)-i*7/(3*11) .3 .18];
  subplot('position', pos);
  gcont(-1e-8*breg(i,:,:), cintb);
  dc2(bcoef(i,:,:));
  color_shade(squeeze(bcoef(i,:,:).^2), 0.25, [0.9 0.9 0.9]);
  if i == 1; 
    title('LP9:  80-270m HC');
  end
end

subplot(3,3,1);
  title('LP9:  0-80m Heat Content');
subplot(3,3,2);
  title('LP9:  80-270m Heat Content');

subplot(3,3,7);
  xlabel(['Contour Interval:  ' num2str(10*cintt) 'x10^7 J m^-^2 std^-^1']);
subplot(3,3,8);
  xlabel(['Contour Interval:  ' num2str(10*cintb) 'x10^7 J m^-^2 std^-^1']);

cd ~/Thesis/Chap4
%  print -dps2 HC0-80m_80-270m_lag-4to2.ps