Documentation of plot_tendencies

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

Help text

  Get pcs

Cross-Reference Information

Listing of script plot_tendencies

clear

biff2 = 1;
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
if biff2 == 1;
  load BP10-60yr_L1-7_EOF.mat; tit = 'Bandpass Filtered Data ( 10 - 60 Years )';
elseif biff2 == 2;
  load HP8_L1-7_EOF.mat; tit = 'Highpass Filtered Data ( < 8 Years )';
elseif biff2 == 3;
  load RAW_L1-7_EOF.mat; tit = 'Unfiltered Data';
end

lags = [-5:5];
tim = 101:550;
lev = 5;
lims = [105 305 -63 63];


%  Get udtdx terms

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Lagged_regressions
   [ubtp, uptb, lat, lon, depth_out] = ...
        Lagged_regressions/get_udtdx.html">get_udtdx(pcs, lims, lags, lev, tim);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Lagged_regressions
   [ubtp, uptb, lat, lon, depth_out] = ...
        Lagged_regressions/get_vdtdy.html">get_vdtdy(pcs, lims, lags, lev, tim);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Lagged_regressions
   [ubtp, uptb, lat, lon, depth_out] = ...
        Lagged_regressions/get_wdtdz.html">get_wdtdz(pcs, lims, lags, lev, tim);

%  Get udHCdx terms

lev = 1:7;
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Lagged_regressions
   [ubtp, uptb, lat, lon, depth_out] = ...
           get_ubar_gradtprime(pcs, lims, lags, tim, lev);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Lagged_regressions
   [ubtp, uptb, lat, lon, depth_out] = ...
           get_uprime_gradtbar(pcs, lims, lags, tim, lev);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Lagged_regressions
   [ubtp, uptb, lat, lon, depth_out] = ...
           get_wbar_dtprimedz(pcs, lims, lags, tim, lev);


%  Plot data

default_global;  FRAME = [110 299 -60 60]; 

tem1 = -1*ubtp; tstring = [tit ':  L1-7 -- U'' dHCbar dx']; fnum = 0; ptit = 'BP_upHCb';
tem1 = -1*uptb; tstring = [tit ':  L1-7 -- V'' dHCbar dy']; fnum = 2; ptit = 'BP_vpHCb';
cint = 1;  clev = [-10:cint:-cint  cint:cint:10];

figure(1+fnum); figure_orient(1);
for i = 1:6;
  tem = squeeze(tem1(i,:,:));
  subplot(3,2,i);
    gcont(tem, clev);
    dc2(tem);
    title(['Lag:  ' num2str(lags(i))]);
  if i >= 5;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(3,2,3);
  ylabel(tstring);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Lagged_regressions/Plot_Tendencies
eval(['print -dps2 ' ptit '_L1-7_neglag.ps']);

figure(2+fnum); figure_orient(2);
for i = 7:12;
  tem = squeeze(tem1(i,:,:));
  subplot(3,2,i-6);
    gcont(tem, clev);
    dc2(tem);
    title(['Lag:  ' num2str(lags(i))]);
  if i >= 5;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(3,2,3);
  ylabel(tstring);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Lagged_regressions/Plot_Tendencies
eval(['print -dps2 ' ptit '_L1-7_poslag.ps']);