/home66/gary/public_html/cloudy/c13_branch/source/iso.h File Reference

#include "two_photon.h"
#include "transition.h"
#include "freebound.h"

Include dependency graph for iso.h:

Data Structures
class	t_isoCTRL
class	extra_tr
class	t_iso_sp
Defines
#define	KILL_BELOW_PLASMA(E_) ( (rfield.lgPlasNu && ((E_)<rfield.plsfrq) ) ? 0.:1. )
#define	N_(A_) (iso_sp[ipISO][nelem].st[A_].n())
#define	L_(A_) (iso_sp[ipISO][nelem].st[A_].l())
#define	S_(A_) (iso_sp[ipISO][nelem].st[A_].S())
#define	J_(A_) (iso_sp[ipISO][nelem].st[A_].j())
#define	IPRAD 0
#define	IPCOLLIS 1
#define	RREC_MAXN 40
#define	LIKE_RREC_MAXN(A_) ( A_ == ipHELIUM ? 40 : 20 )
#define	N_ISO_TE_RECOMB 41
#define	SumUpToThisN 1000
#define	RECOMBMAGIC (110807)
Enumerations
enum	{ ipSINGLET = 1, ipDOUBLET = 2, ipTRIPLET = 3, ipMULTIPLET_END, ipMULTIPLET_BEGIN = ipSINGLET }
Functions
void	iso_cascade (long ipISO, long nelem)
void	iso_charge_transfer_update (long nelem)
void	iso_collapsed_bnl_print (long ipISO, long nelem)
void	iso_collapsed_bnl_set (long ipISO, long nelem)
void	iso_collapsed_Aul_update (long ipISO, long nelem)
void	iso_collapsed_lifetimes_update (long ipISO, long nelem)
void	iso_collide (long ipISO, long nelem)
void	iso_collisional_ionization (long ipISO, long nelem)
void	iso_continuum_lower (long ipISO, long nelem)
void	iso_cool (long ipISO, long nelem)
void	iso_create (void)
double	iso_cross_section (double ERyd, double EthRyd, long n, long l, long S, long globalZ, long globalISO)
void	iso_departure_coefficients (long ipISO, long nelem)
double	iso_dielec_recomb_rate (long ipISO, long nelem, long ipLo)
void	iso_error_generation (long ipISO, long nelem)
long	iso_get_total_num_levels (long ipISO, long nmaxResolved, long numCollapsed)
void	IonHydro ()
void	iso_ionize_recombine (long ipISO, long nelem)
void	iso_level (const long ipISO, const long nelem, double &renorm)
void	iso_photo (long ipISO, long nelem)
void	iso_prt_pops (long ipISO, long nelem, bool lgPrtDeparCoef)
void	iso_put_error (long ipISO, long nelem, long ipHi, long ipLo, long whichData, realnum errorOpt, realnum errorPess)
void	iso_radiative_recomb (long ipISO, long nelem)
void	iso_radiative_recomb_effective (long ipISO, long nelem)
double	iso_recomb_check (long ipISO, long nelem, long level, double temperature)
void	iso_recomb_auxiliary_free (void)
void	iso_recomb_malloc (void)
void	iso_recomb_setup (long ipISO)
double	iso_RRCoef_Te (long ipISO, long nelem, long n)
void	iso_satellite_update (long nelem)
double	iso_state_lifetime (long ipISO, long nelem, long n, long l)
void	iso_solve (long ipISO, long nelem, double &maxerr)
void	iso_suprathermal (long ipISO, long nelem)
void	iso_update_num_levels (long ipISO, long nelem)
void	iso_update_rates (void)
void	iso_collapsed_update (void)
void	iso_set_ion_rates (long ipISO, long nelem)
void	iso_renorm (long nelem, long ipISO, double &renorm)
Variables
long int	max_num_levels
const int	ipH1s = 0
const int	ipH2s = 1
const int	ipH2p = 2
const int	ipH3s = 3
const int	ipH3p = 4
const int	ipH3d = 5
const int	ipH4s = 6
const int	ipH4p = 7
const int	ipH4d = 8
const int	ipH4f = 9
const int	ipHe1s1S = 0
const int	ipHe2s3S = 1
const int	ipHe2s1S = 2
const int	ipHe2p3P0 = 3
const int	ipHe2p3P1 = 4
const int	ipHe2p3P2 = 5
const int	ipHe2p1P = 6
const int	ipHe3s3S = 7
const int	ipHe3s1S = 8
const int	ipHe3p3P = 9
const int	ipHe3d3D = 10
const int	ipHe3d1D = 11
const int	ipHe3p1P = 12
const int	ipH_LIKE = 0
const int	ipHE_LIKE = 1
const int	ipLI_LIKE = 2
const int	ipBE_LIKE = 3
const int	ipB_LIKE = 4
const int	ipC_LIKE = 5
const int	ipN_LIKE = 6
const int	ipO_LIKE = 7
const int	ipF_LIKE = 8
const int	ipNE_LIKE = 9
const int	ipNA_LIKE = 10
const int	ipMG_LIKE = 11
const int	ipAL_LIKE = 12
const int	ipSI_LIKE = 13
const int	ipP_LIKE = 14
const int	ipS_LIKE = 15
const int	ipCL_LIKE = 16
const int	ipAR_LIKE = 17
t_isoCTRL	iso_ctrl
t_iso_sp	iso_sp [NISO][LIMELM]

