Abstract:Context. The study of linear waves and instabilities is necessary to understand the physical evolution of an atmosphere, and can provide physical interpretation of the complex flows found in simulations performed using Global Circulation Models (GCM). In particular, the acceleration of superrotating flow at the equator of hot Jupiters has mostly been studied under several simplifying assumptions, the relaxing of which may impact final results. Aims. We develop and benchmark a publicly available algorithm to identify the eigenmodes of an atmosphere around any initial steady state. We also solve for linear steady states. Methods. We linearise the hydrodynamical equations of a planetary atmosphere in a steady state with arbitrary velocities and thermal profile. We then discretise the linearised equations on an appropriate staggered grid, and solve for eigenvectors and linear steady solutions with the use of a parallel library for linear algebra: ScaLAPACK. We also implement a posteriori calculation of an energy equation in order to obtain more information on the underlying physics of the mode. Results. Our code is benchmarked against classical wave and instability test cases in multiple geometries. The steady linear circulation calculations also reproduce expected results for the atmosphere of hot Jupiters. We finally show the robustness of our energy equation, and its power to obtain physical insight into the modes. Conclusions. We have developed and benchmarked a code for the study of linear processes in planetary atmospheres, with an arbitrary steady state. The calculation of an a posteriori energy equation provides both increased robustness and physical meaning to the obtained eigenmodes. This code can be applied to various problems, and notably to further study the initial spin up of superrotation of GCM simulations of hot Jupiter.

VENUS-LEVIS and its spline-Fourier interpolation of 3D toroidal magnetic field representation for guiding-centre and full-orbit simulations of charged energetic particles

Numerical integration of particle orbits in discontinuous fields using VENUS-LEVIS and SPEC

Higher order Larmor radius corrections to guiding-centre equations and application to fast ion equilibrium distributions

Symplectic Simulations of Radial Diffusion of Fast Alpha Particles in the Presence of Low-Frequency Modes in Rippled Tokamaks

Impact of RMP magnetic field simulation models on fast ion losses

Monte Carlo Orbit-Following Simulations Including the Finite Larmor Radius Effect Based on a Phase-Space Coordinate Transform Method

Variational Symplectic Algorithm for Guiding Center Dynamics and Its Application in Tokamak Geometry

Implementation of Field-Aligned Coordinates in a Semi-Lagrangian Gyrokinetic Code for Tokamak Turbulence Simulation

Energetic particle tracing in optimized quasi-symmetric stellarator equilibria

FP3D: A code for calculating 3D magnetic field and particle motion

Variational Symplectic Integrator for the Guiding Center Motion of Charged Particles for Long Time Simulations in General Magnetic Fields

Energetic Particle Tracing in Optimized Quasisymmetric Stellarator Equilibria

Nonlinear magnetohydrodynamics simulation using high-order finite elements

Eigenvectors, Circulation and Linear Instabilities for Planetary Science in 3 Dimensions (ECLIPS3D)

A well-balanced and exactly divergence-free staggered semi-implicit hybrid finite volume/finite element scheme for the incompressible MHD equations

Variational Symplectic Integrator for Long-Time Simulations of the Guiding-Center Motion of Charged Particles in General Magnetic Fields.

Automation of the guiding center expansion

Guiding center and gyrokinetic orbit theory for large electric field gradients and strong shear flows

Construction and analysis of guiding center distributions for tokamak plasmas with ambient radial electric field

High-Precision Nonsingular Integrator of Guiding-Center Orbit Close to Magnetic Axis in Tokamak Equilibrium Configuration