Documentation of plot_zerolag_terms

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

Help text

load LP9_detrend_L1-7_EOF_yr101-1000.mat

Cross-Reference Information

Listing of script plot_zerolag_terms

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load HP10_hctend_L1-7_yr101-1000_lreg.mat
uptb17 = uptb;
ubtp17 = ubtp; 
vptb17 = vptb; 
vbtp17 = vbtp; 
wptb17 = wptb; 
wbtp17 = wbtp;
load HP10_hctend_L1-3_yr101-1000_lreg.mat
uptb13 = uptb;
ubtp13 = ubtp; 
vptb13 = vptb; 
vbtp13 = vbtp; 
wptb13 = wptb; 
wbtp13 = wbtp;
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP9_hctend_ALL.mat;  pn = 1;
load HP10_hctend_ALL.mat; pn = -1;
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

[lat, lon, depth, lm] = getll('temp', lims);

hflx = getnc('heat', lims, 1, tim);
[lath, lonh, depth] = getll('heat', lims);
hflx_reg = regress_eof(hflx, pcs, nfrm);

[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'];

%  Plot upper level heat content terms

get_global; FRAME = [110 300 -30 30];

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

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

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

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

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

subplot(7,2,13);
  xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2 std^-^1']);
subplot(7,2,12);
  xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2 std^-^1']);

subplot(7,2,13);
  XAX = lonh; YAX = lath;
  gcont(-pn*hflx_reg(lind,:,:), clev);
    XAX = lon; YAX = lat;
    dc2(lm);
  ylabel('Net HFLX')

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