19 static const double dicoef[2][NDIM] = {
20 {1.58e-3,8.38e-3,1.54e-2,3.75e-2,0.117,0.254,0.291,0.150,0.140,0.100,0.200,0.240,
21 0.260,0.190,0.120,0.350,0.066,0.10,0.13,0.23,0.14,0.11,0.041,0.747,0.519,0.},
22 {.456,.323,.310,.411,.359,.0975,.229,4.20,3.30,5.30,1.50,0.700,.600,.5,1.,0.,7.8,
23 6.3,5.5,3.6,4.9,1.6,4.2,.284,.279,0.}
25 static const double dite[2][NDIM] = {
26 {6.00e4,1.94e5,3.31e5,4.32e5,6.28e5,7.50e5,7.73e5,2.62e5,2.50e5,2.57e5,2.84e5,
27 8.69e5,4.21e5,4.57e5,2.85e5,8.18e5,1.51e6,1.30e6,1.19e6,1.09e6,9.62e5,7.23e5,
28 4.23e5,5.84e7,6.01e7,0.},
29 {8.97e4,1.71e5,2.73e5,3.49e5,5.29e5,4.69e5,6.54e5,1.32e6,1.33e6,1.41e6,1.52e6,
30 1.51e6,1.82e6,1.84e6,2.31e6,0.,9.98e6,9.98e6,1.00e7,1.10e7,8.34e6,1.01e7,1.07e7,
33 static const double ditcrt[NDIM] = {1e11,1e11,1e11,1e11,1e11,1e11,1e11,
34 1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,
35 1e11,1e11,1e11,1e11,1e11,1e11,1e11};
36 static const double aa[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
37 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
38 static const double bb[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
39 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
40 static const double cc[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
41 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
42 static const double dd[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
43 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
47 static const double ff[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
48 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
49 static const double fyield[NDIM+1] = {.34,.34,.35,.35,.36,.37,.37,.38,.39,.40,
50 .41,.42,.43,.44,.45,.46,.47,.47,.48,.48,.49,.49,.11,.75,0.,0.,0.};
52 long int i, limit, limit2;
81 enum {DEBUG_LOC=
false};
85 long int iplow ,
iphi , ipop , ns , ion;
102 ion_recomb(
false,(
const double*)dicoef,(
const double*)dite,ditcrt,aa,bb,cc,dd,ff,
ipIRON);
140 for( i=limit; i < limit2; i++ )
158 fprintf(
ioQQQ,
" Fe" );
159 for( i=0; i < 15; i++ )
163 fprintf(
ioQQQ,
"\n" );
168 fprintf(
ioQQQ,
" IonIron-Abund:" );
169 for( i=0; i < 27; i++ )
173 fprintf(
ioQQQ,
"\n" );