20 static const double dicoef[2][NDIM] = {
21 {2.54e-3,6.15e-3,1.62e-3,4.78e-2,3.22e-2,0.}, {4.42e-2,5.88e-2,0.343,0.362,0.315,0.}
23 static const double dite[2][NDIM] = {
24 {1.57e5,1.41e5,8.19e4,3.44e6,4.06e6,0.}, {3.74e5,1.41e5,1.59e5,5.87e5,8.31e5,0.}
26 static const double ditcrt[NDIM] = {1.2e4,1.2e4,1.1e4,4.4e5,7.0e5,1e20};
28 static const double aa[NDIM] = {.0108,1.8267,2.3196,0.,0.,0.};
29 static const double bb[NDIM] = {-0.1075,4.1012,10.7328,0.,0.,0.};
30 static const double cc[NDIM] = {.2810,4.8443,6.8830,0.,0.,0.};
31 static const double dd[NDIM] = {-0.0193,.2261,-0.1824,0.,0.,0.};
33 static const double ff[NDIM] = {0.000,0.5960,0.4101,0.1,0.1,0.};
41 fprintf(
ioQQQ,
" IonCarbo called.\n" );
81 ion_recomb(
false,(
const double*)dicoef,(
const double*)dite,ditcrt,aa,bb,cc,dd,ff,
ipCARBON);
115 fprintf(
ioQQQ,
" IonCarbo returns; fracs=" );
116 for(
int i=0; i < 7; i++ )
121 fprintf(
ioQQQ,
"\n" );
139 ion_recomb(
false,(
double*)dicoef,(
double*)dite,ditcrt,aa,bb,cc,