/home66/gary/public_html/cloudy/c10_branch/source/iso_radiative_recomb.cpp File Reference

#include "cddefines.h"
#include "dynamics.h"
#include "atmdat.h"
#include "conv.h"
#include "cosmology.h"
#include "dense.h"
#include "elementnames.h"
#include "helike.h"
#include "helike_recom.h"
#include "hydrogenic.h"
#include "ionbal.h"
#include "iso.h"
#include "opacity.h"
#include "phycon.h"
#include "physconst.h"
#include "prt.h"
#include "save.h"
#include "thermal.h"
#include "thirdparty.h"
#include "trace.h"
#include "rt.h"
#include "taulines.h"

Include dependency graph for iso_radiative_recomb.cpp:

Go to the source code of this file.

Functions
STATIC double	TempInterp (double TempArray, double ValueArray, long NumElements)
STATIC double	iso_recomb_integrand (double EE)
STATIC void	iso_put_recomb_error (long ipISO, long nelem)
double	iso_radrecomb_from_cross_section (long ipISO, double temp, long nelem, long ipLo)
double	iso_cross_section (double EgammaRyd, double EthRyd, long n, long l, long S, long globalZ, long globalISO)
void	iso_radiative_recomb (long ipISO, long nelem)
void	iso_radiative_recomb_effective (long ipISO, long nelem)
double	iso_RRCoef_Te (long ipISO, long nelem, long n)
double	iso_recomb_check (long ipISO, long nelem, long level, double temperature)
void	iso_recomb_malloc (void)
void	iso_recomb_auxiliary_free (void)
void	iso_recomb_setup (long ipISO)
double	iso_dielec_recomb_rate (long ipISO, long nelem, long ipLo)
Variables
static double ****	RRCoef
static long **	NumLevRecomb
static double ***	TotalRecomb
static double	TeRRCoef [N_ISO_TE_RECOMB]
static double	kTRyd
static double	global_EthRyd
static long int	ipLev
static long int	globalZ
static long int	globalISO
static long int	globalN
static long int	globalL
static long int	globalS

Function Documentation

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 133 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:

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 1125 of file iso_radiative_recomb.cpp.

References t_phycon::alogte, ASSERT, DEBUG_ENTRY, t_iso::DielecRecombVsTemp, hunt_bisect(), ipHe1s1S, ipHE_LIKE, iso, t_iso::lgDielRecom, NUM_DR_TEMPS, t_iso::numLevels_max, phycon, and pow().

Referenced by iso_radiative_recomb().

Here is the call graph for this function:

STATIC void iso_put_recomb_error	(	long	ipISO,
		long	nelem
	)

Definition at line 536 of file iso_radiative_recomb.cpp.

References ASSERT, ipHE_LIKE, ipHELIUM, IPRAD, iso, iso_put_error(), L_, MIN2, t_iso::nCollapsed_max, t_iso::numLevels_max, t_iso::QuantumNumbers2Index, and S_.

Referenced by iso_radiative_recomb().

Here is the call graph for this function:

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 150 of file iso_radiative_recomb.cpp.

References ASSERT, cdEXIT, t_elementnames::chElementSym, t_iso::chISO, t_iso::ConOpacRatio, conv, cosmology, DEBUG_ENTRY, t_iso::DielecRecomb, t_ionbal::DR_Badnell_rate_coef, dynamics, elementnames, t_iso::Error, t_iso::ErrorFactor, fp_equal(), fudge(), Singleton< T >::Inst(), ionbal, ioQQQ, t_save::ioRecom, ipH1s, ipHELIUM, t_iso::ipIsoLevNIonCon, t_trace::ipIsoTrace, IPRAD, ipRecEsc, ipRecNetEsc, ipRecRad, iso, iso_dielec_recomb_rate(), iso_put_error(), iso_put_recomb_error(), iso_radrecomb_from_cross_section(), iso_RRCoef_Te(), iteration, t_opac::lgCaseB, t_iso::lgDielRecom, t_cosmology::lgDo, t_trace::lgHBug, t_trace::lgHeBug, t_save::lgioRecom, t_trace::lgIsoTraceFull, t_iso::lgLevelsLowered, t_iso::lgMustReeval, t_iso::lgNoRecombInterp, t_iso::lgRandErrGen, t_iso::lgTopoff, t_trace::lgTrace, LIMELM, MAX2, MIN2, MyGaussRand(), N_, t_iso::n_HighestResolved_local, t_dynamics::n_initial_relax, t_iso::nCollapsed_local, t_iso::nCollapsed_max, NHYDRO_MAX_LEVEL, NISO, t_conv::nTotalIoniz, t_iso::numLevels_local, t_iso::numLevels_max, opac, phycon, PrintEfmt, t_iso::RadRec_caseB, t_iso::RadRec_effec, t_iso::RadRecomb, 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 578 of file iso_radiative_recomb.cpp.

References ASSERT, t_iso::BranchRatio, t_iso::CascadeProb, DEBUG_ENTRY, dprintf(), EN1RYD, t_iso::Error, ioQQQ, IPRAD, ipRecRad, iso, L_, t_iso::lgRandErrGen, MAX2, N_, t_iso::numLevels_local, t_iso::numLevels_max, phycon, pow(), prt_wl(), t_iso::RadEffec, t_iso::RadRecomb, RYDLAM, S_, t_iso::SigmaAtot, t_iso::SigmaCascadeProb, t_iso::SigmaRadEffec, StatesElemNEW, t_phycon::te, wavelength, and t_iso::xIsoLevNIonRyd.

Referenced by iso_update_rates().

Here is the call graph for this function:

double iso_radrecomb_from_cross_section	(	long	ipISO,
		double	temp,
		long	nelem,
		long	ipLo
	)

Definition at line 47 of file iso_radiative_recomb.cpp.

References DEBUG_ENTRY, global_EthRyd, globalISO, globalL, globalN, globalS, globalZ, Singleton< T >::Inst(), ipH_LIKE, ipHE_LIKE, ipLev, iso, iso_recomb_integrand(), kTRyd, L_, MAX2, MILNE_CONST, MIN2, N_, pow(), qg32(), S_, SMALLDOUBLE, StatesElemNEW, TE1RYD, TotalInsanity(), and t_iso::xIsoLevNIonRyd.

Referenced by iso_radiative_recomb(), iso_recomb_check(), and iso_recomb_setup().

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 801 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: