A three dimensional bounce-averaged Fokker-Planck (FP) numerical code has been newly developed based on fully implicit iterative solving method, and relativistic effect is also included in the code. The code has been tested against various benchmark cases: Ohmic conductivity in the presence of weak Ohmic electric field, runaway losses of electrons in the presence of strong Ohmic electric field, lower hybrid current drive and electron cyclotron current drive via two- or three-dimensional simulation. All the test cases run fast and correctly during calculations. As a result, the code provides a ...
