# 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

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

# Parameters governing the snapshots
Snapshots:
  basename:            stellar_evolution # Common part of the name of output files
  time_first:          0.                # Time of the first output (in internal units)
  delta_time:          3.e-5             # Time difference between consecutive outputs (internal units)
  compression:         4
  recording_triggers_part: [1.0227e-4, 1.0227e-5]   # Recording starts 100M and 10M years before a snapshot
  recording_triggers_bpart: [1.0227e-4, 1.0227e-5]  # Recording starts 100M and 10M years before a snapshot

# Parameters governing the conserved quantities statistics
Statistics:
  time_first:          0.
  delta_time:          1.e-5  # non cosmology time between statistics output
  scale_factor_first:  0.5

# Parameters for the hydrodynamics scheme
SPH:
  resolution_eta:        1.2348   # Target smoothing length in units of the mean inter-particle separation 
  CFL_condition:         0.1      # Courant-Friedrich-Levy condition for time integration.
  minimal_temperature:   10.      # Kelvin

# Properties of the stars
Stars:
  overwrite_birth_time: 1
  birth_time:           0.   # Give the star in the ICs a decent birth time
  luminosity_filename:   ./photometry

# Parameters related to the initial conditions
InitialConditions:
  file_name:  ./stellar_evolution.hdf5       # The file to read
  periodic:   1
  
Scheduler:
  max_top_level_cells: 8
  tasks_per_cell:      200

# Parameters for the EAGLE "equation of state"
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

# Metallicites read in for the gas and star
EAGLEChemistry:              
  init_abundance_metal:      0.01
  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

# Standard EAGLE cooling options
EAGLECooling:
  dir_name:                ./coolingtables/  # Location of the Wiersma+08 cooling tables
  H_reion_z:               11.5              # Redshift of Hydrogen re-ionization
  H_reion_eV_p_H:          2.0
  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
  He_reion_eV_p_H:         2.0               # Energy inject by Helium re-ionization in electron-volt per Hydrogen atom

# PS2020 cooling parameters
PS2020Cooling:
  dir_name:                ./UV_dust1_CR1_G1_shield1.hdf5 # Location of the cooling tables
  H_reion_z:               7.5               # Redshift of Hydrogen re-ionization (Planck 2018)
  H_reion_eV_p_H:          2.0
  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
  He_reion_eV_p_H:         2.0               # Energy inject by Helium re-ionization in electron-volt per Hydrogen atom
  delta_logTEOS_subgrid_properties: 0.3      # delta log T above the EOS below which the subgrid properties use Teq assumption
  rapid_cooling_threshold:          0.333333 # Switch to rapid cooling regime for dt / t_cool above this threshold.

# Standard COLIBRE cooling options
COLIBRECooling:
  dir_name:                ./cooling_files/ # Location of the cooling tables
  filebase_cool:           cooling_rates
  H_reion_z:               7.5              # Redshift of Hydrogen re-ionization
  H_reion_eV_p_H:          2.0
  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
  He_reion_eV_p_H:         2.0               # Energy inject by Helium re-ionization in electron-volt per Hydrogen atom
  delta_logTEOS_subgrid_properties: 0.3      # delta log T above the EOS below which the subgrid properties use Teq assumption
  rapid_cooling_threshold:          0.333333 # Switch to rapid cooling regime for dt / t_cool above this threshold.

