Documentation of w_intbar_dtprimedz

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

Help text

  First, get ubar and vbar

Cross-Reference Information

Listing of script w_intbar_dtprimedz

clear
tim = 101:550;
ctlim = [108 302 -62 62];
lev = 1:7;
cd /home/disk/hayes2/dvimont/csiro/data

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

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

%  Now get tprime

cd /home/disk/hayes2/dvimont/csiro/data
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(temp);

depth = depth/100;
[temp, tbar] = remove_mean(temp);
[ntim, nlev, nlat, nlon] = size(temp);

%  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

j = sqrt(-1);
timeseries = sqrt(2)*pcs(:,num)./std(pcs(:,num));
temp = reshape(temp, ntim, nlev*nlat*nlon);
clear temtim treg
for i = 1:nfrm
  wgt = conj(exp(j * ((i-1) * pi/(lg2*nfrm) + lg) ));
  temtim(:,i) = squeeze(real(wgt .* timeseries));
  treg(i, :) = temtim(:,i)' * temp ./ ntim;
end
treg = reshape(treg, 6, nlev, nlat, nlon);
temp = reshape(temp, ntim, nlev, nlat, nlon);

%  Get div(ubar);

clear divu lat2 lon2
for j = 1:nlev;
  [divu(j,:,:), lat2, lon2] = sph_div1(squeeze(ubar(j,:,:)), ...
                                       squeeze(vbar(j,:,:)), latu, lonu, 0, 1);
end

%  Get vvel (vvel @ lev x = sum(vvel(1:x,:,:))

[nlev, nlat, nlon] = size(divu);
d1 = depth(1);
depthw = [];
for i = 1:nlev;
  depthw = [depthw; depth(i)+d1];
  d1 = depth(i+1) - depth(i) - d1;
end
dz = diff([0; depthw]);
dz = dz(1:7) * ones(1, nlat*nlon);
dz = reshape(dz, nlev, nlat, nlon);
vvel = divu(:,:,:) .* dz;

%  Get dtdz

[tem, nlev, nlat, nlon] = size(treg);
temp2  = reshape(shiftdim(treg, 1), nlev, nlat*nlon*tem);
dz = diff(depth(1:7)) * ones(1, nlat*nlon*tem);
dzdepth = depth(1:6) + 0.5*diff(depth(1:7));
dtdz = diff(temp2) ./ dz;
for i = 1:nlat*nlon*tem;
  dtdz(:,i) = interp1(dzdepth, dtdz(:,i), depthw(1:6));
end
dtdz = shiftdim(reshape(dtdz, (nlev-1), nlat, nlon, tem), 3);

%  Plot the data

default_global; XAX = lon2; YAX = lat2;
lev = 1:5;

figure(1); figure_orient;
cint = 2; clev = [-20:cint:-cint cint:cint:20];
for i = 1:nfrm
  if length(lev) > 1;
    tem = 1e9*squeeze(sum(vvel(lev,:,:))) .* squeeze(dtdz(i,max(lev),:,:));
  else
    tem = 1e9*squeeze(vvel(lev,:,:)) .* squeeze(dtdz(i,lev,:,:));
  end
  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 ':  Vertical Advelction of T'' by  Wbar;  '...
          'Depth = ' num2str(depthw(max(lev))) 'm']);

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