47 static bool lgPntEval =
false;
58 fprintf(
ioQQQ,
" ContCreateMesh called, not evaluating.\n" );
72 fprintf(
ioQQQ,
" ContCreateMesh called first time.\n" );
125 fprintf(
ioQQQ,
" Currently the arrays that hold interpolated tables can only hold %i points.\n",
NCELL);
126 fprintf(
ioQQQ,
" This continuum mesh really needs to have %li points.\n",
rfield.
nupper);
127 fprintf(
ioQQQ,
" Please increase the value of NCELL in rfield.h and recompile.\n Sorry.");
268 ASSERT( fenlo>0. && fenhi>0. && resolv>0. );
271 nbin = (
long int)(log(10.)*log10(fenhi/fenlo)/resolv + 1);
282 fprintf(
ioQQQ,
" FILL improper bounds.\n" );
283 fprintf(
ioQQQ,
" ipnt=%3ld fenlo=%11.4e filbnd(ipnt)=%11.4e\n",
310 fprintf(
ioQQQ,
" Fill would need %ld cells to get to an energy of %.3e\n",
312 fprintf(
ioQQQ,
" This is a major logical error in fill.\n");
318 for( i=0; i < nbin; i++ )
321 aaa = pow( 10. , bbb );
335 " FILL range%2ld from%10.3e to%10.3eR in%4ld cell; ener res=%10.3e WIDTOT=%10.3e\n",
343 fprintf(
ioQQQ,
" The requested range was%10.3e%10.3e The requested resolution was%10.3e\n",
344 fenlo, fenhi, resolv );
374 fprintf(
ioQQQ,
" ChckFill middle test low fail\n" );
384 fprintf(
ioQQQ,
" ChckFill middle test high fail\n" );
393 fprintf(
ioQQQ,
" ChckFill low test low fail\n" );
399 fprintf(
ioQQQ,
" ChckFill low test high fail\n" );
409 fprintf(
ioQQQ,
" ChckFill high test low fail\n" );
418 fprintf(
ioQQQ,
" ChckFill high test high fail\n" );
466 double TeLowestFineOpacity = 1e4;
556 for( i = 1; i <=
LIMELM; i++ )
641 fprintf(
ioQQQ,
" read_continuum_mesh opening continuum_mesh.ini:");
643 ioDATA =
open_data(
"continuum_mesh.ini",
"r" );
648 fprintf(
ioQQQ,
" read_continuum_mesh could not read first line of continuum_mesh.ini.\n");
654 while(
read_whole_line( chLine , (
int)
sizeof(chLine) , ioDATA ) != NULL )
658 if( chLine[0] !=
'#')
678 if( fseek( ioDATA , 0 , SEEK_SET ) != 0 )
680 fprintf(
ioQQQ,
" read_continuum_mesh could not rewind continuum_mesh.ini.\n");
687 fprintf(
ioQQQ,
" read_continuum_mesh could not read first line of continuum_mesh.ini.\n");
698 if( ( i1 != 1 ) || ( i2 != 9 ) || ( i3 != 29 ) )
701 " read_continuum_mesh: the version of continuum_mesh.ini is not the current version.\n" );
703 " I expected to find the number 01 09 29 and got %li %li %li instead.\n" ,
705 fprintf(
ioQQQ,
"Here is the line image:\n==%s==\n", chLine );
712 while(
read_whole_line( chLine , (
int)
sizeof(chLine) , ioDATA ) != NULL )
715 if( chLine[0] !=
'#')
736 " read_continuum_mesh: A continuum resolution was zero - this is not allowed.\n" );
752 " read_continuum_mesh: The continuum definition array energies must be in increasing order.\n" );
759 " read_continuum_mesh: The last continuum array energies must be zero.\n" );