Documentation of tendency_terms

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

Help text

  Load PCS:

Cross-Reference Information

Listing of script tendency_terms

clear

cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
load LP10_L1-7_CEOF.mat; tit = 'Lowpass Filtered Data ( > 10 Years )';
%load HP8_L1-7_CEOF.mat; tit = 'Highpass Filtered Data ( < 8 Years )';
%load RAW_L1-7_CEOF.mat; tit = 'Unfiltered Data';

cd ~/matlab/CSIRO/Heat

nfrm = 24;
tim = 101:550;
lev = 1:7;

biff = 2;

%  First, get WSTP values

global RADUS DEGREE RADIAN

if biff == 1;
  lims = [117.5 169 7.5 24];
elseif biff == 2;
  lims = [179 271 -5 5];
end

%  Get ubar_dHCprimedx and vbar_dHCprimedy

   [utem, vtem, lat_out, lon_out, depth_out] = ...
           get_ubar_gradtprime(pcs, lims, nfrm, tim, lev);
   [tem, nlat, nlon] = size(utem);
   dx = mean(diff(lon_out)); dy = mean(diff(lat_out));
   lon_out = [lon_out; max(lon_out)+dx] - dx/2;
   lat_out = [lat_out; max(lat_out)+dy] - dy/2;
   A = RADUS * ( diff(sin(RADIAN * lat_out)) ) * ( RADIAN * diff(lon_out) )'; 
   A = reshape(A, 1, nlat*nlon);
   ubtp_wstp = -1 * A * reshape(utem, nfrm, nlat*nlon)' ./ sum(A);
   vbtp_wstp = -1 * A * reshape(vtem, nfrm, nlat*nlon)' ./ sum(A);

%  Get uprime_dHCbardx and vprime_dHCbardy

   [utem, vtem, lat_out, lon_out, depth_out] = ...
           get_uprime_gradtbar(pcs, lims, nfrm, tim, lev);
   [tem, nlat, nlon] = size(utem);
   dx = mean(diff(lon_out)); dy = mean(diff(lat_out));
   lon_out = [lon_out; max(lon_out)+dx] - dx/2;
   lat_out = [lat_out; max(lat_out)+dy] - dy/2;
   A = RADUS * ( diff(sin(RADIAN * lat_out)) ) * ( RADIAN * diff(lon_out) )'; 
   A = reshape(A, 1, nlat*nlon);
   uptb_wstp = -1 * A * reshape(utem, nfrm, nlat*nlon)' ./ sum(A);
   vptb_wstp = -1 * A * reshape(vtem, nfrm, nlat*nlon)' ./ sum(A);

%  Get wbar_dHCprimedz

   [wtem, lat_out, lon_out] = ...
           get_wbar_dtprimedz(pcs, lims, nfrm, tim, lev);
   [tem, nlat, nlon] = size(wtem);
   dx = mean(diff(lon_out)); dy = mean(diff(lat_out));
   lon_out = [lon_out; max(lon_out)+dx] - dx/2;
   lat_out = [lat_out; max(lat_out)+dy] - dy/2;
   A = RADUS * ( diff(sin(RADIAN * lat_out)) ) * ( RADIAN * diff(lon_out) )'; 
   A = reshape(A, 1, nlat*nlon);
   wbtp_wstp = -1 * A * reshape(wtem, nfrm, nlat*nlon)' ./ sum(A);

%  Get wprime_dHCbardz

   [wtem, lat_out, lon_out] = ...
           get_wprime_dtbardz(pcs, lims, nfrm, tim, lev);
   [tem, nlat, nlon] = size(wtem);
   dx = mean(diff(lon_out)); dy = mean(diff(lat_out));
   lon_out = [lon_out; max(lon_out)+dx] - dx/2;
   lat_out = [lat_out; max(lat_out)+dy] - dy/2;
   A = RADUS * ( diff(sin(RADIAN * lat_out)) ) * ( RADIAN * diff(lon_out) )'; 
   A = reshape(A, 1, nlat*nlon);
   wptb_wstp = -1 * A * reshape(wtem, nfrm, nlat*nlon)' ./ sum(A);


%  Plot the data

ph = 180/24 * [0:48];
figure(1); figure_orient;
sp(biff);
  h1 = plot(ph, [ubtp_wstp -1*ubtp_wstp ubtp_wstp(1)], 'x-k', ...
       ph, [vbtp_wstp -1*vbtp_wstp vbtp_wstp(1)], 'o-k', ...
       ph, [uptb_wstp -1*uptb_wstp uptb_wstp(1)], '*-k', ...
       ph, [vptb_wstp -1*vptb_wstp vptb_wstp(1)], 's--k', ...
       ph, [wbtp_wstp -1*wbtp_wstp wbtp_wstp(1)], 'v-k', ...;
       ph, [wptb_wstp -1*wptb_wstp wptb_wstp(1)], '^-k');
  grid on
  axis([0 360 -8 8]);
  set(gca, 'XTick', [0:30:360]);
  legend(h1, 'Ubar dHC''dx', 'Vbar dHC''dy', 'U'' dHCbardx', 'V'' dHCbardy', ...
             'Wbar dHC''dz', 'W'' dHCbardz');

end;