Discovered that I was never calculating the gradient of the primitive variables...  Uncommented the appropriate code.
Changes to the preprocessing (now update nonphysical cells) and error checking before applying any residuals to the linear system