PDR sims affected by induced processes
tests reveal unexplained behavior in PDR sims with no ionizing radiation. Modify pdf_leiden_f4 to stop in zone 2. In its default state we get this hydrogen ionization for these two zones:
{{{ HI HII 1.00E+00 5.18E-08 1.00E+00 3.51E-05 }}}
the ionization goes up by assert dex.
Add the command "no induced processes" and we get {{{ HI HII 1.00E+00 5.18E-08 1.00E+00 5.19E-08 }}}
the answer does not change and is equal to the first zone, with induced processes. It is counterintuitive that induced processes would make the ionization higher in the second zone, which should be more shielded.