Define the viscosity beta as a constant instead of a hardcoded value

Set #define const_viscosity_beta in const.h instead of the hardcoded value of 3 in hydro_iact.h.

Note: default cosmology alpha = 0.8, beta = 3; and planetary alpha = 1.5, beta = 4.