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))