ion_recomb_Badnell.cpp File Reference

#include "cddefines.h"
#include "phycon.h"
#include "elementnames.h"
#include "atmdat.h"
#include "atmdat_adfa.h"
#include "iso.h"
#include "ionbal.h"
#include "save.h"
#include "freebound.h"
#include "dense.h"
#include "ran.h"
#include "thirdparty.h"

Include dependency graph for ion_recomb_Badnell.cpp:

Functions
STATIC double	CollisSuppres (long int atomic_number, long int ionic_charge, double eden, double T)
STATIC double	Badnell_DR_rate_eval (int nAtomicNumberCScale, int n_core_e_before_recomb)
STATIC double	Badnell_RR_rate_eval (int nAtomicNumberCScale, int n_core_e_before_recomb)
void	Badnell_rec_init (void)
void	ion_recom_calculate (void)

Variables
static const int	MAX_FIT_PAR_DR = 9
static multi_arr< double, 3 >	DRFitParPart1
static multi_arr< double, 3 >	DRFitParPart2
static multi_arr< int, 2 >	nDRFitPar
static const int	MAX_FIT_PAR_RR = 6
static multi_arr< double, 3 >	RRFitPar
static multi_arr< bool, 2 >	lgDRBadnellDefined
static multi_arr< bool, 2 >	lgDRBadnellDefinedPart2
static multi_arr< bool, 2 >	lgRRBadnellDefined
static multi_arr< bool, 2 >	lgDR_BadWeb_exist
static bool	lgMustAllocateRec =true
static double	RecNoise [LIMELM]
static double	DR_Badnell_rate_coef_mean_ion [LIMELM]
static char	chDRDataSource [LIMELM][LIMELM][10]
static char	chRRDataSource [LIMELM][LIMELM][10]

Function Documentation

STATIC double Badnell_DR_rate_eval	(	int	nAtomicNumberCScale,
		int	n_core_e_before_recomb
	)

Badnell_DR_rate_eval This code is written by Terry Yun, 2005 
It interpolates the rate coefficients in a given temperature.
It receives atomic number on Physics scale, with H = 1, 
and the number of core electrons before recombination, and returns the rate coefficient*
It returns
'-2': initial <= final
init < 0 or init >302 or final < 0 or final > 302
'-1': the transition is not defined
'99': unknown invalid entries                         

Parameters

z_val	atomic number on C scale - He is 1
n_val	number of core electrons before capture of free electron

References ASSERT, chDRDataSource, DEBUG_ENTRY, ipIRON, LIMELM, phycon, sexp(), t_phycon::te, and t_phycon::te32.

Referenced by Badnell_rec_init(), and ion_recom_calculate().

Here is the call graph for this function:

void Badnell_rec_init

(

void

)

Badnell_rec_init This code is written by Terry Yun, 2005 * It reads rate coefficient fits into 3D arrays and output array.out for testing * The testing can be commented out

References multi_arr< T, d, ALLOC, lgBC >::alloc(), ASSERT, Badnell_DR_rate_eval(), Badnell_RR_rate_eval(), cdEXIT, multi_arr< T, d, ALLOC, lgBC >::clone(), DEBUG_ENTRY, dense, freeBound::DielecRecombVsTemp, DR_Badnell_rate_coef_mean_ion, EXIT_FAILURE, exp10(), t_iso_sp::fb, FFmtRead(), fprintf(), t_ionbal::guess_noise, ionbal, ioQQQ, ipH_LIKE, ipHE_LIKE, ipHELIUM, ipHYDROGEN, iso_sp, L_, t_dense::lgElmtOn, lgMustAllocateRec, LIMELM, MAX2, MAX_FIT_PAR_DR, MAX_FIT_PAR_RR, t_iso_sp::nCollapsed_max, NISO, t_ran::normal(), NUM_DR_TEMPS, t_iso_sp::numLevels_max, open_data(), pow3(), t_iso_sp::QN2Index(), ran, read_whole_line(), RecNoise, multi_arr< T, d, ALLOC, lgBC >::reserve(), S, S_, TorF(), TotalInsanity(), and Triangle2().

Referenced by cloudy(), and ParseCompile().

Here is the call graph for this function:

STATIC double Badnell_RR_rate_eval	(	int	nAtomicNumberCScale,
		int	n_core_e_before_recomb
	)

Badnell_RR_rate_eval

Parameters

z_val	atomic number on C scale - He - 1
n_val	number of core electrons before capture of free electron

References ASSERT, chRRDataSource, DEBUG_ENTRY, ipIRON, LIMELM, phycon, pow(), and t_phycon::te.

Referenced by Badnell_rec_init(), and ion_recom_calculate().

Here is the call graph for this function:

STATIC double CollisSuppres	(	long int	atomic_number,
		long int	ionic_charge,
		double	eden,
		double	T
	)

References ASSERT, DEBUG_ENTRY, fprintf(), ioQQQ, pow(), pow2(), powi(), and TotalInsanity().

Referenced by ion_recom_calculate().

Here is the call graph for this function:

void ion_recom_calculate

(

void

)

    ion_recom_calculate called by conv_base to calculate radiative and dielectronic

recombination rate coefficients

References ASSERT, atmdat_dielrec_fe(), Badnell_DR_rate_eval(), Badnell_RR_rate_eval(), chDRDataSource, t_elementnames::chElementNameShort, chRRDataSource, CollisSuppres(), DEBUG_ENTRY, dense, t_ionbal::DR_Badnell_rate_coef, DR_Badnell_rate_coef_mean_ion, t_ionbal::DR_Badnell_suppress_fact, t_ionbal::DR_mean_scale, t_dense::eden, elementnames, fp_equal(), fprintf(), Singleton< t_ADfA >::Inst(), ionbal, ioQQQ, ipHELIUM, ipHYDROGEN, ipIRON, t_save::ipSDSFile, t_ionbal::lgDRsup, t_ionbal::lgRecom_Badnell_print, t_save::lgSDSOn, LIMELM, MIN2, phycon, powpq(), t_ADfA::rad_rec(), RecNoise, t_ionbal::RR_Badnell_rate_coef, t_ionbal::RR_rate_coef_used, t_ionbal::RR_Verner_rate_coef, save, sexp(), t_phycon::te, and t_phycon::te_eV.

Referenced by ConvBase().

Here is the call graph for this function:

Variable Documentation

char chDRDataSource[LIMELM][LIMELM][10]

static

Referenced by Badnell_DR_rate_eval(), and ion_recom_calculate().

char chRRDataSource[LIMELM][LIMELM][10]

static

Referenced by Badnell_RR_rate_eval(), and ion_recom_calculate().

double DR_Badnell_rate_coef_mean_ion[LIMELM]

static

Referenced by Badnell_rec_init(), and ion_recom_calculate().

multi_arr<double,3> DRFitParPart1

static

multi_arr<double,3> DRFitParPart2

static

multi_arr<bool,2> lgDR_BadWeb_exist

static

multi_arr<bool,2> lgDRBadnellDefined

static

multi_arr<bool,2> lgDRBadnellDefinedPart2

static

bool lgMustAllocateRec =true

static

Referenced by Badnell_rec_init().

multi_arr<bool,2> lgRRBadnellDefined

static

const int MAX_FIT_PAR_DR = 9

static

Referenced by Badnell_rec_init().

const int MAX_FIT_PAR_RR = 6

static

Referenced by Badnell_rec_init().

multi_arr<int,2> nDRFitPar

static

double RecNoise[LIMELM]

static

Referenced by Badnell_rec_init(), and ion_recom_calculate().

multi_arr<double,3> RRFitPar

static