Additive operator splitting methods for solving systems of nonlinear finite difference equations