void SetLTEQuantities(void)
{
register int n;
bool_t Debeye = TRUE;
Atom *atom;
for (n = 0; n < atmos.Natom; n++) {
atom = &atmos.atoms[n];
/* --- Get LTE populations for each atom -- -------------- */
LTEpops(atom, Debeye);
if (atom->active) {
/* --- Read the collisional data (in MULTI's GENCOL format).
After this we can close the input file for the active
atom. -- -------------- */
CollisionRate(atom, atom->fp_input);
/* --- Compute the fixed rates and store in Cij -- ------------ */
if (atom->Nfixed > 0) FixedRate(atom);
}
}
}
void SetLTEQuantities(void)
{
register int n;
bool_t Debeye = TRUE;
Atom *atom;
for (n = 0; n < atmos.Natom; n++) {
atom = &atmos.atoms[n];
/* --- Get LTE populations for each atom -- -------------- */
LTEpops(atom, Debeye);
if (atom->active) {
/* --- Read the collisional data (in MULTI's GENCOL format).
After this we can close the input file for the active
atom. -- -------------- */
// Tiago: this must work with atom as string, must use
// atom->offset_coll, and somehow read the file again! (filename gone)
CollisionRate(atom, atom->offset_coll);
/* --- Compute the fixed rates and store in Cij -- ------------ */
if (atom->Nfixed > 0) FixedRate(atom);
}
}
}
