diff --git a/src/chemistry/EAGLE/chemistry.h b/src/chemistry/EAGLE/chemistry.h index a470eef3fafc32c99fe3a853dcf46051a3086441..f3076750478aa669bfcb557676c1eb37b72eee21 100644 --- a/src/chemistry/EAGLE/chemistry.h +++ b/src/chemistry/EAGLE/chemistry.h @@ -197,6 +197,16 @@ __attribute__((always_inline)) INLINE static void chemistry_first_init_spart( sp->chemistry_data.metal_mass_fraction[elem] = data->initial_metal_mass_fraction[elem]; } + + /* Initialize mass fractions for total metals and each metal individually */ + if (data->initial_metal_mass_fraction_total != -1) { + sp->chemistry_data.smoothed_metal_mass_fraction_total = + data->initial_metal_mass_fraction_total; + + for (int elem = 0; elem < chemistry_element_count; ++elem) + sp->chemistry_data.smoothed_metal_mass_fraction[elem] = + data->initial_metal_mass_fraction[elem]; + } } /**