52 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
61 for( ipLo=0; ipLo < ipHi; ipLo++ )
66 for( ipHi=2; ipHi <
iso.
nLyman[ipISO]; ipHi++ )
95 for( i=0; i <
nUTA; i++ )
282 for( nelem=0; nelem <
LIMELM; nelem++ )
390 ASSERT( ipISO<NISO && nelem < LIMELM );
430 3.7e4*BalmerTauOn*AbunRatio + 1e-20);
433 3.7e4*BalmerTauOn*AbunRatio + 1e-20);
436 3.7e4*BalmerTauOn*AbunRatio + 1e-20);
442 balc =
rfield.
qhtot*2.1e-19*BalmerTauOn*AbunRatio + 1e-20;
445 (
realnum)(
MIN2(2e5, balc*3.7e4*BalmerTauOn+1e-20));
448 (
realnum)(
MIN2(2e5, balc*3.7e4*BalmerTauOn+1e-20));
451 (
realnum)(
MIN2(2e5, balc*3.7e4*BalmerTauOn+1e-20));
481 long ipISO =
ipH_LIKE, ipNS, ipNPlusOneP;
485 if( (
N_(ipLo) + 1 ) <=
492 ipNPlusOneP = ipNS + 1;
497 ipNPlusOneP = ipNS + 1;
570 for( nelem=1; nelem <
LIMELM; nelem++ )
621 for( nelem=0; nelem <
LIMELM; nelem++ )
636 "RT_tau_init insanity for h line, Z=%li lo=%li hi=%li tot=%g in=%g \n",
647 fprintf(
ioQQQ,
" stopping due to insanity in RT_tau_init\n");
676 fprintf(
ioQQQ,
"\n\n up Transitions[ipISO][hi][lo].TauTot array as set in RT_tau_init ipZTrace (nelem)= %ld\n",
680 fprintf(
ioQQQ,
" %3ld", ipHi );
681 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
689 fprintf(
ioQQQ,
"\n" );
692 fprintf(
ioQQQ,
"\n\n TauIn array\n" );
695 fprintf(
ioQQQ,
" %3ld", ipHi );
696 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
704 fprintf(
ioQQQ,
"\n" );
707 fprintf(
ioQQQ,
"\n\n Aul As array\n" );
710 fprintf(
ioQQQ,
" %3ld", ipHi );
711 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
719 fprintf(
ioQQQ,
"\n" );
722 fprintf(
ioQQQ,
"\n\n Aul*esc array\n" );
725 fprintf(
ioQQQ,
" %3ld", ipHi );
726 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
737 fprintf(
ioQQQ,
"\n" );
740 fprintf(
ioQQQ,
"\n\n H opac array\n" );
743 fprintf(
ioQQQ,
" %3ld", ipHi );
744 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
752 fprintf(
ioQQQ,
"\n" );
758 fprintf(
ioQQQ,
" RT_tau_init called.\n" );