Add a new diagnostic algorithm for viscous dissipation which works in 2D and 3D, but assumes a constant viscosity tensor. Modify energy stats tests to use this new algorithm.