26 fprintf(
ioQQQ,
" Too many spectra entered. Increase LIMSPC\n" );
63 fprintf(
ioQQQ,
"Too many continuum points were entered.\n" );
65 "The current logic limits the number of possible points to the value of NCELL, which is %i.\n",
NCELL );
67 "Increase the value of NCELL in rfield.h.\nSorry.\n" );
91 if( p.
strcmp(
"CONT") != 0 )
111 fprintf(
ioQQQ,
"There must be at least 2 pairs to interpolate,\nSorry\n" );
132 "An energy of zero was entered for element%3ld in INTERPOLATE and is not allowed.\nSorry\n",
141 for(
long i=0; i < npairs; i++ )
150 for(
long i=0; i < npairs; i++ )
159 for(
long i=0; i < npairs; i++ )
163 fprintf(
ioQQQ,
"An energy of zero was entered for element%3ld in INTERPOLATE and is not allowed.\nSorry\n",
173 enum {DEBUG_LOC=
false};
176 for(
long i=0; i < npairs; i++ )
178 fprintf(
ioQQQ,
"%.4e\t%.3e\n",
186 for(
long i=0; i < npairs-1; i++ )
191 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",
207 "\n NOTE The incident continuum was not defined over the entire energy range. Some energies are set to zero.\n" );
209 " NOTE You may be making a BIG mistake.\n\n" );
232 for( n=0; n < npairs; n++ )
241 for(
long i=0; i < npairs; i++ )
250 fprintf(
ioQQQ,
" Table for this continuum;\ni\tTNU\tTFAC\tTSLOP, npairs=%li\n",
252 for(
long i=0; i < npairs-1; i++ )
254 fprintf(
ioQQQ,
"%li\t%.4e\t%.4e\t%.4e\n",
258 fprintf(
ioQQQ,
"%li\t%.4e\t%.4e\n",
264 double cmin = log10( FLT_MIN );
265 double cmax = log10( FLT_MAX );
267 for(
long i=0; i < npairs; i++ )
273 " The log of the flux specified in interpolate pair %li is not within dynamic range of this CPU - please rescale.\n",i);
275 " The frequency is %f and the log of the flux is %f.\n\n",
283 fprintf(
ioQQQ,
"\n NOTE The log of the flux given in an interpolate command is outside the range of this cpu.\n");
284 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");
285 fprintf(
ioQQQ,
" NOTE Note that the interpolate command only is used for the shape of the continuum.\n");
286 fprintf(
ioQQQ,
" NOTE The order of magnitude of the flux is not used in any way.\n");
287 fprintf(
ioQQQ,
" NOTE For safety this could be of order unity.\n\n");
297 for(
long i=npairs; i <
NCELL; i++ )