cloudy  trunk
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Data Structures | Macros | Functions | Variables
atoms.h File Reference
This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Data Structures

struct  t_atoms

Macros

#define N_OI_LEVELS   6
#define LIMLEVELN   20L

Functions

void AtomSeqBeryllium (double cs12, double cs13, double cs23, transition *t, double a30)
void AtomSeqBoron (transition *t21, transition *t31, transition *t41, transition *t32, transition *t42, transition *t52, double cs51, double cs43, double cs53, double cd54, double pump_rate, const char *chLabel)
void atom_level2 (transition *t)
void atom_level3 (transition *t10, transition *t21, transition *t20)
double atom_pop2 (double omega, double g1, double g2, double a21, double bltz, double abund)
double atom_pop3 (double g1, double g2, double g3, double o12, double o13, double o23, double a21, double a31, double a32, double Tex12, double Tex23, realnum *pop2, double abund, double gam2, double r12, double r13)
void atom_pop5 (double g[], double ex[], double cs12, double cs13, double cs14, double cs15, double cs23, double cs24, double cs25, double cs34, double cs35, double cs45, double a21, double a31, double a41, double a51, double a32, double a42, double a52, double a43, double a53, double a54, double p[], realnum abund)
void atom_levelN (long int nlev, 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 create[], const double destroy[], bool lgCollRateDone, double *cooltl, double *coolder, const char *chLabel, int *nNegPop, bool *lgZeroPop, bool lgDeBug)
void atom_oi_calc (double *coloi)

Variables

EXTERN struct t_atoms atoms

Macro Definition Documentation

#define LIMLEVELN   20L
#define N_OI_LEVELS   6

number of levels in OI atom

Definition at line 221 of file atoms.h.

Referenced by IterStart(), and zero().


Function Documentation

void atom_level2 ( transition t)
void atom_level3 ( transition t10,
transition t21,
transition t20 
)
void atom_levelN ( long int  nlev,
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  create[],
const double  destroy[],
bool  lgCollRateDone,
double *  cooltl,
double *  coolder,
const char *  chLabel,
int *  nNegPop,
bool *  lgZeroPop,
bool  lgDeBug 
)

atom_levelN - compute populations of arbitrary n-level atom

Parameters:
nlevnlev is the number of levels to compute
abundABUND is total abundance of species, used for nth equation
g[]G(ndim) is stat weight of levels
ex[]EX(ndim) is excitation potential of levels, either wn or deg K 0 for first one, NOT d(ENER), but energy rel to ground
chExUnitsthis is 'K' for above ex[] as Kelvin deg, is 'w' for wavenumbers
pops[]populations of each level as deduced here
depart[]departure coefficient derived here
AulEscpnet transition rate, A * esc prob, s-1
col_strcol str rom up to low
AulDestAulDest(ilo,ihi) is destruction rate, from up to low, A * dest prob, [s-1], asserts confirm that ihi,lo is zero
AulPumpAulPump(lo, hi) is pumping rate, A * occ num, (hi,lo) must be zero, [s-1]
CollRatecollision rates, evaluated here and returned for cooling by calling function, unless following flag is true. If true then calling function has already filled in these rates. CollRate[i][j] is rate from i to j
createthis is an additional creation rate, normally zero, units cm-3 s-1
destroy[]this is an additional destruction rate to continuum, normally zero, units s-1
lgCollRateDoneflag saying whether CollRate already done, or we need to do it here
cooltltotal cooling, set here but nothing done with it
coolderderivative of cooling, set here but nothing done with it
chLabelstring used to identify calling program in case of error
lgNegPoplgNegPop flag indicating what we have done positive if negative populations occurred zero if normal calculation done negative if too cold (for some atoms other routine will be called in this case)
lgZeroPoptrue if populations are zero, either due to zero abundance of very low temperature
lgDeBugoption to print matrices for debugging
Postcondition:
atoms.PopLevels[n], atoms.DepLTELevels[n] are set lines added to ots array

Definition at line 13 of file atom_leveln.cpp.

References abund, amat, ASSERT, BOLTZMANN, cdEXIT, t_dense::cdsqte, auto_vec< T >::data(), multi_arr< T, d, ALLOC, lgBC >::data(), DEBUG_ENTRY, dense, getrf_wrapper(), getrs_wrapper(), t_thermal::halfte, ioQQQ, ipiv, t_trace::lgTrace, t_trace::lgTrLevN, MAX2, phycon, sexp(), T1CM, t_phycon::te, t_phycon::te_wn, thermal, TotalInsanity(), trace, t_thermal::tsq1, and multi_arr< T, d, ALLOC, lgBC >::zero().

Referenced by atmol_popsolve(), AtomSeqBoron(), CO_PopsEmisCool(), Fe11Lev5(), Fe13Lev5(), Fe2_cooling(), Fe3Lev14(), Fe4Lev12(), Fe7Lev8(), and H2_Level_low_matrix().

Here is the call graph for this function:

void atom_oi_calc ( double *  coloi)

atom_oi drive the solution of OI level populations, Ly-beta pumping

Parameters:
coloi
Todo:
2 this is not used, should it be?
Todo:
2 following needed to get badbugs/bug8.in to work
Todo:
2 following needed to get badbugs/bug5.in to work

