UAm

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

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

UAm(:,:,1) = 0;

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

UAm(1,1,:) = 1;