Documentation of plot_lp_pac_terms

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

Help text

%%%%%%%%%%%%%%%%  Plot dynamic heat content terms

Cross-Reference Information

Listing of script plot_lp_pac_terms

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP9_hctend_ALL.mat; pn = 1;
load LP9_hflx_terms.mat;
prec = getflx('rnd', lims, tim);
uptb47 = uptb17 - uptb13;
ubtp47 = ubtp17 - ubtp13; 
vptb47 = vptb17 - vptb13; 
vbtp47 = vbtp17 - vbtp13; 
wptb47 = wptb17 - wptb13; 
wbtp47 = wbtp17 - wbtp13;
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP9_detrend_L1-7_EOF_yr101-1000.mat
prec_reg = regress_eof(prec, pcs, nfrm);
[lat, lon, depth, lm] = getll('temp', lims);
[lat1, lon1] = getll('u', lims);
[lat2, lon2] = getll('wl', lims);
order = ['ubtp'; 'vbtp'; 'wbtp'; 'uptb'; 'vptb'; 'wptb'];
llord = num2str([1; 1; 2; 1; 1; 2]);
tit = ['Ubar T''_X'; 'Vbar T''_Y'; 'Wbar T''_Z'; ...
       'U'' Tbar_X'; 'V'' Tbar_Y'; 'W'' Tbar_Z'];

get_global; FRAME = [105 300 -60 60];
figure(1); figure_orient;

%  Use lags  +/- 3 and 0?

lag_plot = -3;
lind = find(lags == lag_plot);

cint = 0.5; clev = [-10:cint:-cint cint:cint:10];

topbot = '47';
for i = 1:3;
  eval(['tem = pn*' order(i,:) topbot '(lind, :, :);']);
  eval(['XAX = lon' llord(i) '; YAX = lat' llord(i) ';']);
  subplot(3,2,2*i-1);
    gcont(tem, clev);
    XAX = lon; YAX = lat;
    dc2(lm);
    yl(i) = ylabel(tit(i,:));
end
for i = 1:3;
  eval(['tem = pn*' order(i+3,:) topbot '(lind, :, :);']);
  eval(['XAX = lon' llord(i+3) '; YAX = lat' llord(i+3) ';']);
  subplot(3,2,2*i);
    gcont(tem, clev);
    XAX = lon; YAX = lat;
    dc2(lm);
    yl(i+3) = ylabel(tit(i+3,:));
end

subplot(3,2,1);
  title(['LP9:  80-270m HC Tendency Terms; Lag = ' num2str(lag_plot)]);
subplot(3,2,2);
  title(['LP9:  80-270m HC Tendency Terms; Lag = ' num2str(lag_plot)]);

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

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap3/Plots
 print -dps2 Lag-3_80-270m_tend_terms.ps

%%%%%%%%%%%%%%%%%  Plot Heat Flux

figure(2); figure_orient(1);
get_global; FRAME = [105 300 -60 60];

order = ['evp'; 'sgn'; 'hfl'; 'rgn'];
tit = ['LH Flux'; 'SW Flux'; 'SH Flux'; 'LW Flux'];

lag_plot = 3;
lind = find(lags == lag_plot);

cint = 0.5; clev = [-10:cint:-cint cint:cint:10];
for i = 1:4;
  eval(['tem = pn*' order(i,:) '_reg(lind, :, :);']);
  XAX = lonh; YAX = lath;
  subplot(3,2,i+2);
    gcont(tem, clev);
    XAX = lon; YAX = lat;
    dc2(lm);
    yl(i) = ylabel(tit(i,:));
end
subplot(3,2,1);
  tem = -1*hflx_reg(lind,:,:);
  XAX = lonh; YAX = lath;
  gcont(tem, clev);
  XAX = lon; YAX = lat;
  dc2(lm);
  ylabel('Net HFLX');

cint = 0.1; clev = [-1:cint:-cint cint:cint:1];
subplot(3,2,2);
  tem = -1*prec_reg(lind,:,:);
  XAX = lonh; YAX = lath;
  gcont(tem, clev);
  XAX = lon; YAX = lat;
  dc2(lm);
  ylabel('Precipitation');

subplot(3,2,1);
  title(['LP9:  HFLX and PREC Regressions; Lag = ' num2str(lag_plot)]);
subplot(3,2,2);
  title(['LP9:  HFLX and PREC Regressions; Lag = ' num2str(lag_plot)]);

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