molecular chemistry is now fully consistent with the ionization balance. This is done by scaling chemical reactions including atoms and positive ions to have rates dependent on a combined species which is the sum of the ions within the chemical network solver. The resultant effects of the chemistry on the ionization balance are allowed for by having the molecular network calculate the net source and sink rates for these ions from molecular chemistry, which are then included in the ionization balance solver.

The nonlinear system for the molecular chemistry is solved using an adaptive timestep implicit solver: in typical usage, this is set up to run the chemical balance to full late time equilibrium. The molecular chemistry and ionization solvers have been adapted to allow the solution of time-dependent and steady-state advective flows, using the approach described by @HenneyEtAl05. in the existing equilibrium solvers to adapt them to find the solution of a backward Euler implicit system for the time advance of the state. This allows the code to take advantage of the existing equilibrium solvers with minimal change, exploiting these to provide an implicit solution for the time advance, which is necessary given the wide range of physical timescales 