Documentation of LP9_spatial_tendterms

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

Help text

prec = getflx('rnd', lims, tim);

Cross-Reference Information

Listing of script LP9_spatial_tendterms

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP9_hctend_ALL.mat; pn = 1;
load LP9_hflx_terms.mat;

uptb17 = uptb47 + uptb13;
ubtp17 = ubtp47 + ubtp13; 
vptb17 = vptb47 + vptb13; 
vbtp17 = vbtp47 + vbtp13; 
wptb17 = wptb47 + wptb13; 
wbtp17 = wbtp47 + 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);
[lat3, lon3] = getll('heat', lims);

order = ['wbtp'; 'uptb'; 'wptb'];
llord = num2str([2; 1; 2]);
tit = ['WT''_z'; 'u''T_X'; 'w''T_Z']



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

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

lag_plot = -4:3:2;

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

topbot = '13'; tit2 = ['0-80m'];
figure(1); fl(1); clf
for j = 1:3;
  for i = 1:3;
    lind = find(lags == lag_plot(i));
    eval(['tem = pn*' order(j,:) topbot '(lind, :, :);']);
    eval(['XAX = lon' llord(j) '; YAX = lat' llord(j) ';']);
    sptalk(3,3,3*(i-1)+j);
    set(gca, 'FontSize', 10);
      gcont(tem, cint);
      XAX = lon; YAX = lat;
      dc2(lm);
      if j == 1;
        yl(i) = ylabel(['Lag = ' num2str(lags(lind))]);
      end
      if i == 1;
        title(['HP10:  ' tit2 ' ' tit(j,:)]);
      elseif i == 3;
        xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2 std^-^2']);
      end
    set(gca, 'FontSize', 9);
  end
end


topbot = '47'; tit2 = ['80-270m'];
figure(2); fl(1); clf
for j = 1:3;
  for i = 1:3;
    lind = find(lags == lag_plot(i));
    eval(['tem = pn*' order(j,:) topbot '(lind, :, :);']);
    eval(['XAX = lon' llord(j) '; YAX = lat' llord(j) ';']);
    sptalk(3,3,3*(i-1)+j);
    set(gca, 'FontSize', 10);
      gcont(tem, cint);
      XAX = lon; YAX = lat;
      dc2(lm);
      if j == 1;
        yl(i) = ylabel(['Lag = ' num2str(lags(lind))]);
      end
      if i == 1;
        title(['HP10:  ' tit2 ' ' tit(j,:)]);
      elseif i == 3;
        xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2 std^-^2']);
      end
    set(gca, 'FontSize', 9);
  end
end


cd ~/Thesis/Talk
figure(1);print -dps2 LP9_0-80m_spatial_terms.ps
figure(2); print -dps2 LP9_80-270m_spatial_terms.ps





%  Plot upper level lag 0 terms

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

lag_plot = 0; lind = find(lags == lag_plot);
cint = 0.5; clev = -8:cint:8;
topbot = '13'; tit2 = ['0-80m'];

figure(1); fo(1); clf
for i = 1:5;
  if i == 5;
    tem = -1*hflx_reg(lind,:,:); XAX = lon3; YAX = lat3;
  elseif i == 2;
    tem = pn*vbtp13(lind,:,:); XAX = lon1; YAX = lat1;
  elseif i == 1;
    tem = pn*wbtp13(lind,:,:); XAX = lon3; YAX = lat3;
  else
    eval(['tem = pn*' order(i-1,:) topbot '(lind, :, :);']);
    eval(['XAX = lon' llord(i-1) '; YAX = lat' llord(i-1) ';']);
  end
  sptalk(5,2,2*i-1);
  set(gca, 'FontSize', 10);
      gcont(tem, cint);
      XAX = lon; YAX = lat;
      h=dc2(lm, 0.5, 1);
  if i == 5;
    yl(i) = ylabel(['SFC HFLX']);
  elseif i == 2;
    yl(i) = ylabel([tit2 ' VT''_z']);
  elseif i == 1
    yl(i) = ylabel([tit2 ' ' tit(i,:)]);
  else
    yl(i) = ylabel([tit2 ' ' tit(i-1,:)]);
  end
  set(gca, 'FontSize', 9);
end

sptalk(5,2,2);
  tem = pn*wbtp47(lind,:,:); XAX = lon2; YAX = lat2;
  gcont(tem, cint);
  XAX = lon; YAX = lat;
  h=dc2(lm, 0.5, 1);
  yl(6) = ylabel(['80-270m ' tit(1,:)]);
  set(gca, 'FontSize', 9);

sptalk(5,2,1); 
t(1) = title('LP9 0-80m HC terms, Lag 0');
sptalk(5,2,2); 
t(2) = title('LP9 80-170m HC terms, Lag 0');
xl(1) = xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2 std^-^2']);
sptalk(5,2,9);
xl(2) = xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2 std^-^2']);

set(t, 'fontsize', 10);
set(yl, 'fontsize', 11);
set(xl, 'fontsize', 10);


lag_plot = -4:3:2;
topbot = '47'; tit2 = ['80-270m'];
figure(2); fo(1); clf
for j = 2:3;
  for i = 1:3;
    lind = find(lags == lag_plot(i));
    eval(['tem = pn*' order(j,:) topbot '(lind, :, :);']);
    eval(['XAX = lon' llord(j) '; YAX = lat' llord(j) ';']);
    sptalk(5,2,(2*(i-1)+j-1));
    set(gca, 'FontSize', 10);
%      tem(abs(tem)<cint) = NaN;
      gcont(tem, cint);
%      mcolor(tem, cint, 'mollweid');
      XAX = lon; YAX = lat;
      dc2(lm);
%      drawmapm
      if j == 2;
        yl(i) = ylabel(['Lag = ' num2str(lags(lind))]);
      end
      if i == 1;
        title(['LP9:  ' tit2 ' ' tit(j,:)]);
      elseif i == 3;
        xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2 std^-^2']);
      end
    set(gca, 'FontSize', 9);
  end
end
set(yl, 'fontsize', 11);

cd ~/Thesis/Talk
figure(1);  print -dps2 LP9_0-80m_spatial_terms.ps
figure(2); print -dps2 LP9_80-270m_spatial_terms.ps