cleaned up commented out code
bug fix in velocity gradient calculation
added print statements
commments for debugging
BUG FIX!!! viscous fluxes now working with ghost cells