Detailed Description

information for isoelectronic sequences

Definition in file iso.h.

Define Documentation

#define IPCOLLIS 1

Definition at line 87 of file iso.h.

Referenced by iso_collide(), iso_collisional_ionization(), and iso_level().

#define IPRAD 0

Definition at line 86 of file iso.h.

Referenced by ForbiddenAuls(), he_1trans(), helike_transprob(), hydro_transprob(), iso_cascade(), iso_ionize_recombine(), iso_level(), iso_put_recomb_error(), iso_radiative_recomb(), iso_radiative_recomb_effective(), and lines_helium().

#define J_ ( A_ ) (iso_sp[ipISO][nelem].st[A_].j())

Definition at line 23 of file iso.h.

Referenced by DGEMM(), DGER(), DGETF2(), DGETRF(), and DTRSM().

#define KILL_BELOW_PLASMA ( E_ ) ( (rfield.lgPlasNu && ((E_)<rfield.plsfrq) ) ? 0.:1. )

This macro is used to zero any radiative process with photon energy below the plasma frequency. The energy must be in Rydbergs!

Definition at line 17 of file iso.h.

Referenced by iso_level().

#define L_ ( A_ ) (iso_sp[ipISO][nelem].st[A_].l())

Definition at line 21 of file iso.h.

Referenced by Badnell_rec_init(), DGEMM(), ForbiddenAuls(), helike_energy(), helike_transprob(), hydro_transprob(), HydroCSInterp(), iso_cascade(), iso_collapsed_Aul_update(), iso_collide(), iso_put_recomb_error(), iso_radiative_recomb_effective(), iso_radrecomb_from_cross_section(), lines_helium(), Opacity_iso_photo_cs(), and SanityCheckBegin().

#define LIKE_RREC_MAXN ( A_ ) ( A_ == ipHELIUM ? 40 : 20 )

Ions of the sequences will go up to this n, h-like He will get same as iso roots.

Definition at line 98 of file iso.h.

Referenced by iso_recomb_malloc(), iso_recomb_setup(), and ParseCompile().

#define N_ ( A_ ) (iso_sp[ipISO][nelem].st[A_].n())

these macros are just an easy way to return the quantum numbers of a given level.

Definition at line 20 of file iso.h.

Referenced by ForbiddenAuls(), t_ADfA::h_coll_str(), he_1trans(), helike_energy(), helike_transprob(), hydro_transprob(), HydroCSInterp(), iso_cascade(), iso_collide(), iso_collisional_ionization(), iso_create(), iso_level(), iso_radiative_recomb(), iso_radiative_recomb_effective(), iso_radrecomb_from_cross_section(), lines_hydro(), Opacity_iso_photo_cs(), RT_tau_init(), and SanityCheckBegin().