Definition at line 23 of file atom_oi.cpp.

References atoms, t_emission::Aul, DEBUG_ENTRY, dense, t_DoppVel::doppler, DoppVel, t_transition::Emis, esc_CRDwing_1side(), t_atoms::esch31, t_transition::Hi, ioQQQ, ipH1s, ipH3p, ipH_LIKE, ipHYDROGEN, ipOXYGEN, ipT1039, ipT1304, ipT4368, ipT8446, ipTO1025, ipTOI11, ipTOI13, ipTOI29, ipTOI46, iteration, t_rfield::lgInducProcess, t_trace::lgTr8446, t_trace::lgTrace, t_transition::Lo, MAX2, MIN2, oi_level_pops(), t_emission::Pelec_esc, t_emission::Pesc, t_atoms::pmph31, t_atoms::pmpo15, t_atoms::pmpo51, t_quantumState::Pop, t_atoms::popoi, t_emission::PopOpc, rfield, StatesElem, t_emission::TauIn, TauLines, t_emission::TauTot, trace, Transitions, and t_dense::xIonDense.

Referenced by CoolOxyg(), and RT_line_all().

Here is the call graph for this function:

double atom_pop2 ( double  omega,
double  g1,
double  g2,
double  a21,
double  bltz,
double  abund 
)

atom_pop2 do level population for simple two level atom, no radiative transfer

Parameters:
omega
g1
g2
a21
bltz
abund

Definition at line 9 of file atom_pop2.cpp.

References ASSERT, t_dense::cdsqte, DEBUG_ENTRY, dense, phycon, and t_phycon::teinv.

Referenced by CoolArgo(), CoolChro(), CoolIron(), CoolMagn(), CoolNeon(), CoolScan(), CoolSili(), CoolSulf(), and IonOxyge().

double atom_pop3 ( double  g1,
double  g2,
double  g3,
double  o12,
double  o13,
double  o23,
double  a21,
double  a31,
double  a32,
double  Tex12,
double  Tex23,
realnum pop2,
double  abund,
double  gam2,
double  r12,
double  r13 
)

atom_pop3 return value is population for 3-level atom, cm^-3

Parameters:
g1statictical weights of level 1
g2statictical weights of level 2
g3statictical weights of level 3
o12collision strengths between three levels
o13collision strengths between three levels
o23collision strengths between three levels
a21transition probabilities between three levels
a31transition probabilities between three levels
a32transition probabilities between three levels
Tex12excitation energy in Kelvin
Tex23excitation energy in Kelvin
*pop2returned population of level 2, cm^-3
abundincoming total abundance of ion
gam2possible photodestruction of level 2, normally 0
r12excitation rates (s-1) due to "other" processes
r13excitation rates (s-1) due to "other" processes

Definition at line 10 of file atom_pop3.cpp.

References ASSERT, atoms, t_atoms::c12, t_atoms::c13, t_dense::cdsqte, DEBUG_ENTRY, dense, phycon, SMALLFLOAT, and t_phycon::te.

Referenced by CoolAlum(), CoolArgo(), CoolCalc(), CoolCarb(), CoolChlo(), CoolChro(), CoolIron(), CoolNeon(), CoolNitr(), CoolOxyg(), CoolPhos(), CoolScan(), CoolSili(), CoolSodi(), CoolSulf(), CoolTita(), and CoolVana().

void atom_pop5 ( double  g[],
double  ex[],
double  cs12,
double  cs13,
double  cs14,
double  cs15,
double  cs23,
double  cs24,
double  cs25,
double  cs34,
double  cs35,
double  cs45,
double  a21,
double  a31,
double  a41,
double  a51,
double  a32,
double  a42,
double  a52,
double  a43,
double  a53,
double  a54,
double  p[],
realnum  abund 
)

atom_pop5 do populations and cooling for five level atom

Parameters:
g[]
ex[]
cs12
cs13
cs14
cs15
cs23
cs24
cs25
cs34
cs35
cs45
a21
a31
a41
a51
a32
a42
a52
a43
a53
a54
p[]
abund
Todo:
2 p(5) was very slightly negative (1e-40) for SII in dqher.in

Definition at line 12 of file atom_pop5.cpp.

References abund, amat, cdEXIT, t_dense::cdsqte, DEBUG_ENTRY, dense, getrf_wrapper(), getrs_wrapper(), ioQQQ, ipiv, MAX2, phycon, sexp(), T1CM, and t_phycon::te.

Referenced by CoolArgo(), CoolCalc(), CoolChlo(), CoolNitr(), CoolOxyg(), and CoolSulf().

Here is the call graph for this function:

void AtomSeqBeryllium ( double  cs12,
double  cs13,
double  cs23,
transition t,
double  a30 
)
void AtomSeqBoron ( transition t21,
transition t31,
transition t41,
transition t32,
transition t42,
transition t52,
double  cs51,
double  cs43,
double  cs53,
double  cd54,
double  pump_rate,
const char *  chLabel 
)

Variable Documentation

EXTERN struct t_atoms atoms

Generated for cloudy by doxygen 1.8.1.1