VAp

Reconstructs \(v\) from \(A_+\).

Developer documentation: this item describes internal implementation details.

---

Discussion

Reconstructs \(v\) from \(A_+\).

These reconstruction coefficients map \(A_+\) onto the \(v\) state variable. In the historical notation of Early et al. (2021), they are the row 2, column 1 entries of \(S\) for the primary internal-gravity-wave and geostrophic solutions in equation C4.

For \(k^2+l^2>0, j>0\) this is written as,

\[\textrm{VAp} \equiv \frac{l \omega + i k f_0}{\omega K}\]

in the manuscript. In code this is computed with,

alpha = atan2(L,K);
fOmega = f./omega;
VAp = (sin(alpha)+sqrt(-1)*fOmega.*cos(alpha));

There are no \(k^2+l^2>0, j=0\) wave solutions for a rigid lid,

UAp(:,:,1) = 0;

The inertial solutions occupy the \(k^2+l^2=0\) portion of the matrix,

VAp(1,1,:) = sqrt(-1);