#include "cddefines.h"
#include "atmdat_adfa.h"
#include "conv.h"
#include "cosmology.h"
#include "elementnames.h"
#include "helike_recom.h"
#include "hydrogenic.h"
#include "ionbal.h"
#include "iso.h"
#include "opacity.h"
#include "phycon.h"
#include "prt.h"
#include "save.h"
#include "thirdparty.h"
#include "trace.h"
#include "rt.h"
#include "freebound.h"
#include "dense.h"
#include "integrate.h"
#include "container_classes.h"

Include dependency graph for iso_radiative_recomb.cpp:

Functions
STATIC double	TempInterp (double TempArray, double ValueArray, long NumElements, double temp)

STATIC double	iso_recomb_integrand (double EE)

STATIC void	iso_put_recomb_error (long ipISO, long nelem)

STATIC 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, double temp, long n)

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

void	iso_recomb_alloc ()

void	iso_recomb_auxiliary_free ()

void	iso_recomb_setup (long ipISO)

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

Variables
static multi_arr< double, 4 >	RRCoef

static multi_arr< long, 2 >	NumLevRecomb

static multi_arr< double, 3 >	TotalRecomb

static double	TeRRCoef [N_ISO_TE_RECOMB]

static double	kTRyd

static double	global_EthRyd

static long int	globalZ

static long int	globalISO

static long int	globalN

static long int	globalL

static long int	globalS

Function Documentation

◆ iso_cross_section()

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

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

Referenced by iso_recomb_integrand().

Here is the call graph for this function:

◆ iso_dielec_recomb_rate()

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

References ASSERT, DEBUG_ENTRY, freeBound::DielecRecombVsTemp, exp10(), hunt_bisect(), ipHe1s1S, ipHE_LIKE, iso_ctrl, iso_sp, NUM_DR_TEMPS, and phycon.

Referenced by iso_radiative_recomb().

Here is the call graph for this function:

◆ iso_put_recomb_error()

STATIC void iso_put_recomb_error	(	long	ipISO,
		long	nelem )

References ASSERT, ipHE_LIKE, ipHELIUM, IPRAD, iso_put_error(), iso_sp, min(), QN2ind(), and STATIC.

Referenced by iso_radiative_recomb().

Here is the call graph for this function:

◆ 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

References ASSERT, cdEXIT, conv, cosmology, DEBUG_ENTRY, dense, elementnames, EXIT_FAILURE, fp_equal(), fprintf(), t_ADfA::H_rad_rec(), Singleton< t_ADfA >::Inst(), ionbal, ioQQQ, ipH1s, ipHE_LIKE, ipHELIUM, 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, LIMELM, MAX2, MIN2, N_, NHYDRO_MAX_LEVEL, NISO, opac, phycon, PrintEfmt, Recomb_Seaton59(), RT_recom_effic(), save, ShowMe(), SumUpToThisN, and trace.

Referenced by iso_update_rates().

Here is the call graph for this function:

◆ iso_radiative_recomb_effective()

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

References ASSERT, DEBUG_ENTRY, dprintf(), fprintf(), ioQQQ, IPRAD, ipRecRad, iso_ctrl, iso_sp, L_, MAX2, N_, phycon, pow2, t_wavl::prt_wl(), S_, t_vac(), and t_wavl::wavlVac().

Referenced by iso_update_rates().

Here is the call graph for this function:

◆ iso_radrecomb_from_cross_section()

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

References DEBUG_ENTRY, global_EthRyd, globalISO, globalL, globalN, globalS, globalZ, t_ADfA::H_rad_rec(), Singleton< t_ADfA >::Inst(), ipH_LIKE, ipHE_LIKE, iso_recomb_integrand(), iso_sp, kTRyd, L_, MAX2, MIN2, N_, powpq(), qg32(), S_, SMALLDOUBLE, STATIC, and TotalInsanity().

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

Here is the call graph for this function:

◆ iso_recomb_alloc()

void iso_recomb_alloc ( )

iso_recomb_alloc - allocate space needed for iso recombination tables.

References DEBUG_ENTRY, dense, iso_get_total_num_levels(), iso_sp, LIKE_RREC_MAXN(), LIMELM, MAX2, N_ISO_TE_RECOMB, NISO, NumLevRecomb, RRCoef, RREC_MAXN, TeRRCoef, and TotalRecomb.

Referenced by iso_create().

Here is the call graph for this function:

◆ iso_recomb_auxiliary_free()

void iso_recomb_auxiliary_free ( )

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

References DEBUG_ENTRY, and NumLevRecomb.

Referenced by iso_create().

◆ iso_recomb_check()

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