parameter_example.yml 27.6 KB
Newer Older
1
# Define the system of units to use internally.
2
InternalUnitSystem:
3
4
5
6
7
  UnitMass_in_cgs:     1   # Grams
  UnitLength_in_cgs:   1   # Centimeters
  UnitVelocity_in_cgs: 1   # Centimeters per second
  UnitCurrent_in_cgs:  1   # Amperes
  UnitTemp_in_cgs:     1   # Kelvin
Matthieu Schaller's avatar
Matthieu Schaller committed
8

9
10
11
12
# Values of some physical constants
PhysicalConstants:
  G:            6.67408e-8 # (Optional) Overwrite the value of Newton's constant used internally by the code.

13
14
15
16
17
18
19
20
21
22
23
24
# Cosmological parameters
Cosmology:
  h:              0.6777        # Reduced Hubble constant
  a_begin:        0.0078125     # Initial scale-factor of the simulation
  a_end:          1.0           # Final scale factor of the simulation
  Omega_m:        0.307         # Matter density parameter
  Omega_lambda:   0.693         # Dark-energy density parameter
  Omega_b:        0.0455        # Baryon density parameter
  Omega_r:        0.            # (Optional) Radiation density parameter
  w_0:            -1.0          # (Optional) Dark-energy equation-of-state parameter at z=0.
  w_a:            0.            # (Optional) Dark-energy equation-of-state time evolution parameter.

25
26
27
28
29
30
# Parameters for the hydrodynamics scheme
SPH:
  resolution_eta:        1.2348   # Target smoothing length in units of the mean inter-particle separation (1.2348 == 48Ngbs with the cubic spline kernel).
  CFL_condition:         0.1      # Courant-Friedrich-Levy condition for time integration.
  h_tolerance:           1e-4     # (Optional) Relative accuracy of the Netwon-Raphson scheme for the smoothing lengths.
  h_max:                 10.      # (Optional) Maximal allowed smoothing length in internal units. Defaults to FLT_MAX if unspecified.
31
  h_min_ratio:           0.       # (Optional) Minimal allowed smoothing length in units of the softening. Defaults to 0 if unspecified.
32
33
  max_volume_change:     1.4      # (Optional) Maximal allowed change of kernel volume over one time-step.
  max_ghost_iterations:  30       # (Optional) Maximal number of iterations allowed to converge towards the smoothing length.
34
35
  initial_temperature:   0        # (Optional) Initial temperature (in internal units) to set the gas particles at start-up. Value is ignored if set to 0.
  minimal_temperature:   0        # (Optional) Minimal temperature (in internal units) allowed for the gas particles. Value is ignored if set to 0.
36
37
  H_mass_fraction:       0.755    # (Optional) Hydrogen mass fraction used for initial conversion from temp to internal energy. Default value is derived from the physical constants.
  H_ionization_temperature: 1e4   # (Optional) Temperature of the transition from neutral to ionized Hydrogen for primoridal gas.
Josh Borrow's avatar
Josh Borrow committed
38
  viscosity_alpha:       0.8      # (Optional) Override for the initial value of the artificial viscosity. In schemes that have a fixed AV, this remains as alpha throughout the run.
Josh Borrow's avatar
Josh Borrow committed
39
40
41
42
43
44
45
46
  viscosity_alpha_max:   2.0      # (Optional) Maximal value for the artificial viscosity in schemes that allow alpha to vary.
  viscosity_alpha_min:   0.1      # (Optional) Minimal value for the artificial viscosity in schemes that allow alpha to vary.
  viscosity_length:      0.1      # (Optional) Decay length for the artificial viscosity in schemes that allow alpha to vary.
  diffusion_alpha:       0.0      # (Optional) Override the initial value for the thermal diffusion coefficient in schemes with thermal diffusion.
  diffusion_beta:        0.01     # (Optional) Override the decay/rise rate tuning parameter for the thermal diffusion.
  diffusion_alpha_max:   1.0      # (Optional) Override the maximal thermal diffusion coefficient that is allowed for a given particle.
  diffusion_alpha_min:   0.0      # (Optional) Override the minimal thermal diffusion coefficient that is allowed for a given particle.

