Go to the documentation of this file.
22 const char chDetail[] )
40 if( strcmp( chMode ,
"pres" )==0 )
47 " PROBLEM ConvFail %li, pressure not converged; itr %li, zone %.2f Te:%.3e Hden:%.4e curr Pres:%.4e Corr Pres:%.4e Pra/gas:%.3e\n",
62 "\n PROBLEM continued, pressure not converged; we are stuck at the sonic point.\n\n");
69 else if( strcmp( chMode,
"eden" ) == 0 )
77 " PROBLEM ConvFail %li, eden not converged itr %li zone %li fnzone %.2f correct=%.3e "
91 fprintf(
ioQQQ,
" Temperature failure also." );
97 fprintf(
ioQQQ,
" Ionization failure also." );
103 fprintf(
ioQQQ,
" Electron density oscillating." );
109 fprintf(
ioQQQ,
" Heating-cooling oscillating." );
112 fprintf(
ioQQQ,
" \n");
115 else if( strcmp( chMode,
"ioni" ) == 0 )
121 fprintf(
ioQQQ,
" PROBLEM ConvFail %li, %s ionization not converged iteration %li zone %li fnzone %.2f \n",
130 else if( strcmp( chMode,
"pops" ) == 0 )
137 fprintf(
ioQQQ,
" PROBLEM ConvFail %li, %s population not converged iteration %li zone %li fnzone %.2f \n",
146 else if( strcmp( chMode,
"grai" ) == 0 )
152 fprintf(
ioQQQ,
" PROBLEM ConvFail %ld, a grain failure occurred iteration %li zone %li fnzone %.2f \n",
161 else if( strcmp( chMode,
"temp" ) == 0 )
168 " PROBLEM ConvFail %ld, Temp not converged itr %li zone %li fnzone %.2f Te=%.4e DTe=%.2e"
169 " Htot=%.3e Ctot=%.3e rel err=%.3e rel tol:%.3e\n",
183 fprintf(
ioQQQ,
" Temperature and d(Cool-Heat)/dT were both oscillating.\n" );
188 fprintf(
ioQQQ,
" Temperature was oscillating.\n" );
193 fprintf(
ioQQQ,
" d(Cool-Heat)/dT was oscillating.\n" );
199 fprintf(
ioQQQ,
" Solution not converged due to %10.10s\n",
206 fprintf(
ioQQQ,
" ConvFail called with insane mode %s detail %s\n",
232 fprintf(
ioQQQ,
" Stop due to excessive convergence failures - there have been %ld so far. \n",
234 fprintf(
ioQQQ,
" This limit can be reset with the FAILURES command.\n" );
240 fprintf(
ioQQQ,
"\n This problem may be solved by adding cosmic rays.\n");
241 fprintf(
ioQQQ,
"\n The gas was cold and neutral.\n");
242 fprintf(
ioQQQ,
"\n The chemistry is not designed to work without a source of ionization.\n");
243 fprintf(
ioQQQ,
" >>> Add galactic background cosmic rays with the COSMIC RAYS BACKBOUND command and try again.\n\n" );
249 fprintf(
ioQQQ,
" These were all pressure failures - we may be near an unstable point in the cooling curve. \n");
250 fprintf(
ioQQQ,
" The PUNCH PRESSURE HISTORY command will show the n-T-P curve, and may be interesting.\n\n");
273 fprintf(
ioQQQ,
" ConvFail sets lgAbort since nTotalFailures=%ld is >= LimFail=%ld\n",
276 fprintf(
ioQQQ,
" This limit can be reset with the FAILURES command.\n");
Generated for cloudy by
1.8.3.1