The TERPSICHORE code is a global 3D linear ideal MHD stability solver that employs hybrid finite elements for the radial discretisation and Fourier decomposition in the angular variables. A transformation from the angular coordinates in ANIMEC to Boozer magnetic coordinates is performed because this frame is more appropriate to follow guiding centre particle orbits and evaluate MHD stability. The normal and binormal components of the displacement vector are the outputs from which the electrostatic potential and the vector potential can be calculated.
The main reference is:
- TERPSICHORE: A Three-Dimensional Ideal Magnetohydrodynamic Stability Program
Anderson, D.V., Cooper, W.A., Gruber, R., Merazzi, S., Schwenn, U. (1990). TERPSICHORE: A Three-Dimensional Ideal Magnetohydrodynamic Stability Program. In: Devreese, J.T., Van Camp, P.E. (eds) Scientific Computing on Supercomputers II. Springer, Boston, MA. https://doi.org/10.1007/978-1-4613-0659-7_8
The license and README are available at https://gitlab.epfl.ch/spc/public/terpsichore
