Correct units for the eddy viscosity
Added diffusion stats for dns mode
Create 2nd and 4th order gradient computations
All standard stats use the correct order interpolation
Merge remote-tracking branch 'chiel/stats' into thijs