We derive an expression for the energy-momentum tensor in the discrete lattice formulation of pure glue QCD. The resulting expression satisfies the continuity equation for energy conservation up to numerical errors with a symmetric procedure for the time discretization. In the case of the momentum conservation equation, we obtain an expression that is of higher accuracy in lattice spacing ( O ( a 2 ) ) than the naive discretization where fields in the continuum expressions are replaced by discretized counterparts. The improvements are verified by performing numerical tests on the derived expressions using classical real-time lattice gauge theory simulations. We demonstrate substantial reductions in relative error of one to several orders of magnitude compared to a naive discretization for both energy and momentum conservation equations. We expect our formulation to have applications in the area of pre-equilibrium dynamics in ultrarelativistic heavy ion collisions, in particular for the extraction of transport coefficients such as shear viscosity.