We develop a unified and easy to use framework to study robust fully discrete numerical methods for nonlinear degenerate diffusion equations ∂tu − L[ϕ(u)] = f (x, t) in R N × (0, T), where L is a general symmetric Lévy type diffusion operator. Included are both local and nonlocal problems with e.g. L = ∆ or L = −(−∆) α 2 , α ∈ (0, 2), and porous medium, fast diffusion, and Stefan type nonlinearities ϕ. By robust methods we mean that they converge even for nonsmooth solutions and under very weak assumptions on the data. We show that they are L p-stable for p ∈ [1, ∞], compact, and convergent in C([0, T ]; L p loc (R N)) for p ∈ [1, ∞). The first part of this project is given in [36] and contains the unified and easy to use theoretical framework. This paper is devoted to schemes and testing. We study many different problems and many different concrete discretizations, proving that the results of [36] apply and testing the schemes numerically. Our examples include fractional diffusions of different orders, and Stefan problems, porous medium, and fast diffusion nonlinearities. Most of the convergence results and many schemes are completely new for nonlocal versions of the equation, including results on high order methods, the powers of the discrete Laplacian method, and discretizations of fast diffusions. Some of the results and schemes are new even for linear and local problems.