#define N_ISO_TE_RECOMB 41

Definition at line 100 of file iso.h.

Referenced by iso_recomb_malloc(), iso_recomb_setup(), and iso_RRCoef_Te().

#define RECOMBMAGIC (110807)

the magic number for the table of recombination coefficients, YYMMDD

Definition at line 106 of file iso.h.

Referenced by iso_recomb_setup().

#define RREC_MAXN 40

this is the number of levels used with the atom xx-like levels large command

Definition at line 95 of file iso.h.

Referenced by iso_recomb_malloc(), iso_recomb_setup(), ParseAtomISO(), and ParseCompile().

#define S_ ( A_ ) (iso_sp[ipISO][nelem].st[A_].S())

Definition at line 22 of file iso.h.

Referenced by Badnell_rec_init(), ForbiddenAuls(), helike_energy(), helike_transprob(), iso_cascade(), iso_put_recomb_error(), iso_radiative_recomb_effective(), iso_radrecomb_from_cross_section(), and Opacity_iso_photo_cs().

#define SumUpToThisN 1000

This is the n to go up to when calculating total recombination. Any change here will not be reflected in total recomb until "compile xxlike" is run

Definition at line 104 of file iso.h.

Referenced by iso_radiative_recomb(), and iso_recomb_setup().

Enumeration Type Documentation

anonymous enum

Enumerator:

ipSINGLET
ipDOUBLET
ipTRIPLET
ipMULTIPLET_END
ipMULTIPLET_BEGIN

Definition at line 81 of file iso.h.

Function Documentation

void IonHydro ( )

IonHydro this controls hydrogen atomic and molecular crosstalk

Definition at line 149 of file iso_solve.cpp.

References atmdat, TransitionProxy::Coll(), colliders, CollisionProxy::ColUL(), t_secondaries::csupra, DEBUG_ENTRY, dense, t_dense::eden, t_dense::EdenHCorr, t_iso_sp::fb, findspecieslocal(), fixit(), fnzone, t_hmi::H2_Solomon_dissoc_rate_BD96_H2g, t_hmi::H2_Solomon_dissoc_rate_TH85_H2g, t_hmi::H2_Solomon_dissoc_rate_used_H2g, t_hmi::H2_total, t_hydro::H_ion_frac_collis, t_hydro::H_ion_frac_photo, t_atmdat::HCharExcIonTotal, t_atmdat::HIonFrac, hmi, hydro, ion_solver(), ioQQQ, ipH1s, ipH2p, ipH_LIKE, ipHYDROGEN, iso_renorm(), iso_sp, t_hydro::lgHiPop2, t_dense::lgSetIoniz, t_trace::lgTrace, MAX2, t_hydro::pop2mx, PrintE82(), t_secondaries::sec2total, secondaries, t_dense::SetIoniz, SMALLDOUBLE, SMALLFLOAT, t_iso_sp::st, trace, t_iso_sp::trans(), and t_dense::xIonDense.

Referenced by ion_wrapper().

Here is the call graph for this function:

void iso_cascade	(	long	ipISO,
		long	nelem
	)

iso_cascade - calculate cascade probabilities, branching ratios, and associated errors

Parameters:

	ipISO
	nelem

Cascade probabilities are as defined in Robbins 68, generalized here for cascade probability for any iso sequence. >>refer He triplets Robbins, R.R. 1968, ApJ 151, 497R >>refer He triplets Robbins, R.R. 1968a, ApJ 151, 511R

Definition at line 1064 of file iso_create.cpp.

References ASSERT, EmissionProxy::Aul(), t_iso_sp::BranchRatio, t_iso_sp::CascadeProb, DEBUG_ENTRY, TransitionProxy::Emis(), ex, t_iso_sp::ex, t_iso_sp::fb, ioQQQ, ipH_LIKE, ipHe2p3P0, ipHe2p3P1, ipHe2p3P2, ipHE_LIKE, ipHELIUM, IPRAD, iso_ctrl, iso_sp, L_, t_opac::lgCaseB, t_isoCTRL::lgRandErrGen, MALLOC, N_, t_iso_sp::n_HighestResolved_max, t_iso_sp::nCollapsed_max, t_iso_sp::numLevels_max, opac, pow(), S, S_, t_isoCTRL::SmallA, and t_iso_sp::trans().

