88 for( nelem=ipISO; nelem<
LIMELM; ++nelem )
107 # define ERR_CHK 1.001
121 " DISASTER radius_increment finds inconsistent populations, Pop2Ion zn %.2f ipISO %li nelem %li %.4e %.4e solver:%s\n",
127 fprintf(
ioQQQ,
" level solver %s found pop_ion_ov_neut of %.5e",
130 fprintf(
ioQQQ,
" ion_solver found pop_ion_ov_neut of %.5e\n",
133 "simple %.3e Test shows Pop2Ion (%.6e) > xIonDense[nelem]/[nelem+1] (%.6e) \n",
157 fprintf(
ioQQQ,
" DISASTER PopLo fnz %.2f ipISO %li nelem %li pop_ion_ov_neut %.2e 2pops %.6e %.6e solver:%s\n",
181 " radius_increment called; radius=%10.3e rinner=%10.3e DRAD=%10.3e drNext=%10.3e ROUTER=%10.3e DEPTH=%10.3e\n",
281 fprintf(
ioQQQ,
" Energy, flux, OTS:\n" );
288 fprintf(
ioQQQ,
"\n" );
297 # if !defined(NDEBUG)
303 fprintf(
ioQQQ,
" radius_increment finds negative intensity in flux.\n" );
304 fprintf(
ioQQQ,
" Intensity, frequency, pointer=%11.3e%11.3e%6ld\n",
310 fprintf(
ioQQQ,
" radius_increment finds negative intensity in otscon.\n" );
311 fprintf(
ioQQQ,
" Intensity, frequency, pointer=%11.3e%11.3e%6ld\n",
317 fprintf(
ioQQQ,
" radius_increment finds negative tmn.\n" );
318 fprintf(
ioQQQ,
" value, frequency, pointer=%11.3e%11.3e%6ld %4.4s\n",
324 fprintf(
ioQQQ,
" radius_increment finds negative intensity in otslin.\n" );
325 fprintf(
ioQQQ,
" Intensity, frequency, pointer=%11.3e%11.3e%6ld %4.4s\n",
331 fprintf(
ioQQQ,
" radius_increment finds negative intensity in outlin.\n" );
332 fprintf(
ioQQQ,
" Intensity, frequency, pointer=%11.3e%11.3e%6ld %4.4s\n",
338 fprintf(
ioQQQ,
" radius_increment finds negative intensity in ConInterOut.\n" );
339 fprintf(
ioQQQ,
" Intensity, frequency, pointer=%11.3e%11.3e%6ld %4.4s\n",
353 fprintf(
ioQQQ,
" Insanity has occurred, this is zone%4ld\n",
537 if( fabs( tau_effec - tau_true ) / t > 0.01 &&
545 " PROBLEM radius_increment Lyman continuum insanity zone %li, effective tau=%g, atomic tau=%g simple tau=%g\n",
628 double Reflec_Diffuse_Cont;
660 aver(
"zone",1.,1.,
" ");
709 for( ion=0; ion<nelem+2; ++ion )
790 for( i=0; i < 5; i++ )
797 for( i=0; i < 3; i++ )
804 for( i=0; i < 4; i++ )
820 for( nelem=0; nelem <
LIMELM; nelem++ )
832 for( nd=0; nd <
gv.
nBin; nd++ )
898 fprintf(
ioQQQ,
" radius_increment returns\n" );
918 fprintf( ioFile,
"%10.2e %10.2e \n",
rfield.
anu[i],