Only compute the viscosity if a strain rate was given. This avoids
accessing empty arrays and also matches the documented behavior.
Add a testcase.
Add an entry to changes.h.
Merge pull request #45 from bangerth/why-evaluate-viscosity-in-adiabatic-conditions

Why evaluate viscosity in adiabatic conditions