# Define some meta-data about the simulation
MetaData:
  run_name: CosmoBox00

# Define the system of units to use internally. 
InternalUnitSystem:
  UnitMass_in_cgs:     1.98841e43    # 10^10 M_sun in grams
  UnitLength_in_cgs:   3.08567758e24 # Mpc in centimeters
  UnitVelocity_in_cgs: 1e5           # km/s in centimeters per second
  UnitCurrent_in_cgs:  1             # Amperes
  UnitTemp_in_cgs:     1             # Kelvin


# Cosmological parameters
Cosmology:
  h:              0.67321       # Reduced qHubble constant
  a_begin:        0.01282051    # Initial scale-factor of the simulation (z=77)
  a_end:          1.0           # Final scale factor of the simulation
  Omega_cdm:      0.26641139    # Cold Dark Matter density parameter  (value without neutrinos : 0.26499)
  Omega_lambda:   0.6842        # Dark-energy density parameter
  Omega_b:        0.049389      # Baryon density parameter 

# Parameters governing the time integration
TimeIntegration:
  dt_min:     1e-10 # 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).
  
# Parameters governing the snapshots
Snapshots:
  basename:            snap/snapshot # Common part of the name of output files
  output_list_on:      1
  output_list:         ./output_list.txt
  compression:         4

# Parameters governing the conserved quantities statistics
Statistics:
  delta_time:           1.01
  scale_factor_first:   0.02

# Parameters for the self-gravity scheme
Gravity:
  eta:                         0.025     # Constant dimensionless multiplier for time integration.
  MAC:                         adaptive  # Use the geometric opening angle condition
  theta_cr:                    0.7       # Opening angle (Multipole acceptance criterion)
  epsilon_fmm:                 0.001     # Adaptive opening angle
  mesh_side_length:            256
  comoving_DM_softening:       0.0066   # Comoving softening for DM
  max_physical_DM_softening:   0.00174  # Physical softening for DM
  comoving_baryon_softening:       0.0066   # Comoving softening for DM
  max_physical_baryon_softening:   0.00174  # Physical softening for DM
    
    
# 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).
  h_min_ratio:                       0.01     # Minimal smoothing length in units of softening.
  h_max:                             0.5      # Maximal smoothing length in co-moving internal units.
  CFL_condition:                     0.2      # Courant-Friedrich-Levy condition for time integration.
  minimal_temperature:               10.0      # (internal units)

# Cooling with Grackle 3.0
GrackleCooling:
  cloudy_table: CloudyData_UVB=HM2012.h5       # Name of the Cloudy Table (available on the grackle bitbucket repository)
  with_UV_background: 0                        # Enable or not the UV background
  redshift: -1                                 # Redshift to use (-1 means time based redshift)
  with_metal_cooling: 0                        # Enable or not the metal cooling
  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: 0                     # (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
  HydrogenFractionByMass : 0.76                # Hydrogen fraction by mass (default is 0.76)
  use_radiative_transfer : 0                   # Arrays of ionization and heating rates are provided
  H2_three_body_rate : 1                       # (optional) H2 formation three body rate (see grackle documentation)
  H2_cie_cooling : 0                           # Enable/disable H2 collision-induced emission cooling from Ripamonti & Abel (2004)
  cmb_temperature_floor : 1                    # Enable/disable an effective CMB temperature floor
  initial_nHII_to_nH_ratio:   2.393761E-004    # initial nHII   to nH ratio (number density ratio). Value is ignored if set to -1.
  initial_nHeI_to_nH_ratio:   7.900000E-002    # initial nHeI   to nH ratio (number density ratio). Value is ignored if set to -1.  
  initial_nHeII_to_nH_ratio:  5.692753E-020    # initial nHeII  to nH ratio (number density ratio). Value is ignored if set to -1.
  initial_nHeIII_to_nH_ratio: 0                # initial nHeIII to nH ratio (number density ratio). Value is ignored if set to -1.
  initial_nDI_to_nH_ratio:    2.529538E-005    # initial nDI    to nH ratio (number density ratio). Value is ignored if set to -1.
  initial_nDII_to_nH_ratio:   4.284496E-009    # initial nDII   to nH ratio (number density ratio). Value is ignored if set to -1.
  initial_nHM_to_nH_ratio:    6.509693E-012    # initial nHM    to nH ratio (number density ratio). Value is ignored if set to -1.
  initial_nH2I_to_nH_ratio:   1.891128E-006    # initial nH2I   to nH ratio (number density ratio). Value is ignored if set to -1.
  initial_nH2II_to_nH_ratio:  4.153181E-014    # initial nH2II  to nH ratio (number density ratio). Value is ignored if set to -1.
  initial_nHDI_to_nH_ratio:   3.356904E-010    # initial nHDI   to nH ratio (number density ratio). Value is ignored if set to -1.
  maximal_density_Hpcm3: -1                    # Maximal density (in hydrogen atoms/cm^3) for cooling. Higher densities are floored to this value to ensure grackle works properly when interpolating beyond the cloudy_table maximal density. A value < 0 deactivates this parameter.




Scheduler:
  max_top_level_cells:   12
  #cell_split_size:       200
  
Restarts:
  onexit:       1
  delta_hours:  6.0
  max_run_time: 71.5                 # Three days minus fergie time
  #resubmit_on_exit:   1
  #resubmit_command:   ./resub.sh

# Parameters related to the initial conditions
InitialConditions:
  file_name:  homogeneousCube_32.hdf5
  periodic:   1
  cleanup_h_factors: 0               # Remove the h-factors inherited from Gadget
  cleanup_velocity_factors: 0        # Remove the sqrt(a) factor in the velocities inherited from Gadget
  