Referenced by iso_collapsed_update(), and iso_create().

Here is the call graph for this function:

void iso_charge_transfer_update ( long nelem )

iso_charge_transfer_update - update rate coefficients for CT of H and He with everything else

Definition at line 20 of file iso_ionize_recombine.cpp.

References atmdat, DEBUG_ENTRY, dense, t_atmdat::HCharExcIonOf, t_atmdat::HCharExcIonTotal, t_atmdat::HCharExcRecTo, t_atmdat::HCharExcRecTotal, t_atmdat::HeCharExcIonOf, t_atmdat::HeCharExcIonTotal, t_atmdat::HeCharExcRecTo, t_atmdat::HeCharExcRecTotal, ipHELIUM, ipHYDROGEN, ipLITHIUM, LIMELM, and t_dense::xIonDense.

Referenced by ion_solver().

void iso_collapsed_Aul_update	(	long	ipISO,
		long	nelem
	)

iso_collapsed_Aul_update - update decays from collapsed levels

Parameters:

	ipISO
	nelem

Definition at line 1723 of file iso_create.cpp.

References ASSERT, EmissionProxy::Aul(), t_iso_sp::CachedAs, DEBUG_ENTRY, TransitionProxy::Emis(), ipH_LIKE, ipHE_LIKE, iso_sp, L_, t_iso_sp::n_HighestResolved_max, t_iso_sp::nCollapsed_max, t_iso_sp::numLevels_max, t_iso_sp::QuantumNumbers2Index, t_iso_sp::st, TotalInsanity(), and t_iso_sp::trans().

Referenced by iso_collapsed_update(), and iso_zero().

Here is the call graph for this function:

void iso_collapsed_bnl_print	(	long	ipISO,
		long	nelem
	)

iso_collapsed_bnl_print - print departure coefficients for collapsed levels

Parameters:

	ipISO
	nelem

Definition at line 1677 of file iso_create.cpp.

References t_elementnames::chElementSym, t_isoCTRL::chISO, DEBUG_ENTRY, elementnames, ioQQQ, ipH_LIKE, ipHE_LIKE, iso_ctrl, iso_sp, t_iso_sp::n_HighestResolved_max, and t_iso_sp::nCollapsed_max.

void iso_collapsed_bnl_set	(	long	ipISO,
		long	nelem
	)

iso_collapsed_bnl_set - set departure coefficients for collapsed levels

Parameters:

	ipISO
	nelem

these are just sanity checks, the interpolated value should be between values at interpolation points

Definition at line 1482 of file iso_create.cpp.

References ASSERT, t_iso_sp::bnl_effective, DEBUG_ENTRY, dense, t_dense::eden, hunt_bisect(), ipH_LIKE, ipHE_LIKE, ipHELIUM, ipHYDROGEN, iso_sp, MAX2, MAX4, MIN2, MIN4, t_iso_sp::n_HighestResolved_max, t_iso_sp::nCollapsed_max, phycon, t_phycon::te, and TotalInsanity().

Referenced by iso_collapsed_update(), and iso_zero().

Here is the call graph for this function:

void iso_collapsed_lifetimes_update	(	long	ipISO,
		long	nelem
	)

iso_collapsed_lifetimes_update - update lifetimes of collapsed levels

Parameters:

	ipISO
	nelem

Definition at line 1770 of file iso_create.cpp.

References ASSERT, EmissionProxy::Aul(), EmissionProxy::dampXvel(), DEBUG_ENTRY, TransitionProxy::Emis(), TransitionProxy::EnergyWN(), iso_ctrl, iso_sp, t_iso_sp::numLevels_max, PI4, t_isoCTRL::SmallA, SMALLFLOAT, t_iso_sp::st, and t_iso_sp::trans().

Referenced by iso_collapsed_update(), and iso_zero().

Here is the call graph for this function:

void iso_collapsed_update ( void )

Definition at line 27 of file iso_solve.cpp.

