Skip to content
Snippets Groups Projects
Commit 43c7bfba authored by Folkert Nobels's avatar Folkert Nobels
Browse files

Resolve misleading rho_crit and update parameter file

parent de5de7c1
No related branches found
No related tags found
1 merge request!705Star formation following Schaye08
......@@ -294,11 +294,10 @@ SchayeSF:
SchmidtLawExponent: 1.4 # The power law of the Kennicutt-Schmidt law
SchmidtLawHighDensExponent: 2.0 # The high density exponent for the Kennicutt-Schmidt law
SchmidtLawHighDens_thresh_HpCM3: 1e3
Schaye2004: 1 # whether to use the metallicity dependent critical star formation of Schaye (2004) (1) or not (0).
thresh_norm_HpCM3: 0.1 # Critical sf normalization to use (is not a normalization when Schaye2004=0, than it is the value.
thresh_norm_HpCM3: 0.1 # Critical sf normalization to use
thresh_max_norm_HpCM3: 10.0 # Maximum norm of the critical sf density
MetDep_Z0: 0.002 # Scale metallicity to use for the equation (not used for Schaye2004=0)
MetDep_SFthresh_Slope: -0.64 # Scaling of the critical density with the metallicity (not used for Schaye2004=0)
MetDep_Z0: 0.002 # Scale metallicity to use for the equation
MetDep_SFthresh_Slope: -0.64 # Scaling of the critical density with the metallicity
thresh_MaxPhysDensOn: 0 # Default is 0.
thresh_MaxOverDens_HpCM3: 1e5 # Density at which the SF law changes
EOS_Jeans_GammaEffective: 1.33333 # The polytropic index
......
......@@ -144,7 +144,7 @@ INLINE static int star_formation_potential_to_become_star(
/* Read the critical overdensity factor and the critical density of
* the universe to determine the critical density to form stars*/
const double rho_crit = cosmo->critical_density * starform->min_over_den;
const double rho_crit_times_min_over_den = cosmo->critical_density * starform->min_over_den;
/* double tempp = cooling_get_temperature() */
double tempp;
......@@ -155,7 +155,7 @@ INLINE static int star_formation_potential_to_become_star(
* threshold is reached or if the metallicity dependent
* threshold is reached, after this we calculate if the
* temperature is appropriate */
if (p->rho > rho_crit) {
if (p->rho > rho_crit_times_min_over_den) {
/* In this case there are actually multiple possibilities
* because we also need to check if the physical density exceeded
* the appropriate limit */
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment