“…It is possible to further reduce the model by averaging out, in addition to the cyclotron motion, the trapped particle bounce (or banana) motion and considering adiabatic passing particles with kinetic trapped ions (from zero to supra-thermal, although non-relativistic, energies), thus focusing on slow phenomena on a time-scale of the toroidal precession of trapped particles with thermal velocity, which is around 10 −2 s. In this work we use the TERESA code [6,7,10,14,19,31,32,40] which is based on this reduced gyrokinetic model and which is less computationally intensive than standard gyrokinetic codes. This axisymmetric, electrostatic code solves the Vlasov equation coupled to the quasi-neutrality constraint, in a 4 dimensional space: 2 spatial variables (α, ψ) and 2 adiabatic invariants κ and E. α = ϕ − qθ is the precession angle with ϕ the toroidal angle, θ the poloidal angle and q the safety factor, taken independent of the radius.…”