15 long int np, nd, ns, nc;
27 eff_charge = nelem + 1 - ipISO;
31 a = sqrt( 1.8887E8 * (nelem+1.) / pow((
double)
dense.xNucleiTotal, 0.333) );
37 a = 2.6E7 * eff_charge * eff_charge * pow(
phycon.te/
dense.eden, 0.25);
45 a = 3171. * pow( (
double)eff_charge, 0.8 ) * pow(
dense.eden + (
double)
dense.xIonDense[
ipHYDROGEN][1]
50 nc =
MIN3(np, nd, ns);
54 if( nc <= sp->n_HighestResolved_max)
96 if( !
conv.nTotalIoniz )
110 for(
long ipLo=0; ipLo < ipHi; ++ipLo )
116 fprintf(
ioQQQ,
" iso_continuum_lower: ipISO %li nelem %li nc %li (np:%li,nd:%li,ns:%li) numLevels %li nCollapsed %li n_HighestResolved %li \n",
CollisionProxy Coll() const
TransitionProxy trans(const long ipHi, const long ipLo)
long int n_HighestResolved_local
long int n_HighestResolved_max
long int nCollapsed_local
void CollisionZero(const CollisionProxy &t)
t_iso_sp iso_sp[NISO][LIMELM]
long iso_get_total_num_levels(long ipISO, long nmaxResolved, long numCollapsed)
void iso_continuum_lower(long ipISO, long nelem)