transformToRadialWavenumber

transforms in the spectral domain from (j,kl) to (j,kRadial)

---

Declaration

 [varargout] = transformToRadialWavenumber(varargin) 

Parameters

• varargin variables with dimensions \((j,kl)\)

Returns

• varargout variables with dimensions \((kRadial)\) or \((kRadial,j)\)

Discussion

Sums all the variance/energy in radial bins kRadial.

The following example takes the total energy of the geostrophic part of flow, converts it to a one-dimensional spectrum in \(k\), and then plots it with pcolor. The next plot then sums over all wavenumber, and produces plots the total energy spectrum as a function of vertical mode \(j\) only.

  figure
  tiledlayout('flow')
  Ekj = wvt.transformToRadialWavenumber( wvt.A0_TE_factor .* abs(wvt.A0).^2 );
  nexttile, pcolor(wvt.kRadial,wvt.j,Ekj), shading flat
  nexttile, plot(wvt.j,sum(Ekj,2))