  We show for a variety of classes of conservative PDEs that discrete gradient
methods designed to have a conserved quantity (here called energy) also have a
time-discrete conservation law. The discrete conservation law has the same
conserved density as the continuous conservation law, while its flux is found
by replacing all derivatives of the conserved density appearing in the
continuous flux by discrete gradients.