47

48
49
# Parameters for the self-gravity scheme
Gravity:
50
  mesh_side_length:       32        # Number of cells along each axis for the periodic gravity mesh.
51
  eta:          0.025               # Constant dimensionless multiplier for time integration.
52
  theta:        0.7                 # Opening angle (Multipole acceptance criterion).
53
54
  comoving_softening:     0.0026994 # Comoving softening length (in internal units).
  max_physical_softening: 0.0007    # Physical softening length (in internal units).
55
  rebuild_frequency:      0.01      # (Optional) Frequency of the gravity-tree rebuild in units of the number of g-particles (this is the default value).
56
57
58
  a_smooth:     1.25                # (Optional) Smoothing scale in top-level cell sizes to smooth the long-range forces over (this is the default value).
  r_cut_max:    4.5                 # (Optional) Cut-off in number of top-level cells beyond which no FMM forces are computed (this is the default value).
  r_cut_min:    0.1                 # (Optional) Cut-off in number of top-level cells below which no truncation of FMM forces are performed (this is the default value).
59

60
# Parameters for the task scheduling
Matthieu Schaller's avatar
Matthieu Schaller committed
61
Scheduler:
62
63
64
65
66
67
  nr_queues:                 0         # (Optional) The number of task queues to use. Use 0  to let the system decide.
  cell_max_size:             8000000   # (Optional) Maximal number of interactions per task if we force the split (this is the default value).
  cell_sub_size_pair_hydro:  256000000 # (Optional) Maximal number of interactions per sub-pair hydro task  (this is the default value).
  cell_sub_size_self_hydro:  32000     # (Optional) Maximal number of interactions per sub-self hydro task  (this is the default value).
  cell_sub_size_pair_grav:   256000000 # (Optional) Maximal number of interactions per sub-pair gravity task  (this is the default value).
  cell_sub_size_self_grav:   32000     # (Optional) Maximal number of interactions per sub-self gravity task  (this is the default value).
68
69
  cell_sub_size_pair_stars:  256000000 # (Optional) Maximal number of interactions per sub-pair stars task  (this is the default value).
  cell_sub_size_self_stars:  32000     # (Optional) Maximal number of interactions per sub-self stars task  (this is the default value).
70
  cell_split_size:           400       # (Optional) Maximal number of particles per cell (this is the default value).
71
  cell_subdepth_diff_grav:   4         # (Optional) Maximal depth difference between leaves and a cell that gravity tasks can be pushed down to (this is the default value).
72
73
74
  cell_extra_parts:          0         # (Optional) Number of spare parts per top-level allocated at rebuild time for on-the-fly creation.
  cell_extra_gparts:         0         # (Optional) Number of spare gparts per top-level allocated at rebuild time for on-the-fly creation.
  cell_extra_sparts:         400       # (Optional) Number of spare sparts per top-level allocated at rebuild time for on-the-fly creation.
75
76
  max_top_level_cells:       12        # (Optional) Maximal number of top-level cells in any dimension. The number of top-level cells will be the cube of this (this is the default value).
  tasks_per_cell:            0         # (Optional) The average number of tasks per cell. If not large enough the simulation will fail (means guess...).
77
  links_per_tasks:           10        # (Optional) The average number of links per tasks (before adding the communication tasks). If not large enough the simulation will fail (means guess...). Defaults to 10.
78
  mpi_message_limit:         4096      # (Optional) Maximum MPI task message size to send non-buffered, KB.
Matthieu Schaller's avatar
Matthieu Schaller committed
79

80
# Parameters governing the time integration (Set dt_min and dt_max to the same value for a fixed time-step run.)
Matthieu Schaller's avatar
Matthieu Schaller committed
81
TimeIntegration:
82
83
84
85
86
  time_begin:        0.    # The starting time of the simulation (in internal units).
  time_end:          1.    # The end time of the simulation (in internal units).
  dt_min:            1e-6  # The minimal time-step size of the simulation (in internal units).
  dt_max:            1e-2  # The maximal time-step size of the simulation (in internal units).
  max_dt_RMS_factor: 0.25  # (Optional) Dimensionless factor for the maximal displacement allowed based on the RMS velocities.
