26 bool lgHi_Down =
false;
31 bool lgLo_Down =
false;
84 trimhi =
MIN2( trimhi , 1e-20f );
96 trimhi =
MIN2( trimhi , 1e-17f );
101 trimhi =
MIN2( trimhi , 1e-12f );
113 trimhi =
MIN2(trimhi, 1e-20f);
189 "PROBLEM DISASTER the density of ion %li of element %s is too small to be computed on this cpu.\n",
193 "Turn off the element with the command ELEMENT XXX OFF or do not consider "
194 "gas with low density, the current hydrogen density is %.2e cm-3.\n",
197 "The outer radius of the cloud is %.2e cm - should a stopping "
201 "The abort flag is being set - the calculation is stopping.\n");
233 lgHi_Up_enable =
true;
236 lgHi_Up_enable =
false;
239 lgHi_Up_enable =
false;
242 lgHi_Up_enable =
false;
245 lgHi_Up_enable =
false;
260 abundnew > abundold*1.01 )
289 else if(
nzone < 10 &&
356 enum {DEBUG_LOC=
false};
357 if( DEBUG_LOC && nelem ==
ipHELIUM )
359 if( lgHi_Down ||lgHi_Up ||lgLo_Up ||lgLo_Down )
361 fprintf(
ioQQQ,
"DEBUG TrimZone\t%li\t",
nzone );
364 fprintf(
ioQQQ,
"high dn %li to %li",
370 fprintf(
ioQQQ,
"high up %li to %li",
376 fprintf(
ioQQQ,
"low up %li to %li",
382 fprintf(
ioQQQ,
"low dn %li to %li",
386 fprintf(
ioQQQ,
"\n" );
397 if( lgHi_Down || lgHi_Up || lgLo_Up || lgLo_Down )
402 enum {DEBUG_LOC=
false};
405 fprintf(
ioQQQ,
"DEBUG ion_trim zone\t%.2f \t%li\t",
fnzone, nelem);
407 fprintf(
ioQQQ,
"\tHi_Down");
409 fprintf(
ioQQQ,
"\tHi_Up");
411 fprintf(
ioQQQ,
"\tLo_Up");
413 fprintf(
ioQQQ,
"\tLo_Down");