Summary: The boundary value problem for the singularly perturbed reaction-diffusion parabolic equation in a ball in the case of spherical symmetry is considered. The derivatives with respect to the radial variable appearing in the equation are written in divergent form. The third kind boundary condition, which admits the Dirichlet and Neumann conditions, is specified on the boundary of the domain. The Laplace operator in the differential equation involves a perturbation parameter, \(\varepsilon^2\) where, \(\varepsilon\) takes arbitrary values in the half-open interval \((0, 1]\). When \(\varepsilon \to 0\), the solution of such a problem has a parabolic boundary layer in a neighborhood of the boundary. Using the integro-interpolational method and the condensing grid technique, conservative finite difference schemes on flux grids are constructed that converge \(\epsilon\)-uniformly at a rate of \(O(N^{-2}\ln^2N + N_0^{-1})\), where \(N + 1\) and \(N_0 + 1\) are the numbers of the mesh points in the radial and time variables, respectively.