38 fprintf(
ioQQQ,
" Too many spectra entered. Increase LIMSPC\n" );
60 while( !lgDONE && !*lgEOF )
69 while( !lgEOL && npairs<
NCELL )
80 fprintf(
ioQQQ,
"Too many continuum points were entered.\n" );
82 "The current logic limits the number of possible points to the value of NCELL, which is %i.\n",
NCELL );
84 "Increase the value of NCELL in rfield.h.\nSorry.\n" );
102 cap4( chLab4 , chCard );
107 fprintf(
ioQQQ,
" * ");
110 while( chCard[k]!=
'\0' )
112 fprintf(
ioQQQ,
"%c",chCard[k]);
117 fprintf(
ioQQQ,
"%c",
' ');
120 fprintf(
ioQQQ,
"*\n");
127 if( strncmp(chCard,
"CONT",4) != 0 )
134 if( chCard[0] ==
' ' )
148 fprintf(
ioQQQ,
"There must be at least 2 pairs to interpolate,\nSorry\n" );
154 else if( npairs >=
NCELL - 2 )
156 fprintf(
ioQQQ,
" Too many continuum points entered.\n" );
158 "The current logic limits the number of possible points to the value of NCELL, which is %i.\nSorry.\n",
NCELL );
159 fprintf(
ioQQQ,
" Increase NCELL (in cddefines.h) to more than the number of continuum points.\n" );
180 "An energy of zero was entered for element%3ld in INTERPOLATE and is not allowed.\nSorry\n",
189 for( i=0; i < (npairs + 1); i++ )
199 for( i=0; i < (npairs + 1); i++ )
208 for( i=0; i < (npairs + 1); i++ )
212 fprintf(
ioQQQ,
"An energy of zero was entered for element%3ld in INTERPOLATE and is not allowed.\nSorry\n",
222 enum {DEBUG_LOC=
false};
225 for( i=0; i < npairs; i++ )
227 fprintf(
ioQQQ,
"%.4e\t%.3e\n",
235 for( i=0; i < npairs; i++ )
240 fprintf(
ioQQQ,
"The energies MUST be in increasing order. Energy #%3ld=%10.2e Ryd was greater than or equal to the next one.\nSorry.\n",
255 if( npairs + 2 <
NCELL )
260 for( i=npairs + 1; i <
NCELL; i++ )
271 "\n NOTE The incident continuum was not defined over the entire energy range. Some energies are set to zero.\n" );
273 " NOTE You may be making a BIG mistake.\n\n" );
297 for( i=0; i <= npairs; i++ )
304 if( cmax - cmin > 74. )
306 fprintf(
ioQQQ,
"The dynamic range of the specified continuum is too large.\nSorry.\n" );
314 for( i=0; i <= npairs; i++ )
320 else if( cmax > 37. )
324 for( i=0; i <= npairs; i++ )
333 fprintf(
ioQQQ,
" Table for this continuum;\ni\tTNU\tTFAC\tTSLOP, npairs=%li\n",
335 for( i=0; i < npairs; i++ )
337 fprintf(
ioQQQ,
"%li\t%.4e\t%.4e\t%.4e\n",
342 fprintf(
ioQQQ,
"%li\t%.4e\t%.4e\n",
361 fprintf(
ioQQQ,
"ParseInput finds insane i after rfield.tNuRyd loop\n");
368 for( i=0; i <= npairs; i++ )
375 cmin = log10( FLT_MIN );
376 cmax = log10( FLT_MAX );
378 for( i=0; i <= npairs; i++ )
384 " The log of the flux specified in interpolate pair %li is not within dynamic range of this CPU - please rescale.\n",i);
386 " The frequency is %f and the log of the flux is %f.\n\n",
394 fprintf(
ioQQQ,
"\n NOTE The log of the flux given in an interpolate command is outside the range of this cpu.\n");
395 fprintf(
ioQQQ,
" NOTE I will try to renormalize it to be within the range of this cpu, but if I crash, this is a likely reason.\n");
396 fprintf(
ioQQQ,
" NOTE Note that the interpolate command only is used for the shape of the continuum.\n");
397 fprintf(
ioQQQ,
" NOTE The order of magnitude of the flux is not used in any way.\n");
398 fprintf(
ioQQQ,
" NOTE For safety this could be of order unity.\n\n");
416 fprintf(
ioQQQ,
" Too many continua entered; increase LIMSPC\n" );