36 fprintf(
ioQQQ,
" UPDATE estimating new optical depths\n" );
39 fprintf(
ioQQQ,
" New Hydrogen outward optical depths:\n" );
42 fprintf(
ioQQQ,
"%3ld", ipHi );
43 for( ipLo=0; ipLo < ipHi; ipLo++ )
46 fprintf(
ioQQQ,
"%10.2e", 1e-30 );
48 fprintf(
ioQQQ,
"%10.2e",
51 fprintf(
ioQQQ,
"\n" );
61 for( i=0; i<
NCOLD; ++i )
88 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
98 for( ipLo=0; ipLo < ipHi; ipLo++ )
100 enum {DEBUG_LOC=
false};
103 if( ipISO==1 && nelem==1 && ipHi==
iso.
nLyaLevel[ipISO] && ipLo==0 )
104 fprintf(
ioQQQ,
"DEBUG rt before loop %li %li %li %li tot %.3e in %.3e\n",
105 ipISO, nelem, ipHi , ipLo ,
114 if( ipISO==1 && nelem==1 && ipHi==
iso.
nLyaLevel[ipISO] && ipLo==0 )
115 fprintf(
ioQQQ,
"DEBUG rt after loop %li %li %li %li tot %.3e in %.3e\n",
116 ipISO, nelem, ipHi , ipLo ,
138 for( nelem=0; nelem <
LIMELM; nelem++ )
168 for( nelem=1; nelem <
LIMELM; nelem++ )
249 for( i=0; i <
nUTA; i++ )