87

88
89
# Parameters governing the snapshots
Snapshots:
90
  basename:   output      # Common part of the name of output files
91
92
  scale_factor_first: 0.1 # (Optional) Scale-factor of the first snapshot if cosmological time-integration.
  time_first: 0.          # (Optional) Time of the first output if non-cosmological time-integration (in internal units)
93
  delta_time: 0.01        # Time difference between consecutive outputs (in internal units)
94
  invoke_stf: 0           # (Optional) Call VELOCIraptor every time a snapshot is written irrespective of the VELOCIraptor output strategy.
95
  compression: 0          # (Optional) Set the level of compression of the HDF5 datasets [0-9]. 0 does no compression.
96
  int_time_label_on:   0  # (Optional) Enable to label the snapshots using the time rounded to an integer (in internal units)
97
98
99
100
101
  UnitMass_in_cgs:     1  # (Optional) Unit system for the outputs (Grams)
  UnitLength_in_cgs:   1  # (Optional) Unit system for the outputs (Centimeters)
  UnitVelocity_in_cgs: 1  # (Optional) Unit system for the outputs (Centimeters per second)
  UnitCurrent_in_cgs:  1  # (Optional) Unit system for the outputs (Amperes)
  UnitTemp_in_cgs:     1  # (Optional) Unit system for the outputs (Kelvin)
102
103
  output_list_on:      0  # (Optional) Enable the output list
  output_list:         snaplist.txt # (Optional) File containing the output times (see documentation in "Parameter File" section)
lhausamm's avatar
lhausamm committed
104
105
106

# Parameters governing the logger snapshot system
Logger:
107
108
109
110
  delta_step:           10     # Update the particle log every this many updates
  initial_buffer_size:  1      # buffer size in GB
  buffer_scale:		10     # (Optional) When buffer size is too small, update it with required memory times buffer_scale
  basename:             index  # Common part of the filenames
111
  
112
113
# Parameters governing the conserved quantities statistics
Statistics:
114
115
116
  delta_time:           1e-2     # Time between statistics output
  scale_factor_first:     0.1    # (Optional) Scale-factor of the first statistics dump if cosmological time-integration.
  time_first:             0.     # (Optional) Time of the first stats output if non-cosmological time-integration (in internal units)
Matthieu Schaller's avatar
Matthieu Schaller committed
117
  energy_file_name:    energy    # (Optional) File name for energy output
118
  timestep_file_name:  timesteps # (Optional) File name for timing information output. Note: No underscores "_" allowed in file name
119
120
  output_list_on:      0   	 # (Optional) Enable the output list
  output_list:         statlist.txt # (Optional) File containing the output times (see documentation in "Parameter File" section)
121

122
# Parameters related to the initial conditions
Matthieu Schaller's avatar
Matthieu Schaller committed
123
InitialConditions:
124
  file_name:  SedovBlast/sedov.hdf5 # The file to read
Peter W. Draper's avatar
typo    
Peter W. Draper committed
125
  periodic:                    1    # Are we running with periodic ICs?
126
  generate_gas_in_ics:         0    # (Optional) Generate gas particles from the DM-only ICs (e.g. from panphasia).
127
  cleanup_h_factors:           0    # (Optional) Clean up the h-factors used in the ICs (e.g. in Gadget files).
128
  cleanup_velocity_factors:    0    # (Optional) Clean up the scale-factors used in the definition of the velocity variable in the ICs (e.g. in Gadget files).
129
130
  cleanup_smoothing_lengths:   0    # (Optional) Clean the values of the smoothing lengths that are read in to remove stupid values. Set to 1 to activate.
  smoothing_length_scaling:    1.   # (Optional) A scaling factor to apply to all smoothing lengths in the ICs.
131
  shift:      [0.0,0.0,0.0]         # (Optional) A shift to apply to all particles read from the ICs (in internal units).
132
  replicate:  2                     # (Optional) Replicate all particles along each axis a given integer number of times. Default 1.
Matthieu Schaller's avatar
Matthieu Schaller committed
133

