InternalModesWKBHydrostatic

Compute hydrostatic WKB mode approximations from a spectrally resolved stratification.


Declaration

classdef InternalModesWKBHydrostatic < InternalModesSpectral

Overview

InternalModesWKBHydrostatic uses the spectral initialization path of InternalModesSpectral to resolve \(N^2(z)\), then applies a hydrostatic WKB approximation for the fixed-\(\omega\) problem. This class is useful as an asymptotic comparison tool rather than as the primary production solver.

In this approximation, the vertical phase coordinate is built from the positive part of \(N(z) - \omega\), and the modal depth is estimated from

\[h_j = \frac{1}{g}\left(\frac{d}{j\pi}\right)^2,\]

where \(d\) is the accumulated hydrostatic WKB phase over the oscillatory region.

im = InternalModesWKBHydrostatic(rho=rho, zIn=zIn, zOut=zOut, latitude=latitude);
[F, G, h, k] = im.modesAtFrequency(5*im.f0);

Topics