Documentation of plot_pac_regressions

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

Help text

11111111111111111111111111111111111111  Plot the data

Cross-Reference Information

Listing of script plot_pac_regressions

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

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

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

    pos = [.13+.2811*(j-4) (9.25/11)-(i)*(9.25-2.2)/(3*11) ...
           .2128 .18]; 

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

subplot(3,3,1);
  title('HP10:  0-80m Heat Content');
subplot(3,3,2);
  title('HP10:  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/Chap3
%  print -dps2 HC0-80m_80-270m_lag-1to1.ps




%2222222222222222222222222222222222222222222  Plot the data

default_global;  FRAME = [105 300 -60 60];
[lat1, lon1] = getll('taux', lims);

for i = 1:3;
  subplot(3,2,2*i-1);
  xdat = squeeze(treg(i,:,:)); ydat = squeeze(breg(i,:,:));
  XAX = lon1; YAX = lat1;
  gquiv(xdat, ydat, -0.08, 2, 'dynes cm^-^2 std^-^1');
  XAX = lon; YAX = lat;
  dc2(lm);
%  color_shade(squeeze(tcoef(i,:,:).^2), 0.25, [0.9 0.9 0.9]);
  ylabel(['Lag = ' num2str(lags(i))]);
  if i ~= 3;
    xlabel('');
  end
end

[lat1, lon1] = getll('z250', [-0.1 360 -5 90]);
FRAME = [0 360 15 90];
cint = 2; clev = [-20:cint:-cint cint:cint:20];
for i = 1:3;
  subplot(3,2,2*i);
  XAX = lon1; YAX = lat1;
  mcont(squeeze(treg(i,:,:)), clev);
  axis(1.6 * [-1 1 -1 1])
end

subplot(3,2,1);
  title('HP10:  Wind Stress Anomalies');
subplot(3,2,2);
  title('HP10:  250mb Height Anomalies');

subplot(3,2,6);
  xlabel(['Contour Interval:  ' num2str(cint) 'm std^-^1']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap2/Plots
%  print -dps2 HC0-80m_80-270m_lag-1to1.ps