fixed bug in compute_max_principle_viscosity()