Add functions for checking mhd state and flux for invalid values.
Check time step and cell area in mhd solve function.
Use mhd state and flux checking functions in hll_athena.