Commit 4b7ef88f authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Merge branch 'Ploeckinger_cooling' into 'master'

Add Ploeckinger+20 cooling

See merge request !1113
parents a214d648 5f2ab077
...@@ -1534,6 +1534,19 @@ case "$with_subgrid" in ...@@ -1534,6 +1534,19 @@ case "$with_subgrid" in
with_subgrid_sink=none with_subgrid_sink=none
enable_fof=yes enable_fof=yes
;; ;;
EAGLE-XL)
with_subgrid_cooling=COLIBRE
with_subgrid_chemistry=EAGLE
with_subgrid_tracers=EAGLE
with_subgrid_entropy_floor=EAGLE
with_subgrid_stars=EAGLE
with_subgrid_star_formation=EAGLE
with_subgrid_feedback=EAGLE
with_subgrid_black_holes=EAGLE
with_subgrid_task_order=EAGLE
with_subgrid_sink=none
enable_fof=yes
;;
*) *)
AC_MSG_ERROR([Unknown subgrid choice: $with_subgrid]) AC_MSG_ERROR([Unknown subgrid choice: $with_subgrid])
;; ;;
...@@ -1860,6 +1873,9 @@ case "$with_cooling" in ...@@ -1860,6 +1873,9 @@ case "$with_cooling" in
EAGLE) EAGLE)
AC_DEFINE([COOLING_EAGLE], [1], [Cooling following the EAGLE model]) AC_DEFINE([COOLING_EAGLE], [1], [Cooling following the EAGLE model])
;; ;;
COLIBRE)
AC_DEFINE([COOLING_COLIBRE], [1], [Cooling following the COLIBRE model])
;;
*) *)
AC_MSG_ERROR([Unknown cooling function: $with_cooling]) AC_MSG_ERROR([Unknown cooling function: $with_cooling])
;; ;;
...@@ -2269,6 +2285,9 @@ AM_CONDITIONAL([HAVEQLACOOLING], [test $with_cooling = "QLA"]) ...@@ -2269,6 +2285,9 @@ AM_CONDITIONAL([HAVEQLACOOLING], [test $with_cooling = "QLA"])
# Check if using EAGLE cooling # Check if using EAGLE cooling
AM_CONDITIONAL([HAVEEAGLECOOLING], [test $with_cooling = "EAGLE"]) AM_CONDITIONAL([HAVEEAGLECOOLING], [test $with_cooling = "EAGLE"])
# Check if using COLIBRE cooling
AM_CONDITIONAL([HAVECOLIBRECOOLING], [test $with_cooling = "COLIBRE"])
# Check if using EAGLE feedback # Check if using EAGLE feedback
AM_CONDITIONAL([HAVEEAGLEFEEDBACK], [test $with_feedback = "EAGLE"]) AM_CONDITIONAL([HAVEEAGLEFEEDBACK], [test $with_feedback = "EAGLE"])
......
#!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/CoolingTables/COLIBRE/UV_dust1_CR1_G1_shield1.hdf5
...@@ -123,6 +123,17 @@ EAGLECooling: ...@@ -123,6 +123,17 @@ EAGLECooling:
He_reion_z_sigma: 0.5 He_reion_z_sigma: 0.5
He_reion_eV_p_H: 2.0 He_reion_eV_p_H: 2.0
# COLIBRE cooling parameters
COLIBRECooling:
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.
# EAGLE star formation parameters # EAGLE star formation parameters
EAGLEStarFormation: EAGLEStarFormation:
EOS_density_norm_H_p_cm3: 0.1 # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3. EOS_density_norm_H_p_cm3: 0.1 # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3.
...@@ -133,7 +144,7 @@ EAGLEStarFormation: ...@@ -133,7 +144,7 @@ EAGLEStarFormation:
min_over_density: 57.7 # The over-density above which star-formation is allowed. min_over_density: 57.7 # The over-density above which star-formation is allowed.
KS_high_density_threshold_H_p_cm3: 1e3 # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3. KS_high_density_threshold_H_p_cm3: 1e3 # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3.
KS_high_density_exponent: 2.0 # Slope of the Kennicut-Schmidt law above the high-density threshold. KS_high_density_exponent: 2.0 # Slope of the Kennicut-Schmidt law above the high-density threshold.
EOS_entropy_margin_dex: 0.5 # Logarithm base 10 of the maximal entropy above the EOS at which stars can form. EOS_entropy_margin_dex: 0.3 # Logarithm base 10 of the maximal entropy above the EOS at which stars can form.
threshold_norm_H_p_cm3: 0.1 # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3. threshold_norm_H_p_cm3: 0.1 # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3.
threshold_Z0: 0.002 # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation. threshold_Z0: 0.002 # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation.
threshold_slope: -0.64 # Slope of the metal-dependant star formation threshold threshold_slope: -0.64 # Slope of the metal-dependant star formation threshold
...@@ -141,13 +152,13 @@ EAGLEStarFormation: ...@@ -141,13 +152,13 @@ EAGLEStarFormation:
# Parameters for the EAGLE "equation of state" # Parameters for the EAGLE "equation of state"
EAGLEEntropyFloor: 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_density_threshold_H_p_cm3: 1e-4 # 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_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_temperature_norm_K: 800 # 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 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_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_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_temperature_norm_K: 10. # Temperature of the EAGLE Cool limiter entropy floor at the density threshold expressed in Kelvin. (NOTE: This is below the min T and hence this floor does nothing)
Cool_gamma_effective: 1. # Slope the of the EAGLE Cool limiter entropy floor Cool_gamma_effective: 1. # Slope the of the EAGLE Cool limiter entropy floor
# EAGLE feedback model # EAGLE feedback model
......
...@@ -166,6 +166,13 @@ Iterate_cm_flag=0 #do not interate to determine centre-of-mass ...@@ -166,6 +166,13 @@ Iterate_cm_flag=0 #do not interate to determine centre-of-mass
Sort_by_binding_energy=1 #sort particles by binding energy Sort_by_binding_energy=1 #sort particles by binding energy
Reference_frame_for_properties=2 #use the minimum potential as reference frame about which to calculate properties Reference_frame_for_properties=2 #use the minimum potential as reference frame about which to calculate properties
################################
#Extra fields
################################
Gas_internal_property_names=ElementMassFractions,SpeciesFractions,SpeciesFractions,SpeciesFractions,
Gas_internal_property_index_in_file=0,0,1,2,
Gas_internal_property_calculation_type =averagemassweighted,averagemassweighted,averagemassweighted,averagemassweighted,
################################ ################################
#output related #output related
################################ ################################
......
...@@ -123,6 +123,17 @@ EAGLECooling: ...@@ -123,6 +123,17 @@ EAGLECooling:
He_reion_z_sigma: 0.5 He_reion_z_sigma: 0.5
He_reion_eV_p_H: 2.0 He_reion_eV_p_H: 2.0
# COLIBRE cooling parameters
COLIBRECooling:
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.
# EAGLE star formation parameters # EAGLE star formation parameters
EAGLEStarFormation: EAGLEStarFormation:
EOS_density_norm_H_p_cm3: 0.1 # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3. EOS_density_norm_H_p_cm3: 0.1 # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3.
...@@ -133,7 +144,7 @@ EAGLEStarFormation: ...@@ -133,7 +144,7 @@ EAGLEStarFormation:
min_over_density: 57.7 # The over-density above which star-formation is allowed. min_over_density: 57.7 # The over-density above which star-formation is allowed.
KS_high_density_threshold_H_p_cm3: 1e3 # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3. KS_high_density_threshold_H_p_cm3: 1e3 # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3.
KS_high_density_exponent: 2.0 # Slope of the Kennicut-Schmidt law above the high-density threshold. KS_high_density_exponent: 2.0 # Slope of the Kennicut-Schmidt law above the high-density threshold.
EOS_entropy_margin_dex: 0.5 # Logarithm base 10 of the maximal entropy above the EOS at which stars can form. EOS_entropy_margin_dex: 0.3 # Logarithm base 10 of the maximal entropy above the EOS at which stars can form.
threshold_norm_H_p_cm3: 0.1 # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3. threshold_norm_H_p_cm3: 0.1 # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3.
threshold_Z0: 0.002 # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation. threshold_Z0: 0.002 # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation.
threshold_slope: -0.64 # Slope of the metal-dependant star formation threshold threshold_slope: -0.64 # Slope of the metal-dependant star formation threshold
...@@ -141,13 +152,13 @@ EAGLEStarFormation: ...@@ -141,13 +152,13 @@ EAGLEStarFormation:
# Parameters for the EAGLE "equation of state" # Parameters for the EAGLE "equation of state"
EAGLEEntropyFloor: 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_density_threshold_H_p_cm3: 1e-4 # 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_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_temperature_norm_K: 800 # 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 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_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_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_temperature_norm_K: 10. # Temperature of the EAGLE Cool limiter entropy floor at the density threshold expressed in Kelvin. (NOTE: This is below the min T and hence this floor does nothing)
Cool_gamma_effective: 1. # Slope the of the EAGLE Cool limiter entropy floor Cool_gamma_effective: 1. # Slope the of the EAGLE Cool limiter entropy floor
# EAGLE feedback model # EAGLE feedback model
......
...@@ -166,6 +166,13 @@ Iterate_cm_flag=0 #do not interate to determine centre-of-mass ...@@ -166,6 +166,13 @@ Iterate_cm_flag=0 #do not interate to determine centre-of-mass
Sort_by_binding_energy=1 #sort particles by binding energy Sort_by_binding_energy=1 #sort particles by binding energy
Reference_frame_for_properties=2 #use the minimum potential as reference frame about which to calculate properties Reference_frame_for_properties=2 #use the minimum potential as reference frame about which to calculate properties
################################
#Extra fields
################################
Gas_internal_property_names=ElementMassFractions,SpeciesFractions,SpeciesFractions,SpeciesFractions,
Gas_internal_property_index_in_file=0,0,1,2,
Gas_internal_property_calculation_type =averagemassweighted,averagemassweighted,averagemassweighted,averagemassweighted,
################################ ################################
#output related #output related
################################ ################################
......
...@@ -121,6 +121,17 @@ EAGLECooling: ...@@ -121,6 +121,17 @@ EAGLECooling:
He_reion_z_sigma: 0.5 He_reion_z_sigma: 0.5
He_reion_eV_p_H: 2.0 He_reion_eV_p_H: 2.0
# COLIBRE cooling parameters
COLIBRECooling:
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.
# EAGLE star formation parameters # EAGLE star formation parameters
EAGLEStarFormation: EAGLEStarFormation:
EOS_density_norm_H_p_cm3: 0.1 # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3. EOS_density_norm_H_p_cm3: 0.1 # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3.
...@@ -131,7 +142,7 @@ EAGLEStarFormation: ...@@ -131,7 +142,7 @@ EAGLEStarFormation:
min_over_density: 57.7 # The over-density above which star-formation is allowed. min_over_density: 57.7 # The over-density above which star-formation is allowed.
KS_high_density_threshold_H_p_cm3: 1e3 # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3. KS_high_density_threshold_H_p_cm3: 1e3 # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3.
KS_high_density_exponent: 2.0 # Slope of the Kennicut-Schmidt law above the high-density threshold. KS_high_density_exponent: 2.0 # Slope of the Kennicut-Schmidt law above the high-density threshold.
EOS_entropy_margin_dex: 0.5 # Logarithm base 10 of the maximal entropy above the EOS at which stars can form. EOS_entropy_margin_dex: 0.3 # Logarithm base 10 of the maximal entropy above the EOS at which stars can form.
threshold_norm_H_p_cm3: 0.1 # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3. threshold_norm_H_p_cm3: 0.1 # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3.
threshold_Z0: 0.002 # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation. threshold_Z0: 0.002 # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation.
threshold_slope: -0.64 # Slope of the metal-dependant star formation threshold threshold_slope: -0.64 # Slope of the metal-dependant star formation threshold
...@@ -139,13 +150,13 @@ EAGLEStarFormation: ...@@ -139,13 +150,13 @@ EAGLEStarFormation:
# Parameters for the EAGLE "equation of state" # Parameters for the EAGLE "equation of state"
EAGLEEntropyFloor: 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_density_threshold_H_p_cm3: 1e-4 # 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_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_temperature_norm_K: 800 # 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 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_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_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_temperature_norm_K: 10. # Temperature of the EAGLE Cool limiter entropy floor at the density threshold expressed in Kelvin. (NOTE: This is below the min T and hence this floor does nothing)
Cool_gamma_effective: 1. # Slope the of the EAGLE Cool limiter entropy floor Cool_gamma_effective: 1. # Slope the of the EAGLE Cool limiter entropy floor
# EAGLE feedback model # EAGLE feedback model
......
...@@ -166,6 +166,13 @@ Iterate_cm_flag=0 #do not interate to determine centre-of-mass ...@@ -166,6 +166,13 @@ Iterate_cm_flag=0 #do not interate to determine centre-of-mass
Sort_by_binding_energy=1 #sort particles by binding energy Sort_by_binding_energy=1 #sort particles by binding energy
Reference_frame_for_properties=2 #use the minimum potential as reference frame about which to calculate properties Reference_frame_for_properties=2 #use the minimum potential as reference frame about which to calculate properties
################################
#Extra fields
################################
Gas_internal_property_names=ElementMassFractions,SpeciesFractions,SpeciesFractions,SpeciesFractions,
Gas_internal_property_index_in_file=0,0,1,2,
Gas_internal_property_calculation_type =averagemassweighted,averagemassweighted,averagemassweighted,averagemassweighted,
################################ ################################
#output related #output related
################################ ################################
......
...@@ -120,6 +120,17 @@ EAGLECooling: ...@@ -120,6 +120,17 @@ EAGLECooling:
He_reion_z_sigma: 0.5 He_reion_z_sigma: 0.5
He_reion_eV_p_H: 2.0 He_reion_eV_p_H: 2.0
# COLIBRE cooling parameters
COLIBRECooling:
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.
# EAGLE star formation parameters # EAGLE star formation parameters
EAGLEStarFormation: EAGLEStarFormation:
EOS_density_norm_H_p_cm3: 0.1 # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3. EOS_density_norm_H_p_cm3: 0.1 # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3.
...@@ -130,7 +141,7 @@ EAGLEStarFormation: ...@@ -130,7 +141,7 @@ EAGLEStarFormation:
min_over_density: 57.7 # The over-density above which star-formation is allowed. min_over_density: 57.7 # The over-density above which star-formation is allowed.
KS_high_density_threshold_H_p_cm3: 1e3 # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3. KS_high_density_threshold_H_p_cm3: 1e3 # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3.
KS_high_density_exponent: 2.0 # Slope of the Kennicut-Schmidt law above the high-density threshold. KS_high_density_exponent: 2.0 # Slope of the Kennicut-Schmidt law above the high-density threshold.
EOS_entropy_margin_dex: 0.5 # Logarithm base 10 of the maximal entropy above the EOS at which stars can form. EOS_entropy_margin_dex: 0.3 # Logarithm base 10 of the maximal entropy above the EOS at which stars can form.
threshold_norm_H_p_cm3: 0.1 # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3. threshold_norm_H_p_cm3: 0.1 # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3.
threshold_Z0: 0.002 # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation. threshold_Z0: 0.002 # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation.
threshold_slope: -0.64 # Slope of the metal-dependant star formation threshold threshold_slope: -0.64 # Slope of the metal-dependant star formation threshold
...@@ -138,13 +149,13 @@ EAGLEStarFormation: ...@@ -138,13 +149,13 @@ EAGLEStarFormation:
# Parameters for the EAGLE "equation of state" # Parameters for the EAGLE "equation of state"
EAGLEEntropyFloor: 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_density_threshold_H_p_cm3: 1e-4 # 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_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_temperature_norm_K: 800 # 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 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_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_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_temperature_norm_K: 10. # Temperature of the EAGLE Cool limiter entropy floor at the density threshold expressed in Kelvin. (NOTE: This is below the min T and hence this floor does nothing)
Cool_gamma_effective: 1. # Slope the of the EAGLE Cool limiter entropy floor Cool_gamma_effective: 1. # Slope the of the EAGLE Cool limiter entropy floor
# EAGLE feedback model # EAGLE feedback model
......
...@@ -166,6 +166,13 @@ Iterate_cm_flag=0 #do not interate to determine centre-of-mass ...@@ -166,6 +166,13 @@ Iterate_cm_flag=0 #do not interate to determine centre-of-mass
Sort_by_binding_energy=1 #sort particles by binding energy Sort_by_binding_energy=1 #sort particles by binding energy
Reference_frame_for_properties=2 #use the minimum potential as reference frame about which to calculate properties Reference_frame_for_properties=2 #use the minimum potential as reference frame about which to calculate properties
################################
#Extra fields
################################
Gas_internal_property_names=ElementMassFractions,SpeciesFractions,SpeciesFractions,SpeciesFractions,
Gas_internal_property_index_in_file=0,0,1,2,
Gas_internal_property_calculation_type =averagemassweighted,averagemassweighted,averagemassweighted,averagemassweighted,
################################ ################################
#output related #output related
################################ ################################
......
...@@ -3,6 +3,10 @@ the EAGLE suite of simulations. The cosmology, resolution ...@@ -3,6 +3,10 @@ the EAGLE suite of simulations. The cosmology, resolution
and phases are the same as used in the original suite. The only and phases are the same as used in the original suite. The only
difference is the file format, adapted for SWIFT. difference is the file format, adapted for SWIFT.
For standard runs, the code should be configured with
--with-subgrid=EAGLE-XL --with-hydro=sphenix --with-kernel=wendland-C2
Compared to the original EAGLE runs of Schaye et al. 2015), Compared to the original EAGLE runs of Schaye et al. 2015),
the following changes have been made (at standard resolution): the following changes have been made (at standard resolution):
...@@ -17,15 +21,33 @@ the following changes have been made (at standard resolution): ...@@ -17,15 +21,33 @@ the following changes have been made (at standard resolution):
- The metallicity-dependent density threshold for star formation - The metallicity-dependent density threshold for star formation
uses the smoothed metallicities and not the raw metallicities uses the smoothed metallicities and not the raw metallicities
any more. any more.
- The cooling uses the Ploeckinger+20 tables including an updated
UV background, newer version of CLOUDY and self-shielding
(was using the Wiersma+09 tables before).
- The redshift of H reionization has been lowered to 7.5 (from 11.5). - The redshift of H reionization has been lowered to 7.5 (from 11.5).
- The entropy floor between n_H = 10^-5 and 10^-1 cm^3 at 8000K has
been removed as the new cooling tables handle this correctly. The
gamma=4/3 floor has been extended to lower densities
(i.e. 800K at n_H = 10^-4 cm^-3) as a fail-safe.
- Particles can be star-forming if they are within 0.3 dex of the
entropy floor (was 0.5 dex). These particles also get their
subgrid properties (rho, T) computed.
- The minimal mass of SNII stars has been raised to 8 Msun (from 6). - The minimal mass of SNII stars has been raised to 8 Msun (from 6).
- The SNII feedback delay is done by sampling the stellar age - The SNII feedback delay is done by sampling the stellar age
distribution and not using a fixed delay of 30 Myr any more. distribution and not using a fixed delay of 30 Myr any more.
- The halo mass for BH seeding has been lowered to 10^10 M_sun
(from 10^10/h Msun).
- The black hole accretion rate is now limited to 100% of the - The black hole accretion rate is now limited to 100% of the
Eddington rate (from 100/h = 150 %). Eddington rate (from 100/h = 150 %).
- The circular velocity threshold for BH mergers is measured - The circular velocity threshold for BH mergers is measured
at the actual distance linking the BHs not at the kernel support at the actual distance linking the BHs not at the kernel support
length any more. length any more.
- The BHs can reposition only onto gas particles and not onto any
type as before.
- The maximal relative velocity between a BH and gas particle
eligible to be a reposition target has been increased to 0.5 times
the local sound-speed around the BH (was 0.25 * c_sound).
The scripts in this directory download the tables required to The scripts in this directory download the tables required to
run the EAGLE model. Plotting scripts are also provided run the EAGLE model. Plotting scripts are also provided
...@@ -38,5 +60,5 @@ cmake -DVR_USE_GAS=ON -DVR_USE_STAR=ON -DV_USE_BH=ON ...@@ -38,5 +60,5 @@ cmake -DVR_USE_GAS=ON -DVR_USE_STAR=ON -DV_USE_BH=ON
and run using and run using
stf -C vrconfig_3dfof_subhalos_SO_hydro.cfg -i eagle_0035 -o halos_0035 -I 2 stf -C vrconfig_3dfof_subhalos_SO_hydro.cfg -i eagle_0036 -o halos_0036 -I 2
#!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/CoolingTables/COLIBRE/UV_dust1_CR1_G1_shield1.hdf5
...@@ -99,6 +99,7 @@ EAGLEChemistry: # Solar abundances ...@@ -99,6 +99,7 @@ EAGLEChemistry: # Solar abundances
init_abundance_Silicon: 6.825874e-4 init_abundance_Silicon: 6.825874e-4
init_abundance_Iron: 1.1032152e-3 init_abundance_Iron: 1.1032152e-3
# EAGLE cooling parameters
EAGLECooling: EAGLECooling:
dir_name: ./coolingtables/ dir_name: ./coolingtables/
H_reion_z: 11.5 H_reion_z: 11.5
...@@ -107,6 +108,17 @@ EAGLECooling: ...@@ -107,6 +108,17 @@ EAGLECooling:
He_reion_z_sigma: 0.5 He_reion_z_sigma: 0.5
He_reion_eV_p_H: 2.0 He_reion_eV_p_H: 2.0
# 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.3
rapid_cooling_threshold: 0.333333
# EAGLE star formation parameters # EAGLE star formation parameters
EAGLEStarFormation: EAGLEStarFormation:
EOS_density_norm_H_p_cm3: 0.1 # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3. EOS_density_norm_H_p_cm3: 0.1 # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3.
...@@ -117,7 +129,7 @@ EAGLEStarFormation: ...@@ -117,7 +129,7 @@ EAGLEStarFormation:
min_over_density: 57.7 # The over-density above which star-formation is allowed. min_over_density: 57.7 # The over-density above which star-formation is allowed.
KS_high_density_threshold_H_p_cm3: 1e3 # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3. KS_high_density_threshold_H_p_cm3: 1e3 # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3.
KS_high_density_exponent: 2.0 # Slope of the Kennicut-Schmidt law above the high-density threshold. KS_high_density_exponent: 2.0 # Slope of the Kennicut-Schmidt law above the high-density threshold.
EOS_entropy_margin_dex: 0.5 # Logarithm base 10 of the maximal entropy above the EOS at which stars can form. EOS_entropy_margin_dex: 0.3 # Logarithm base 10 of the maximal entropy above the EOS at which stars can form.
threshold_norm_H_p_cm3: 0.1 # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3. threshold_norm_H_p_cm3: 0.1 # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3.
threshold_Z0: 0.002 # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation. threshold_Z0: 0.002 # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation.
threshold_slope: -0.64 # Slope of the metal-dependant star formation threshold threshold_slope: -0.64 # Slope of the metal-dependant star formation threshold
......
...@@ -107,6 +107,7 @@ EAGLEChemistry: # Solar abundances ...@@ -107,6 +107,7 @@ EAGLEChemistry: # Solar abundances
init_abundance_Silicon: 6.825874e-4 init_abundance_Silicon: 6.825874e-4
init_abundance_Iron: 1.1032152e-3 init_abundance_Iron: 1.1032152e-3
# EAGLE cooling parameters
EAGLECooling: EAGLECooling:
dir_name: ./coolingtables/ dir_name: ./coolingtables/
H_reion_z: 11.5 H_reion_z: 11.5
...@@ -115,6 +116,17 @@ EAGLECooling: ...@@ -115,6 +116,17 @@ EAGLECooling:
He_reion_z_sigma: 0.5 He_reion_z_sigma: 0.5
He_reion_eV_p_H: 2.0 He_reion_eV_p_H: 2.0
# 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.3
rapid_cooling_threshold: 0.333333
# EAGLE star formation parameters # EAGLE star formation parameters
EAGLEStarFormation: EAGLEStarFormation:
EOS_density_norm_H_p_cm3: 0.1 # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3. EOS_density_norm_H_p_cm3: 0.1 # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3.
...@@ -125,7 +137,7 @@ EAGLEStarFormation: ...@@ -125,7 +137,7 @@ EAGLEStarFormation:
min_over_density: 57.7 # The over-density above which star-formation is allowed. min_over_density: 57.7 # The over-density above which star-formation is allowed.
KS_high_density_threshold_H_p_cm3: 1e3 # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3. KS_high_density_threshold_H_p_cm3: 1e3 # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3.
KS_high_density_exponent: 2.0 # Slope of the Kennicut-Schmidt law above the high-density threshold. KS_high_density_exponent: 2.0 # Slope of the Kennicut-Schmidt law above the high-density threshold.
EOS_entropy_margin_dex: 0.5 # Logarithm base 10 of the maximal entropy above the EOS at which stars can form. EOS_entropy_margin_dex: 0.3 # Logarithm base 10 of the maximal entropy above the EOS at which stars can form.
threshold_norm_H_p_cm3: 0.1 # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3. threshold_norm_H_p_cm3: 0.1 # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3.
threshold_Z0: 0.002 # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation. threshold_Z0: 0.002 # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation.
threshold_slope: -0.64 # Slope of the metal-dependant star formation threshold threshold_slope: -0.64 # Slope of the metal-dependant star formation threshold
......
...@@ -98,6 +98,7 @@ EAGLEChemistry: # Solar abundances ...@@ -98,6 +98,7 @@ EAGLEChemistry: # Solar abundances
init_abundance_Silicon: 6.825874e-4 init_abundance_Silicon: 6.825874e-4
init_abundance_Iron: 1.1032152e-3 init_abundance_Iron: 1.1032152e-3
# EAGLE cooling parameters
EAGLECooling: EAGLECooling:
dir_name: ./coolingtables/ dir_name: ./coolingtables/
H_reion_z: 11.5 H_reion_z: 11.5
...@@ -106,6 +107,17 @@ EAGLECooling: ...@@ -106,6 +107,17 @@ EAGLECooling:
He_reion_z_sigma: 0.5 He_reion_z_sigma: 0.5
He_reion_eV_p_H: 2.0 He_reion_eV_p_H: 2.0
# 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.3
rapid_cooling_threshold: 0.333333
# EAGLE star formation parameters # EAGLE star formation parameters
EAGLEStarFormation: EAGLEStarFormation:
EOS_density_norm_H_p_cm3: 0.1 # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3. EOS_density_norm_H_p_cm3: 0.1 # Physical density used for the normalisation of the EOS assumed for the star-forming gas in Hydrogen atoms per cm^3.
...@@ -116,7 +128,7 @@ EAGLEStarFormation: ...@@ -116,7 +128,7 @@ EAGLEStarFormation:
min_over_density: 57.7 # The over-density above which star-formation is allowed. min_over_density: 57.7 # The over-density above which star-formation is allowed.
KS_high_density_threshold_H_p_cm3: 1e3 # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3. KS_high_density_threshold_H_p_cm3: 1e3 # Hydrogen number density above which the Kennicut-Schmidt law changes slope in Hydrogen atoms per cm^3.
KS_high_density_exponent: 2.0 # Slope of the Kennicut-Schmidt law above the high-density threshold. KS_high_density_exponent: 2.0 # Slope of the Kennicut-Schmidt law above the high-density threshold.
EOS_entropy_margin_dex: 0.5 # Logarithm base 10 of the maximal entropy above the EOS at which stars can form. EOS_entropy_margin_dex: 0.3 # Logarithm base 10 of the maximal entropy above the EOS at which stars can form.
threshold_norm_H_p_cm3: 0.1 # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3. threshold_norm_H_p_cm3: 0.1 # Normalisation of the metal-dependant density threshold for star formation in Hydrogen atoms per cm^3.
threshold_Z0: 0.002 # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation. threshold_Z0: 0.002 # Reference metallicity (metal mass fraction) for the metal-dependant threshold for star formation.
threshold_slope: -0.64 # Slope of the metal-dependant star formation threshold threshold_slope: -0.64 # Slope of the metal-dependant star formation threshold
......