The drift-kinetic code VENUS-LEVIS is designed for simulating a wide variety of physical phenomena related to fast particle motion in electromagnetic fields. The code employs a 4th order Runge-Kutta method to solve the single particle equations of motion, either in the guiding-centre approximation or following the full particle orbits. The formulation is independent of coordinate choice and handles 3D time-varying electromagnetic fields. The interaction with the background plasma as well