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

An efficient numerical code FP3D has been developed to calculate particle orbits and evaluate particle confinement in 3D magnetic fields including stellarators and tokamaks with 3D fields. The magnetic field is either calculated from coils directly or obtained from equilibrium codes. FP3D has been verified with the 3D equilibrium code Variational Moments Equilibrium Code [Hirshman and Whitson, Phys. Fluids 26, 3553–3568 (1983)] for magnetic field calculation and with the drift-kinetic code Stellarator Fokker–Planck Iterative Neoclassical Conservative Solver [Landreman et al., Phys. Plasmas 21(4), 042503 (2014)] for neoclassical transport. The code has been applied successfully to the National Compact Stellarator Experiment stellarator [Nelson et al., Fusion Eng. Des. 66, 169–174 (2003)] for the calculation of neoclassical transport coefficient with the 3D magnetic field obtained directly from coils. FP3D is also used to calculate ripple losses in the tokamak Experimental Advanced Superconducting Tokamak [Y. Wan, Plasma Sci. Technol. 8(3) (2006)].