Fixed an error in computing viscosity which was causing poor convergence.