[BUGFIX] Correcting the timestep sent to grackle with RT

#391 Merged at 37f9795
Repository
jwise77
Branch
week-of-code
Repository
enzo
Branch
week-of-code
Author
  1. John Wise
Reviewers
Description
  • Using grackle data when computing the RT timestep.
  • [BUGFIX] In the intermediate RT grackle call, the RT timestep is the proper value instead of the hydro timestep.
  • Commit status

Comments (5)

  1. Britton Smith

    Great work, John! Nice catch on using the Grackle rates for the timestep as well. So glad to see this solved.

  2. John Wise author

    Here is a before/after image of the ionization front propagation with this fix. Now it's nearly identical to Enzo's internal solver.

    It lags behind the analytical solution for known reasons (see Wise & Abel 2011).compare.png

  3. John Regan

    Nice! I noticed this bug sometime ago when checking my fork of Grackle and updated it. I assumed that I introduced though - I didn't realise it was in the original code. Good to get it sorted.