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”, J. Phys. Oceanogr.41, 412-428. [PDF]
NEW VERSION: Fourier-Galerkin method is faster and more accurate.wave_analysis_FG.m : code to solve demonstration problem
nash_data.txt : observational data for demonstration problem
vTG_FG.m : the main subroutine
BaryL.m : auxiliary subroutine
vTG_FGprep.m : auxiliary subroutine (call first)
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.