Three dimensional toroidal magnetohydrodynamic particle code

The magnetohydrodynamic particle code has been developed to three dimensions in a cylindrical coordinate system in order to describe the plasma in a torus. To keep the noise level down, the finite differences are defined halfway between grid points and the magentic force is defined in a non-conservative manner. Two practical examples of using such a code for physics applications are reported: simulations of high amplitude Global Alfven Eigenmodes and stabilization of flute modes by a hot electron ring