Create visc.cpp

Viscous flux caluclation