12 static const double dicoef[2][NDIM] = {
13 {1.10e-3,5.87e-3,5.03e-3,5.43e-3,8.86e-3,1.68e-2,2.49e-2,3.13e-2,
14 4.25e-2,6.18e-2,1.38e-2,.327,.189,0.},
15 {0.,.753,.188,.450,0.,1.80,1.88,2.01,1.22,.303,1.42,.306,.286,0.}
17 static const double dite[2][NDIM] = {
18 {7.70e4,9.63e4,8.75e4,1.05e6,1.14e6,4.85e5,4.15e5,3.66e5,3.63e5,
19 3.88e5,2.51e5,1.88e7,1.99e7,0.},
20 {0.,6.46e4,4.71e4,7.98e5,0.,1.03e6,1.91e6,2.11e6,2.14e6,1.12e6,
21 3.93e6,3.60e6,4.14e6,0.}
23 static const double ditcrt[NDIM] = {1.1e4,1.1e4,1.1e4,1.7e5,9.5e4,8.0e4,
24 7.4e4,6.8e4,6.6e4,6.5e4,4.5e4,3.7e6,6.3e6,1e20};
25 static const double aa[NDIM] = {-0.0219,3.2163,0.1203,0.,0.,0.,0.,0.,
27 static const double bb[NDIM] = {0.4364,-12.0571,-2.6900,0.,0.,0.,0.,
28 0.,0.,0.,0.,0.,0.,0.};
29 static const double cc[NDIM] = {0.0684,16.2118,19.1943,0.,0.,0.,0.,
30 0.,0.,0.,0.,0.,0.,0.};
31 static const double dd[NDIM] = {-0.0032,-0.5886,-0.1479,0.,0.,0.,0.,
32 0.,0.,0.,0.,0.,0.,0.};
33 static const double ff[NDIM] = {0.1342,0.5613,0.1118,0.1,0.,0.,0.,0.,
36 bool lgPrtDebug=
false;
49 static long nzUsed = -1;
50 static double OldRate = 0.;
55 if(
nzone > 1 && OldRate > 0. )
76 ion_recomb(
false,(
const double*)dicoef,(
const double*)dite,ditcrt,aa,bb,cc,dd,ff,
ipSILICON);
89 for(
int i=0; i <= 10; i++ )
93 fprintf(
ioQQQ,
"\n" );
96 for( i=0; i < 2; i++ )
100 fprintf(
ioQQQ,
"\t%.3e",
co.hevmol[ipATSI]/
SDIV(
co.hevmol[ipSIP]));
107 fprintf(
ioQQQ,
"\n" );