UAp

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

Developer documentation: this item describes internal implementation details.

---

Discussion

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

These reconstruction coefficients map \(A_+\) onto the \(u\) state variable. In the historical notation of Early et al. (2021), they are the row 1, 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{UAp} \equiv \frac{k \omega - i l f_0}{\omega K}\]

in the manuscript. In code this is computed with,

alpha = atan2(L,K);
fOmega = f./omega;
UAp = (cos(alpha)-sqrt(-1)*fOmega.*sin(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,

UAp(1,1,:) = 1;