tune down accuracy for solving velocity in VISH2+1 to stable numerical simulation with bulk viscosity
change to read in bulk viscous delta f coefficient from a table (for different PCE EOS)