115 #if (COOLING == MINEq)
137 #elif COOLING == TABULATED
142 #elif COOLING == SNEq
152 #elif COOLING == H2_COOL
157 v[nv] =
MAX(v[nv], 0.0);
158 v[nv] =
MIN(v[nv], 1.0);
165 double hn = v[
X_HII];
187 #elif COOLING == MINEq
193 for (i = 0; i <
NIONS; i++) {
208 mmw2 += elem_ab[
el_H]*(1.0 - v[
X_HI])*2.;
215 for (i = 0; i <
NIONS; i++) {
221 print(
">>> Error! MMW2 NaN!\n");
222 for (i = 0; i <
NIONS; i++) {
const double rad_rec_z[31]
double MeanMolecularWeight(double *v)
const double elem_mass[8]
long int g_stepNumber
Gives the current integration step number.
#define CONST_AZ
Mean atomic weight of heavy elements.
void print(const char *fmt,...)
#define CONST_AHe
Atomic weight of Helium.
#define QUIT_PLUTO(e_code)