114 double factor, ConvLTEPOP, edenCorr;
119 {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0},
120 {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0} };
126 {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0},
127 {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0} };
146 " iso_collide called %s nelem %li - no reeval Boltz fac, LTE dens\n",
157 " iso_collide called %s nelem %li - will reeval Boltz fac, LTE dens\n",
170 for( ipLo=0; ipLo<ipHi; ipLo++ )
212 if(
iso.
PopLTE[ipISO][nelem][ipLo] <= 0. )
228 if( (TeUsedForCS[ipISO][nelem] == 0.) ||
229 ( TeUsedForCS[ipISO][nelem]/
phycon.
te > 1.15 ) ||
230 ( TeUsedForCS[ipISO][nelem]/
phycon.
te < 0.85 ) )
233 TeUsedForCS[ipISO][nelem] =
phycon.
te;
237 for( ipLo=0; ipLo < ipHi; ipLo++ )
271 if(
N_(ipHi) !=
N_(ipLo) )
311 Transitions[ipISO][nelem][ipHi][ipLo].Coll.cs = 0.;
326 for( ipLo=0; ipLo<ipHi; ipLo++ )
370 Transitions[ipISO][nelem][ipHi][ipLo].Coll.ColUL *=
383 Transitions[ipISO][nelem][ipHi][ipLo].Coll.ColUL *=
392 fprintf(
ioQQQ,
" iso_collide: %s Z=%li de-excitation rates coefficients\n",
iso.
chISO[ipISO], nelem + 1 );
395 fprintf(
ioQQQ,
" %li\t", ipHi );
396 for( ipLo=0; ipLo < ipHi; ipLo++ )
400 fprintf(
ioQQQ,
"\n" );
403 fprintf(
ioQQQ,
" iso_collide: %s Z=%li collisional ionization coefficients\n",
iso.
chISO[ipISO], nelem + 1 );
408 fprintf(
ioQQQ,
"\n" );
410 fprintf(
ioQQQ,
" iso_collide: %s Z=%li continuum boltzmann factor\n",
iso.
chISO[ipISO], nelem + 1 );
415 fprintf(
ioQQQ,
"\n" );
417 fprintf(
ioQQQ,
" iso_collide: %s Z=%li continuum boltzmann factor\n",
iso.
chISO[ipISO], nelem + 1 );
422 fprintf(
ioQQQ,
"\n" );