134
135
# Parameters controlling restarts
Restarts:
136
137
138
139
140
141
142
143
144
  enable:             1          # (Optional) whether to enable dumping restarts at fixed intervals.
  save:               1          # (Optional) whether to save copies of the previous set of restart files (named .prev)
  onexit:             0          # (Optional) whether to dump restarts on exit (*needs enable*)
  subdir:             restart    # (Optional) name of subdirectory for restart files.
  basename:           swift      # (Optional) prefix used in naming restart files.
  delta_hours:        6.0        # (Optional) decimal hours between dumps of restart files.
  stop_steps:         100        # (Optional) how many steps to process before checking if the <subdir>/stop file exists. When present the application will attempt to exit early, dumping restart files first.
  max_run_time:       24.0       # (optional) Maximal wall-clock time in hours. The application will exit when this limit is reached.
  resubmit_on_exit:   0          # (Optional) whether to run a command when exiting after the time limit has been reached.
145
  resubmit_command:   ./resub.sh # (Optional) Command to run when time limit is reached. Compulsory if resubmit_on_exit is switched on. Note potentially unsafe.
146

147
# Parameters governing domain decomposition
Matthieu Schaller's avatar
Matthieu Schaller committed
148
DomainDecomposition:
149
150
151
  initial_type:     memory    # (Optional) The initial decomposition strategy: "grid",
                              #            "region", "memory", or "vectorized".
  initial_grid: [10,10,10]    # (Optional) Grid sizes if the "grid" strategy is chosen.
152

153
154
155
156
157
158
159
160
161
162
163
  repartition_type: fullcosts # (Optional) The re-decomposition strategy, one of:
                              # "none", "fullcosts", "edgecosts", "memory" or
                              # "timecosts".
  trigger:          0.05      # (Optional) Fractional (<1) CPU time difference between MPI ranks required to trigger a
                              # new decomposition, or number of steps (>1) between decompositions
  minfrac:          0.9       # (Optional) Fractional of all particles that should be updated in previous step when
                              # using CPU time trigger
  usemetis:         0         # Use serial METIS when ParMETIS is also available.
  adaptive:         1         # Use adaptive repartition when ParMETIS is available, otherwise simple refinement.
  itr:              100       # When adaptive defines the ratio of inter node communication time to data redistribution time, in the range 0.00001 to 10000000.0.
                              # Lower values give less data movement during redistributions, at the cost of global balance which may require more communication.
164
  use_fixed_costs:  0         # If 1 then use any compiled in fixed costs for
165
166
                              # task weights in first repartition, if 0 only use task timings, if > 1 only use
                              # fixed costs, unless none are available.
167

168
169
170
171
# Structure finding options (requires velociraptor)
StructureFinding:
  config_file_name:     stf_input.cfg # Name of the STF config file.
  basename:             ./stf         # Common part of the name of output files.
172
173
174
175
176
  scale_factor_first:   0.92          # (Optional) Scale-factor of the first snaphot (cosmological run)
  time_first:           0.01          # (Optional) Time of the first structure finding output (in internal units).
  delta_time:           1.10          # (Optional) Time difference between consecutive structure finding outputs (in internal units) in simulation time intervals.
  output_list_on:       0   	      # (Optional) Enable the output list
  output_list:          stflist.txt   # (Optional) File containing the output times (see documentation in "Parameter File" section)
177

178
179
180
# Parameters related to the equation of state ------------------------------------------

EoS:
181
182
  isothermal_internal_energy: 20.26784  # Thermal energy per unit mass for the case of isothermal equation of state (in internal units).

183
184
185
186
  planetary_use_Til:    1   # (Optional) Whether to prepare the Tillotson EOS
  planetary_use_HM80:   0   # (Optional) Whether to prepare the Hubbard & MacFarlane (1980) EOS
  planetary_use_ANEOS:  0   # (Optional) Whether to prepare the ANEOS EOS
  planetary_use_SESAME: 0   # (Optional) Whether to prepare the SESAME EOS
187
                            # (Optional) Table file paths
