# Define the system of units to use internally. 
InternalUnitSystem:
  UnitMass_in_cgs:     1.99e33   # g
  UnitLength_in_cgs:   3.086e21   # cm
  UnitVelocity_in_cgs: 1.e5   # Centimeters per second
  UnitCurrent_in_cgs:  1   # Amperes
  UnitTemp_in_cgs:     1   # Kelvin

# Parameters governing the time integration
TimeIntegration:
  time_begin: 0.    # The starting time of the simulation (in internal units).
  time_end:   1.e-3  # The end time of the simulation (in internal units).
  dt_min:     1.e-13  # The minimal time-step size of the simulation (in internal units).
  dt_max:     1.e-3  # The maximal time-step size of the simulation (in internal units).

# Parameters governing the snapshots
Snapshots:
  basename:            sedov # Common part of the name of output files
  time_first:          0.    # Time of the first output (in internal units)
  delta_time:          1.e-5  # Time difference between consecutive outputs (in internal units)

# Parameters governing the conserved quantities statistics
Statistics:
  delta_time:          1e-5 # Time between statistics output

# 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.
  minimal_temperature:   1.e2
  max_ghost_iterations:  1000
  h_max:                 0.04     # Maximum smoothing length. Since the density in the post-shock region becomes very low, we need to make sure that this does not exceed 1/3 of the periodic box (including kernel gamma ~ 2)

# Parameters related to the initial conditions
InitialConditions:
  file_name:  ./sedov.hdf5          # The file to read
  periodic:   1

EAGLECooling:
  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
  dir_name: ./coolingtables/

EAGLEChemistry:
  init_abundance_metal:        0.014        # Mass fraction in *all* metals
  init_abundance_Hydrogen:     0.70649785   # Mass fraction in Hydrogen
  init_abundance_Helium:       0.28055534   # Mass fraction in Helium
  init_abundance_Carbon:       2.0665436e-3 # Mass fraction in Carbon
  init_abundance_Nitrogen:     8.3562563e-4 # Mass fraction in Nitrogen
  init_abundance_Oxygen:       5.4926244e-3 # Mass fraction in Oxygen
  init_abundance_Neon:         1.4144605e-3 # Mass fraction in Neon
  init_abundance_Magnesium:    5.907064e-4  # Mass fraction in Magnesium
  init_abundance_Silicon:      6.825874e-4  # Mass fraction in Silicon
  init_abundance_Iron:         1.1032152e-3 # Mass fraction in Iron

EoS:
  isothermal_internal_energy: 1240419161676.6465

LambdaTTableCooling:
  file_name: cooling.dat

ConstCooling:
  cooling_rate: 1.e6
  min_energy: 0.01
  cooling_tstep_mult: 2

LambdaCooling:
  lambda_nH2_cgs: 2.e-23
  cooling_tstep_mult: 0.1
  rapid_cooling: 1

Scheduler:
  max_top_level_cells: 12
  cell_max_size: 8000000
  cell_sub_size_pair_hydro: 256000000
  cell_sub_size_self_hydro: 32000
  cell_sub_size_pair_stars: 256000000
  cell_sub_size_self_stars: 32000
  cell_sub_size_pair_grav: 256000000
  cell_sub_size_self_grav: 32000
  cell_split_size: 400
  cell_subdepth_diff_grav: 4
  cell_extra_parts: 0
  cell_extra_sparts: 100
  cell_extra_gparts: 0
  cell_extra_bparts: 0
  cell_extra_sinks: 0
  engine_max_parts_per_ghost: 1000
  engine_max_sparts_per_ghost: 1000
  engine_max_parts_per_cooling: 10000
  nr_queues: 4
  dependency_graph_frequency: 0
  task_level_output_frequency: 0
  tasks_per_cell: 100
  links_per_tasks: 25
  mpi_message_limit: 4