Hello.
I am thinking that for the adiabatic and non-adiabatic run of GYRE, the eigenfunctions are internally available.
Then it is worthwhile to get the quasi-adiabatic work integral, and also to compare that with the full non-adiabatic mode stability.
If this would be done internally, then there would be no need to store all eigenfunctions, and do them on a post-processing phase.
Therefore, I want to kindly ask for this additional feature in the next releases of the code.
Cheers.
Ehsan.