Jacob Kegerreis 's avatar
Jacob Kegerreis committed
188
189
190
191
192
193
194
  planetary_HM80_HHe_table_file:        ./EoSTables/planetary_HM80_HHe.txt
  planetary_HM80_ice_table_file:        ./EoSTables/planetary_HM80_ice.txt
  planetary_HM80_rock_table_file:       ./EoSTables/planetary_HM80_rock.txt
  planetary_SESAME_iron_table_file:     ./EoSTables/planetary_SESAME_iron_2140.txt
  planetary_SESAME_basalt_table_file:   ./EoSTables/planetary_SESAME_basalt_7530.txt
  planetary_SESAME_water_table_file:    ./EoSTables/planetary_SESAME_water_7154.txt
  planetary_SS08_water_table_file:      ./EoSTables/planetary_SS08_water.txt
195

196
# Parameters related to external potentials --------------------------------------------
197

198
# Point mass external potentials
199
PointMassPotential:
200
  useabspos:       0        # 0 -> positions based on centre, 1 -> absolute positions 
201
202
203
204
  position:        [50.,50.0,50.]      # location of external point mass (internal units)
  mass:            1e10                # mass of external point mass (internal units)
  timestep_mult:   0.03                # Dimensionless pre-factor for the time-step condition
  softening:       0.05                # For point-mass-softened option
205

206
# Isothermal potential parameters
207
IsothermalPotential:
208
  useabspos:       0        # 0 -> positions based on centre, 1 -> absolute positions 
209
  position:        [100.,100.,100.]    # Location of centre of isothermal potential with respect to centre of the box (internal units)
210
211
  vrot:            200.     # Rotation speed of isothermal potential (internal units)
  timestep_mult:   0.03     # Dimensionless pre-factor for the time-step condition
212
  epsilon:         0.1      # Softening size (internal units)
213
214
215
216
217
  
# Hernquist potential parameters
HernquistPotential:
  useabspos:       0        # 0 -> positions based on centre, 1 -> absolute positions 
  position:        [100.,100.,100.]    # Location of centre of isothermal potential with respect to centre of the box (if 0) otherwise absolute (if 1) (internal units)
218
219
220
  idealizeddisk:   0        # (Optional) Whether to run with idealizeddisk or without, 0 used the mass and scalelength as mandatory parameters, while 1 uses more advanced disk dependent paramters
  mass:            1e10     # (Optional 0) default parameter, Mass of the Hernquist potential
  scalelength:     10.0     # (Optional 0) default parameter, Scale length of the potential
221
                            # If multiple X200 values are given, only one is used, in the order M200 > V200 > R200.
222
  M200:            3e11     # (Optional 1a) M200 of the galaxy+halo (when used V200 and R200 are not used)
Folkert Nobels's avatar
Folkert Nobels committed
223
224
  V200:            100.     # (Optional 1b) V200 of the galaxy+halo (when used M200 and R200 are not used, if M200 is given M200 is used)
  R200:            10.      # (Optional 1c) R200 of the galaxy+halo (when used M200 and V200 are not used, if M200 or V200 are given they are used)
225
226
  h:               0.704    # (Optional 1) reduced Hubble constant
  concentration:   7.1      # (Optional 1) concentration of the Halo
227
228
  diskfraction:              0.0434370991372   # (Optional 1) Disk mass fraction (equal to MD in MakeNewDisk and GalIC)
  bulgefraction:              0.00705852860979  # (Optional 1) Bulge mass fraction (equal to MB in MakeNewDisk and GalIC)
229
230
231
232
233
234
235
236
237
238
239
  timestep_mult:   0.01     # Dimensionless pre-factor for the time-step condition, basically determines the fraction of the orbital time we use to do the time integration
  epsilon:         0.1      # Softening size (internal units)
 
# Isothermal potential parameters
NFWPotential:
  useabspos:          0
  position:           [0.0,0.0,0.0]      # Location of centre of isothermal potential with respect to centre of the box (internal units) if useabspos=0 otherwise with respect to the 0,0,0, coordinates.
  concentration:      8.       # Concentration of the halo
  M_200:              2.0e+12  # Mass of the halo (M_200 in internal units)
  critical_density:   127.4    # Critical density (internal units).
  timestep_mult:      0.01     # Dimensionless pre-factor for the time-step condition, basically determines fraction of orbital time we need to do an integration step
