20 #define NVAR_COOLING (NVAR+1) 27 double CompEquil (
double,
double,
double *);
29 double GetMaxRate (
double *,
double *,
double);
30 void Jacobian (
double *,
double *,
double **);
31 void NormalizeIons (
double *);
33 #if COOLING == POWER_LAW 36 void Radiat (
double *,
double *);
42 double SolveODE_ROS34 (
double *,
double *,
double *,
double,
double);
48 #if COOLING == H2_COOL 53 #define X_H2 (NFLX + 1) 54 #define X_HII (NFLX + 2) 61 #elif COOLING == MINEq 76 #define C_EXPAND(a,b,c,d,e) 78 #define C_EXPAND(a,b,c,d,e) ,a 80 #define C_EXPAND(a,b,c,d,e) ,a,b 82 #define C_EXPAND(a,b,c,d,e) ,a,b,c 84 #define C_EXPAND(a,b,c,d,e) ,a,b,c,d 86 #define C_EXPAND(a,b,c,d,e) ,a,b,c,d,e 90 #define N_EXPAND(a,b,c,d,e) 92 #define N_EXPAND(a,b,c,d,e) ,a 94 #define N_EXPAND(a,b,c,d,e) ,a,b 96 #define N_EXPAND(a,b,c,d,e) ,a,b,c 98 #define N_EXPAND(a,b,c,d,e) ,a,b,c,d 100 #define N_EXPAND(a,b,c,d,e) ,a,b,c,d,e 104 #define O_EXPAND(a,b,c,d,e) 106 #define O_EXPAND(a,b,c,d,e) ,a 108 #define O_EXPAND(a,b,c,d,e) ,a,b 110 #define O_EXPAND(a,b,c,d,e) ,a,b,c 112 #define O_EXPAND(a,b,c,d,e) ,a,b,c,d 114 #define O_EXPAND(a,b,c,d,e) ,a,b,c,d,e 118 #define Ne_EXPAND(a,b,c,d,e) 120 #define Ne_EXPAND(a,b,c,d,e) ,a 122 #define Ne_EXPAND(a,b,c,d,e) ,a,b 124 #define Ne_EXPAND(a,b,c,d,e) ,a,b,c 126 #define Ne_EXPAND(a,b,c,d,e) ,a,b,c,d 128 #define Ne_EXPAND(a,b,c,d,e) ,a,b,c,d,e 132 #define S_EXPAND(a,b,c,d,e) 134 #define S_EXPAND(a,b,c,d,e) ,a 136 #define S_EXPAND(a,b,c,d,e) ,a,b 138 #define S_EXPAND(a,b,c,d,e) ,a,b,c 140 #define S_EXPAND(a,b,c,d,e) ,a,b,c,d 142 #define S_EXPAND(a,b,c,d,e) ,a,b,c,d,e 146 #define Fe_EXPAND(a,b,c) 148 #define Fe_EXPAND(a,b,c) ,a 150 #define Fe_EXPAND(a,b,c) ,a,b 152 #define Fe_EXPAND(a,b,c) ,a,b,c 165 X_HI = NFLX, X_HeI, X_HeII
166 C_EXPAND(X_CI, X_CII, X_CIII, X_CIV, X_CV)
167 N_EXPAND(X_NI, X_NII, X_NIII, X_NIV, X_NV)
168 O_EXPAND(X_OI, X_OII, X_OIII, X_OIV, X_OV)
169 Ne_EXPAND(X_NeI, X_NeII, X_NeIII, X_NeIV, X_NeV)
170 S_EXPAND(X_SI, X_SII, X_SIII, X_SIV, X_SV)
171 Fe_EXPAND(X_FeI, X_FeII, X_FeIII)
174 #define NIONS (3+C_IONS+N_IONS+O_IONS+Ne_IONS+S_IONS+Fe_IONS) 177 double find_N_rho ();
179 #elif COOLING == POWER_LAW 183 #elif COOLING == TABULATED double SolveODE_RKF12(double *, double *, double *, double, intList *)
Definition: cooling_ode_solver.c:265
double SolveODE_RKF23(double *, double *, double *, double, intList *)
Definition: cooling_ode_solver.c:327
void Jacobian(double *, double *, double **)
Definition: jacobian.c:4
void PowerLawCooling(Data_Arr, double, timeStep *, Grid *)
Definition: cooling.c:50
double SolveODE_RK4(double *, double *, double *, double, intList *)
Definition: cooling_ode_solver.c:221
void Radiat(double *, double *)
Definition: radiat.c:113
Definition: structs.h:503
Definition: structs.h:124
void CoolingSource(const Data *, double, timeStep *, Grid *)
Definition: cooling_source.c:44
double H_MassFrac(void)
Definition: radiat.c:738
Definition: structs.h:514
double SolveODE_CK45(double *, double *, double *, double, double, intList *)
Definition: cooling_ode_solver.c:26
double GetMaxRate(double *, double *, double)
Definition: maxrate.c:4
void Numerical_Jacobian(double *, double **)
Definition: cooling_source.c:229
void H2RateTables(double, double *)
Definition: radiat.c:25
double CompEquil(double, double, double *)
Definition: comp_equil.c:44
Definition: structs.h:363