Header file to implement viscosity in the solver