70 for( i=0; i < ip; i++ )
104 fprintf( io,
"#>>>> Temperature not converged.\n" );
108 fprintf( io,
"#>>>> Electron density not converged.\n" );
112 fprintf( io,
"#>>>> Ionization not converged.\n" );
116 fprintf( io,
"#>>>> Pressure not converged.\n" );
119 if( strcmp(chJob,
"EACH") == 0 )
122 fprintf( io,
"%.5e\t%.4e\t%.4e",
126 double debug_ctot = 0.;
128 for(
int i=0 ; i <
LIMELM ; i++)
221 enum{ DEBUG_COOLING =
false };
227 "fractional diff: %.4e\n",
230 debug_ctot / cool_total - 1. );
235 if( fabs( (debug_ctot - cool_total)/cool_total ) > 1e-10 )
237 fprintf(
ioQQQ ,
"PROBLEM with the SAVE EACH COOLING output\n" );
238 fprintf(
ioQQQ ,
"PROBLEM One or more coolants have been lost, the sum of the reported cooling is %.4e\n", debug_ctot );
239 fprintf(
ioQQQ ,
"PROBLEM The total cooling is %.4e\n", cool_total );
240 fprintf(
ioQQQ ,
"PROBLEM The difference is %.4e\n", cool_total - debug_ctot );
244 else if( strcmp(chJob,
"COOL") == 0 )
249 fprintf( io,
"%.5e\t%.4e\t%.4e\t%.4e",
261 for( is=0; is < ip; is++ )
265 fprintf( io,
"\t%s %.1f\t%.7f",
268 sign(csav[index[is]],sgnsav[index[is]]) );
void CoolSave(FILE *io, const char chJob[])
double heavycollcool[LIMELM]
NORETURN void TotalInsanity(void)
char chClntLab[NCOLNT][NCOLNT_LAB_LEN+1]
sys_float safe_div(sys_float x, sys_float y, sys_float res_0by0)
#define DEBUG_ENTRY(funcname)
double elementcool[LIMELM+1]
int fprintf(const Output &stream, const char *format,...)
void spsort(realnum x[], long int n, long int iperm[], int kflag, int *ier)