One considers an elliptic system of variational inequalities with diagonal linear operators. The unknown vector functions on the boundary of the considered domain belong to a convex set of an arbitrary form. One investigates the regularity of the solution of such a system.