Skip to content
Snippets Groups Projects
Commit d3b1994c authored by Mladen Ivkovic's avatar Mladen Ivkovic Committed by Matthieu Schaller
Browse files

GEARRT tiny fix.

parent f200b0dc
No related branches found
No related tags found
5 merge requests!1715Update planetary strength after planetary plus's master rebase,!1693More threapool plotting tweaks,!1668before Mag.Egy in all the flavors,!1662Initial sync from previous months,!1648GEARRT tiny fix.
...@@ -30,7 +30,7 @@ Compiling for GEAR RT ...@@ -30,7 +30,7 @@ Compiling for GEAR RT
you to select a hydro Riemann solver, e.g ``--with-riemann-solver=hllc``. you to select a hydro Riemann solver, e.g ``--with-riemann-solver=hllc``.
- The thermochemistry requires the `grackle <https://github.com/grackle-project/grackle>`_ - The thermochemistry requires the `grackle <https://github.com/grackle-project/grackle>`_
library. Grackle is a chemistry and cooling library presented in library version 3.2. Grackle is a chemistry and cooling library presented in
`B. Smith et al. 2017 <https://ui.adsabs.harvard.edu/abs/2017MNRAS.466.2217S>`_. `B. Smith et al. 2017 <https://ui.adsabs.harvard.edu/abs/2017MNRAS.466.2217S>`_.
Please note that the current implementation is not (yet) as Please note that the current implementation is not (yet) as
advanced as the :ref:`GEAR subgrid model grackle cooling <gear_grackle_cooling>`, advanced as the :ref:`GEAR subgrid model grackle cooling <gear_grackle_cooling>`,
......
...@@ -239,17 +239,17 @@ rt_tchem_set_particle_radiation_field_for_test( ...@@ -239,17 +239,17 @@ rt_tchem_set_particle_radiation_field_for_test(
const double time_to_cgs = units_cgs_conversion_factor(us, UNIT_CONV_TIME); const double time_to_cgs = units_cgs_conversion_factor(us, UNIT_CONV_TIME);
const double t_Myr = time * time_to_cgs / (3600. * 24. * 365. * 1e6); const double t_Myr = time * time_to_cgs / (3600. * 24. * 365. * 1e6);
double fixed_fluxes[RT_NGROUPS]; /* NOTE: this assumes that the test is set up with 3 photon groups. */
for (int g = 0; g < RT_NGROUPS; g++) fixed_fluxes[g] = 0.; double fixed_fluxes[3];
for (int g = 0; g < 3; g++) fixed_fluxes[g] = 0.;
if (t_Myr < 0.5) { if (t_Myr < 0.5) {
/* Be vocal, just in case somebody forgets you exist. */ /* Be vocal, just in case somebody forgets you exist. */
if (p->id == 1) message("Setting fixed radiation field."); if (p->id == 1) message("Setting fixed radiation field.");
/* Set fixed radiation fields, in cgs*/ /* Set fixed radiation fields, in cgs*/
fixed_fluxes[0] = 0.; fixed_fluxes[0] = 1.350e01;
fixed_fluxes[1] = 1.350e01; fixed_fluxes[1] = 2.779e01;
fixed_fluxes[2] = 2.779e01; fixed_fluxes[2] = 6.152e00;
fixed_fluxes[3] = 6.152e00;
} }
const double flux_to_cgs = const double flux_to_cgs =
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment