Replacing MMS for 3D Navier-Stokes with Sutherland viscosity

The previous code had the wrong thermal conductivity k (it was constant). It's fixed.