The drift-kinetic code VENUS-LEVIS was designed to simulate a wide variety of physical phenomena related to fast particles in electromagnetic fields. The code uses 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 as