A full wave code for ion cyclotron waves in toroidal plasmas

The code TORIC solves the finite Larmor radius wave equations in the ion cyclotron frequency range in arbitrary axisymmetric toroidal geometry. The model used describes the compressional and torsional Alfven waves (or, depending on the parallel phase velocity, the kinetic counterpart of the latter), and ion Bernstein waves excited by mode conversion near the first ion cyclotron harmonic. In the ion response the broadening of the absorption regions due to the finite width of the cyclotron resonance of individual ions in toroidal geometry is taken into account. The parallel component of the wave electric field is evaluated on the same footing as the transverse ones; the response of the electrons includes Landau damping, Transit Time damping and the mixed term. The numerical approach uses a spectral representation of the solution in the poloidal angle θ, and cubic finite elements in the radial variable ψ. Great flexibility is provided in the way ion Bernstein waves excited by mode conversion are damped when their wavelength becomes comparable with the ion Larmor radius, in the regularization of Alfven resonances, and in the treatment of the outer plasma layers. As an option, we have also implemented the Order Reduction Algorithm, which provides a particularly fast, yet accurate evaluation of the power deposition profiles in toroidal geometry. Thee present report describes the model and its numerical implementation, and provides the information needed to use the code. A few examples illustrating applications of TORIC are also included. (orig.)