Commit 377e9bcb authored by Folkert Nobels's avatar Folkert Nobels
Browse files

Add temperature constant exp10(Tdex) to be used in the criteria

parent 95a69be1
...@@ -70,6 +70,9 @@ struct star_formation { ...@@ -70,6 +70,9 @@ struct star_formation {
/*! Dalla Vecchia & Schaye temperature criteria */ /*! Dalla Vecchia & Schaye temperature criteria */
double temperature_margin_threshold_dex; double temperature_margin_threshold_dex;
/*! 10^Tdex of Dalla Vecchia & SChaye temperature criteria */
double ten_to_temperature_margin_threshold_dex;
/*! gas fraction */ /*! gas fraction */
double fgas; double fgas;
...@@ -259,7 +262,7 @@ INLINE static int star_formation_is_star_forming( ...@@ -259,7 +262,7 @@ INLINE static int star_formation_is_star_forming(
/* Check the Scahye & Dalla Vecchia 2012 EOS-based temperature critrion */ /* Check the Scahye & Dalla Vecchia 2012 EOS-based temperature critrion */
return (temperature < return (temperature <
temperature_eos * exp10(starform->temperature_margin_threshold_dex)); temperature_eos * starform->ten_to_temperature_margin_threshold_dex);
} }
/** /**
...@@ -554,6 +557,9 @@ INLINE static void starformation_init_backend( ...@@ -554,6 +557,9 @@ INLINE static void starformation_init_backend(
starform->temperature_margin_threshold_dex = parser_get_opt_param_double( starform->temperature_margin_threshold_dex = parser_get_opt_param_double(
parameter_file, "EAGLEStarFormation:KS_temperature_margin_dex", FLT_MAX); parameter_file, "EAGLEStarFormation:KS_temperature_margin_dex", FLT_MAX);
starform->ten_to_temperature_margin_threshold_dex = exp10(
starform->temperature_margin_threshold_dex);
/* Read the normalization of the metallicity dependent critical /* Read the normalization of the metallicity dependent critical
* density*/ * density*/
starform->density_threshold_HpCM3 = parser_get_param_double( starform->density_threshold_HpCM3 = parser_get_param_double(
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment