modesAtQuadraturePoints
Return resolved modes on mode-adapted quadrature points.
Declaration
[F,G,h,zq] = modesAtQuadraturePoints(self,options)
Parameters
selfInternalModesSpectral instanceoptions.nPointsnumber of quadrature points and modes requestedoptions.omegafixed frequency in radians per second; specify eitheromegaorkoptions.kfixed horizontal wavenumber; specify eitheromegaork
Returns
Fhorizontal-velocity mode matrix onzqGvertical-velocity mode matrix onzqhequivalent-depth row vectorzqquadrature-point depth locations
Discussion
This helper chooses quadrature points from the roots of a resolved higher mode, temporarily evaluates the requested modes on those points, and restores the original output grid before returning.