# CHIMES cooling parameters
CHIMESCooling: 
  data_path:                  ./chimes-data     # Path to chimes-data repository
  EqmAbundanceTable:          colibre_HHe              # Path to equilibrium abundance tables
  PhotoIonTable_UVB:          SP20_cross_sections      # Path to cross-sections tables; UVB 
  PhotoIonTable_ISRF:         cross_sections_B87.hdf5  # Path to cross-sections tables; ISRF
  UV_field:             COLIBRE           # Flag specifying radiation field model 
  Shielding_model:      COLIBRE           # Flag specifying local shielding model 
  UVB_z_dependence:     COLIBRE           # Flag to switch on redshift-dependent UVB 
  shielding_length_factor:    0.5         # Multiplicative factor to scale shielding length
  max_shielding_length_kpc:   100.0       # Max shielding length in code units
  rad_field_norm_factor:      0.01        # Multiplicative factor to scale ISRF
  init_abundance_mode:     read           # Flag to specify how to set initial CHIMES abundances
  colibre_metal_depletion:    0           # Flag to specify metal depletion model
  turbulent_velocity_dispersion_km_p_s:  6.0 # 1D velocity dispersion, in km/s.
  colibre_use_turbulent_jeans_length:  1  # Integer flag, 0 - use only thermal Jeans length, 1 - use maximum of thermal or turbulent Jeans length.
  colibre_saturate_radiation_field:    1  # Integer flag, 0 - no max rad field, 1 - impose max rad field.
  colibre_col_dens_saturate_rad:  1.3e22  # Column density at which the Colibre rad field saturates, in cm^-2. 
  colibre_saturate_cr_rate:   1           # Integer flag, 0 - no max CR rate, 1 - impose max CR rate.
  colibre_col_dens_saturate_cr:  3.63e20  # Column density at which the CR rate saturates, in cm^-2. 
  colibre_cr_plaw_index:      1.0         # Power law index of the scaling between CR rate and column density. 
  relativeTolerance:          1e-3        # Relative tolerance for CHIMES intergration
  absoluteTolerance:          1e-10       # Absolute tolerance for CHIMES chemistry integration
  explicitTolerance:          0.1         # Tolerance below which we use the explicit solution in CHIMES
  scale_metal_tolerances:     1           # Flag to scale absolute tolerances by element abundance
  Tmol_K:                     1.0e5       # Maximum temperature for the molecular network
  ChemistryEqmMode:           0           # Flag to use pre-computed equilibrium abundances throughout
  ThermEvolOn:                1           # Flag to switch on thermal evolution in CHIMES
  chimes_debug:               0           # Flag to include extra debug outputs in CHIMES
  cosmic_ray_rate_cgs:        1.8e-16     # HI ionisation rate from cosmic rays in units of HI per second.
  max_dust_boost_factor:      100.0       # Max boost factor applied to all reactions on dust grains
  dust_boost_nH_min_cgs:      1.0         # Density below which the dust boost factor is equal to unity, in units of cm^-3
  dust_boost_nH_max_cgs:      100.0       # Density below which the dust boost factor is equal to unity, in units of cm^-3
  delta_logTEOS_subgrid_properties:  0.3  # delta logT above the EOS below which set CHIMES to eqm and use subgrid properties
  use_colibre_subgrid_EOS:    1           # Flag to set subgrid rho and T to thermal and pressure eqm on the EOS
  set_FB_particles_to_eqm:    1           # Flag to set particles to eqm if they have been heated by feedback 
  use_hybrid_cooling:         1           # Use hybrid cooling, to read in eqm cooling rates from metals not included in CHIMES
  rapid_cooling_threshold:    0.333333    # Threshold in dt / t_cool to switch between rapid- and slow-cooling modes
  colibre_table_path:         ./cooling_files  # Path to COLIBRE cooling tables (Ploeckinger & Schaye 2020)
  colibre_table_filebase:     UVB_dust1_CR1_G1_shield1_v1  # Filebase of COLIBRE cooling tables split in redshift (Ploeckinger & Schaye 2020)
  UVB_cutoff_z:               127.0            # Redshift above which the UVB is switched off 
  ISRF_low_dens_cutoff_z:     7.5              # Redshift above which the ISRF is cut off at low densities (UV_field_flag == 2 only)
  S_over_Si_in_solar:         1.0         # S / Si relative to Solar 
  Ca_over_Si_in_solar:        1.0         # Ca / Si relative to Solar 
  IncludeCarbon:              0           # Include Carbon in the CHIMES network
  IncludeNitrogen:            0           # Include Nitrogen in the CHIMES network
  IncludeOxygen:              0           # Include Oxygen in the CHIMES network
  IncludeNeon:                0           # Include Neon in the CHIMES network
  IncludeMagnesium:           0           # Include Magnesium in the CHIMES network
  IncludeSilicon:             0           # Include Silicon in the CHIMES network
  IncludeSulphur:             0           # Include Sulphur in the CHIMES network
  IncludeCalcium:             0           # Include Calcium in the CHIMES network
  IncludeIron:                0           # Include Iron in the CHIMES network

