Conventional explicit time-domain methods used to solve Maxwell's equations are reliable and robust but are conditionally stable and often require the fulfillment of the condition CFL ≤ 1. While this is acceptable for many applications, in some instances where the Maxwell's equations are solved alongside systems with slower propagation velocities, explicit methods prove costly. This is the case for non-relativistic electromagnetic Particle-In-Cell methods which are required to study plasma thrusters. Several algorithms have been proposed to retain a nearly explicit formulation using large time steps to achieve higher CFL values. Among these is the semi-Lagrangian Constrained Interpolation Profile method. While the ability of this method to handle CFL > 1 has been demonstrated for planar 2D-3D cases, this has not been done for 2D cases with cylindrical symmetry. In this paper, a procedure is presented to compute the electromagnetic wave propagation in 2D domains with cylindrical symmetry using the Constrained Interpolation Profile (CIP) method. The CIP scheme is extended for CFL ≥ 1 cases, and a ghost node method is proposed to deal with the axis singularity and with the wall boundary condition. The results are compared to the fields of a Hertzian dipole and with a coaxial cable, and they show a good agreement.