The viscous terms are evaluated in two steps. First, we obtain the spatial derivatives of the primitive variables using the discontinuous Galerkin approach. Then, we repeat the process for each of the viscous fluxes using these derivatives. If the user supplies Dirichlet boundary conditions for the momentum and energy variables we set these terms explicitly after the fluxes have been evaluated and then project the result using the orthogonal basis.