Changes to the chemistry / energy solver

#289 Merged at 1d9d8bc
Repository
jwise77
Branch
week-of-code
Repository
enzo
Branch
week-of-code
Author
  1. John Wise
Reviewers
Description
  • Adding a timestep limiter based on cooling time. New parameters: UseCoolingTimestep, CoolingTimestepSafetyFactor (default = 0.1).

  • Use RadiationFieldRedshift for Compton cooling in non-cosmological runs.

  • Adding H2 self-shielding according to the Wolcott-Green et al. (2011) Sobolev-like, spherically averaged prescription. Turned on with RadiationShield = 2.

  • Restricting high-density dn_H/dt timestep limiter to cells without ionizing radiation.

Comments (8)

    1. John Wise author

      I think the failures were caused by moving timestep limiter to inside the high-density if-block. In hindsight, this isn't necessary, and I've reverted this change. Let's see what the tests reveal.

  1. James Larrue

    There is my one comment on the documentation and the conflict reported by bitbucket. Once those are addressed, I will approve.

  2. Brian OShea

    All issues have been dealt with and it passes the test suite. I'm approving and merging.