initial solvers fail badly in atomic region with SNR abundances
attached sim fails with non-conservation in C08 and trunk r3667. Robin reports does fine in newmole. This might be an interesting sim to add to the test suite after newmole comes to trunk
CoolHeatError: Te: 3.0133e+03 C: 2.3546e-16 H: 1.8927e-17 (C-H)/H: 1.1440e+01 CoolHeatError: Te: 2.9530e+03 C: 2.2965e-16 H: 1.8987e-17 (C-H)/H: 1.1095e+01 CoolHeatError: Te: 2.8939e+03 C: 2.2428e-16 H: 1.9051e-17 (C-H)/H: 1.0773e+01
PROBLEM DISASTER PressureTotal: the chemical species are not conserved. The sum of the densities of the atoms and ions for Carbon is 8.146e+02 which is greater than the gas-phase abundance of the element, 1.245e+02. Something that cannot happen, has happened. This is TotalInsanity, I live in ../service.cpp.