I'm currently trying to reproduce some g-mode analysis that involves a formula with eigenvalues of Laplace's tidal equation. This is beyond my usual area of expertise but, given a vague notion that this is something GYRE needs for the Traditional Approximation of Rotation (TAR), I went digging around inside GYRE to see how it evaluates the eigenvalues.
As far as I can tell, this is done through approximations in terms of Chebyshev polynomials, the coefficients of which are stored in data/tar. Fortunately, the segment of code, located in src/tar, include nifty Python routines for plotting evaluations of these functions. So I tried to reproduce Fig. 1 of Lee & Saio (1997).
At risk of offending the copyright gods, here's their figure:
Whereas here's my output from going to src/tar and running
Code: Select all
python3 plot_tar_fit.py -10.0 10.0 201 ../../data/tar/tar_fit.m-2.*.h5 LS97_Fig1.png