43 fprintf(
ioQQQ,
" PROBLEM DISASTER - the kinetic temperature, %.3eK,"
44 " is above the upper limit of the code, %.3eK.\n",
46 fprintf(
ioQQQ,
" This calculation is aborting.\n Sorry.\n");
54 fprintf(
ioQQQ,
" PROBLEM DISASTER - the kinetic temperature, %.3eK,"
55 " is below the lower limit of the code, %.3eK.\n",
57 fprintf(
ioQQQ,
" Consider setting a lowest temperature with the SET TEMPERATURE FLOOR command.\n");
58 fprintf(
ioQQQ,
" This calculation is aborting.\n Sorry.\n");
65 fprintf(
ioQQQ,
"temp_change: temp change floor hit, TempNew=%.3e TeFloor=%.3e, "
66 "setting constant temperature, nTotalIoniz=%li\n",
100 fprintf(
ioQQQ,
" PROBLEM DISASTER - the kinetic temperature, %.3eK,"
101 " is above the upper limit of the code, %.3eK.\n",
103 fprintf(
ioQQQ,
" This calculation is aborting.\n Sorry.\n");
111 fprintf(
ioQQQ,
" PROBLEM DISASTER - the kinetic temperature, %.3eK,"
112 " is below the lower limit of the code, %.3eK.\n",
114 fprintf(
ioQQQ,
" Consider setting a lowest temperature with the SET TEMPERATURE FLOOR command.\n");
115 fprintf(
ioQQQ,
" This calculation is aborting.\n Sorry.\n");
134 static double tgffused=-1.;
135 static double ttused = 0.;
136 static bool lgZLogSet =
false;
150 fprintf(
ioQQQ,
"tfidle called with a zero or negative electron density,%10.2e\n",
158 fprintf(
ioQQQ,
"tfidle called with a negative electron temperature,%10.2e\n",
166 for(
long nelem=0; nelem<
LIMELM; ++nelem )
195 for(
int i=1; i < 7; i++ )
359 ASSERT( massAMU < 10000. );
412 fixit(
"turbulence was included here for molecules but not ions. Now neither. Resolve.");
413 return (
realnum)sqrt(8.*BOLTZMANN/PI/ATOMIC_MASS_UNIT*
phycon.
te/massAMU);
long int ipEnergyBremsThin
NORETURN void TotalInsanity(void)
double widflx(size_t i) const
sys_float sexp(sys_float x)
bool lgTemperatureConstant
double anu(size_t i) const
void TempChange(double TempNew, bool lgForceUpdate)
STATIC void tfidle(bool lgForceUpdate)
void vexp(const double x[], double y[], long nlo, long nhi)
long int nflux_with_check
const double TEMP_LIMIT_LOW
double xIonDense[LIMELM][LIMELM+1]
bool fp_equal(sys_float x, sys_float y, int n=3)
bool lgBallistic(void) const
const double TEMP_LIMIT_HIGH
realnum GetDopplerWidth(realnum massAMU)
realnum GetAveVelocity(realnum massAMU)
double anumin(size_t i) const
#define DEBUG_ENTRY(funcname)
int fprintf(const Output &stream, const char *format,...)
double pow(double x, int i)
bool lgStatic(void) const
void vexpm1(const double x[], double y[], long nlo, long nhi)