changes to numerics_viscous.cpp