From e303f7a11787cc4047fd2e2b87efd493bca9b56b Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <schaller@strw.leidenuniv.nl>
Date: Sat, 25 Jul 2020 14:29:06 +0200
Subject: [PATCH] Add the collection of the BH subgrid masses to the default VR
configs in the EAGLE ICs examples
---
.../vrconfig_3dfof_subhalos_SO_hydro.cfg | 10 ++++++++++
.../vrconfig_3dfof_subhalos_SO_hydro.cfg | 10 ++++++++++
.../vrconfig_3dfof_subhalos_SO_hydro.cfg | 10 ++++++++++
.../vrconfig_3dfof_subhalos_SO_hydro.cfg | 10 ++++++++++
examples/EAGLE_ICs/README | 17 +++++++++++++----
5 files changed, 53 insertions(+), 4 deletions(-)
diff --git a/examples/EAGLE_ICs/EAGLE_12/vrconfig_3dfof_subhalos_SO_hydro.cfg b/examples/EAGLE_ICs/EAGLE_12/vrconfig_3dfof_subhalos_SO_hydro.cfg
index be762a68c8..ad342fa599 100644
--- a/examples/EAGLE_ICs/EAGLE_12/vrconfig_3dfof_subhalos_SO_hydro.cfg
+++ b/examples/EAGLE_ICs/EAGLE_12/vrconfig_3dfof_subhalos_SO_hydro.cfg
@@ -169,9 +169,19 @@ Reference_frame_for_properties=2 #use the minimum potential as reference frame a
################################
#Extra fields
################################
+
+# Also compute the H abundance as well as the HI, HII and H_2 fractions
Gas_internal_property_names=ElementMassFractions,SpeciesFractions,SpeciesFractions,SpeciesFractions,
Gas_internal_property_index_in_file=0,0,1,2,
+Gas_internal_property_input_output_unit_conversion_factors=1.0,1.0,1.0,1.0
Gas_internal_property_calculation_type =averagemassweighted,averagemassweighted,averagemassweighted,averagemassweighted,
+Gas_internal_property_output_units=unitless,unitless,unitless,unitless,
+
+# Collect the BH subgrid masses and compute the max, min, average and total mass in apertures
+BH_internal_property_names=SubgridMasses,SubgridMasses,SubgridMasses,SubgridMasses,
+BH_internal_property_input_output_unit_conversion_factors=1.0e10,1.0e10,1.0e10,1.0e10,
+BH_internal_property_calculation_type=max,min,average,aperture_total,
+BH_internal_property_output_units=solar_mass,solar_mass,solar_mass,solar_mass,
################################
#output related
diff --git a/examples/EAGLE_ICs/EAGLE_25/vrconfig_3dfof_subhalos_SO_hydro.cfg b/examples/EAGLE_ICs/EAGLE_25/vrconfig_3dfof_subhalos_SO_hydro.cfg
index be762a68c8..ad342fa599 100644
--- a/examples/EAGLE_ICs/EAGLE_25/vrconfig_3dfof_subhalos_SO_hydro.cfg
+++ b/examples/EAGLE_ICs/EAGLE_25/vrconfig_3dfof_subhalos_SO_hydro.cfg
@@ -169,9 +169,19 @@ Reference_frame_for_properties=2 #use the minimum potential as reference frame a
################################
#Extra fields
################################
+
+# Also compute the H abundance as well as the HI, HII and H_2 fractions
Gas_internal_property_names=ElementMassFractions,SpeciesFractions,SpeciesFractions,SpeciesFractions,
Gas_internal_property_index_in_file=0,0,1,2,
+Gas_internal_property_input_output_unit_conversion_factors=1.0,1.0,1.0,1.0
Gas_internal_property_calculation_type =averagemassweighted,averagemassweighted,averagemassweighted,averagemassweighted,
+Gas_internal_property_output_units=unitless,unitless,unitless,unitless,
+
+# Collect the BH subgrid masses and compute the max, min, average and total mass in apertures
+BH_internal_property_names=SubgridMasses,SubgridMasses,SubgridMasses,SubgridMasses,
+BH_internal_property_input_output_unit_conversion_factors=1.0e10,1.0e10,1.0e10,1.0e10,
+BH_internal_property_calculation_type=max,min,average,aperture_total,
+BH_internal_property_output_units=solar_mass,solar_mass,solar_mass,solar_mass,
################################
#output related
diff --git a/examples/EAGLE_ICs/EAGLE_50/vrconfig_3dfof_subhalos_SO_hydro.cfg b/examples/EAGLE_ICs/EAGLE_50/vrconfig_3dfof_subhalos_SO_hydro.cfg
index be762a68c8..ad342fa599 100644
--- a/examples/EAGLE_ICs/EAGLE_50/vrconfig_3dfof_subhalos_SO_hydro.cfg
+++ b/examples/EAGLE_ICs/EAGLE_50/vrconfig_3dfof_subhalos_SO_hydro.cfg
@@ -169,9 +169,19 @@ Reference_frame_for_properties=2 #use the minimum potential as reference frame a
################################
#Extra fields
################################
+
+# Also compute the H abundance as well as the HI, HII and H_2 fractions
Gas_internal_property_names=ElementMassFractions,SpeciesFractions,SpeciesFractions,SpeciesFractions,
Gas_internal_property_index_in_file=0,0,1,2,
+Gas_internal_property_input_output_unit_conversion_factors=1.0,1.0,1.0,1.0
Gas_internal_property_calculation_type =averagemassweighted,averagemassweighted,averagemassweighted,averagemassweighted,
+Gas_internal_property_output_units=unitless,unitless,unitless,unitless,
+
+# Collect the BH subgrid masses and compute the max, min, average and total mass in apertures
+BH_internal_property_names=SubgridMasses,SubgridMasses,SubgridMasses,SubgridMasses,
+BH_internal_property_input_output_unit_conversion_factors=1.0e10,1.0e10,1.0e10,1.0e10,
+BH_internal_property_calculation_type=max,min,average,aperture_total,
+BH_internal_property_output_units=solar_mass,solar_mass,solar_mass,solar_mass,
################################
#output related
diff --git a/examples/EAGLE_ICs/EAGLE_50_low_res/vrconfig_3dfof_subhalos_SO_hydro.cfg b/examples/EAGLE_ICs/EAGLE_50_low_res/vrconfig_3dfof_subhalos_SO_hydro.cfg
index be762a68c8..ad342fa599 100644
--- a/examples/EAGLE_ICs/EAGLE_50_low_res/vrconfig_3dfof_subhalos_SO_hydro.cfg
+++ b/examples/EAGLE_ICs/EAGLE_50_low_res/vrconfig_3dfof_subhalos_SO_hydro.cfg
@@ -169,9 +169,19 @@ Reference_frame_for_properties=2 #use the minimum potential as reference frame a
################################
#Extra fields
################################
+
+# Also compute the H abundance as well as the HI, HII and H_2 fractions
Gas_internal_property_names=ElementMassFractions,SpeciesFractions,SpeciesFractions,SpeciesFractions,
Gas_internal_property_index_in_file=0,0,1,2,
+Gas_internal_property_input_output_unit_conversion_factors=1.0,1.0,1.0,1.0
Gas_internal_property_calculation_type =averagemassweighted,averagemassweighted,averagemassweighted,averagemassweighted,
+Gas_internal_property_output_units=unitless,unitless,unitless,unitless,
+
+# Collect the BH subgrid masses and compute the max, min, average and total mass in apertures
+BH_internal_property_names=SubgridMasses,SubgridMasses,SubgridMasses,SubgridMasses,
+BH_internal_property_input_output_unit_conversion_factors=1.0e10,1.0e10,1.0e10,1.0e10,
+BH_internal_property_calculation_type=max,min,average,aperture_total,
+BH_internal_property_output_units=solar_mass,solar_mass,solar_mass,solar_mass,
################################
#output related
diff --git a/examples/EAGLE_ICs/README b/examples/EAGLE_ICs/README
index a80e1a17a1..16dfaccd07 100644
--- a/examples/EAGLE_ICs/README
+++ b/examples/EAGLE_ICs/README
@@ -7,17 +7,23 @@ For standard runs, the code should be configured with
--with-subgrid=EAGLE-XL --with-hydro=sphenix --with-kernel=wendland-C2
+(and any other optimization, library or compiler parameters specific
+to your system.)
+
Compared to the original EAGLE runs of Schaye et al. 2015),
the following changes have been made (at standard resolution):
+ - The flavour of SPH has been changed from the Pressure-Entropy
+ "ANARCHY" implementation to the Density-energy based "SPHENIX"
+ implementation of Borrow+20.
- The dark matter softening lengths have been increased to 1.3
pkpc and 3.32 ckpc. The comoving baryon softening lengths have
been changed to 1.79 ckpc. This follows the recommendations of
Ludlow et al. 2019. Old values were 0.7 pkpc and 2.69 ckpc for
all the particle species.
- - SPH particles with a mass larger than 7*10^6 Msun (~4x the initial
- gas particle mass) are now split into 2 equal mass particles
- within the smoothing length of the original particle.
+ - SPH particles reaching a mass larger than 7*10^6 Msun (~4x the
+ initial gas particle mass) are now split into 2 equal mass
+ particles within the smoothing length of the original particle.
- The metallicity-dependent density threshold for star formation
uses the smoothed metallicities and not the raw metallicities
any more.
@@ -31,7 +37,7 @@ the following changes have been made (at standard resolution):
(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.
+ subgrid properties (rho, T as well as the HI and H_2 frac) computed.
- 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
distribution and not using a fixed delay of 30 Myr any more.
@@ -53,6 +59,9 @@ The scripts in this directory download the tables required to
run the EAGLE model. Plotting scripts are also provided
for basic quantities.
+To use the cooling model based on the Wiersma+09 tables, replace
+EAGLE-XL by EAGLE in the configuration command line.
+
VELOCIraptor can be run on the output. The code is compiled
using
--
GitLab