The Elastodynamic Finite Integration Technique (EFIT), originally developed by Fellinger et al.,1–3 represents a stable and efficient numerical code to model elastic wave propagation in linearly-elastic isotropic and anisotropic, homogeneous and heterogeneous as well as dissipative and nondissipative media. In previous works, the FIT discretization of the basic equations of linear elasticity, Hooke's law and Cauchy's equation of motion, was exclusively carried out in Cartesian coordinates. For problems in cylindrical geometries it is more suitable to use cylindrical coordinates. By that, axisymmetric problems can be treated in a two-dimensional staggered grid in the r,z-plane. The paper presents an EFIT version for axisymmetric problems in cylindrical coordinates called Cylindrical EFIT (CEFIT). After demonstrating the accuracy of the numerical code by a comparison between simulation results and analytical solutions, different examples of application are given. These examples include modeling of sound fields of ultrasonic transducers, thermoelastic laser sources, geophysical borehole probes, impact-echo measurements in layered media, and load simulations of the European Spallation Source (ESS) mercury target.