Documentation of get_uptb

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

Function Synopsis

[udtdx, vdtdy, lat_out, lon_out, depth_out] = ...

Help text

  [uptb, vptb, lat_out, lon_out, depth_out] = ...
          get_uprime_gradtbar(pcs, lims, nfrm, tim, lev);

Cross-Reference Information

Listing of function get_uptb

function [udtdx, vdtdy, lat_out, lon_out, depth_out] = ...
          get_uprime_gradtbar(pcs, lims, nfrm, tim, lev);

%lims = [117.5 169 7.5 24];
if nargin == 1;
  lims = [110 300 -60 60]; nfrm = 6; tim = 101:550; lev = 1;
elseif nargin == 2;
  nfrm = 6; tim = 101:550; lev = 1;
elseif nargin == 3;
  tim = 101:550; lev = 1;
elseif nargin == 4;
  lev = 1;
end

if length(lev) > 1; 
  error('Single levels only')
end;

cdtem = ['cd ' eval('pwd')];

cd /home/disk/hayes2/dvimont/csiro/data

%  First get tbar

temp = getnc('temp', lims, lev, tim);
[latt, lont, deptht] = getll('temp', lims);
temp = squeeze(mean(temp));

%  Now get u' and v'

[u, ue] = getnc('u', 'ue', lims, lev, tim);
[latu, lonu, depthu] = getll('u', lims);
[ntim, nlat, nlon] = size(ue);
ue = (ue(:,1:(nlat-1),:) + ue(:,2:nlat,:))/2;
u = u + ue;
clear ue;

[u, climu] = remove_mean(u/100);

[v, ve] = getnc('v', 've', lims, lev, tim);
[ntim, nlat, nlon] = size(ve);
ve = (ve(:,:,1:(nlon-1)) + ve(:,:,2:nlon))/2;
v = v + ve;
clear ve;

[v, climv] = remove_mean(v/100);

[ntim, nlat, nlon] = size(u);

%  PCS should be sent in the command line

%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';

%  Get regressions

if isreal(pcs);
  ureg = regress_eof(u, pcs, nfrm);
  vreg = regress_eof(v, pcs, nfrm);
else
  ureg = regress_ceof(u, pcs, nfrm);
  vreg = regress_ceof(v, pcs, nfrm);
end

%  Get dtdx and dtdy

[nlat, nlon] = size(temp);
global DEGREE RADIAN RADUS

ty = (temp(:,1:(nlon-1)) + temp(:,2:nlon))/2;
tx = (temp(1:(nlat-1),:) + temp(2:nlat,:))/2;
[dtdy temlat] = sph_grady1(ty, RADIAN*latt, RADIAN*lonu, 0);
[dtdx temlon] = sph_gradx1(tx, RADIAN*latu, RADIAN*lont, 0);

%  Get ubar*dtdx and vbar*dtdy, and sum vertically

for i = 1:nfrm;
  udtdx(i,:,:) = squeeze(ureg(i,:,:)) .* dtdx;
  vdtdy(i,:,:) = squeeze(vreg(i,:,:)) .* dtdy;
end

lat_out = latu;
lon_out = lonu;
depth_out = depthu(lev);

eval(cdtem);