VAm

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 2 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{VAm} \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;
VAm = (sin(alpha)-sqrt(-1)*fOmega.*cos(alpha));

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

VAm(:,:,1) = 0;

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

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