References conv, dense, t_dense::IonHigh, t_dense::IonLow, ipH_LIKE, ipHELIUM, ipHYDROGEN, iso_cascade(), iso_collapsed_Aul_update(), iso_collapsed_bnl_set(), iso_collapsed_lifetimes_update(), t_dense::lgElmtOn, MIN2, NISO, and t_conv::nTotalIoniz.

Referenced by ConvBase().

Here is the call graph for this function:

void iso_collide	(	long	ipISO,
		long	nelem
	)

iso_collide - calculate collision data for ipISO, nelem

Parameters:

	ipISO
	nelem

Definition at line 120 of file iso_collide.cpp.

References t_rfield::anu, ASSERT, ATOMIC_MASS_UNIT, t_dense::AtomicWeight, BIGDOUBLE, t_isoCTRL::chISO, CollisionProxy::col_str(), TransitionProxy::Coll(), COLL_CONST, ColliderMass, colliders, t_ionbal::CollIonRate_Ground, CollisionProxy::ColUL(), conv, DEBUG_ENTRY, dense, t_dense::density_low_limit, dsexp(), t_dense::EdenHCorr, ELECTRON_MASS, TransitionProxy::Emis(), EN1RYD, TransitionProxy::EnergyWN(), t_iso_sp::fb, fixit(), fp_equal(), GetGF(), EmissionProxy::gf(), Heavy, HeCSInterp(), HION_LTE_POP, HydroCSInterp(), HydroEinstA(), ionbal, ioQQQ, ipALPHA, IPCOLLIS, ipELECTRON, ipH_LIKE, ipHE_LIKE, t_Heavy::ipHeavy, t_trace::ipIsoTrace, ipNCOLLIDER, iso_collisional_ionization(), iso_ctrl, iso_put_error(), iso_sp, iso_suprathermal(), L_, t_opac::lgCaseB_HummerStorey, t_isoCTRL::lgColl_excite, t_isoCTRL::lgColl_l_mixing, t_trace::lgHBug, t_trace::lgHeBug, t_trace::lgIsoTraceFull, t_iso_sp::lgPopLTE_OK, t_trace::lgTrace, LIMELM, MAX_DENSITY, MAX_POP_LTE, MIN2, N_, t_iso_sp::n_HighestResolved_max, NISO, t_isoCTRL::nLyaLevel, t_conv::nTotalIoniz, t_iso_sp::numLevels_local, t_iso_sp::numLevels_max, opac, phycon, pow(), PrintEfmt, CollisionProxy::rate_coef_ul(), CollisionProxy::rate_coef_ul_set(), rfield, SMALLDOUBLE, t_phycon::sqrte, t_iso_sp::st, t_isoCTRL::stat_ion, t_phycon::te, t_phycon::te32, t_phycon::te_ryd, TotalInsanity(), trace, and t_iso_sp::trans().

Referenced by iso_update_rates(), and SaveLineData().

Here is the call graph for this function:

void iso_collisional_ionization	(	long	ipISO,
		long	nelem
	)

iso_collisional_ionization - calculate collisional ionization rate for ipISO, nelem

Parameters:

	ipISO
	nelem

Definition at line 25 of file iso_collide.cpp.

References ASSERT, t_ADfA::coll_ion_wrapper(), DEBUG_ENTRY, t_iso_sp::fb, Hion_coll_ioniz_ratecoef(), hydro_vs_ioniz(), Singleton< t_ADfA >::Inst(), IPCOLLIS, iso_ctrl, iso_put_error(), iso_sp, t_isoCTRL::lgColl_ionize, t_iso_sp::lgLevelsLowered, N_, NISO, t_iso_sp::numLevels_max, phycon, and t_phycon::te.

Referenced by iso_collide().

Here is the call graph for this function:

void iso_continuum_lower	(	long	ipISO,
		long	nelem
	)

iso_continuum_lower - limit max prin. quan. no. due to continuum lowering processes

Parameters:

	ipISO
	nelem

Definition at line 12 of file iso_continuum_lower.cpp.

