Some code for bulk viscosity