Documentation of HP10_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 HP10_spatial_tendterms

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load HP10_hctend_ALL.mat; pn = -1;
load HP10_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 HP10_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); figure_landscape(1);
get_global; FRAME = [105 295 -60 60];

lag_plot = -1:1;
cint = 0.5; clev = -8:cint:8;

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);
%      tem(abs(tem)<cint) = NaN;
      gcont(tem, cint);
%      mcolor(tem, cint, 'mollweid');
      XAX = lon; YAX = lat;
      h=dc2(lm, 0.5, 1);
%      drawmapm
      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);
    caxis([-3 3]);
  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);
%      tem(abs(tem)<cint) = NaN;
      gcont(tem, cint);
%      mcolor(tem, cint, 'mollweid');
      XAX = lon; YAX = lat;
      dc2(lm);
%      drawmapm
      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 HP10_0-80m_spatial_terms.ps
figure(2); print -dps2 HP10_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:4;
  if i ~= 4;
    eval(['tem = pn*' order(i,:) topbot '(lind, :, :);']);
    eval(['XAX = lon' llord(i) '; YAX = lat' llord(i) ';']);
  else
    tem = hflx_reg(lind,:,:); XAX = lon3; YAX = lat3;
  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 ~= 4;
    yl(i) = ylabel([tit2 ' ' tit(i,:)]);
  else
    yl(i) = ylabel(['SFC HFLX']);
  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(5) = ylabel(['80-270m ' tit(1,:)]);
  set(gca, 'FontSize', 9);

sptalk(5,2,1); 
t(1) = title('HP10 0-80m HC terms, Lag 0');
sptalk(5,2,2); 
t(2) = title('HP10 80-170m HC terms, Lag 0');
xl(1) = xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2 std^-^2']);
sptalk(5,2,7);
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 = -1:1;
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(['HP10:  ' 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 HP10_0-80m_spatial_terms.ps
figure(2); print -dps2 HP10_80-270m_spatial_terms.ps