References ASSERT, TransitionProxy::Coll(), CollisionZero(), conv, dense, t_dense::eden, INT16_MAX, ioQQQ, ipHELIUM, ipHYDROGEN, iso_ctrl, iso_get_total_num_levels(), iso_sp, t_trace::lgHBug, t_trace::lgHeBug, t_iso_sp::lgLevelsEverLowered, t_iso_sp::lgLevelsLowered, t_iso_sp::lgMustReeval, t_trace::lgTrace, LIMELM, MAX2, MAX_DENSITY, MIN2, MIN3, t_iso_sp::n_HighestResolved_local, t_iso_sp::n_HighestResolved_max, t_iso_sp::nCollapsed_local, t_iso_sp::nCollapsed_max, t_isoCTRL::nLyman, t_isoCTRL::nLyman_malloc, t_conv::nTotalIoniz, t_iso_sp::numLevels_local, t_iso_sp::numLevels_max, phycon, pow(), t_phycon::te, trace, t_iso_sp::trans(), t_dense::xIonDense, and t_dense::xNucleiTotal.

Referenced by ContSetIntensity(), and iso_update_rates().

Here is the call graph for this function:

void iso_cool	(	long	ipISO,
		long	nelem
	)

iso_cool compute net heating/cooling due to hydrogenc atom species

Parameters:

	ipISO	the isoelectronic sequence, 0 for H
	nelem	is element, so 0 for H itself

Referenced by CoolEvaluate().

void iso_create ( void )

iso_create create storage space data for iso sequences, 1 one time per coreload

Todo:: 2 this should be changed for j-resolved levels

Definition at line 39 of file iso_create.cpp.

References abscf(), TransitionProxy::AddLine2Stack(), ASSERT, EmissionProxy::Aul(), EmissionProxy::dampXvel(), DEBUG_ENTRY, dense, DoFSMixing(), TransitionProxy::Emis(), TransitionProxy::EnergyWN(), EVRYD, t_iso_sp::ex, ExtraLymanLines, t_iso_sp::fb, FillExtraLymanLine(), g, GetGF(), EmissionProxy::gf(), HeCollidSetup(), helike_energy(), helike_transprob(), HelikeTransProbSetup(), hydro_transprob(), Singleton< t_ADfA >::Inst(), TransitionProxy::ipEmis(), ipExtraLymanLines, ipH1s, ipH2p, ipH2s, ipH3d, ipH3p, ipH3s, ipH_LIKE, ipHe1s1S, ipHe2s3S, ipHE_LIKE, ipHELIUM, ipHYDROGEN, ipLY_A, t_isoCTRL::ipLyaRedist, t_isoCTRL::ipResoRedist, t_isoCTRL::ipSubRedist, EmissionProxy::iRedisFun(), iso_allocate(), iso_assign_quantum_numbers(), iso_cascade(), iso_ctrl, iso_recomb_auxiliary_free(), iso_recomb_malloc(), iso_recomb_setup(), iso_satellite(), iso_satellite_update(), iso_sp, iso_zero(), t_dense::lgElmtOn, t_isoCTRL::lgFSM, LIMELM, N_, NISO, t_isoCTRL::nLyaLevel, t_isoCTRL::nLyman_malloc, t_iso_sp::numLevels_max, EmissionProxy::opacity(), t_ADfA::ph1(), PI4, POW2, RefIndex(), RYD_INF, t_isoCTRL::SmallA, SMALLFLOAT, t_iso_sp::st, t_isoCTRL::stat_ion, TotalInsanity(), t_iso_sp::trans(), WAVNRYD, and TransitionProxy::WLAng().

Referenced by ContCreatePointers().

Here is the call graph for this function:

double iso_cross_section	(	double	ERyd,
		double	EthRyd,
		long	n,
		long	l,
		long	S,
		long	globalZ,
		long	globalISO
	)

iso_cross_section get cross section for a particular level of an iso sequence ion

Parameters:

	ERyd
	EthRyd
	n
	l
	S
	Z
	ipISO

Definition at line 132 of file iso_radiative_recomb.cpp.

References cross_section(), DEBUG_ENTRY, H_cross_section(), He_cross_section(), ipH_LIKE, ipHE_LIKE, and TotalInsanity().

