We present a variant restricted Additive Schwarz preconditioner and apply Partial-Newton-Krylov-Schwarz algorithm to solve nonlinear algebraic equations of two-dimensional three-temperature systems. Iteration and CPU time for convergence are decreased. Numerical results show efficiency of the method.