# Properties of the EAGLE feedback and enrichment model. This matches the EAGLE-Ref (Schaye+2015) model.
EAGLEFeedback:
  use_SNII_feedback:                    0               # Global switch for SNII thermal (stochastic) feedback.
  use_SNIa_feedback:                    0               # Global switch for SNIa thermal (continuous) feedback.
  use_AGB_enrichment:                   1               # Global switch for enrichement from AGB stars.
  use_SNII_enrichment:                  1               # Global switch for enrichement from SNII stars.
  use_SNIa_enrichment:                  1               # Global switch for enrichement from SNIa stars.
  filename:                             ./yieldtables/  # Path to the directory containing the EAGLE yield tables.
  IMF_min_mass_Msun:                    0.1             # Minimal stellar mass considered for the Chabrier IMF in solar masses.
  IMF_max_mass_Msun:                  100.0             # Maximal stellar mass considered for the Chabrier IMF in solar masses.
  SNII_min_mass_Msun:                   6.0             # Minimal mass considered for SNII stars in solar masses.
  SNII_max_mass_Msun:                 100.0             # Maximal mass considered for SNII stars in solar masses.
  SNII_feedback_model:                  Random          # Feedback modes: Random, Isotropic, MinimumDistance, MinimumDensity
  SNII_sampled_delay:                   0               # Sample the SNII lifetimes to do feedback.
  SNII_wind_delay_Gyr:                  0.03            # Time in Gyr between a star's birth and the SNII thermal feedback event when not sampling.
  SNII_delta_T_K:                       3.16228e7       # Change in temperature to apply to the gas particle in a SNII thermal feedback event in Kelvin.
  SNII_energy_erg:                      1.0e51          # Energy of one SNII explosion in ergs.
  SNII_energy_fraction_function:        EAGLE
  SNII_energy_fraction_min:             0.3             # Minimal fraction of energy applied in a SNII feedback event.
  SNII_energy_fraction_max:             3.0             # Maximal fraction of energy applied in a SNII feedback event.
  SNII_energy_fraction_Z_0:             0.0012663729    # Pivot point for the metallicity dependance of the SNII energy fraction (metal mass fraction).
  SNII_energy_fraction_n_0_H_p_cm3:     0.67            # Pivot point for the birth density dependance of the SNII energy fraction in cm^-3.
  SNII_energy_fraction_n_Z:             0.8686          # Power-law for the metallicity dependance of the SNII energy fraction.
  SNII_energy_fraction_n_n:             0.8686          # Power-law for the birth density dependance of the SNII energy fraction.
  SNII_energy_fraction_use_birth_density: 0             # Are we using the density and metallicity at birth to compute f_E or at feedback time?
  SNII_energy_fraction_use_birth_metallicity: 0         # Are we using the density and metallicity at birth to compute f_E or at feedback time?
  SNIa_DTD:                             Exponential     # Use the EAGLE-Ref SNIa DTD.
  SNIa_DTD_delay_Gyr:                   0.04            # Age of the most massive SNIa in Gyr.
  SNIa_DTD_exp_timescale_Gyr:           2.0             # Time-scale of the exponential decay of the SNIa rates in Gyr.
  SNIa_DTD_exp_norm_p_Msun:             0.002           # Normalisation of the SNIa rates in inverse solar masses.
  SNIa_energy_erg:                     1.0e51           # Energy of one SNIa explosion in ergs.
  AGB_ejecta_velocity_km_p_s:          10.0             # Velocity of the AGB ejectas in km/s.
  stellar_evolution_age_cut_Gyr:       0.1              # Stellar age in Gyr above which the enrichment is down-sampled.
  stellar_evolution_sampling_age_fraction: 0.05         # Time-step of the stars in units of the star's current age when the star is older than the age cut.
  SNII_yield_factor_Hydrogen:           1.0             # (Optional) Correction factor to apply to the Hydrogen yield from the SNII channel.
  SNII_yield_factor_Helium:             1.0             # (Optional) Correction factor to apply to the Helium yield from the SNII channel.
  SNII_yield_factor_Carbon:             0.5             # (Optional) Correction factor to apply to the Carbon yield from the SNII channel.
  SNII_yield_factor_Nitrogen:           1.0             # (Optional) Correction factor to apply to the Nitrogen yield from the SNII channel.
  SNII_yield_factor_Oxygen:             1.0             # (Optional) Correction factor to apply to the Oxygen yield from the SNII channel.
  SNII_yield_factor_Neon:               1.0             # (Optional) Correction factor to apply to the Neon yield from the SNII channel.
  SNII_yield_factor_Magnesium:          2.0             # (Optional) Correction factor to apply to the Magnesium yield from the SNII channel.
  SNII_yield_factor_Silicon:            1.0             # (Optional) Correction factor to apply to the Silicon yield from the SNII channel.
  SNII_yield_factor_Iron:               0.5             # (Optional) Correction factor to apply to the Iron yield from the SNII channel.

