small_cosmo_volume.yml 6.76 KB
Newer Older
1
2
# Define the system of units to use internally. 
InternalUnitSystem:
3
  UnitMass_in_cgs:     1.98841e43    # 10^10 M_sun
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
  UnitLength_in_cgs:   3.08567758e24 # 1 Mpc
  UnitVelocity_in_cgs: 1e5           # 1 km/s
  UnitCurrent_in_cgs:  1             # Amperes
  UnitTemp_in_cgs:     1             # Kelvin

Cosmology:                      # WMAP9 cosmology
  Omega_m:        0.276
  Omega_lambda:   0.724
  Omega_b:        0.0455
  h:              0.703
  a_begin:        0.019607843	# z_ini = 50.
  a_end:          1.0		# z_end = 0.

# Parameters governing the time integration
TimeIntegration:
19
  dt_min:     1e-8
20
21
22
23
  dt_max:     1e-2 

# Parameters for the self-gravity scheme
Gravity:
24
  eta:          0.025
25
  MAC:          adaptive
26
27
  theta_cr:     0.7
  epsilon_fmm:  0.001
28
29
30
31
  comoving_DM_softening:         0.0889     # 1/25th of the mean inter-particle separation: 88.9 kpc
  max_physical_DM_softening:     0.0889     # 1/25th of the mean inter-particle separation: 88.9 kpc
  comoving_baryon_softening:     0.0889     # 1/25th of the mean inter-particle separation: 88.9 kpc
  max_physical_baryon_softening: 0.0889     # 1/25th of the mean inter-particle separation: 88.9 kpc
32
33
34
35
36
  mesh_side_length:       64

# Parameters of the hydro scheme
SPH:
  resolution_eta:      1.2348   # "48 Ngb" with the cubic spline kernel
37
  h_min_ratio:         0.1
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
  CFL_condition:       0.1
  initial_temperature: 7075.    # (1 + z_ini)^2 * 2.72K
  minimal_temperature: 100.

# Parameters governing the snapshots
Snapshots:
  basename:            snap
  delta_time:          1.02
  scale_factor_first:  0.02
  
# Parameters governing the conserved quantities statistics
Statistics:
  delta_time:          1.02
  scale_factor_first:  0.02
  
Scheduler:
  max_top_level_cells: 8
  cell_split_size:     50
56
57
  cell_extra_gparts:         100         # (Optional) Number of spare gparts per top-level allocated at rebuild time for on-the-fly creation.
  cell_extra_sparts:         100       # (Optional) Number of spare sparts per top-level allocated at rebuild time for on-the-fly creation.
58
59
60
61
  
# Parameters related to the initial conditions
InitialConditions:
  file_name:  small_cosmo_volume.hdf5
62
  periodic:                    1
63
64
  cleanup_h_factors:           1    
  cleanup_velocity_factors:    1  
65
66
  generate_gas_in_ics:         1    # Generate gas particles from the DM-only ICs
  cleanup_smoothing_lengths:   1    # Since we generate gas, make use of the (expensive) cleaning-up procedure.
67
68
69

# Constant lambda cooling function
LambdaCooling:
70
  lambda_nH2_cgs:              1e-26 # Cooling rate divided by square Hydrogen number density (in cgs units [erg * s^-1 * cm^3])
71
72

# EAGLE cooling function
73
EAGLECooling:
74
  dir_name:                 ./coolingtables/
75
  H_reion_z:               11.5
76
  H_reion_eV_p_H:          2.0
77
78
  He_reion_z_centre:       3.5
  He_reion_z_sigma:        0.5
79
  He_reion_eV_p_H:         2.0
80

Matthieu Schaller's avatar
Matthieu Schaller committed
81
82
83
84
85
86
87
88
89
90
91
# COLIBRE cooling parameters
COLIBRECooling:
  dir_name:                ./UV_dust1_CR1_G1_shield1.hdf5
  H_reion_z:               11.5
  H_reion_eV_p_H:          2.0
  He_reion_z_centre:       3.5
  He_reion_z_sigma:        0.5
  He_reion_eV_p_H:         2.0
  delta_logTEOS_subgrid_properties: 0.5
  rapid_cooling_threshold:          0.333333

92
93
# Impose primoridal metallicity
EAGLEChemistry:
94
95
96
97
98
99
100
101
102
103
  init_abundance_metal:     0.0
  init_abundance_Hydrogen:  0.752
  init_abundance_Helium:    0.248
  init_abundance_Carbon:    0.0
  init_abundance_Nitrogen:  0.0
  init_abundance_Oxygen:    0.0
  init_abundance_Neon:      0.0
  init_abundance_Magnesium: 0.0
  init_abundance_Silicon:   0.0
  init_abundance_Iron:      0.0
104

105
106
107
108
109
110
111
112
113
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
  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
114
115
116
117
118

GEARPressureFloor:
  jeans_factor: 10.       # Number of particles required to suppose a resolved clump and avoid the pressure floor.


119
120
# Cooling with Grackle 3.0
GrackleCooling:
121
122
123
124
125
126
127
128
129
130
131
132
  cloudy_table: CloudyData_UVB=HM2012.h5       # Name of the Cloudy Table (available on the grackle bitbucket repository)
  with_UV_background: 1                        # Enable or not the UV background
  redshift: 0                                  # Redshift to use (-1 means time based redshift)
  with_metal_cooling: 1                        # Enable or not the metal cooling
  provide_volumetric_heating_rates: 0          # (optional) User provide volumetric heating rates
  provide_specific_heating_rates: 0            # (optional) User provide specific heating rates
  max_steps: 10000                             # (optional) Max number of step when computing the initial composition
  convergence_limit: 1e-2                      # (optional) Convergence threshold (relative) for initial composition
  thermal_time_myr: 5                          # (optional) Time (in Myr) for adiabatic cooling after a feedback event.
  self_shielding_method: -1                    # (optional) Grackle (1->3 for Grackle's ones, 0 for none and -1 for GEAR)
  self_shielding_threshold_atom_per_cm3: 0.007 # Required only with GEAR's self shielding. Density threshold of the self shielding

133

134
# GEAR chemistry model (Revaz and Jablonka 2018)
135
GEARChemistry:
136
137
  initial_metallicity: 1         # Initial metallicity of the gas (mass fraction)
  scale_initial_metallicity: 1   # Should we scale the initial metallicity with the solar one?
Loic Hausammann's avatar
Loic Hausammann committed
138

139
140
141
142
143
144
145
146
147
148
149
150
  # GEAR star formation model (Revaz and Jablonka 2018)
GEARStarFormation:
  star_formation_efficiency: 0.01   # star formation efficiency (c_*)
  maximal_temperature:  3e4         # Upper limit to the temperature of a star forming particle
  n_stars_per_particle: 4
  min_mass_frac: 0.5

# GEAR feedback model
GEARFeedback:
  supernovae_energy_erg: 0.1e51                            # Energy released by a single supernovae.
  yields_table: chemistry-AGB+OMgSFeZnSrYBaEu-16072013.h5  # Table containing the yields.
  discrete_yields: 0                                       # Should we use discrete yields or the IMF integrated one?