240

241
# Disk-patch potential parameters
242
243
244
DiscPatchPotential:
  surface_density: 10.      # Surface density of the disc (internal units)
  scale_height:    100.     # Scale height of the disc (internal units)
245
246
247
  z_disc:          400.     # Position of the disc along the z-axis (internal units)
  z_trunc:         300.     # (Optional) Distance from the disc along z-axis above which the potential gets truncated.
  z_max:           380.     # (Optional) Distance from the disc along z-axis above which the potential is set to 0.
248
  timestep_mult:   0.03     # Dimensionless pre-factor for the time-step condition
249
  growth_time:     5.       # (Optional) Time for the disc to grow to its final size (multiple of the dynamical time)
250

251
252
253
254
255
# Sine Wave potential
SineWavePotential:
  amplitude:        10.     # Amplitude of the sine wave (internal units)
  timestep_limit:   1.      # Time-step dimensionless pre-factor.
  growth_time:      0.      # (Optional) Time for the potential to grow to its final size.
256

257
258
259
260
261
262
263
# Parameters related to entropy floors    ----------------------------------------------

EAGLEEntropyFloor:
  Jeans_density_threshold_H_p_cm3: 0.1       # Physical density above which the EAGLE Jeans limiter entropy floor kicks in expressed in Hydrogen atoms per cm^3.
  Jeans_over_density_threshold:    10.       # Overdensity above which the EAGLE Jeans limiter entropy floor can kick in.
  Jeans_temperature_norm_K:        8000      # Temperature of the EAGLE Jeans limiter entropy floor at the density threshold expressed in Kelvin.
  Jeans_gamma_effective:           1.3333333 # Slope the of the EAGLE Jeans limiter entropy floor
264
265
266
267
  Cool_density_threshold_H_p_cm3:  1e-5      # Physical density above which the EAGLE Cool limiter entropy floor kicks in expressed in Hydrogen atoms per cm^3.
  Cool_over_density_threshold:     10.       # Overdensity above which the EAGLE Cool limiter entropy floor can kick in.
  Cool_temperature_norm_K:         8000      # Temperature of the EAGLE Cool limiter entropy floor at the density threshold expressed in Kelvin.
  Cool_gamma_effective:            1.        # Slope the of the EAGLE Cool limiter entropy floor
268
  
269
270
271
272
273
274
275
276
277
278
# Parameters related to cooling function  ----------------------------------------------

# Constant du/dt cooling function
ConstCooling:
  cooling_rate: 1.          # Cooling rate (du/dt) (internal units)
  min_energy:   1.          # Minimal internal energy per unit mass (internal units)
  cooling_tstep_mult: 1.    # Dimensionless pre-factor for the time-step condition

# Constant lambda cooling function
LambdaCooling:
279
  lambda_nH2_cgs:              1e-22 # Cooling rate divided by square Hydrogen number density (in cgs units [erg * s^-1 * cm^3])
280
  cooling_tstep_mult:          1.0   # (Optional) Dimensionless pre-factor for the time-step condition.
281

282
# Parameters of the EAGLE cooling model (Wiersma+08 cooling tables).
283
EAGLECooling:
284
  dir_name:                  ./coolingtables/  # Location of the Wiersma+08 cooling tables
285
  H_reion_z:                 8.5               # Redshift of Hydrogen re-ionization
286
287
  He_reion_z_centre:         3.5               # Redshift of the centre of the Helium re-ionization Gaussian
  He_reion_z_sigma:          0.5               # Spread in redshift of the  Helium re-ionization Gaussian
288
289
290
291
  He_reion_eV_p_H:           2.0               # Energy inject by Helium re-ionization in electron-volt per Hydrogen atom
  newton_integration:        0                 # (Optional) Set to 1 to use the Newton-Raphson method to solve the xplicit cooling problem.
  Ca_over_Si_in_solar:       1.                # (Optional) Ratio of Ca/Si to use in units of solar. If set to 1, the code uses [Ca/Si] = 0, i.e. Ca/Si = 0.0941736.
  S_over_Si_in_solar:        1.                # (Optional) Ratio of S/Si to use in units of solar. If set to 1, the code uses [S/Si] = 0, i.e. S/Si = 0.6054160.
