More flexible energy scaling function
Merge request to add the option for different energy scaling functions. Three options are now available, set by the parameter EAGLEFeedback:SNII_energy_fraction_function:
-
EAGLE: as in the original EAGLE simulations. This gives iso-fth contours that are diagonal lines in the log-Z vs. log-n plane. -
Separable: a transition from fE_min to fE_max that is separable in density and metallicity. This gives "boxy" iso-fth contours which curve upwards relative to the EAGLE scaling. -
Independent: a model where the energy scalings with metallicity and density are fully independent of each other. This creates a "four-box" pattern, with constant f_th towards each of the corners of the log-Z/log-n diagram and smooth transitions between them. This model has an additional parameter,EAGLEFeedback:SNII_energy_fraction_delta_E_n, which sets the maximum energy difference between low and high density at fixed metallicity.
The change only applies to the EAGLE_thermal model, the EAGLE_kinetic model is not changed @jborrow, please also have a look and check whether you can spot any issues.