#include "cddefines.h"
#include "atmdat.h"
#include "conv.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 "punch.h"
#include "thermal.h"
#include "thirdparty.h"
#include "trace.h"
#include "rt.h"

Include dependency graph for iso_radiative_recomb.cpp:

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

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

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]

Function Documentation

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 985 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, and phycon.

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

References ASSERT, cdEXIT, t_elementnames::chElementSym, t_iso::chISO, t_iso::ConOpacRatio, conv, DEBUG_ENTRY, t_iso::DielecRecomb, t_ionbal::DR_Badnell_rate_coef, elementnames, t_iso::ErrorFactor, hydro, Singleton< T >::Inst(), ionbal, ioQQQ, t_punch::ioRecom, ipH1s, ipHELIUM, t_iso::ipIsoLevNIonCon, t_trace::ipIsoTrace, IPRAD, ipRecEsc, ipRecNetEsc, ipRecRad, iso, iso_dielec_recomb_rate(), iso_radrecomb_from_cross_section(), iso_RRCoef_Te(), t_opac::lgCaseB, t_iso::lgDielRecom, t_punch::lgioRecom, t_trace::lgIsoTraceFull, t_iso::lgLevelsLowered, t_iso::lgNoRecombInterp, t_iso::lgRandErrGen, t_hydro::lgReevalRecom, t_iso::lgTopoff, t_trace::lgTrace, LIMELM, MAX2, MIN2, N_, 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(), punch, t_iso::RadRec_caseB, t_iso::RadRec_effec, t_iso::RadRecomb, Recomb_Seaton59(), RT_recom_effic(), ShowMe(), SumUpToThisN, t_opac::TauAbsGeo, t_phycon::te, and trace.

Referenced by iso_solve().

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 440 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, prt_wl(), t_iso::RadEffec, t_iso::RadRecomb, RYDLAM, S_, t_iso::SigmaAtot, t_iso::SigmaCascadeProb, t_iso::SigmaRadEffec, StatesElem, t_phycon::te, wavelength, and t_iso::xIsoLevNIonRyd.

Referenced by iso_solve().

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

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

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 661 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

ipISO
nelem	the chemical element, 1 for He
level	the level, 0 for ground
temperature	the temperature to be used

Definition at line 572 of file iso_radiative_recomb.cpp.

References DEBUG_ENTRY, iso_radrecomb_from_cross_section(), iso_RRCoef_Te(), MAX2, phycon, t_phycon::te, and TempChange().

Referenced by SanityCheckBegin().

Here is the call graph for this function:

void iso_recomb_malloc ( void )

iso_recomb_malloc - malloc space needed for iso recombination tables.

Definition at line 601 of file iso_radiative_recomb.cpp.

References DEBUG_ENTRY, dense, iso, iso_get_total_num_levels(), t_dense::lgElmtOn, LIKE_RREC_MAXN, LIMELM, MALLOC, MAX2, N_ISO_TE_RECOMB, NISO, t_iso::numLevels_max, NumLevRecomb, RRCoef, RREC_MAXN, TeRRCoef, and TotalRecomb.

Referenced by iso_create().

Here is the call graph for this function:

void iso_recomb_setup ( long ipISO )

iso_recomb_setup - read in or compile iso recombination tables.

Parameters

ipISO

Establish radiative recombination rate coefficients - RRC

Definition at line 674 of file iso_radiative_recomb.cpp.

References AS_LOCAL_ONLY, ASSERT, cdEXIT, t_elementnames::chElementSym, t_iso::chISO, chLine_LENGTH, DEBUG_ENTRY, dense, elementnames, FFmtRead(), INPUT_LINE_LENGTH, Singleton< T >::Inst(), ioQQQ, ipHe1s1S, iso, iso_radrecomb_from_cross_section(), t_iso::lgCompileRecomb, t_dense::lgElmtOn, t_iso::lgNoRecombInterp, t_trace::lgTrace, LIKE_RREC_MAXN, LIMELM, t_iso::n_HighestResolved_max, N_ISO_TE_RECOMB, t_iso::nCollapsed_max, NHYDRO_MAX_LEVEL, NISO, t_iso::numLevels_max, NumLevRecomb, open_data(), read_whole_line(), Recomb_Seaton59(), RECOMBMAGIC, RRCoef, RREC_MAXN, SumUpToThisN, TeRRCoef, TotalRecomb, and trace.

Referenced by iso_create().

Here is the call graph for this function:

double iso_RRCoef_Te	(	long	ipISO,
		long	nelem,
		long	n
	)

iso_RRCoef_Te - interpolate iso recomb coeff as function of temperature

Parameters

ipISO
nelem
n

Definition at line 546 of file iso_radiative_recomb.cpp.

References ASSERT, DEBUG_ENTRY, iso, t_iso::lgNoRecombInterp, N_ISO_TE_RECOMB, t_iso::nCollapsed_max, t_iso::numLevels_max, RRCoef, TempInterp(), TeRRCoef, and TotalRecomb.

Referenced by iso_radiative_recomb(), and iso_recomb_check().

Here is the call graph for this function:

STATIC double TempInterp	(	double *	TempArray,
		double *	ValueArray,
		long	NumElements
	)

Todo:: 2 use a canned interpolation routine, no need for special one here

Definition at line 1059 of file iso_radiative_recomb.cpp.

References t_phycon::alogte, ASSERT, cdEXIT, DEBUG_ENTRY, hunt_bisect(), ioQQQ, lagrange(), max(), min(), phycon, and t_phycon::te.

Referenced by iso_RRCoef_Te().

Here is the call graph for this function:

Variable Documentation

long** NumLevRecomb

static

Definition at line 31 of file iso_radiative_recomb.cpp.

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

double**** RRCoef

static

Definition at line 30 of file iso_radiative_recomb.cpp.

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

double TeRRCoef[N_ISO_TE_RECOMB]

static

Definition at line 35 of file iso_radiative_recomb.cpp.

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

double*** TotalRecomb

static

Definition at line 32 of file iso_radiative_recomb.cpp.

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

Functions

Variables

Function Documentation

Variable Documentation