This paper describes a numerical treatment (ring average) to relax the time step in finite-difference schemes when using spherical and cylindrical coordinates with axis singularities. The ring average is used to develop a high-resolution thermosphere–ionosphere coupled community model. The technique is a significant improvement in space weather modeling capability, and it can also be adapted to more general finite-difference solvers for hyperbolic equations in spherical and polar geometries.