bug fixes
started change to gradient calculation on faces - need to update viscous flux routines
added gradient calculation on face
removed old gradient calculation function
bug fixes in viscous gradient calculation on face
added function to zero a tensor
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