An energy-correction method is proposed as an addition to mainstream integrators for equations of motion of systems of classical spins. This solves the problem of nonconservation of energy in long computations and makes mainstream integrators competitive with symplectic integrators for spin systems that for different-site interactions conserve the energy explicitly. The proposed method is promising for spin systems with single-site interactions for which symplectic integrators do not conserve energy and thus have no edge against mainstream integrators. From the energy balance in the spin system with a phenomenological damping and Langevin fields, a formula for the dynamical spin temperature in the presence of single-site anisotropy is obtained.