Short-N characteristic radiative transfer method in 2D cylindrical RZ coordinates