Attempted code to use floating point exception handling, but it didn't work on Mac OS, so it's commented out for future potential use on Linux.  Also commented out code for checking indiv values of pGrid physical quanitities in the main
More code to check the density in every cell for NaNs