Matlab tools to solve the viscous Taylor Goldstein equation for both instabilities and waves

wave_analysis.m : code to solve
demonstration problem

nash_data.txt : observational data for
demonstration problem

SSF.m : the main subroutine

ddz.m , ddz2.m , ddz4.m :
auxiliary subroutines

Reference:

Smyth, W.D., J.N. Moum and J.D. Nash, 2011: “Narrowband, high-frequency oscillations at the equator. Part II: Properties of shear instabilities",-----------------------------------------------------------------------------------------

NEW VERSION: Fourier-Galerkin method is faster and more accurate.

wave_analysis_FG.m : code to solve demonstration problemnash_data.txt : observational data for demonstration problem

vTG_FG.m : the main subroutine

BaryL.m : auxiliary subroutine

vTG_FGprep.m : auxiliary subroutine (call first)

Reference:

Lian, Q. , Z. Liu, W.D. Smyth 2019: ``Numerical computation of instabilities and internal waves from in situ measurements via the viscous Taylor-Goldstein problem'', J Atm. Oc. Tech. (in press) [PDF]

-----------------------------------------------------------------------------------------

Watch this space for more efficient methods currently in the works.