Skip to content
Snippets Groups Projects
Commit 6485b6a9 authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Do not ignore the BH reads from the ICs.

parent 4440b85f
No related branches found
No related tags found
1 merge request!543Fof
......@@ -69,7 +69,7 @@ FOF:
linking_length_ratio: 0.2 # Linking length in units of the main inter-particle separation.
black_hole_seed_halo_mass_Msun: 1.5e10 # Minimal halo mass in which to seed a black hole (in solar masses).
absolute_linking_length: -1. # (Optional) Absolute linking length (in internal units)
run_freq: 100 # (Optional) The no. of steps between each FOF search.
run_freq: 10 # (Optional) The no. of steps between each FOF search.
group_id_default: 2147483647 # (Optional) Sets the group ID of particles in groups below the minimum size. Defaults to 2^31 - 1 if unspecified. Has to be positive.
group_id_offset: 1 # (Optional) Sets the offset of group ID labeling. Defaults to 1 if unspecified.
......@@ -80,6 +80,10 @@ InitialConditions:
cleanup_h_factors: 1 # Remove the h-factors inherited from Gadget
cleanup_velocity_factors: 1 # Remove the sqrt(a) factor in the velocities inherited from Gadget
# EAGLE black hole seeding parameters
EAGLEBlackHoles:
seed_halo_mass: 1.5e10 # The minimum group mass required for black hole seeding
EAGLEChemistry: # Solar abundances
init_abundance_metal: 0.014
init_abundance_Hydrogen: 0.70649785
......@@ -171,4 +175,3 @@ EAGLEAGN:
coupling_efficiency: 0.15 # Fraction of the radiated energy that couples to the gas in feedback events.
AGN_delta_T_K: 3.16228e8 # Change in temperature to apply to the gas particle in an AGN feedback event in Kelvin.
AGN_num_ngb_to_heat: 1. # Target number of gas neighbours to heat in an AGN feedback event.
......@@ -444,8 +444,6 @@ void read_ic_single(const char* fileName,
for (int ptype = 0; ptype < swift_type_count; ++ptype)
N[ptype] = (numParticles[ptype]) + (numParticles_highWord[ptype] << 32);
N[5] = 0;
/* Get the box size if not cubic */
dim[0] = boxSize[0];
dim[1] = (boxSize[1] < 0) ? boxSize[0] : boxSize[1];
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment