Discrete nonnegativity principles are established for finite element approximations of nonlinear parabolic PDE systems with mixed boundary conditions. Previous results of the authors are extended such that diagonal dominance ...