33 double DepartCoefInv ,
62 for( ion=0; ion<nelem+1; ++ion )
78 for( i=loop_min; i < loop_max; i++ )
99 static double *TotBremsAllIons;
100 static bool lgFirstTime=
true;
101 double BremsThisEner,bfac, bhfac,sumion[LIMELM+1];
102 long int ion_lo , ion_hi;
114 for(ion=1; ion<=LIMELM; ++ion )
125 sumion[ion] *=
POW2((
double)ion)*bfac;
132 while( sumion[ion_lo]==0 && ion_lo<LIMELM-1 )
135 while( sumion[ion_hi]==0 && ion_hi>0 )
147 BremsThisEner = bhfac *
rfield.
gff[ion][i]*
157 for(ion=ion_lo; ion<=ion_hi; ++ion )
159 BremsThisEner += sumion[ion]*
rfield.
gff[ion][i];
161 TotBremsAllIons[i] = BremsThisEner;
230 for( i=0; i < limit; i++ )
258 hydro.lgHLyaMased =
true;
346 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
431 double opacity_this_species =
435 double opacity_fraction = 1. - opacity_this_species /
opac.
opacity_abs[ip-1];
436 if(opacity_fraction < 0)
437 opacity_fraction = 0.;
500 fprintf(
ioQQQ,
" OpacityAddTotal returns; grd rec eff (opac) for Hn=1,4%10.2e%10.2e%10.2e%10.2e%10.2e%10.2e%10.2e%10.2e%10.2e%10.2e HeI,II:%10.2e%10.2e\n",
517 enum {DEBUG_LOC=
false};
519 if( DEBUG_LOC && (
nzone>=378) )
525 fprintf(
ioQQQ,
"rtotsbugggg\t%li\t%.3e\t%.3e\t%.3e\t%.3e\n",