13 static const double dicoef[2][NDIM] = {
14 {2.54e-3,6.15e-3,1.62e-3,4.78e-2,0.}, {4.42e-2,5.88e-2,0.343,0.362,0.}
16 static const double dite[2][NDIM] = {
17 {1.57e5,1.41e5,8.19e4,3.44e6,0.}, {3.74e5,1.41e5,1.59e5,5.87e5,0.}
19 static const double ditcrt[NDIM] = {1.2e4,1.2e4,1.1e4,4.4e5,1e20};
20 static const double aa[NDIM] = {1.8267,2.3196,0.,0.,0.};
21 static const double bb[NDIM] = {4.1012,10.7328,0.,0.,0.};
22 static const double cc[NDIM] = {4.8443,6.8830,0.,0.,0.};
23 static const double dd[NDIM] = {.2261,-0.1824,0.,0.,0.};
24 static const double ff[NDIM] = {0.5960,0.4101,0.1,0.1,0.};
56 ion_recomb(
false,(
const double*)dicoef,(
const double*)dite,ditcrt,aa,bb,cc,dd,ff,
ipBORON);
63 fprintf(
ioQQQ,
" Boroni returns; frac=" );
64 for(
int i=0; i <
ipBORON+2; i++ )
69 fprintf(
ioQQQ,
"\n" );