13 static const double dicoef[2][NDIM] = {
14 {1.0e-3,2.6e-3,6.0e-3,1.1e-2,8.0e-3,1.0e-2,3.2e-2,1.2e-2,2.2e-1,1.8e-1,0.},
15 {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}
17 static const double dite[2][NDIM] = {
18 {3.6e5,3.8e5,3.4e5,3.0e5,2.7e5,2.8e5,3.2e5,1.9e5,1.2e7,1.3e7,0.},
19 {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}
21 static const double ditcrt[NDIM] = {0.,1.3e5,4.4e5,7.2e4,2.8e5,4.0e5,
22 1.3e6,1.6e6,1.9e6,3.4e4,0.};
23 static const double aa[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
24 static const double bb[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
25 static const double cc[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
26 static const double dd[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
27 static const double ff[NDIM] = {0.1,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
51 ion_recomb(
false,(
const double*)dicoef,(
const double*)dite,ditcrt,aa,bb,cc,dd,ff,
ipSODIUM);
58 fprintf(
ioQQQ,
" IonSodiu returns; frac=" );
59 for(
int i=0; i < 10; i++ )
64 fprintf(
ioQQQ,
"\n" );