21 static const double dicoef[2][NDIM] = {
22 {1.11e-3,5.07e-3,1.48e-2,1.84e-2,4.13e-3,1.06e-1,6.23e-2,0.},
23 {.0925,.181,.305,.1,.162,.34,.304,0.}
25 static const double dite[2][NDIM] = {
26 {1.75e5,1.98e5,2.41e5,2.12e5,1.25e5,6.25e6,7.01e6,0.},
27 {1.45e5,3.35e5,2.83e5,2.83e5,2.27e5,1.12e6,1.47e6,0.}
29 static const double ditcrt[NDIM] = {2.7e4,2.2e4,2.4e4,2.5e4,1.6e4,1.0e6,1.5e6,1e20};
30 static const double aa[NDIM] = {0.,-0.0036,0.,0.0061,-2.8425,0.,0.,0.};
31 static const double bb[NDIM] = {0.0238,0.7519,21.8790,0.2269,0.2283,0.,0.,0.};
32 static const double cc[NDIM] = {0.0659,1.5252,16.2730,32.1419,40.4072,0.,0.,0.};
33 static const double dd[NDIM] = {0.0349,-0.0838,-0.7020,1.9939,-3.4956,0.,0.,0.};
34 static const double ff[NDIM] = {0.5334,0.2769,1.1899,-0.0646,1.7558,0.,0.,0.};
38 double aeff, save_rec;
69 ion_recomb(
false,(
const double*)dicoef,(
const double*)dite,ditcrt,aa,bb,cc,dd,ff,
ipOXYGEN);
102 enum {DEBUG_LOC=
false};
106 fprintf(
ioQQQ,
"pop rel %.1e rate %.1e grnd rate %.1e\n",
194 co.hevmol[ipATO]/
co.hevmol[ipOP];
200 if(0 &&
nzone > 100 )
206 fprintf(
ioQQQ,
"DEBUG O\t%.3e\t%.3e\tH\t%.3e\t%.3e\n",
225 fprintf(
ioQQQ,
" IonOxyge returns; frac=" );
226 for(
int i=1; i <= 9; i++ )
231 fprintf(
ioQQQ,
"\n" );