Referenced by iso_recomb_integrand().

Here is the call graph for this function:

void iso_departure_coefficients	(	long	ipISO,
		long	nelem
	)

iso_departure_coefficients - calculate departure coefficients

Parameters:

	ipISO
	nelem

Definition at line 363 of file iso_solve.cpp.

References DEBUG_ENTRY, dense, t_dense::eden, t_iso_sp::fb, iso_sp, t_iso_sp::numLevels_local, t_iso_sp::numLevels_max, safe_div(), SMALLFLOAT, and t_dense::xIonDense.

Referenced by ion_solver().

Here is the call graph for this function:

double iso_dielec_recomb_rate	(	long	ipISO,
		long	nelem,
		long	ipLo
	)

iso_dielec_recomb_rate - get state-specific dielectronic recombination rate

Parameters:

	ipISO
	nelem
	ipLo

Definition at line 1149 of file iso_radiative_recomb.cpp.

References t_phycon::alogte, ASSERT, DEBUG_ENTRY, freeBound::DielecRecombVsTemp, t_iso_sp::fb, hunt_bisect(), ipHe1s1S, ipHE_LIKE, iso_ctrl, iso_sp, t_isoCTRL::lgDielRecom, NUM_DR_TEMPS, phycon, and pow().

Referenced by iso_radiative_recomb().

Here is the call graph for this function:

void iso_error_generation	(	long	ipISO,
		long	nelem
	)

iso_error_generation generate gaussian errors

Parameters:

	ipISO
	nelem

Definition at line 39 of file iso_error.cpp.

References ASSERT, DEBUG_ENTRY, t_iso_sp::ex, ex, iso_sp, t_iso_sp::lgErrGenDone, MyGaussRand(), and t_iso_sp::numLevels_max.

Referenced by iso_update_rates().

Here is the call graph for this function:

long iso_get_total_num_levels	(	long	ipISO,
		long	nmaxResolved,
		long	numCollapsed
	)

iso_get_total_num_levels - get total number of levels with the given number of resolved and collapsed

Parameters:

	ipISO
	nmaxResolved
	numCollapsed

Definition at line 1428 of file iso_create.cpp.

References DEBUG_ENTRY, ipH_LIKE, ipHE_LIKE, and TotalInsanity().

Referenced by iso_continuum_lower(), iso_recomb_malloc(), and iso_update_num_levels().

Here is the call graph for this function:

void iso_ionize_recombine	(	long	ipISO,
		long	nelem
	)

iso_ionize_recombine evaluate state specific creation and destruction processes

Parameters:

	ipISO
	nelem

Referenced by iso_update_rates().

void iso_level	(	const long	ipISO,
		const long	nelem,
		double &	renorm
	)

iso_level solve for iso-sequence ionization balance

Parameters:

	ipISO
	nelem

Referenced by iso_solve().

void iso_photo	(	long	ipISO,
		long	nelem
	)

iso_photo do photoionization rates for element nelem on the ipISO isoelectronic sequence

Parameters:

	ipISO
	nelem

Referenced by iso_update_rates().

void iso_prt_pops	(	long	ipISO,
		long	nelem,
		bool	lgPrtDeparCoef
	)

iso_prt_pops routine to print level pops or departure coefficients for iso sequences

Parameters:

	ipISO
	nelem
	lgPrtDeparCoef

Definition at line 386 of file iso_solve.cpp.

References ASSERT, t_elementnames::chElementSym, t_isoCTRL::chISO, DEBUG_ENTRY, elementnames, ioQQQ, ipH_LIKE, ipHe2p3P0, ipHe2p3P1, ipHe2p3P2, ipHE_LIKE, iso_ctrl, iso_sp, ITEM_TO_PRINT, t_iso_sp::nCollapsed_local, NISO, t_iso_sp::numLevels_local, t_iso_sp::QuantumNumbers2Index, and t_iso_sp::st.

Referenced by PrtZone().

void iso_put_error	(	long	ipISO,
		long	nelem,
		long	ipHi,
		long	ipLo,
		long	whichData,
		realnum	errorOpt,
		realnum	errorPess
	)

