81 ASSERT( ipLoEnr >= 0 && ipHiEnr >= 0 );
94 RateInducRecCool = 0.;
110 RateInducRec += prod;
112 RateInducRecCool += prod*(
rfield.
anu[i-1] - emin);
118 for( i=ipLoEnr; i < limit; i++ )
130 RateInducRec += prod;
132 RateInducRecCool += prod*(
rfield.
anu[i] - emin);
149 for( i=ilo-1; i < iup; i++ )
158 RateInducRec += prod;
160 RateInducRecCool += prod*(
rfield.
anu[i] - emin);
180 *rcool = RateInducRecCool*
EN1RYD;
181 *ainduc = RateInducRec;
204 fprintf(
ioQQQ,
" GammaPrtShells nz\t%.2f \t%.2li %.2li ",
214 fprintf(
ioQQQ,
"\ttot\t%.2e", sum);
272 if( ipLoEnr >=
rfield.
nflux || ipLoEnr >= ipHiEnr )
277 fprintf( ioFILE,
" GammaPrt %.2f from ",
fnzone);
279 fprintf( ioFILE,
" to ");
281 fprintf( ioFILE,
"R rates >");
282 fprintf( ioFILE,
PrintEfmt(
"%9.2e",threshold));
283 fprintf( ioFILE,
" of total=");
284 fprintf( ioFILE,
PrintEfmt(
"%9.2e",total));
285 fprintf( ioFILE,
" (frac inc, otslin, otscon, ConInterOut, outlin ConOTS_local_OTS_rate ) chL, C\n");
287 if( threshold <= 0. || total <= 0. )
301 if( phisig > threshold || phisig < 0.)
306 fprintf( ioFILE,
"[%5ld]" , i );
308 fprintf( ioFILE,
PrintEfmt(
"%9.2e",phisig/total));
309 fprintf( ioFILE,
"%5.2f%5.2f%5.2f%5.2f%5.2f%5.2f %4.4s %4.4s %.2e \n",
322 for( i=k; i < j; i++ )
325 if( phisig > threshold || phisig < 0.)
331 fprintf( ioFILE,
"%5ld", i );
333 fprintf(ioFILE,
PrintEfmt(
"%9.2e",phisig/total));
334 fprintf( ioFILE,
"%5.2f%5.2f%5.2f%5.2f%5.2f%5.2f %4.4s %4.4s %.2e \n",
406 ipOffset = ipOpac - ipLoEnr;
420 for( i=ipLoEnr; i < limit; i++ )
441 for( i=ilo-1; i < iup; i++ )
533 for( i=ipLoEnr; i < limit; i++ )
554 for( i=ilo-1; i < iup; i++ )
563 gamPL_v += GamHi*csThresh;
620 if( ipLoEnr >=
rfield.
nflux || ipLoEnr >= ipHiEnr )
643 RateInducRecCool = 0.;
659 RateInducRec += prod;
661 RateInducRecCool += prod*(
rfield.
anu[ipLoEnr-1] - emin);
664 for( i=ipLoEnr; i < limit; i++ )
673 RateInducRec += prod;
675 RateInducRecCool += prod*(
rfield.
anu[i] - emin);
687 for( i=ilo-1; i < iup; i++ )
695 RateInducRec += prod;
697 RateInducRecCool += prod*(
rfield.
anu[i] - emin);
699 RateInducRec *= csThresh;
700 RateInducRecCool *= csThresh;
703 bnfunPL_v = (g + GamHi)*csThresh;
719 *rcool = RateInducRecCool*
EN1RYD;
720 *ainduc = RateInducRec;
744 long int nshell , ns;
752 fprintf(ioFILE ,
"GammaPrtRate: %li %li",ion , nelem );
753 for( ns=nshell-1; ns>=0; --ns )
760 fprintf(ioFILE ,
"\n");
774 fprintf(ioFILE ,
"\n");