incorporated viscous flux jacobian calculation into main program
changed error message