iso_put_error put an error bar on a piece of data, to be used with Gaussian random noise gen

Parameters:

	ipISO
	nelem
	ipHi
	ipLo
	whichData
	errorOpt
	errorPess

Referenced by ForbiddenAuls(), he_1trans(), helike_transprob(), hydro_transprob(), iso_collide(), iso_collisional_ionization(), iso_put_recomb_error(), and iso_radiative_recomb().

void iso_radiative_recomb	(	long	ipISO,
		long	nelem
	)

iso_radiative_recomb - get rad recomb rate coefficients for iso sequences.

Parameters:

	ipISO
	nelem

Definition at line 149 of file iso_radiative_recomb.cpp.

References ASSERT, cdEXIT, t_elementnames::chElementSym, t_isoCTRL::chISO, conv, cosmology, DEBUG_ENTRY, t_ionbal::DR_Badnell_rate_coef, elementnames, EXIT_FAILURE, t_iso_sp::fb, fp_equal(), t_ADfA::H_rad_rec(), t_iso_sp::HighestLevelOpacStack, Singleton< t_ADfA >::Inst(), ionbal, ioQQQ, t_save::ioRecom, ipH1s, ipHELIUM, t_trace::ipIsoTrace, IPRAD, ipRecEsc, ipRecNetEsc, ipRecRad, iso_ctrl, iso_dielec_recomb_rate(), iso_put_error(), iso_put_recomb_error(), iso_radrecomb_from_cross_section(), iso_RRCoef_Te(), iso_sp, t_opac::lgCaseB, t_isoCTRL::lgDielRecom, t_cosmology::lgDo, t_trace::lgHBug, t_trace::lgHeBug, t_save::lgioRecom, t_trace::lgIsoTraceFull, t_iso_sp::lgMustReeval, t_isoCTRL::lgNoRecombInterp, t_isoCTRL::lgRandErrGen, t_isoCTRL::lgTopoff, t_trace::lgTrace, LIMELM, MAX2, MIN2, N_, t_iso_sp::n_HighestResolved_local, t_iso_sp::nCollapsed_local, NHYDRO_MAX_LEVEL, NISO, t_conv::nTotalIoniz, t_iso_sp::numLevels_local, t_iso_sp::numLevels_max, opac, t_opac::OpacStack, phycon, PrintEfmt, t_iso_sp::RadRec_caseB, t_iso_sp::RadRec_effec, Recomb_Seaton59(), RT_recom_effic(), save, ShowMe(), SumUpToThisN, t_opac::TauAbsGeo, t_phycon::te, and trace.

Referenced by iso_update_rates().

Here is the call graph for this function:

void iso_radiative_recomb_effective	(	long	ipISO,
		long	nelem
	)

iso_radiative_recomb_effective - get effective recomb rate coefficients into each level (including indirect)

Parameters:

	ipISO
	nelem

Definition at line 602 of file iso_radiative_recomb.cpp.

References ASSERT, t_iso_sp::BranchRatio, t_iso_sp::CascadeProb, DEBUG_ENTRY, dprintf(), EN1RYD, ex, t_iso_sp::fb, ioQQQ, IPRAD, ipRecRad, iso_ctrl, iso_sp, L_, t_isoCTRL::lgRandErrGen, MAX2, N_, t_iso_sp::numLevels_local, phycon, pow(), prt_wl(), RYDLAM, S_, t_phycon::te, and wavelength.

Referenced by iso_update_rates().

Here is the call graph for this function:

void iso_recomb_auxiliary_free ( void )

iso_recomb_auxiliary_free - free up some auxiliary space associated with iso recombination tables.

Definition at line 825 of file iso_radiative_recomb.cpp.

References DEBUG_ENTRY, NISO, and NumLevRecomb.

Referenced by iso_create().

double iso_recomb_check	(	long	ipISO,
		long	nelem,
		long	level,
		double	temperature
	)

iso_recomb_check - called by SanityCheck to confirm that recombination coef are ok, return value is relative error between new calculation of recom, and interp value

Parameters: