  We present a strategy for a non-perturbative determination of the finite
renormalization constants of the energy-momentum tensor in the SU(3) Yang-Mills
theory. The computation is performed by imposing on the lattice suitable Ward
Identites at finite temperature in presence of shifted boundary conditions. We
show accurate preliminary numerical data for values of the bare coupling g_0^2
ranging for 0 to 1.