# Solar abundances
COLIBREChemistry:
  init_abundance_metal:      0.01           # Inital fraction of particle mass in *all* metals 
  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
  init_abundance_Europium:   0.0
  init_abundance_Strontium: 0.0    # Initial fraction of particle mass in Strontium
  init_abundance_Barium:    0.0    # Inital fraction of particle mass in Barium
  metal_diffusion_constant:         0.   # No diffusion
  metal_diffusion_timestep_mult:    100. # No diffusion

# COLIBRE feedback model
COLIBREFeedback:
  use_SNII_feedback:                0               # Global switch for SNII thermal (stochastic) feedback.
  use_SNIa_feedback:                1               # Global switch for SNIa thermal (continuous) feedback.
  use_AGB_enrichment:               1               # Global switch for enrichment from AGB stars.
  use_SNII_enrichment:              1               # Global switch for enrichment from SNII stars.
  use_SNIa_enrichment:              1               # Global switch for enrichment from SNIa stars.
  with_r_process_enrichment:        1               # Global switch for r process enrichment from neutron star mergers and rare core-collapse SN (collapsars and common envelop jets SN).
  with_HIIRegions:                  0               # Switch for early feedback (HII regions).
  with_StellarWinds:                0               # Switch for early feedback (Stellar winds).
  with_RadiationPressure:           0               # Switch for early feedback (Radiation pressure).
  filename:                         ./yieldtables/  # Path to the directory containing the EAGLE yield tables.
  IMF_min_mass_Msun:                0.1             # Minimal stellar mass considered for the Chabrier IMF in solar masses.
  IMF_max_mass_Msun:              100.0             # Maximal stellar mass considered for the Chabrier IMF in solar masses.
  SNII_min_mass_Msun:               8.0             # Minimal mass considered for SNII feedback (not SNII enrichment!) in solar masses.
  SNII_max_mass_Msun:             100.0             # Maximal mass considered for SNII feedback (not SNII enrichment!) in solar masses.
  SNII_wind_delay_Gyr:              0.03            # Time in Gyr between a star's birth and the SNII thermal feedback event.
  SNII_delay_thermal_feedback_Gyr:  0.02            # Time in Gyr between a star's birth and the earliest possible SNII thermal feedback event. If SNe occur at earlier times, these events are postponed.
  SNII_delta_T_K:                   3.16228e7       # Change in temperature to apply to the gas particle in a SNII thermal feedback event in Kelvin.
  SNII_energy_erg:                  1.0e51          # Energy of one SNII explosion in ergs.
  SNII_delta_v_km_p_s:              25.0            # Change in the gas-particle velocity that is kicked in a SNII kinetic feedback event
  SNII_f_kinetic:                   0.1             # Fraction of SNII energy injected into the gas in kinetic form (the remaining fraction is injected in thermal form).
  SNII_energy_fraction_min:         0.3             # Maximal fraction of energy applied in a SNII feedback event.
  SNII_energy_fraction_max:         3.0             # Minimal fraction of energy applied in a SNII feedback event.
  SNII_energy_fraction_Z_0:         0.0012663729    # Pivot point for the metallicity dependance of the SNII energy fraction (metal mass fraction).
  SNII_energy_fraction_n_0_H_p_cm3: 0.67            # Pivot point for the birth density dependance of the SNII energy fraction in cm^-3.
  SNII_energy_fraction_n_Z:         0.8686          # Power-law for the metallicity dependance of the SNII energy fraction.
  SNII_energy_fraction_n_n:         0.8686          # Power-law for the birth density dependance of the SNII energy fraction.
  SNII_delta_T_n_0_H_p_cm3:         0.3             # Pivot (physical) density in SNII variable-dT thermal feedback in H/cm3.  
  SNII_delta_T_K_pivot:             3.16228e6       # Pivot temperature in SNII variable-dT thermal feedback in Kelvin. 
  SNII_delta_T_slope:               0.666666666     # Slope in SNII variable-dT thermal feedback.
  SNII_delta_T_K_min:               3.16228e6       # Minimal change in temperature to apply to the gas particle in an SNII thermal feedback event in Kelvin.
  SNII_delta_T_K_max:               1e8             # Maximal change in temperature to apply to the gas particle in na SNII thermal feedback event in Kelvin.
  SNII_energy_fraction_P_0_K_p_cm3: 1.5e4
  SNII_energy_fraction_sigma_P:     0.3
  SNII_enrichment_min_mass_Msun:    8.0
  SNII_enrichment_max_mass_Msun:    40.0
  SNIa_energy_erg:                 1.0e51           # Energy of one SNIa explosion in ergs.
  SNIa_delta_T_K:                   3.16228e7       # Change in temperature to apply to the gas particle in a SNIa thermal feedback event in Kelvin.
  SNIa_energy_fraction:             1.0             # Fraction of energy applied in a SNIa feedback event.
  SNIa_delta_T_n_0_H_p_cm3:         0.3             # Pivot (physical) density in SNIa variable-dT thermal feedback in H/cm3.
  SNIa_delta_T_K_pivot:             3.16228e6       # Pivot temperature in SNIa variable-dT thermal feedback in Kelvin. 
  SNIa_delta_T_slope:               0.666666666     # Slope in SNIa variable-dT thermal feedback.
  SNIa_delta_T_K_min:               3.16228e6       # Minimal change in temperature to apply to the gas particle in an SNIa thermal feedback event in Kelvin.
  SNIa_delta_T_K_max:               1e8             # Maximal change in temperature to apply to the gas particle in na SNIa thermal feedback event in Kelvin.
  AGB_ejecta_velocity_km_p_s:      10.0             # Velocity of the AGB ejectas in km/s.
  stellar_evolution_age_cut_Gyr:    1.0             # Stellar age in Gyr above which the enrichment is down-sampled.
  stellar_evolution_sampling_age_fraction: 0.1      # Fraction of the particle's age that can be used as a maximal time-step size
  SNII_yield_factor_Hydrogen:       1.0             # (Optional) Correction factor to apply to the Hydrogen yield from the SNII channel.
  SNII_yield_factor_Helium:         1.0             # (Optional) Correction factor to apply to the Helium yield from the SNII channel.
  SNII_yield_factor_Carbon:         3.0             # (Optional) Correction factor to apply to the Carbon yield from the SNII channel.
  SNII_yield_factor_Nitrogen:       0.8             # (Optional) Correction factor to apply to the Nitrogen yield from the SNII channel.
  SNII_yield_factor_Oxygen:         1.0             # (Optional) Correction factor to apply to the Oxygen yield from the SNII channel.
  SNII_yield_factor_Neon:           1.0             # (Optional) Correction factor to apply to the Neon yield from the SNII channel.
  SNII_yield_factor_Magnesium:      1.0             # (Optional) Correction factor to apply to the Magnesium yield from the SNII channel.
  SNII_yield_factor_Silicon:        0.8             # (Optional) Correction factor to apply to the Silicon yield from the SNII channel.
  SNII_yield_factor_Iron:           2.0             # (Optional) Correction factor to apply to the Iron yield from the SNII channel.
  num_of_NSM_per_Msun:              1.3e-5          # Number of neutron star mergers per Msolar.
  yield_Eu_from_NSM_event_Msun:     1e-5            # Amount of europium (in units of Msolar) released by NSM.
  num_of_CEJSN_per_Msun:            1.22e-5         # Number of rare core-collapse SN (CEJSN) per Msolar.
  yield_Eu_from_CEJSN_event_Msun:   2e-5            # Amount of europium (in units of Msolar) released by CEJSN(=Common envelop jet SN).
  num_of_collapsar_per_Msun:        2.6e-5          # Number of rare core-collapse SN (collapsars) per Msolar.
  yield_Eu_from_collapsar_event_Msun: 1e-6          # Amount of europium (in units of Msolar) released by collapsars.
  collapsar_min_mass_Msun:          10.0            # Minimal mass considered for r-process enrichment from collapsars in solar masses.
  collapsar_max_mass_Msun:          40.0            # Maximal mass considered for r-process enrichment from collapsars in solar masses.
  HIIregion_ionization_fraction:    1.0             # Ionization fraction for gas particles tagged as HII regions (between 0.5 and 1.)
  HIIregion_temperature_K:          1.e4            # Temperature of gas particles tagged as HII regions

# COLIBRE feedback SNIa delay time distribution struct
SNIaDTD:
  SNIa_efficiency_p_Msun:       0.002       # Normalisation of the SNIa rates in inverse solar masses, used when there is only one functional form in the DTD.
  normalization_timescale_Gyr:  13.6        # Normalization time for DTDs that cannot be integrated to infinity 
  SNIa_delay_time_Gyr:          0.04        # Delay time before the DTD kicks in.
  SNIa_timescale_Gyr:           2.0

Event_logger:
  delta_time_SNIa_Myr:        10
  delta_time_SNII_Myr:        1
  delta_time_r_processes_Myr: 10

DustEvolution:
  pair_to_cooling: 0

XrayEmissivity:
  xray_table_path:    ./X_Ray_tables.hdf5   # Path to the X-ray emissivity tables
