santa_barbara.yml 2.61 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
# Define the system of units to use internally. 
InternalUnitSystem:
  UnitMass_in_cgs:     1.98848e43    # 10^10 Msun 
  UnitLength_in_cgs:   3.08567758e24 # 1 Mpc 
  UnitVelocity_in_cgs: 1e5           # 1 km/s 
  UnitCurrent_in_cgs:  1             # Amperes
  UnitTemp_in_cgs:     1             # Kelvin

# Cosmological parameters
Cosmology:
  h:              0.5        
  a_begin:        0.047619048        # z_ini = 20
  a_end:          1.0                # z_end = 0
  Omega_m:        1.0        
  Omega_lambda:   0.0        
  Omega_b:        0.1        
  
# Parameters governing the time integration
TimeIntegration:
  dt_max:     0.01
  dt_min:     1e-10

Scheduler:
  max_top_level_cells: 16
  cell_split_size:     100

# Parameters governing the snapshots
Snapshots:
  basename:            santabarbara 
  scale_factor_first:  0.05
  delta_time:          1.02

# Parameters governing the conserved quantities statistics
Statistics:
  delta_time:           1.02
  scale_factor_first:   0.05

# Parameters for the self-gravity scheme
Gravity:
  eta:                    0.025  
  theta:                  0.5
  comoving_softening:     0.01    # 10 kpc = 1/25 mean inter-particle separation
  max_physical_softening: 0.00263 # 10 ckpc = 2.63 pkpc at z=2.8 (EAGLE-like evolution of softening).
  mesh_side_length:       128

# Parameters of the hydro scheme
SPH:
  resolution_eta:      1.2348   # "48 Ngb" with the cubic spline kernel
49
  h_min_ratio:         0.1
50
51
52
53
54
55
56
  CFL_condition:       0.1
  initial_temperature: 1200.    # (1 + z_ini)^2 * 2.72K
  minimal_temperature: 100.

# Parameters related to the initial conditions
InitialConditions:
  file_name:  ./SantaBarbara.hdf5
57
  periodic:   1
58
59
60
  cleanup_h_factors: 1              # ICs were generated for Gadget, we need to get rid of h-factors
  cleanup_velocity_factors: 1       # ICs were generated for Gadget, we need to get rid of sqrt(a) factors in the velocity
  generate_gas_in_ics: 1            # Generate gas particles from the DM-only ICs
61
62
  cleanup_smoothing_lengths: 1      # Since we generate gas, make use of the (expensive) cleaning-up procedure.

63
# Impose primoridal metallicity
64
EAGLEChemistry:
65
66
67
68
69
70
71
72
73
74
  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
75

76
EAGLECooling:
77
  dir_name:                ./coolingtables/
78
  H_reion_z:               11.5 
79
80
81
  He_reion_z_centre:       3.5
  He_reion_z_sigma:        0.5
  He_reion_eV_p_H:         2.0
82