292
  
293
# Cooling with Grackle 3.0
lhausamm's avatar
lhausamm committed
294
GrackleCooling:
295
  CloudyTable: CloudyData_UVB=HM2012.h5 # Name of the Cloudy Table (available on the grackle bitbucket repository)
Matthieu Schaller's avatar
Matthieu Schaller committed
296
297
298
299
300
301
302
303
304
  WithUVbackground: 1                   # Enable or not the UV background
  Redshift: 0                           # Redshift to use (-1 means time based redshift)
  WithMetalCooling: 1                   # Enable or not the metal cooling
  ProvideVolumetricHeatingRates: 0      # (optional) User provide volumetric heating rates
  ProvideSpecificHeatingRates: 0        # (optional) User provide specific heating rates
  SelfShieldingMethod: 0                # (optional) Grackle (<= 3) or Gear self shielding method
  OutputMode: 0                         # (optional) Write in output corresponding primordial chemistry mode
  MaxSteps: 10000                       # (optional) Max number of step when computing the initial composition
  ConvergenceLimit: 1e-2                # (optional) Convergence threshold (relative) for initial composition
305
306
307
308
309

# Parameters related to chemistry models  -----------------------------------------------

# EAGLE model
EAGLEChemistry:
310
311
312
313
314
315
316
317
318
319
  init_abundance_metal:     0.           # Inital fraction of particle mass in *all* metals
  init_abundance_Hydrogen:  0.752        # Inital fraction of particle mass in Hydrogen
  init_abundance_Helium:    0.248        # Inital fraction of particle mass in Helium
  init_abundance_Carbon:    0.000        # Inital fraction of particle mass in Carbon
  init_abundance_Nitrogen:  0.000        # Inital fraction of particle mass in Nitrogen
  init_abundance_Oxygen:    0.000        # Inital fraction of particle mass in Oxygen
  init_abundance_Neon:      0.000        # Inital fraction of particle mass in Neon
  init_abundance_Magnesium: 0.000        # Inital fraction of particle mass in Magnesium
  init_abundance_Silicon:   0.000        # Inital fraction of particle mass in Silicon
  init_abundance_Iron:      0.000        # Inital fraction of particle mass in Iron
320

321
322
# Parameters related to star formation models  -----------------------------------------------

323
324
325
326
327
328
329
330
331
332
333
334
# EAGLE star formation model (Schaye and Dalla Vecchia 2008)
EAGLEStarFormation:
  EOS_density_norm_H_p_cm3:          0.1       # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3.
  EOS_temperature_norm_K:            8000      # Temperature om the polytropic EOS assumed for star-forming gas at the density normalisation in Kelvin.
  EOS_gamma_effective:               1.3333333 # Slope the of the polytropic EOS assumed for the star-forming gas.
  gas_fraction:                      0.25      # (Optional) The gas fraction used internally by the model (Defaults to 1).
  KS_normalisation:                  1.515e-4  # The normalization of the Kennicutt-Schmidt law in Msun / kpc^2 / yr.
  KS_exponent:                       1.4       # The exponent of the Kennicutt-Schmidt law.
  KS_min_over_density:               57.7      # The over-density above which star-formation is allowed.
  KS_high_density_threshold_H_p_cm3: 1e3       # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3.
  KS_high_density_exponent:          2.0       # Slope of the Kennicut-Schmidt law above the high-density threshold.
  KS_max_density_threshold_H_p_cm3:  1e5       # (Optional) Density above which a gas particle gets automatically turned into a star in Hydrogen atoms per cm^3 (Defaults to FLT_MAX).
335
  KS_temperature_margin_dex:         0.5       # (Optional) Logarithm base 10 of the maximal temperature difference above the EOS allowed to form stars (Defaults to FLT_MAX).
336
337
338
339
340
  threshold_norm_H_p_cm3:            0.1       # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3.
  threshold_Z0:                      0.002     # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation.
  threshold_slope:                   -0.64     # Slope of the metal-dependant star formation threshold
  threshold_max_density_H_p_cm3:     10.0      # Maximal density of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3.