13 static const double dicoef[2][NDIM] = {
14 {5.5e-3,6.0e-3,7.5e-3,5.6e-3,1.2e-2,1.9e-2,1.6e-2,1.7e-2,4.4e-2,1.7e-2,3.0e-1,3.4e-0,0.},
15 {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}
17 static const double dite[2][NDIM] = {
18 {9.0e4,7.5e4,8.5e5,5.0e5,4.4e5,3.8e5,3.4e5,3.4e5,3.9e5,2.3e5,1.6e7,1.7e7,1e20},
19 {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}
21 static const double ditcrt[NDIM] = {8.0e3,7.0e3,1.2e5,7.0e4,6.5e4,5.5e4,
22 5.5e4,5.5e4,5.8e4,4.0e4,2.9e6,5.5e6,1e20};
23 static const double aa[NDIM] = {0.0219,0.7086,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
24 static const double bb[NDIM] = {-0.4528,-3.1083,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
25 static const double cc[NDIM] = {2.5427,7.0422,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
26 static const double dd[NDIM] = {-0.1678,0.5998,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
27 static const double ff[NDIM] = {0.2276,0.4194,0.1,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
62 fprintf(
ioQQQ,
" IonAlumi returns; frac=" );
63 for(
int i=0; i < 10; i++ )
68 fprintf(
ioQQQ,
"\n" );