Documentation of uprime_gradtbar

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

Help text

  Start with u'dtbardx and v'dtbardy

Cross-Reference Information

Listing of script uprime_gradtbar

clear

tim = 101:550;
ctlim = [110 300 -60 60];
lev = 5;

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

%  First get tbar

filin = 'temp_A_L1-10.nc';
nc = netcdf(filin, 'nowrite');
  depth = nc{'depth'}(:);
  latt = nc{'latitude'}(:);
  lont = nc{'longitude'}(:);
  [xk, yk] = keep_var(ctlim, lont, latt);
  temp = nc{'temp'}(tim, lev, yk, xk);
  mv = nc{'temp'}.missing_value(:);
nc = close(nc);
temp(temp == mv) = NaN;
latt = latt(yk); lont = lont(xk);
temp = squeeze(mean2(temp));

%  Now get u' and v'

filin = 'ul_L1-10.nc';
nc = netcdf(filin, 'nowrite');
  latu = nc{'latitude'}(:);
  lonu = nc{'longitude'}(:);
  [xk, yk] = keep_var(ctlim, lonu, latu);
  u = nc{'u'}(tim,lev,yk,xk);
  mv = nc{'u'}.missing_value(:);
nc = close(nc);
u(u == mv) = NaN;
latu = latu(yk); lonu = lonu(xk);
[u, climu] = remove_mean(u/100);

filin = 'vl_L1-10.nc';
nc = netcdf(filin, 'nowrite');
  v = nc{'v'}(tim,lev,yk,xk);
  mv = nc{'v'}.missing_value(:);
nc = close(nc);
v(v == mv) = NaN;
[v, climv] = remove_mean(v/100);

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

%  Load CPCs

cd /home/disk/hayes2/dvimont/csiro/matlab_data/20Deg_Isotherm
load 20deg_CEOF_LP10_yr101-550.mat; tit = 'Lowpass Filtered Data ( > 10 Years )';
%load 20deg_CEOF_HP8_yr101-550.mat; tit = 'Highpass Filtered Data ( < 8 Years )';
%load 20deg_CEOF_RAW_yr101-550.mat; tit = 'Unfiltered Data';

%  Get regressions

lag = 0; lg = lag*pi/180; lg2 = 1;
num = 1; lind = 1;
nfrm = 6;

%  Store all regressions under one variable

[nlat, nlon] = size(squeeze(climu));
j = sqrt(-1);
timeseries = sqrt(2)*pcs(:,num)./std(pcs(:,num));
u = reshape(u, ntim, nlat*nlon);
v = reshape(v, ntim, nlat*nlon);
clear temtim ureg vreg
for i = 1:nfrm
  wgt = conj(exp(j * ((i-1) * pi/(lg2*nfrm) + lg) ));
  temtim(:,i) = squeeze(real(wgt .* timeseries));
  ureg(i, :) = temtim(:,i)' * u ./ ntim;
  vreg(i, :) = temtim(:,i)' * v ./ ntim;
end
ureg = reshape(ureg, 6, nlat, nlon);
vreg = reshape(vreg, 6, nlat, nlon);

%  Get dtdx and dtdy

[nlat, nlon] = size(temp);
global DEGREE RADIAN RADUS
clear ty tx dtdy dtdx j
for j = 1:nlat;
  ty(j,:) = interp1(lont, temp(j,:), lonu)';
end
for j = 1:nlon;
  tx(:,j) = interp1(latt, temp(:,j), latu);
end  
[dtdy temlat] = sph_grady1(ty, RADIAN*latt, RADIAN*lonu, 1);
[dtdx temlon] = sph_gradx1(tx, RADIAN*latu, RADIAN*lont, 0);


% Plot the data

get_global; FRAME = ctlim; XAX = lonu; YAX = latu;

figure(1); figure_orient;
cint = 1; clev = [-50:cint:-cint cint:cint:50];
for i = 1:nfrm
  tem = -1e9 * squeeze(ureg(i,:,:)) .* dtdx;
  subplot(3,2,i);
    gcont(tem, clev);
    dc
  if i > 4;
    xlabel(['Contour Interval:  ' num2str(cint) ' x 10^-^9 K s^-^1']);
  end
end
subplot(3,2,3);
  ylabel([tit ':  Zonal Advection of Tbar by  U'';  '...
          'Depth = ' num2str(depth(lev)/100) 'm']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Velocity/Plots_Tadv
%print -dps2 LP10_uprime_dtbardx_L5_yr1.ps

figure(2); figure_orient;
for i = 1:nfrm
  tem = -1e9 * squeeze(vreg(i,:,:)) .* dtdy;
  subplot(3,2,i);
    gcont(tem, clev);
    dc
  if i > 4;
    xlabel(['Contour Interval:  ' num2str(cint) ' x 10^-^9 K s^-^1']);
  end
end
subplot(3,2,3);
  ylabel([tit ':  Meridional Advection of T'' by Mean V;  '...
          'Depth = ' num2str(depth(lev)/100) 'm']);

%print -dps2 LP10_vprime_dtbardy_L5_yr1.ps