98 abundan =
dense.xIonDense[ (*t10.
Hi()).nelem() -1][(*t10.
Hi()).IonStg()-1];
104 (*t10.
Lo()).Pop() = 0.;
105 (*t20.
Lo()).Pop() = 0.;
106 (*t30.
Lo()).Pop() = 0.;
107 (*t21.
Lo()).Pop() = 0.;
108 (*t31.
Lo()).Pop() = 0.;
109 (*t41.
Lo()).Pop() = 0.;
118 (*t10.
Hi()).Pop() = 0.;
119 (*t20.
Hi()).Pop() = 0.;
120 (*t30.
Hi()).Pop() = 0.;
121 (*t21.
Hi()).Pop() = 0.;
122 (*t31.
Hi()).Pop() = 0.;
123 (*t41.
Hi()).Pop() = 0.;
172 atoms.PopLevels[i] = 0.;
173 atoms.DepLTELevels[i] = 1.;
205 stat[0] = (*t10.
Lo()).g();
206 stat[1] = (*t10.
Hi()).g();
207 stat[2] = (*t20.
Hi()).g();
208 stat[3] = (*t30.
Hi()).g();
209 stat[4] = (*t41.
Hi()).g();
210 ASSERT( stat[0]>0. && stat[1]>0. &&stat[2]>0. &&stat[3]>0. &&stat[4]>0.);
211 ASSERT( fabs((*t10.
Lo()).g()/2.-1.) < FLT_EPSILON);
212 ASSERT( fabs((*t10.
Hi()).g()/4.-1.) < FLT_EPSILON);
213 ASSERT( fabs((*t20.
Lo()).g()/2.-1.) < FLT_EPSILON);
214 ASSERT( fabs((*t20.
Hi()).g()/2.-1.) < FLT_EPSILON);
215 ASSERT( fabs((*t30.
Lo()).g()/2.-1.) < FLT_EPSILON);
216 ASSERT( fabs((*t30.
Hi()).g()/4.-1.) < FLT_EPSILON);
217 ASSERT( fabs((*t21.
Lo()).g()/4.-1.) < FLT_EPSILON);
218 ASSERT( fabs((*t21.
Hi()).g()/2.-1.) < FLT_EPSILON);
219 ASSERT( fabs((*t31.
Lo()).g()/4.-1.) < FLT_EPSILON);
220 ASSERT( fabs((*t31.
Hi()).g()/4.-1.) < FLT_EPSILON);
221 ASSERT( fabs((*t41.
Lo()).g()/4.-1.) < FLT_EPSILON);
222 ASSERT( fabs((*t41.
Hi()).g()/6.-1.) < FLT_EPSILON);
230 ASSERT( excit[1]>0. &&excit[2]>0. &&excit[3]>0. &&excit[4]>0.);
320 (*t10.
Lo()).Pop() =
pops[0];
321 (*t20.
Lo()).Pop() =
pops[0];
322 (*t30.
Lo()).Pop() =
pops[0];
323 (*t21.
Lo()).Pop() =
pops[1];
324 (*t31.
Lo()).Pop() =
pops[1];
325 (*t41.
Lo()).Pop() =
pops[1];
334 (*t10.
Hi()).Pop() =
pops[1];
335 (*t20.
Hi()).Pop() =
pops[2];
336 (*t30.
Hi()).Pop() =
pops[3];
337 (*t21.
Hi()).Pop() =
pops[2];
338 (*t31.
Hi()).Pop() =
pops[3];
339 (*t41.
Hi()).Pop() =
pops[4];
void atom_levelN(long int nLevelCalled, realnum abund, const double g[], const double ex[], char chExUnits, double pops[], double depart[], double ***AulEscp, double ***col_str, double ***AulDest, double ***AulPump, double ***CollRate, const double source[], const double sink[], bool lgCollRateDone, double *cooltl, double *coolder, const char *chLabel, int *nNegPop, bool *lgZeroPop, bool lgDeBug, bool lgLTE, multi_arr< double, 2 > *Cool, multi_arr< double, 2 > *dCooldT)
void AtomSeqBoron(const TransitionProxy &t10, const TransitionProxy &t20, const TransitionProxy &t30, const TransitionProxy &t21, const TransitionProxy &t31, const TransitionProxy &t41, double cs40, double cs32, double cs42, double cs43, double pump_rate, const char *chLabel)