#include <h2_priv.h>

Collaboration diagram for diatomics:

Public Member Functions
double	Abund () const

void	GetIndices (long &ipHi, long &ipLo, const char *chLine, long &i) const

void	CalcPhotoionizationRate (void)

long	OpacityCreate (vector< double > &stack)

double	GetHeatRate (const diss_tran &tran)

double	GetDissociationRate (const diss_tran &tran)

double	MolDissocOpacity (const diss_tran &tran, const double &Mol_Ene)

double	Cont_Diss_Heat_Rate (void)

void	Mol_Photo_Diss_Rates (void)

void	Read_Mol_Diss_cross_sections (void)

void	SolveExcitedElectronicLevels (void)

void	SolveSomeGroundElectronicLevels (void)

double	GetExcitedElecDensity (void)

realnum	GetXColden (long iVib, long iRot)

long int	getLine (long iElecHi, long iVibHi, long iRotHi, long iElecLo, long iVibLo, long iRotLo, double relint, double absint)

realnum	H2_CollidRateEvalOne (long iVibHi, long iRotHi, long iVibLo, long iRotLo, long ipHi, long ipLo, long nColl, double temp_K)

void	H2_Calc_Average_Rates (void)

void	H2_X_sink_and_source (void)

void	H2_X_coll_rate_evaluate (void)

void	H2_Level_low_matrix (realnum abundance)

void	H2_ReadEnergies ()

void	H2_ReadEnergies (long int nelec, vector< int > &n, vector< int > &v, vector< int > &J, vector< double > &eWN)

void	H2_ReadDissprob (long int nelec)

void	H2_CollidRateEvalAll (void)

void	H2_CollidRateRead (long int nColl)

void	H2_ReadTransprob (long int nelec, TransitionList &trans)

void	H2_Read_hminus_distribution (void)

void	mole_H2_form (void)

void	mole_H2_LTE (void)

void	H2_Solomon_rate (void)

double	gs_rate (void)

void	H2_zero_pops_too_low (void)

void	init (void)

void	H2_ContPoint (void)

double	H2_DR (void)

double	H2_Accel (void)

void	H2_RT_OTS (void)

double	H2_RadPress (void)

void	H2_LinesAdd (void)

void	H2_Reset (void)

double	H2_InterEnergy (void)

void	H2_Colden (const char *chLabel)

void	H2_Cooling (void)

double	LTE_Cooling_per_H2 ()

void	H2_Punch_line_data (FILE *ioPUN, bool lgDoAll)

void	H2_PunchLineStuff (FILE *io, realnum xLimit, long index)

void	H2_RT_diffuse (void)

void	H2_RTMake (linefunc line_one)

void	H2_RT_tau_inc (void)

void	H2_Prt_Zone (void)

void	H2_PrtDepartCoef (void)

void	H2_LineZero (void)

void	H2_RT_tau_reset (void)

void	H2_LevelPops (bool &lgPopsConverged, double &old_value, double &new_value)

void	H2_PunchDo (FILE *io, char chJOB[], const char chTime[], long int ipPun)

void	H2_ParseSave (Parser &p, ostringstream &chHeader)

double	H2_itrzn (void)

void	H2_Prt_column_density (FILE *ioMEAN)

void	set_numLevelsMatrix (long numLevels)

void	H2_ReadDissocEnergies (void)

	diatomics (const string &a, const double &e_star, const double const abund, double(fun)(double))

void	H2_Read_LTE_cooling_per_H2 ()

double	interpolate_LTE_Cooling (double Temp)

Public Attributes
double(*	photoion_opacity_fun )(double energy)

bool	lgREAD_DATA

double	photoionize_rate

double	photo_heat_soft

double	photo_heat_hard

double	photodissoc_BigH2_H2s

double	photodissoc_BigH2_H2g

double	spon_diss_tot

double	Solomon_dissoc_rate_g

double	Solomon_dissoc_rate_s

double	Solomon_elec_decay_g

double	Solomon_elec_decay_s

double	rate_grain_op_conserve

double	rate_grain_J1_to_J0

double	Cont_Dissoc_Rate_H2s

double	Cont_Dissoc_Rate_H2g

multi_arr< double, 3 >	Cont_Dissoc_Rate

double	rel_pop_LTE_g

double	rel_pop_LTE_s

double	average_energy_g

double	average_energy_s

double	HeatDiss

double	HeatDexc

double	HeatDexc_old

double	HeatDexc_deriv

double	HeatChangeOld

double	HeatChange

double	Average_A

double	Average_collH2_deexcit

double	Average_collH_deexcit

double	Average_collH2_excit

double	Average_collH_excit

double	Average_collH_dissoc_g

double	Average_collH_dissoc_s

double	Average_collH2_dissoc_g

double	Average_collH2_dissoc_s

bool	lgEvaluated

long	ip_photo_opac_thresh

long	ip_photo_opac_offset

t_coll_source	coll_source [N_X_COLLIDER]

double	ortho_density

double	para_density

realnum	ortho_density_f

realnum	para_density_f

double	ortho_colden

double	para_colden

double	ortho_para_old

double	ortho_para_older

double	ortho_para_current

double	renorm_max

double	renorm_min

long int	nCall_this_zone

bool	lgEnabled

int	nElecLevelOutput

bool	lgColl_gbar

bool	lgColl_deexec_Calc

bool	lgColl_dissoc_coll

bool	lgH2_grain_deexcitation

bool	lgLTE

bool	lgH2_ortho_para_coll_on

bool	lgH2_He_ORNL

bool	lgH2_ORH2_ORNL

bool	lgH2_PAH2_ORNL

bool	lgH2_NOISE

bool	lgH2_NOISECOSMIC

long int	loop_h2_oscil

long int	nzoneEval

double	xMeanNoise

double	xSTDNoise

double	H2_to_H_limit

realnum	mass_amu

int	nTRACE

int	n_trace_final

int	n_trace_iterations

int	n_trace_full

int	n_trace_matrix

long int	n_elec_states

double	frac_matrix

double	TeUsedBoltz

double	TeUsedColl

molecule *	sp

molecule *	sp_star

qList	states

TransitionList	trans

TransitionList::iterator	rad_end

vector< diss_tran >	Diss_Trans

const double	ENERGY_H2_STAR

Private Attributes
string	label

string	shortlabel

string	path

const double *const	dense_total

char	chH2ColliderLabels [N_X_COLLIDER][chN_X_COLLIDER]

long int	nEner_H2_ground

multi_arr< double, 2 >	pops_per_vib

double	H2_renorm_chemistry

multi_arr< realnum, 2 >	H2_X_coll_rate

double	H2_DissocEnergies [N_ELEC]

long int	nVib_hi [N_ELEC]

valarray< long >	nRot_hi [N_ELEC]

long int	Jlowest [N_ELEC]

long int	nLevels_per_elec [N_ELEC]

double	pops_per_elec [N_ELEC]

multi_arr< realnum, 3 >	CollRateCoeff

multi_arr< realnum, 3 >	CollRateErrFac

vector< CollRateCoeffArray >	RateCoefTable

multi_arr< realnum, 3 >	H2_dissprob

multi_arr< realnum, 3 >	H2_disske

multi_arr< double, 3 >	H2_rad_rate_out

multi_arr< double, 3 >	H2_old_populations

multi_arr< double, 3 >	H2_populations_LTE

multi_arr< bool, 3 >	H2_lgOrtho

long int	nzoneAsEval

long int	iterationAsEval

multi_arr< int, 2 >	H2_ipPhoto

multi_arr< double, 2 >	H2_col_rate_in

multi_arr< double, 2 >	H2_col_rate_out

multi_arr< double, 2 >	H2_rad_rate_in

multi_arr< realnum, 2 >	H2_X_formation

multi_arr< realnum, 2 >	H2_X_Hmin_back

multi_arr< realnum, 2 >	H2_X_colden

multi_arr< realnum, 2 >	H2_X_colden_LTE

multi_arr< double, 2 >	H2_X_rate_from_elec_excited

multi_arr< double, 2 >	H2_X_rate_to_elec_excited

multi_arr< realnum, 2 >	H2_coll_dissoc_rate_coef

multi_arr< realnum, 2 >	H2_coll_dissoc_rate_coef_H2

valarray< realnum >	H2_X_source

valarray< realnum >	H2_X_sink

multi_arr< realnum, 3 >	H2_X_grain_formation_distribution

double	H2_den_s

double	H2_den_g

multi_arr< realnum, 3 >	H2_X_hminus_formation_distribution

valarray< long >	ipVib_H2_energy_sort

valarray< long >	ipElec_H2_energy_sort

valarray< long >	ipRot_H2_energy_sort

multi_arr< long int, 3 >	ipEnergySort

multi_arr< long int, 2 >	ipTransitionSort

long int	nXLevelsMatrix

long int	ndim_allocated

multi_arr< double, 2 >	AulEscp

multi_arr< double, 2 >	AulDest

multi_arr< double, 2 >	AulPump

multi_arr< double, 2 >	CollRate_levn

vector< double >	pops

vector< double >	create

vector< double >	destroy

vector< double >	depart

vector< double >	stat_levn

vector< double >	excit

long int	levelAsEval

bool	lgFirst

long int	nzone_eval

long int	iteration_evaluated

multi_arr< realnum, 6 >	H2_SaveLine

multi_arr< bool, 2 >	lgH2_radiative

long int	nH2_pops

long int	nH2_zone

long int	nzone_nlevel_set

long int	nCall_this_iteration

vector< double >	LTE_Temp

vector< double >	LTE_cool

bool	lgPrtMatrix

bool	lgImgMatrix

Constructor & Destructor Documentation

◆ diatomics()

diatomics::diatomics	(	const string &	a,
		const double &	e_star,
		const double *const	abund,
		double(*	fun )(double) )

explicit

Todo: 3 these should be const since cannot change, are flags

References abund, Average_A, Average_collH2_deexcit, Average_collH2_dissoc_g, Average_collH2_dissoc_s, Average_collH2_excit, Average_collH_deexcit, Average_collH_dissoc_g, Average_collH_dissoc_s, Average_collH_excit, chH2ColliderLabels, DEBUG_ENTRY, dense_total, ENERGY_H2_STAR, fixit, H2_to_H_limit, HeatChange, HeatChangeOld, HeatDexc, HeatDexc_deriv, HeatDexc_old, HeatDiss, iteration_evaluated, iterationAsEval, label, levelAsEval, lgColl_deexec_Calc, lgColl_dissoc_coll, lgColl_gbar, lgEnabled, lgEvaluated, lgFirst, lgH2_grain_deexcitation, lgH2_He_ORNL, lgH2_NOISE, lgH2_NOISECOSMIC, lgH2_ORH2_ORNL, lgH2_ortho_para_coll_on, lgH2_PAH2_ORNL, lgLTE, lgREAD_DATA, loop_h2_oscil, N_ELEC, n_elec_states, n_trace_final, n_trace_full, n_trace_iterations, n_trace_matrix, nCall_this_zone, ndim_allocated, nElecLevelOutput, nH2_pops, nH2_zone, nTRACE, nXLevelsMatrix, nzone_eval, nzone_nlevel_set, nzoneAsEval, nzoneEval, ortho_colden, ortho_para_current, ortho_para_old, ortho_para_older, para_colden, path, photo_heat_hard, photo_heat_soft, photodissoc_BigH2_H2g, photodissoc_BigH2_H2s, photoion_opacity_fun, rate_grain_J1_to_J0, rate_grain_op_conserve, renorm_max, renorm_min, shortlabel, Solomon_dissoc_rate_g, Solomon_dissoc_rate_s, spon_diss_tot, states, TeUsedBoltz, TeUsedColl, toupper(), and trans.

Here is the call graph for this function:

Member Function Documentation

◆ Abund()

double diatomics::Abund ( ) const

inline

References dense_total.

Referenced by H2_PunchDo().

◆ CalcPhotoionizationRate()

void diatomics::CalcPhotoionizationRate ( void )

References DEBUG_ENTRY, GammaK(), t_phoHeat::HeatHiEnr, t_phoHeat::HeatLowEnr, ionbal, ip_photo_opac_offset, ip_photo_opac_thresh, ipHYDROGEN, iteration, iteration_evaluated, nzone, nzone_eval, photo_heat_hard, photo_heat_soft, photoionize_rate, and rfield.

Referenced by H2_LevelPops().

Here is the call graph for this function:

◆ Cont_Diss_Heat_Rate()

double diatomics::Cont_Diss_Heat_Rate ( void )

References DEBUG_ENTRY, Diss_Trans, lgEnabled, Mol_Photo_Diss_Rates(), and mole_global.

Here is the call graph for this function:

◆ GetDissociationRate()

double diatomics::GetDissociationRate ( const diss_tran & tran )

References DEBUG_ENTRY, diss_tran::initial, ipEnergySort, diss_level::j, diss_level::n, diss_tran::rate_coeff, states, and diss_level::v.

◆ GetExcitedElecDensity()

double diatomics::GetExcitedElecDensity ( void )

References n_elec_states, and pops_per_elec.

Referenced by H2_Solomon_rate().

◆ GetHeatRate()

double diatomics::GetHeatRate ( const diss_tran & tran )

References DEBUG_ENTRY, density(), diss_tran::energies, energy(), diss_tran::initial, ipEnergySort, ipoint(), diss_level::j, MAX2, MIN2, MolDissocCrossSection(), diss_level::n, rfield, states, and diss_level::v.

Here is the call graph for this function:

◆ GetIndices()

void diatomics::GetIndices	(	long &	ipHi,
		long &	ipLo,
		const char *	chLine,
		long &	i ) const

References ASSERT, FFmtRead(), ipEnergySort, Jlowest, nRot_hi, and nVib_hi.

Here is the call graph for this function:

◆ getLine()

long int diatomics::getLine	(	long	iElecHi,
		long	iVibHi,
		long	iRotHi,
		long	iElecLo,
		long	iVibLo,
		long	iRotLo,
		double *	relint,
		double *	absint )

References ASSERT, DEBUG_ENTRY, energy(), H2_lgOrtho, H2_SaveLine, ipEnergySort, lgH2_radiative, LineSave, radius, and states.

Here is the call graph for this function:

◆ GetXColden()

realnum diatomics::GetXColden	(	long	iVib,
		long	iRot )

References DEBUG_ENTRY, fprintf(), H2_X_colden, ioQQQ, nRot_hi, and nVib_hi.

Referenced by cdH2_colden().

Here is the call graph for this function:

◆ gs_rate()

double diatomics::gs_rate ( void )

gs_rate evaluate rate between ground and star states

References EmissionProxy::Aul(), DEBUG_ENTRY, H2_den_g, H2_dissprob, H2_rad_rate_out, ipEnergySort, ipTransitionSort, Jlowest, lgH2_radiative, MAX2, n_elec_states, nEner_H2_ground, nRot_hi, nVib_hi, EmissionProxy::Ploss(), SDIV(), states, and trans.

Referenced by H2_LevelPops().

Here is the call graph for this function:

◆ H2_Accel()

double diatomics::H2_Accel ( void )

radiative acceleration due to H2 called in rt_line_driving

References ASSERT, DEBUG_ENTRY, lgEnabled, rad_end, and trans.

◆ H2_Calc_Average_Rates()

void diatomics::H2_Calc_Average_Rates ( void )

References Average_A, Average_collH2_deexcit, Average_collH2_excit, Average_collH_deexcit, Average_collH_excit, CollRateCoeff, ENERGY_H2_STAR, H2_lgOrtho, hmi, ipEnergySort, ipTransitionSort, lgH2_radiative, safe_div(), SDIV(), states, and trans.

Referenced by H2_LevelPops().

Here is the call graph for this function:

◆ H2_Colden()

void diatomics::H2_Colden ( const char * chLabel )

H2_Colden maintain H2 column densities within X

Parameters

*chLabel

References cdEXIT, DEBUG_ENTRY, dense_total, EXIT_FAILURE, fprintf(), H2_populations_LTE, H2_X_colden, H2_X_colden_LTE, ioQQQ, lgEnabled, radius, and states.

Here is the call graph for this function:

◆ H2_CollidRateEvalAll()

void diatomics::H2_CollidRateEvalAll ( void )

H2_CollidRateEvalAll - set H2 collision rates

References ASSERT, CollRateCoeff, DEBUG_ENTRY, energy(), fixit, fprintf(), H2_coll_dissoc_rate_coef, H2_coll_dissoc_rate_coef_H2, H2_CollidRateEvalOne(), H2_DissocEnergies, ioQQQ, ipRot_H2_energy_sort, ipVib_H2_energy_sort, label, lgColl_dissoc_coll, n_trace_full, N_X_COLLIDER, nLevels_per_elec, nTRACE, phycon, sexp(), and states.

Referenced by H2_LevelPops().

Here is the call graph for this function:

◆ H2_CollidRateEvalOne()

realnum diatomics::H2_CollidRateEvalOne	(	long	iVibHi,
		long	iRotHi,
		long	iVibLo,
		long	iRotLo,
		long	ipHi,
		long	ipLo,
		long	nColl,
		double	temp_K )

References CollRateErrFac, DEBUG_ENTRY, GbarRateCoeff(), H2_lgOrtho, InterpCollRate(), lgColl_deexec_Calc, lgColl_gbar, lgH2_NOISE, lgH2_ortho_para_coll_on, RateCoefTable, and states.

Referenced by H2_CollidRateEvalAll().

Here is the call graph for this function:

◆ H2_CollidRateRead()

void diatomics::H2_CollidRateRead ( long int nColl )

read collision rates

Parameters

nColl

References cdEXIT, coll_source, cpu, DEBUG_ENTRY, EXIT_FAILURE, fprintf(), ioQQQ, nLevels_per_elec, open_data(), path, RateCoefTable, read_whole_line(), and ReadCollisionRateTable().

Referenced by init().

Here is the call graph for this function:

◆ H2_ContPoint()

void diatomics::H2_ContPoint ( void )

set the ipCont struc element for the H2 molecule, called by ContCreatePointers

References ASSERT, DEBUG_ENTRY, ipFineCont(), ipLineEnergy(), label, lgEnabled, rad_end, and trans.

Here is the call graph for this function:

◆ H2_Cooling()

void diatomics::H2_Cooling ( void )

H2_cooling evaluate cooling and heating due to H2 molecule

References ASSERT, collider_density, CollRateCoeff, conv, DEBUG_ENTRY, energy(), fnzone, fprintf(), H2_disske, H2_dissprob, H2_renorm_chemistry, HeatDexc, HeatDexc_deriv, HeatDiss, ioQQQ, lgEnabled, n_trace_full, N_X_COLLIDER, nCall_this_iteration, nLevels_per_elec, nTRACE, phycon, POW2, PRT_POPS, safe_div(), states, and thermal.

Referenced by H2_LevelPops().

Here is the call graph for this function:

◆ H2_DR()

double diatomics::H2_DR ( void )

H2_DR choose next zone thickness based on H2 big molecule

References BIGFLOAT.

◆ H2_InterEnergy()

double diatomics::H2_InterEnergy ( void )

internal energy of H2 called in PresTotCurrent

References ipPun.

◆ H2_itrzn()

double diatomics::H2_itrzn ( void )

H2_itrzn - average number of H2 pop evaluations per zone

References lgEnabled, nH2_pops, and nH2_zone.

◆ H2_Level_low_matrix()

void diatomics::H2_Level_low_matrix ( realnum abundance )

References ASSERT, AulDest, AulEscp, AulPump, cdEXIT, CollRate_levn, create, DEBUG_ENTRY, dense_total, depart, destroy, dsexp(), energy(), excit, EXIT_FAILURE, fixit, fprintf(), H2_lgOrtho, H2_X_coll_rate, H2_X_rate_from_elec_excited, H2_X_rate_to_elec_excited, H2_X_sink, H2_X_source, ioQQQ, ipRot_H2_energy_sort, ipTransitionSort, ipVib_H2_energy_sort, iteration, iterationAsEval, levelAsEval, lgFirst, lgH2_radiative, lgImgMatrix, lgPrtMatrix, n_trace_full, n_trace_matrix, ndim_allocated, nLevels_per_elec, nTRACE, nXLevelsMatrix, nzone, nzoneAsEval, phycon, pops, rate_grain_J1_to_J0, rate_grain_op_conserve, stat_levn, states, and trans.

Referenced by H2_LevelPops().

Here is the call graph for this function:

◆ H2_LevelPops()

void diatomics::H2_LevelPops	(	bool &	lgPopsConverged,
		double &	old_value,
		double &	new_value )

do level populations for H2, called by iso_solve

References ASSERT, Average_collH2_dissoc_g, Average_collH2_dissoc_s, Average_collH_dissoc_g, Average_collH_dissoc_s, average_energy_g, average_energy_s, CalcPhotoionizationRate(), conv, DEBUG_ENTRY, dense, dense_total, ENERGY_H2_STAR, exp10(), fixit, fnzone, fp_equal(), fp_equal_tol(), fprintf(), FRAC, frac_matrix, gs_rate(), h2, H2_Calc_Average_Rates(), H2_coll_dissoc_rate_coef, H2_coll_dissoc_rate_coef_H2, H2_CollidRateEvalAll(), H2_Cooling(), H2_den_g, H2_den_s, H2_DISS_ALLISON_DALGARNO, H2_DissocEnergies, H2_Level_low_matrix(), H2_lgOrtho, H2_old_populations, H2_populations_LTE, H2_rad_rate_in, H2_rad_rate_out, H2_renorm_chemistry, H2_Solomon_rate(), H2_to_H_limit, H2_X_coll_rate_evaluate(), H2_X_rate_from_elec_excited, H2_X_rate_to_elec_excited, H2_X_sink_and_source(), H2_zero_pops_too_low(), hd, HeatChange, HeatChangeOld, HeatDexc, HeatDexc_old, hmi, ioQQQ, ipElec_H2_energy_sort, ipEnergySort, ipoint(), ipRot_H2_energy_sort, ipVib_H2_energy_sort, label, lgEnabled, lgEvaluated, lgLTE, LIM_H2_POP_LOOP, loop_h2_oscil, MAX2, mole, mole_H2_form(), mole_H2_LTE(), mole_update_species_cache(), N_ELEC, n_trace_final, n_trace_full, n_trace_iterations, nCall_this_iteration, nCall_this_zone, nH2_pops, nH2_zone, nLevels_per_elec, nTRACE, nXLevelsMatrix, nzone, nzone_nlevel_set, nzoneEval, ortho_density, ortho_density_f, ortho_para_current, ortho_para_old, ortho_para_older, para_density, para_density_f, photodissoc_BigH2_H2g, photodissoc_BigH2_H2s, phycon, pops_per_elec, pops_per_vib, PRT_POPS, rad_end, rel_pop_LTE_g, rel_pop_LTE_s, rfield, safe_div(), SDIV(), set_xIntensity(), sexp(), SMALLFLOAT, Solomon_dissoc_rate_g, Solomon_dissoc_rate_s, SolveExcitedElectronicLevels(), SolveSomeGroundElectronicLevels(), sp, states, TeUsedColl, thermal, TorF(), TotalInsanity(), trace, and trans.

Here is the call graph for this function:

◆ H2_LinesAdd()

void diatomics::H2_LinesAdd ( void )

add in explicit lines from the large H2 molecule, called by lines_molecules

References DEBUG_ENTRY, H2_SaveLine, label, lgEnabled, LineSave, nElecLevelOutput, PutLine(), rad_end, radius, and trans.

Here is the call graph for this function:

◆ H2_LineZero()

void diatomics::H2_LineZero ( void )

initialize optical depths in H2, called from RT_tau_init

References DEBUG_ENTRY, lgEnabled, rad_end, and trans.

◆ H2_ParseSave()

void diatomics::H2_ParseSave	(	Parser &	p,
		ostringstream &	chHeader )

H2_ParseSave parse the save h2 command

References cdEXIT, chH2ColliderLabels, DEBUG_ENTRY, EXIT_FAILURE, fprintf(), Parser::getNumberDefault(), Parser::getNumberDefaultNegImplLog(), ioQQQ, N_X_COLLIDER, nElecLevelOutput, Parser::nMatch(), save, sncatf(), and thresh_punline_h2.

Here is the call graph for this function:

◆ H2_Prt_column_density()

void diatomics::H2_Prt_column_density ( FILE * ioMEAN )

H2_Prt_column_density print H2 info into zone results, called from prtzone for each printed zone

Parameters

*ioMEAN this is stream used for io, is stdout when called by final, is save unit when save output generated

References DEBUG_ENTRY, fprintf(), H2_X_colden, lgEnabled, nCall_this_zone, ortho_colden, para_colden, and SDIV().

Here is the call graph for this function:

◆ H2_Prt_Zone()

void diatomics::H2_Prt_Zone ( void )

H2_Prt_Zone print H2 info into zone results, called from prtzone for each printed zone

References DEBUG_ENTRY, dense_total, fprintf(), ioQQQ, ipEnergySort, label, lgEnabled, nCall_this_zone, ortho_density, para_density, pops_per_vib, PrintEfmt, SDIV(), and states.

Here is the call graph for this function:

◆ H2_PrtDepartCoef()

void diatomics::H2_PrtDepartCoef ( void )

References DEBUG_ENTRY, depart, fprintf(), ioQQQ, ipEnergySort, Jlowest, label, lgEnabled, n_elec_states, nCall_this_zone, nRot_hi, and nVib_hi.

Here is the call graph for this function:

◆ H2_Punch_line_data()

void diatomics::H2_Punch_line_data	(	FILE *	ioPUN,
		bool	lgDoAll )

save H2 line data

Parameters

ioPUN	io unit for save
lgDoAll	save all levels if true, only subset if false

References cdEXIT, DEBUG_ENTRY, EXIT_FAILURE, fprintf(), ioQQQ, lgEnabled, Save1LineData(), and trans.

Here is the call graph for this function:

◆ H2_PunchDo()

void diatomics::H2_PunchDo	(	FILE *	io,
		char	chJOB[],
		const char	chTime[],
		long int	ipPun )

save some properties of the large H2 molecule

Parameters

io
chJOB[]
chTime[]
ipPun

References Abund(), ASSERT, EmissionProxy::Aul(), average_energy_g, average_energy_s, chMolBranch(), colden, CollRateCoeff, column(), DEBUG_ENTRY, dense, dense_total, energy(), findspecieslocal(), fprintf(), frac(), gv, H2_col_rate_in, H2_col_rate_out, H2_den_g, H2_den_s, H2_lgOrtho, H2_old_populations, H2_populations_LTE, H2_rad_rate_in, H2_rad_rate_out, H2_renorm_chemistry, H2_SaveLine, H2_X_colden, H2_X_colden_LTE, HeatDexc, HeatDiss, hmi, hyperfine, ipEnergySort, ipHYDROGEN, ipPun, ipRot_H2_energy_sort, ipTransitionSort, ipVib_H2_energy_sort, Jlowest, label, lgEnabled, lgEvaluated, lgH2_radiative, LineSave, LTE_Cooling_per_H2(), MAX2, MIN2, mole, N_X_COLLIDER, nCall_this_zone, nElecLevelOutput, nLevels_per_elec, nRot_hi, nVib_hi, ortho_colden, ortho_density, para_colden, para_density, phycon, EmissionProxy::Ploss(), pops_per_elec, rad_end, radius, rfield, save, SDIV(), secondaries, SMALLFLOAT, Solomon_dissoc_rate_g, Solomon_dissoc_rate_s, Solomon_elec_decay_g, Solomon_elec_decay_s, sp, sp_star, spsort(), states, thermal, thresh_punline_h2, and trans.

Here is the call graph for this function:

◆ H2_PunchLineStuff()

void diatomics::H2_PunchLineStuff	(	FILE *	io,
		realnum	xLimit,
		long	index )

include H2 lines in punched optical depths, etc, called from SaveLineStuff

Parameters

io
xLimit
index

References DEBUG_ENTRY, dense, GetDopplerWidth(), ipHYDROGEN, lgEnabled, Save1Line(), and trans.

Here is the call graph for this function:

◆ H2_RadPress()

double diatomics::H2_RadPress ( void )

rad pre due to h2 lines called in PresTotCurrent

References ASSERT, DEBUG_ENTRY, fprintf(), GetDopplerWidth(), ioQQQ, lgEnabled, mass_amu, n_trace_full, nCall_this_zone, nTRACE, PressureRadiationLine(), rad_end, SMALLFLOAT, and trans.

Here is the call graph for this function:

◆ H2_Read_hminus_distribution()

void diatomics::H2_Read_hminus_distribution ( void )

H2_Read_hminus_distribution read distribution function for H2 population following formation from H minus

References ASSERT, BadRead(), cdEXIT, cpu, DEBUG_ENTRY, EXIT_FAILURE, exp10(), FFmtRead(), fprintf(), H2_X_hminus_formation_distribution, ioQQQ, Jlowest, nRot_hi, nTE_HMINUS, nVib_hi, open_data(), path, and read_whole_line().

Referenced by init().

Here is the call graph for this function:

◆ H2_Read_LTE_cooling_per_H2()

void diatomics::H2_Read_LTE_cooling_per_H2 ( )

References cdEXIT, cpu, DEBUG_ENTRY, EXIT_SUCCESS, fprintf(), LTE_cool, LTE_Temp, open_data(), path, read_whole_line(), and TotalInsanity().

Here is the call graph for this function:

◆ H2_ReadDissocEnergies()

void diatomics::H2_ReadDissocEnergies ( void )

References ASSERT, cdEXIT, cpu, DEBUG_ENTRY, EXIT_FAILURE, FFmtRead(), fprintf(), H2_DissocEnergies, ioQQQ, N_ELEC, open_data(), path, and read_whole_line().

Referenced by init().

Here is the call graph for this function:

◆ H2_ReadDissprob()

void diatomics::H2_ReadDissprob ( long int nelec )

read dissociation probabilities and kinetic energies for all electronic levels

Parameters

nelec

References ASSERT, cdEXIT, cpu, DEBUG_ENTRY, EXIT_FAILURE, FFmtRead(), fprintf(), H2_disske, H2_dissprob, ioQQQ, Jlowest, N_ELEC, nRot_hi, nVib_hi, open_data(), path, and read_whole_line().

Referenced by init().

Here is the call graph for this function:

◆ H2_ReadEnergies() [1/2]

void diatomics::H2_ReadEnergies ( )

References ASSERT, DEBUG_ENTRY, H2_ReadEnergies(), Jlowest, levels(), MAX2, n_elec_states, nLevels_per_elec, nRot_hi, nVib_hi, sp, and states.

Referenced by H2_ReadEnergies(), and init().

Here is the call graph for this function:

◆ H2_ReadEnergies() [2/2]

void diatomics::H2_ReadEnergies	(	long int	nelec,
		vector< int > &	n,
		vector< int > &	v,
		vector< int > &	J,
		vector< double > &	eWN )

References ASSERT, cdEXIT, cpu, DEBUG_ENTRY, EXIT_FAILURE, FFmtRead(), fprintf(), ioQQQ, Jlowest, MAX2, MIN2, N_ELEC, nLevels_per_elec, nVib_hi, open_data(), path, and read_whole_line().

Here is the call graph for this function:

◆ H2_ReadTransprob()

void diatomics::H2_ReadTransprob	(	long int	nelec,
		TransitionList &	trans )

read transition probabilities

Parameters

nelec

References ASSERT, cdEXIT, cpu, DEBUG_ENTRY, TransitionList::Emis(), EXIT_FAILURE, FFmtRead(), fprintf(), ioQQQ, ipEnergySort, ipTransitionSort, lgH2_radiative, N_ELEC, nRot_hi, nVib_hi, open_data(), path, read_whole_line(), ShowMe(), and states.

Referenced by init().

Here is the call graph for this function:

◆ H2_Reset()

void diatomics::H2_Reset ( void )

H2_Reset called to reset variables that are needed after an iteration

References DEBUG_ENTRY, fnzone, fprintf(), H2_SaveLine, ioQQQ, iteration, iterationAsEval, label, lgEvaluated, lgPrtMatrix, n_elec_states, nCall_this_iteration, nCall_this_zone, nElecLevelOutput, nH2_pops, nH2_zone, nTRACE, nzone_nlevel_set, nzoneAsEval, prt, renorm_max, renorm_min, TeUsedBoltz, and TeUsedColl.

Here is the call graph for this function:

◆ H2_RT_diffuse()

void diatomics::H2_RT_diffuse ( void )

do emission from H2 - called from RT_diffuse

References DEBUG_ENTRY, lgEnabled, nCall_this_zone, rad_end, and trans.

◆ H2_RT_OTS()

void diatomics::H2_RT_OTS ( void )

H2_RT_OTS - add H2 ots fields

References DEBUG_ENTRY, lgEnabled, nCall_this_zone, rad_end, RT_OTS_AddLine(), and trans.

Here is the call graph for this function:

◆ H2_RT_tau_inc()

void diatomics::H2_RT_tau_inc ( void )

increment optical depth for the H2 molecule, called from RT_tau_inc

References ASSERT, DEBUG_ENTRY, GetDopplerWidth(), H2_renorm_chemistry, lgEnabled, mass_amu, MAX2, MIN2, nCall_this_iteration, nzone, rad_end, renorm_max, renorm_min, RT_line_one_tauinc(), and trans.

Here is the call graph for this function:

◆ H2_RT_tau_reset()

void diatomics::H2_RT_tau_reset ( void )

the large H2 molecule, called from RT_tau_reset

References DEBUG_ENTRY, lgEnabled, rad_end, RT_line_one_tau_reset(), and trans.

Here is the call graph for this function:

◆ H2_RTMake()

void diatomics::H2_RTMake ( linefunc line_one )

do RT for H2 lines

References DEBUG_ENTRY, GetDopplerWidth(), lgEnabled, mass_amu, rad_end, and trans.

Here is the call graph for this function:

◆ H2_Solomon_rate()

void diatomics::H2_Solomon_rate ( void )

H2_Solomon_rate find rates between H2s and H2g and other levels, for use in the chemistry

References DEBUG_ENTRY, dense_total, ENERGY_H2_STAR, GetExcitedElecDensity(), H2_den_g, H2_den_s, H2_dissprob, H2_rad_rate_out, hmi, rad_end, SDIV(), SMALLFLOAT, Solomon_dissoc_rate_g, Solomon_dissoc_rate_s, Solomon_elec_decay_g, Solomon_elec_decay_s, and trans.

Referenced by H2_LevelPops().

Here is the call graph for this function:

◆ H2_X_coll_rate_evaluate()

void diatomics::H2_X_coll_rate_evaluate ( void )

References ASSERT, collider_density, CollRateCoeff, DEBUG_ENTRY, molezone::den, dense, findspecieslocal(), fp_equal(), fprintf(), h2, H2_X_coll_rate, hmi, ioQQQ, ipHELIUM, ipHYDROGEN, lgColl_deexec_Calc, n_trace_full, N_X_COLLIDER, nLevels_per_elec, and nTRACE.

Referenced by H2_LevelPops().

Here is the call graph for this function:

◆ H2_X_sink_and_source()

void diatomics::H2_X_sink_and_source ( void )

References ASSERT, collider_density, collider_density_total_not_H2, Cont_Dissoc_Rate, DEBUG_ENTRY, dense, dense_total, energy(), ENERGY_H2_STAR, fixit, H2_coll_dissoc_rate_coef, H2_coll_dissoc_rate_coef_H2, H2_den_s, H2_DISS_ALLISON_DALGARNO, H2_ipPhoto, H2_populations_LTE, H2_X_formation, H2_X_Hmin_back, H2_X_sink, H2_X_source, hmi, ipRot_H2_energy_sort, ipVib_H2_energy_sort, lgColl_deexec_Calc, mole, mole_global, nLevels_per_elec, rfield, SDIV(), sp, sp_star, spon_diss_tot, and states.

Referenced by H2_LevelPops().

Here is the call graph for this function:

◆ H2_zero_pops_too_low()

void diatomics::H2_zero_pops_too_low ( void )

H2_zero_pops_too_low - zero out some H2 variables if we decide not to compute the full sim, called by H2_LevelPops

References DEBUG_ENTRY, H2_old_populations, H2_populations_LTE, HeatDexc, HeatDexc_deriv, HeatDiss, N_ELEC, photodissoc_BigH2_H2g, photodissoc_BigH2_H2s, pops_per_elec, pops_per_vib, rad_end, Solomon_dissoc_rate_g, Solomon_dissoc_rate_s, states, and trans.

Referenced by H2_LevelPops().

◆ init()

void diatomics::init ( void )

create H2 molecules, called by ContCreatePointers

Todo: 1 add this as a Lya excitation process

Todo: 2 put supra thermal excitation into excitation of electronic bands

References abscf(), AllTransitions, ASSERT, TransitionList::begin(), BIGFLOAT, cdEXIT, CollRateCoeff, CollRateErrFac, compareEmis(), DEBUG_ENER, DEBUG_ENTRY, EH2_eval(), energy(), ENERGY_H2_STAR, energy_off, EXIT_FAILURE, EXIT_SUCCESS, exp10f(), findspecies(), fixit, fprintf(), GetGF(), h2, H2_col_rate_in, H2_col_rate_out, H2_coll_dissoc_rate_coef, H2_coll_dissoc_rate_coef_H2, H2_CollidRateRead(), H2_disske, H2_DissocEnergies, H2_dissprob, H2_ipPhoto, H2_lgOrtho, H2_old_populations, H2_populations_LTE, H2_rad_rate_in, H2_rad_rate_out, H2_Read_hminus_distribution(), H2_ReadDissocEnergies(), H2_ReadDissprob(), H2_ReadEnergies(), H2_ReadTransprob(), H2_SaveLine, H2_TOP, H2_vib_dist(), H2_X_colden, H2_X_colden_LTE, H2_X_coll_rate, H2_X_formation, H2_X_grain_formation_distribution, H2_X_Hmin_back, H2_X_hminus_formation_distribution, H2_X_rate_from_elec_excited, H2_X_rate_to_elec_excited, H2_X_sink, H2_X_source, hd, hmi, ioQQQ, ipCRDW, ipElec_H2_energy_sort, ipEnergySort, ipH2, ipH2p, ipH_LIKE, ipHYDROGEN, ipNCOLLIDER, ipoint(), ipRot_H2_energy_sort, ipTransitionSort, ipVib_H2_energy_sort, is_odd(), iso_sp, Jlowest, label, lgEnabled, lgH2_NOISE, lgH2_radiative, lgRadiative(), lgREAD_DATA, mass_amu, mole, mole_global, N_ELEC, n_elec_states, N_X_COLLIDER, nEner_H2_ground, nLevels_per_elec, nRot_hi, nTE_HMINUS, null_mole, nVib_hi, nXLevelsMatrix, opac, pops_per_vib, POW2, rad_end, ran, RateCoefTable, Read_Mol_Diss_cross_sections(), TransitionList::resize(), sexp(), shortlabel, sp, sp_star, states, TransitionList::states(), TotalInsanity(), trans, wn2angVac(), xMeanNoise, and xSTDNoise.

Here is the call graph for this function:

◆ interpolate_LTE_Cooling()

double diatomics::interpolate_LTE_Cooling ( double Temp )

References DEBUG_ENTRY, LTE_cool, and LTE_Temp.

◆ LTE_Cooling_per_H2()

double diatomics::LTE_Cooling_per_H2 ( )

LTE_Cooling_per_H2 compute exact cooling at LTE, currently for transitions within the X band only.

References DEBUG_ENTRY, H2_populations_LTE, iterations, mole_H2_LTE(), and trans.

Referenced by H2_PunchDo().

Here is the call graph for this function:

◆ Mol_Photo_Diss_Rates()

void diatomics::Mol_Photo_Diss_Rates ( void )

References ASSERT, Cont_Dissoc_Rate, Cont_Dissoc_Rate_H2g, Cont_Dissoc_Rate_H2s, DEBUG_ENTRY, Diss_Trans, ENERGY_H2_STAR, GetDissociationRateCoeff(), H2_den_g, H2_den_s, hmi, ipEnergySort, lgEnabled, mole_global, SDIV(), and states.

Referenced by Cont_Diss_Heat_Rate().

Here is the call graph for this function:

◆ MolDissocOpacity()

double diatomics::MolDissocOpacity	(	const diss_tran &	tran,
		const double &	Mol_Ene )

References cross_section(), DEBUG_ENTRY, diss_tran::initial, ipEnergySort, diss_level::j, MolDissocCrossSection(), diss_level::n, states, and diss_level::v.

Here is the call graph for this function:

◆ mole_H2_form()

void diatomics::mole_H2_form ( void )

mole_H2_form find state specific rates grains and H- form H2

References ASSERT, DEBUG_ENTRY, molezone::den, dense, deut, energy(), ENERGY_H2_STAR, findspecieslocal(), fixit, fp_equal_tol(), gv, h2, H2_logte_hminus, H2_X_formation, H2_X_grain_formation_distribution, H2_X_Hmin_back, H2_X_hminus_formation_distribution, hd, hmi, ipEnergySort, ipH2, ipHYDROGEN, Jlowest, mole, nRot_hi, nTE_HMINUS, nVib_hi, phycon, and states.

Referenced by H2_LevelPops().

Here is the call graph for this function:

◆ mole_H2_LTE()

void diatomics::mole_H2_LTE ( void )

mole_H2_LTE sets Boltzmann factors and LTE unit population of large H2 molecular

References ASSERT, DEBUG_ENTRY, dsexp(), fp_equal(), fp_equal_tol(), fprintf(), H2_populations_LTE, ioQQQ, n_trace_full, nTRACE, phycon, states, and TeUsedBoltz.

Referenced by H2_LevelPops(), and LTE_Cooling_per_H2().

Here is the call graph for this function:

◆ OpacityCreate()

long diatomics::OpacityCreate ( vector< double > & stack )

References ASSERT, DEBUG_ENTRY, ip_photo_opac_thresh, NULL, photoion_opacity_fun, and rfield.

◆ Read_Mol_Diss_cross_sections()

void diatomics::Read_Mol_Diss_cross_sections ( void )

References ASSERT, Cont_Dissoc_Rate, cpu, DEBUG_ENTRY, Diss_Trans, energy(), lgEnabled, n_elec_states, nRot_hi, nVib_hi, open_data(), and read_whole_line().

Referenced by init().

Here is the call graph for this function:

◆ set_numLevelsMatrix()

void diatomics::set_numLevelsMatrix ( long numLevels )

References lgREAD_DATA, and nXLevelsMatrix.

Referenced by ParseDatabaseH2().

◆ SolveExcitedElectronicLevels()

void diatomics::SolveExcitedElectronicLevels ( void )

References multi_arr< T, d, ALLOC, lgBC >::alloc(), ASSERT, DEBUG_ENTRY, dense_total, fixit, fprintf(), H2_den_s, H2_dissprob, H2_lgOrtho, H2_old_populations, H2_rad_rate_in, H2_rad_rate_out, H2_X_rate_from_elec_excited, H2_X_rate_to_elec_excited, hmi, ioQQQ, n_elec_states, n_trace_full, nTRACE, pops_per_elec, pops_per_vib, rad_end, SDIV(), secondaries, spon_diss_tot, states, trans, and multi_arr< T, d, ALLOC, lgBC >::zero().

Referenced by H2_LevelPops().

Here is the call graph for this function:

◆ SolveSomeGroundElectronicLevels()

void diatomics::SolveSomeGroundElectronicLevels ( void )

References ASSERT, DEBUG_ENTRY, dsexp(), energy(), fp_equal(), H2_col_rate_in, H2_col_rate_out, H2_lgOrtho, H2_old_populations, H2_rad_rate_in, H2_rad_rate_out, H2_X_coll_rate, H2_X_sink, H2_X_source, hmi, ipElec_H2_energy_sort, ipRot_H2_energy_sort, ipTransitionSort, ipVib_H2_energy_sort, nLevels_per_elec, nXLevelsMatrix, phycon, rate_grain_J1_to_J0, rate_grain_op_conserve, SDIV(), secondaries, states, and trans.

Referenced by H2_LevelPops().

Here is the call graph for this function:

Member Data Documentation

◆ AulDest

multi_arr<double,2> diatomics::AulDest

private

Referenced by H2_Level_low_matrix().

◆ AulEscp

multi_arr<double,2> diatomics::AulEscp

private

Referenced by H2_Level_low_matrix().

◆ AulPump

multi_arr<double,2> diatomics::AulPump

private

Referenced by H2_Level_low_matrix().

◆ Average_A

double diatomics::Average_A

Average Einstein A for H2s to H2g transition

Referenced by diatomics(), and H2_Calc_Average_Rates().

◆ Average_collH2_deexcit

double diatomics::Average_collH2_deexcit

Average noreactive collisional rate for H2s to H2g transition

Referenced by diatomics(), and H2_Calc_Average_Rates().

◆ Average_collH2_dissoc_g

double diatomics::Average_collH2_dissoc_g

Referenced by diatomics(), and H2_LevelPops().

◆ Average_collH2_dissoc_s

double diatomics::Average_collH2_dissoc_s

Referenced by diatomics(), and H2_LevelPops().

◆ Average_collH2_excit

double diatomics::Average_collH2_excit

Referenced by diatomics(), and H2_Calc_Average_Rates().

◆ Average_collH_deexcit

double diatomics::Average_collH_deexcit

Referenced by diatomics(), and H2_Calc_Average_Rates().

◆ Average_collH_dissoc_g

double diatomics::Average_collH_dissoc_g

Average collisional dissociation of H2g and H2s by H and H2

Referenced by diatomics(), and H2_LevelPops().

◆ Average_collH_dissoc_s

double diatomics::Average_collH_dissoc_s

Referenced by diatomics(), and H2_LevelPops().

◆ Average_collH_excit

double diatomics::Average_collH_excit

Referenced by diatomics(), and H2_Calc_Average_Rates().

◆ average_energy_g

double diatomics::average_energy_g

average energy level of H2g and H2s

Referenced by H2_LevelPops(), and H2_PunchDo().

◆ average_energy_s

double diatomics::average_energy_s

Referenced by H2_LevelPops(), and H2_PunchDo().

◆ chH2ColliderLabels

char diatomics::chH2ColliderLabels[N_X_COLLIDER][chN_X_COLLIDER]

private

Referenced by diatomics(), and H2_ParseSave().

◆ coll_source

t_coll_source diatomics::coll_source[N_X_COLLIDER]

Referenced by H2_CollidRateRead(), and ParseDatabaseH2().

◆ CollRate_levn

multi_arr<double,2> diatomics::CollRate_levn

private

Referenced by H2_Level_low_matrix().

◆ CollRateCoeff

multi_arr<realnum,3> diatomics::CollRateCoeff

private

Referenced by H2_Calc_Average_Rates(), H2_CollidRateEvalAll(), H2_Cooling(), H2_PunchDo(), H2_X_coll_rate_evaluate(), and init().

◆ CollRateErrFac

multi_arr<realnum,3> diatomics::CollRateErrFac

private

Referenced by H2_CollidRateEvalOne(), and init().

◆ Cont_Dissoc_Rate

multi_arr<double,3> diatomics::Cont_Dissoc_Rate

Referenced by H2_X_sink_and_source(), Mol_Photo_Diss_Rates(), and Read_Mol_Diss_cross_sections().

◆ Cont_Dissoc_Rate_H2g

double diatomics::Cont_Dissoc_Rate_H2g

Referenced by Mol_Photo_Diss_Rates().

◆ Cont_Dissoc_Rate_H2s

double diatomics::Cont_Dissoc_Rate_H2s

H2 continuum photodissociation rate coefficient (not scaled by density) from P.C. Stancil data

Referenced by Mol_Photo_Diss_Rates().

◆ create

vector<double> diatomics::create

private

Referenced by H2_Level_low_matrix().

◆ dense_total

const double* const diatomics::dense_total

private

Referenced by Abund(), diatomics(), H2_Colden(), H2_Level_low_matrix(), H2_LevelPops(), H2_Prt_Zone(), H2_PunchDo(), H2_Solomon_rate(), H2_X_sink_and_source(), and SolveExcitedElectronicLevels().

◆ depart

vector<double> diatomics::depart

private

Referenced by H2_Level_low_matrix(), and H2_PrtDepartCoef().

◆ destroy

vector<double> diatomics::destroy

private

Referenced by H2_Level_low_matrix().

◆ Diss_Trans

vector< diss_tran > diatomics::Diss_Trans

Referenced by Cont_Diss_Heat_Rate(), Mol_Photo_Diss_Rates(), and Read_Mol_Diss_cross_sections().

◆ ENERGY_H2_STAR

const double diatomics::ENERGY_H2_STAR

this is the energy (in cm-1), above which levels are considered to be H2*, and below which they are H2

Referenced by diatomics(), H2_Calc_Average_Rates(), H2_LevelPops(), H2_Solomon_rate(), H2_X_sink_and_source(), init(), Mol_Photo_Diss_Rates(), and mole_H2_form().

◆ excit

vector<double> diatomics::excit

private

Referenced by H2_Level_low_matrix().

◆ frac_matrix

double diatomics::frac_matrix

Referenced by H2_LevelPops().

◆ H2_col_rate_in

multi_arr<double,2> diatomics::H2_col_rate_in

private

Referenced by H2_PunchDo(), init(), and SolveSomeGroundElectronicLevels().

◆ H2_col_rate_out

multi_arr<double,2> diatomics::H2_col_rate_out

private

Referenced by H2_PunchDo(), init(), and SolveSomeGroundElectronicLevels().

◆ H2_coll_dissoc_rate_coef

multi_arr<realnum,2> diatomics::H2_coll_dissoc_rate_coef

private

save rate coef (cm3 s-1) for collisional dissociation

Referenced by H2_CollidRateEvalAll(), H2_LevelPops(), H2_X_sink_and_source(), and init().

◆ H2_coll_dissoc_rate_coef_H2

multi_arr<realnum,2> diatomics::H2_coll_dissoc_rate_coef_H2

private

save rate coef (cm3 s-1) for collisional dissociation with H2g and H2s

Referenced by H2_CollidRateEvalAll(), H2_LevelPops(), H2_X_sink_and_source(), and init().

◆ H2_den_g

double diatomics::H2_den_g

private

Referenced by gs_rate(), H2_LevelPops(), H2_PunchDo(), H2_Solomon_rate(), and Mol_Photo_Diss_Rates().

◆ H2_den_s

double diatomics::H2_den_s

private

density of H2s and H2g during current iteration

Referenced by H2_LevelPops(), H2_PunchDo(), H2_Solomon_rate(), H2_X_sink_and_source(), Mol_Photo_Diss_Rates(), and SolveExcitedElectronicLevels().

◆ H2_disske

multi_arr<realnum,3> diatomics::H2_disske

private

Referenced by H2_Cooling(), H2_ReadDissprob(), and init().

◆ H2_DissocEnergies

double diatomics::H2_DissocEnergies[N_ELEC]

private

Referenced by H2_CollidRateEvalAll(), H2_LevelPops(), H2_ReadDissocEnergies(), and init().

◆ H2_dissprob

multi_arr<realnum,3> diatomics::H2_dissprob

private

these will mostly become xxx[elec][vib][rot]

Referenced by gs_rate(), H2_Cooling(), H2_ReadDissprob(), H2_Solomon_rate(), init(), and SolveExcitedElectronicLevels().

◆ H2_ipPhoto

multi_arr<int,2> diatomics::H2_ipPhoto

private

Referenced by H2_X_sink_and_source(), and init().

◆ H2_lgOrtho

multi_arr<bool,3> diatomics::H2_lgOrtho

private

this is true if state is para, false if ortho

Referenced by getLine(), H2_Calc_Average_Rates(), H2_CollidRateEvalOne(), H2_Level_low_matrix(), H2_LevelPops(), H2_PunchDo(), init(), SolveExcitedElectronicLevels(), and SolveSomeGroundElectronicLevels().

◆ H2_old_populations

multi_arr<double,3> diatomics::H2_old_populations

private

these will mostly become xxx[elec][vib][rot]

Referenced by H2_LevelPops(), H2_PunchDo(), H2_zero_pops_too_low(), init(), SolveExcitedElectronicLevels(), and SolveSomeGroundElectronicLevels().

◆ H2_populations_LTE

multi_arr<double,3> diatomics::H2_populations_LTE

private

Referenced by H2_Colden(), H2_LevelPops(), H2_PunchDo(), H2_X_sink_and_source(), H2_zero_pops_too_low(), init(), LTE_Cooling_per_H2(), and mole_H2_LTE().

◆ H2_rad_rate_in

multi_arr<double,2> diatomics::H2_rad_rate_in

private

Referenced by H2_LevelPops(), H2_PunchDo(), init(), SolveExcitedElectronicLevels(), and SolveSomeGroundElectronicLevels().

◆ H2_rad_rate_out

multi_arr<double,3> diatomics::H2_rad_rate_out

private

Referenced by gs_rate(), H2_LevelPops(), H2_PunchDo(), H2_Solomon_rate(), init(), SolveExcitedElectronicLevels(), and SolveSomeGroundElectronicLevels().

◆ H2_renorm_chemistry

double diatomics::H2_renorm_chemistry

private

the renorm factor for this H2 to the chemistry - should be unity

Referenced by H2_Cooling(), H2_LevelPops(), H2_PunchDo(), and H2_RT_tau_inc().

◆ H2_SaveLine

multi_arr<realnum,6> diatomics::H2_SaveLine

private

this is array of accumulated line intensities, used for save he lines command

Referenced by getLine(), H2_LinesAdd(), H2_PunchDo(), H2_Reset(), and init().

◆ H2_to_H_limit

double diatomics::H2_to_H_limit

limit to the ratio H2/Htot - if ratio is below this, large atom is not called

Referenced by diatomics(), H2_LevelPops(), and ParseDatabaseH2().

◆ H2_X_colden

multi_arr<realnum,2> diatomics::H2_X_colden

private

column density within X only vib and rot

Referenced by GetXColden(), H2_Colden(), H2_Prt_column_density(), H2_PunchDo(), and init().

◆ H2_X_colden_LTE

multi_arr<realnum,2> diatomics::H2_X_colden_LTE

private

LTE column density within X only vib and rot

Referenced by H2_Colden(), H2_PunchDo(), and init().

◆ H2_X_coll_rate

multi_arr<realnum,2> diatomics::H2_X_coll_rate

private

rate [s-1] for collisions from ihi to ilo

Referenced by H2_Level_low_matrix(), H2_X_coll_rate_evaluate(), init(), and SolveSomeGroundElectronicLevels().

◆ H2_X_formation

multi_arr<realnum,2> diatomics::H2_X_formation

private

formation into specific states within X only vib and rot, includes both H- and H2 routes

Referenced by H2_X_sink_and_source(), init(), and mole_H2_form().

◆ H2_X_grain_formation_distribution

multi_arr<realnum,3> diatomics::H2_X_grain_formation_distribution

private

distribution function for formation on grain surfaces, vib, rot, last dim is grain type

Referenced by init(), and mole_H2_form().

◆ H2_X_Hmin_back

multi_arr<realnum,2> diatomics::H2_X_Hmin_back

private

backwards destruction of v,J levels due to the H- route

Referenced by H2_X_sink_and_source(), init(), and mole_H2_form().

◆ H2_X_hminus_formation_distribution

multi_arr<realnum,3> diatomics::H2_X_hminus_formation_distribution

private

vib, rot, last dim is temperature

Referenced by H2_Read_hminus_distribution(), init(), and mole_H2_form().

◆ H2_X_rate_from_elec_excited

multi_arr<double,2> diatomics::H2_X_rate_from_elec_excited

private

rates [cm-3 s-1] from elec excited states into X only vib and rot

Referenced by H2_Level_low_matrix(), H2_LevelPops(), init(), and SolveExcitedElectronicLevels().

◆ H2_X_rate_to_elec_excited

multi_arr<double,2> diatomics::H2_X_rate_to_elec_excited

private

rates [s-1] to elec excited states from X only vib and rot

Referenced by H2_Level_low_matrix(), H2_LevelPops(), init(), and SolveExcitedElectronicLevels().

◆ H2_X_sink

valarray<realnum> diatomics::H2_X_sink

private

Referenced by H2_Level_low_matrix(), H2_X_sink_and_source(), init(), and SolveSomeGroundElectronicLevels().

◆ H2_X_source

valarray<realnum> diatomics::H2_X_source

private

Referenced by H2_Level_low_matrix(), H2_X_sink_and_source(), init(), and SolveSomeGroundElectronicLevels().

◆ HeatChange

double diatomics::HeatChange

Referenced by diatomics(), and H2_LevelPops().

◆ HeatChangeOld

double diatomics::HeatChangeOld

Referenced by diatomics(), and H2_LevelPops().

◆ HeatDexc

double diatomics::HeatDexc

Referenced by diatomics(), H2_Cooling(), H2_LevelPops(), H2_PunchDo(), and H2_zero_pops_too_low().

◆ HeatDexc_deriv

double diatomics::HeatDexc_deriv

Referenced by diatomics(), H2_Cooling(), and H2_zero_pops_too_low().

◆ HeatDexc_old

double diatomics::HeatDexc_old

Referenced by diatomics(), and H2_LevelPops().

◆ HeatDiss

double diatomics::HeatDiss

Referenced by diatomics(), H2_Cooling(), H2_PunchDo(), and H2_zero_pops_too_low().

◆ ip_photo_opac_offset

long diatomics::ip_photo_opac_offset

Referenced by CalcPhotoionizationRate().

◆ ip_photo_opac_thresh

long diatomics::ip_photo_opac_thresh

index for threshold for photoionization

Referenced by CalcPhotoionizationRate(), and OpacityCreate().

◆ ipElec_H2_energy_sort

valarray<long> diatomics::ipElec_H2_energy_sort

private

Referenced by H2_LevelPops(), init(), and SolveSomeGroundElectronicLevels().

◆ ipEnergySort

multi_arr<long int,3> diatomics::ipEnergySort

private

Referenced by GetDissociationRate(), GetHeatRate(), GetIndices(), getLine(), gs_rate(), H2_Calc_Average_Rates(), H2_LevelPops(), H2_Prt_Zone(), H2_PrtDepartCoef(), H2_PunchDo(), H2_ReadTransprob(), init(), Mol_Photo_Diss_Rates(), MolDissocOpacity(), and mole_H2_form().

◆ ipRot_H2_energy_sort

valarray<long> diatomics::ipRot_H2_energy_sort

private

Referenced by H2_CollidRateEvalAll(), H2_Level_low_matrix(), H2_LevelPops(), H2_PunchDo(), H2_X_sink_and_source(), init(), and SolveSomeGroundElectronicLevels().

◆ ipTransitionSort

multi_arr<long int,2> diatomics::ipTransitionSort

private

Referenced by gs_rate(), H2_Calc_Average_Rates(), H2_Level_low_matrix(), H2_PunchDo(), H2_ReadTransprob(), init(), and SolveSomeGroundElectronicLevels().

◆ ipVib_H2_energy_sort

valarray<long> diatomics::ipVib_H2_energy_sort

private

Referenced by H2_CollidRateEvalAll(), H2_Level_low_matrix(), H2_LevelPops(), H2_PunchDo(), H2_X_sink_and_source(), init(), and SolveSomeGroundElectronicLevels().

◆ iteration_evaluated

long int diatomics::iteration_evaluated

private

Referenced by CalcPhotoionizationRate(), and diatomics().

◆ iterationAsEval

long int diatomics::iterationAsEval

private

Referenced by diatomics(), H2_Level_low_matrix(), and H2_Reset().

◆ Jlowest

long int diatomics::Jlowest[N_ELEC]

private

this gives the first rotational state for each electronic state - J=0 does not exist when Lambda = 1

Referenced by GetIndices(), gs_rate(), H2_PrtDepartCoef(), H2_PunchDo(), H2_Read_hminus_distribution(), H2_ReadDissprob(), H2_ReadEnergies(), H2_ReadEnergies(), init(), and mole_H2_form().

◆ label

string diatomics::label

private

Referenced by diatomics(), H2_CollidRateEvalAll(), H2_ContPoint(), H2_LevelPops(), H2_LinesAdd(), H2_Prt_Zone(), H2_PrtDepartCoef(), H2_PunchDo(), H2_Reset(), and init().

◆ levelAsEval

long int diatomics::levelAsEval

private

Referenced by diatomics(), and H2_Level_low_matrix().

◆ lgColl_deexec_Calc

bool diatomics::lgColl_deexec_Calc

this is option to turn off the calculated collision rates

Referenced by diatomics(), H2_CollidRateEvalOne(), H2_X_coll_rate_evaluate(), H2_X_sink_and_source(), and ParseDatabaseH2().

◆ lgColl_dissoc_coll

bool diatomics::lgColl_dissoc_coll

this is option to turn off guesses of collisional dissociation rates

Referenced by diatomics(), H2_CollidRateEvalAll(), and ParseDatabaseH2().

◆ lgColl_gbar

bool diatomics::lgColl_gbar

this is option to use estimates of the collision rates from g-bar approximations turn mole.lgColl_gbar on/off with atom h2 gbar on off

Referenced by diatomics(), H2_CollidRateEvalOne(), and ParseDatabaseH2().

◆ lgEnabled

◆ lgEvaluated

bool diatomics::lgEvaluated

says whether model has ever been evaluated in this run - if it has not been then use TH85 physics for mole balance and cooling

Referenced by diatomics(), H2_LevelPops(), H2_PunchDo(), and H2_Reset().

◆ lgFirst

bool diatomics::lgFirst

private

Referenced by diatomics(), and H2_Level_low_matrix().

◆ lgH2_grain_deexcitation

bool diatomics::lgH2_grain_deexcitation

Referenced by diatomics(), and ParseDatabaseH2().

◆ lgH2_He_ORNL

bool diatomics::lgH2_He_ORNL

Referenced by diatomics(), and ParseDatabaseH2().

◆ lgH2_NOISE

bool diatomics::lgH2_NOISE

put noise into collision rates

Referenced by diatomics(), H2_CollidRateEvalOne(), init(), and ParseDatabaseH2().

◆ lgH2_NOISECOSMIC

bool diatomics::lgH2_NOISECOSMIC

noise for the CR collisions

Referenced by diatomics(), and ParseDatabaseH2().

◆ lgH2_ORH2_ORNL

bool diatomics::lgH2_ORH2_ORNL

Referenced by diatomics(), and ParseDatabaseH2().

◆ lgH2_ortho_para_coll_on

bool diatomics::lgH2_ortho_para_coll_on

option to turn off ortho-para collisions, command SPECIES H2 COLLISIONS ORTHO PARA OFF

Referenced by diatomics(), H2_CollidRateEvalOne(), and ParseDatabaseH2().

◆ lgH2_PAH2_ORNL

bool diatomics::lgH2_PAH2_ORNL

Referenced by diatomics(), and ParseDatabaseH2().

◆ lgH2_radiative

multi_arr<bool,2> diatomics::lgH2_radiative

private

fully defined array saying whether (true) or not (false) a radiative decay is defined by the standard emission line structure

Referenced by getLine(), gs_rate(), H2_Calc_Average_Rates(), H2_Level_low_matrix(), H2_PunchDo(), H2_ReadTransprob(), and init().

◆ lgImgMatrix

bool diatomics::lgImgMatrix

private

create image for Matrix input to solver

Referenced by H2_Level_low_matrix().

◆ lgLTE

bool diatomics::lgLTE

flag to force LTE level populations, atom H2 LTE

Referenced by diatomics(), H2_LevelPops(), and ParseDatabaseH2().

◆ lgPrtMatrix

bool diatomics::lgPrtMatrix

private

print Matrix input to solver

Referenced by H2_Level_low_matrix(), and H2_Reset().

◆ lgREAD_DATA

bool diatomics::lgREAD_DATA

flag saying whether molecular data have been read in yet

Referenced by diatomics(), init(), ParseDatabaseH2(), and set_numLevelsMatrix().

◆ loop_h2_oscil

long int diatomics::loop_h2_oscil

Referenced by diatomics(), and H2_LevelPops().

◆ LTE_cool

vector<double> diatomics::LTE_cool

private

Referenced by H2_Read_LTE_cooling_per_H2(), and interpolate_LTE_Cooling().

◆ LTE_Temp

vector<double> diatomics::LTE_Temp

private

Referenced by H2_Read_LTE_cooling_per_H2(), and interpolate_LTE_Cooling().

◆ mass_amu

realnum diatomics::mass_amu

Referenced by H2_RadPress(), H2_RT_tau_inc(), H2_RTMake(), and init().

◆ n_elec_states

long int diatomics::n_elec_states

Referenced by diatomics(), GetExcitedElecDensity(), gs_rate(), H2_PrtDepartCoef(), H2_ReadEnergies(), H2_Reset(), init(), ParseDatabaseH2(), Read_Mol_Diss_cross_sections(), and SolveExcitedElectronicLevels().

◆ n_trace_final

int diatomics::n_trace_final

this sets how fine a trace we want for atom h2 trace

Referenced by diatomics(), H2_LevelPops(), and ParseDatabaseH2().

◆ n_trace_full

int diatomics::n_trace_full

Referenced by diatomics(), H2_CollidRateEvalAll(), H2_Cooling(), H2_Level_low_matrix(), H2_LevelPops(), H2_RadPress(), H2_X_coll_rate_evaluate(), mole_H2_LTE(), ParseDatabaseH2(), and SolveExcitedElectronicLevels().

◆ n_trace_iterations

int diatomics::n_trace_iterations

Referenced by diatomics(), H2_LevelPops(), and ParseDatabaseH2().

◆ n_trace_matrix

int diatomics::n_trace_matrix

Referenced by diatomics(), H2_Level_low_matrix(), and ParseDatabaseH2().

◆ nCall_this_iteration

long int diatomics::nCall_this_iteration

private

the number of times the H2 molecules has been called in this iteration. For the very first call we will use lte for the level populations, for later calls use the last solution

Referenced by H2_Cooling(), H2_LevelPops(), H2_Reset(), and H2_RT_tau_inc().

◆ nCall_this_zone

long int diatomics::nCall_this_zone

Referenced by diatomics(), H2_LevelPops(), H2_Prt_column_density(), H2_Prt_Zone(), H2_PrtDepartCoef(), H2_PunchDo(), H2_RadPress(), H2_Reset(), H2_RT_diffuse(), and H2_RT_OTS().

◆ ndim_allocated

long int diatomics::ndim_allocated

private

Referenced by diatomics(), and H2_Level_low_matrix().

◆ nElecLevelOutput

int diatomics::nElecLevelOutput

Referenced by diatomics(), H2_LinesAdd(), H2_ParseSave(), H2_PunchDo(), and H2_Reset().

◆ nEner_H2_ground

long int diatomics::nEner_H2_ground

private

number of levels in H2g

Referenced by gs_rate(), and init().

◆ nH2_pops

long int diatomics::nH2_pops

private

counters used by H2_itrzn to find number of calls of h2 per zone

Referenced by diatomics(), H2_itrzn(), H2_LevelPops(), and H2_Reset().

◆ nH2_zone

long int diatomics::nH2_zone

private

Referenced by diatomics(), H2_itrzn(), H2_LevelPops(), and H2_Reset().

◆ nLevels_per_elec

long int diatomics::nLevels_per_elec[N_ELEC]

private

the number of ro-vib levels in each elec state

Referenced by H2_CollidRateEvalAll(), H2_CollidRateRead(), H2_Cooling(), H2_Level_low_matrix(), H2_LevelPops(), H2_PunchDo(), H2_ReadEnergies(), H2_ReadEnergies(), H2_X_coll_rate_evaluate(), H2_X_sink_and_source(), init(), and SolveSomeGroundElectronicLevels().

◆ nRot_hi

valarray<long> diatomics::nRot_hi[N_ELEC]

private

number of rotation levels within each elec - vib

Referenced by GetIndices(), GetXColden(), gs_rate(), H2_PrtDepartCoef(), H2_PunchDo(), H2_Read_hminus_distribution(), H2_ReadDissprob(), H2_ReadEnergies(), H2_ReadTransprob(), init(), mole_H2_form(), and Read_Mol_Diss_cross_sections().

◆ nTRACE

int diatomics::nTRACE

turn on trace information

Referenced by diatomics(), H2_CollidRateEvalAll(), H2_Cooling(), H2_Level_low_matrix(), H2_LevelPops(), H2_RadPress(), H2_Reset(), H2_X_coll_rate_evaluate(), mole_H2_LTE(), ParseDatabaseH2(), and SolveExcitedElectronicLevels().

◆ nVib_hi

long int diatomics::nVib_hi[N_ELEC]

private

number of vib states within electronic states

Referenced by GetIndices(), GetXColden(), gs_rate(), H2_PrtDepartCoef(), H2_PunchDo(), H2_Read_hminus_distribution(), H2_ReadDissprob(), H2_ReadEnergies(), H2_ReadEnergies(), H2_ReadTransprob(), init(), mole_H2_form(), and Read_Mol_Diss_cross_sections().

◆ nXLevelsMatrix

long int diatomics::nXLevelsMatrix

private

number of levels within X which are done with matrix solver, set with atom h2 matrix command

Referenced by diatomics(), H2_Level_low_matrix(), H2_LevelPops(), init(), set_numLevelsMatrix(), and SolveSomeGroundElectronicLevels().

◆ nzone_eval

long int diatomics::nzone_eval

private

Referenced by CalcPhotoionizationRate(), and diatomics().

◆ nzone_nlevel_set

long int diatomics::nzone_nlevel_set

private

this is used to establish zone number for evaluation of number of levels in matrix

Referenced by diatomics(), H2_LevelPops(), and H2_Reset().

◆ nzoneAsEval

long int diatomics::nzoneAsEval

private

Referenced by diatomics(), H2_Level_low_matrix(), and H2_Reset().

◆ nzoneEval

long int diatomics::nzoneEval

Referenced by diatomics(), and H2_LevelPops().

◆ ortho_colden

double diatomics::ortho_colden

column density in ortho and para H2

Referenced by cdH2_colden(), diatomics(), H2_Prt_column_density(), and H2_PunchDo().

◆ ortho_density

double diatomics::ortho_density

the density (cm-3) of ortho H2

Referenced by H2_LevelPops(), H2_Prt_Zone(), and H2_PunchDo().

◆ ortho_density_f

realnum diatomics::ortho_density_f

Referenced by H2_LevelPops().

◆ ortho_para_current

double diatomics::ortho_para_current

Referenced by diatomics(), and H2_LevelPops().

◆ ortho_para_old

double diatomics::ortho_para_old

Referenced by diatomics(), and H2_LevelPops().

◆ ortho_para_older

double diatomics::ortho_para_older

Referenced by diatomics(), and H2_LevelPops().

◆ para_colden

double diatomics::para_colden

Referenced by cdH2_colden(), diatomics(), H2_Prt_column_density(), and H2_PunchDo().

◆ para_density

double diatomics::para_density

the density (cm-3) of para H2

Referenced by H2_LevelPops(), H2_Prt_Zone(), and H2_PunchDo().

◆ para_density_f

realnum diatomics::para_density_f

Referenced by H2_LevelPops().

◆ path

string diatomics::path

private

Referenced by diatomics(), H2_CollidRateRead(), H2_Read_hminus_distribution(), H2_Read_LTE_cooling_per_H2(), H2_ReadDissocEnergies(), H2_ReadDissprob(), H2_ReadEnergies(), and H2_ReadTransprob().

◆ photo_heat_hard

double diatomics::photo_heat_hard

Referenced by CalcPhotoionizationRate(), and diatomics().

◆ photo_heat_soft

double diatomics::photo_heat_soft

Referenced by CalcPhotoionizationRate(), and diatomics().

◆ photodissoc_BigH2_H2g

double diatomics::photodissoc_BigH2_H2g

Referenced by diatomics(), H2_LevelPops(), and H2_zero_pops_too_low().

◆ photodissoc_BigH2_H2s

double diatomics::photodissoc_BigH2_H2s

Referenced by diatomics(), H2_LevelPops(), and H2_zero_pops_too_low().

◆ photoion_opacity_fun

double(* diatomics::photoion_opacity_fun) (double energy)

Referenced by diatomics(), and OpacityCreate().

◆ photoionize_rate

double diatomics::photoionize_rate

Referenced by CalcPhotoionizationRate().

◆ pops

vector<double> diatomics::pops

private

Referenced by H2_Level_low_matrix().

◆ pops_per_elec

double diatomics::pops_per_elec[N_ELEC]

private

the total population in each elec state

Referenced by GetExcitedElecDensity(), H2_LevelPops(), H2_PunchDo(), H2_zero_pops_too_low(), and SolveExcitedElectronicLevels().

◆ pops_per_vib

multi_arr<double,2> diatomics::pops_per_vib

private

total population in each vib state

Referenced by H2_LevelPops(), H2_Prt_Zone(), H2_zero_pops_too_low(), init(), and SolveExcitedElectronicLevels().

◆ rad_end

TransitionList::iterator diatomics::rad_end

Referenced by H2_Accel(), H2_ContPoint(), H2_LevelPops(), H2_LinesAdd(), H2_LineZero(), H2_PunchDo(), H2_RadPress(), H2_RT_diffuse(), H2_RT_OTS(), H2_RT_tau_inc(), H2_RT_tau_reset(), H2_RTMake(), H2_Solomon_rate(), H2_zero_pops_too_low(), init(), and SolveExcitedElectronicLevels().

◆ rate_grain_J1_to_J0

double diatomics::rate_grain_J1_to_J0

Referenced by diatomics(), H2_Level_low_matrix(), and SolveSomeGroundElectronicLevels().

◆ rate_grain_op_conserve

double diatomics::rate_grain_op_conserve

rate H2 goes from all X into either J=1 (ortho) or (J=0) para on grain surfaces - units s-1

Referenced by diatomics(), H2_Level_low_matrix(), and SolveSomeGroundElectronicLevels().

◆ RateCoefTable

vector<CollRateCoeffArray> diatomics::RateCoefTable

private

Referenced by H2_CollidRateEvalOne(), H2_CollidRateRead(), and init().

◆ rel_pop_LTE_g

double diatomics::rel_pop_LTE_g

LTE pops of g and s used for H- back reactions

Referenced by H2_LevelPops().

◆ rel_pop_LTE_s

double diatomics::rel_pop_LTE_s

Referenced by H2_LevelPops().

◆ renorm_max

double diatomics::renorm_max

Referenced by diatomics(), H2_Reset(), and H2_RT_tau_inc().

◆ renorm_min

double diatomics::renorm_min

Referenced by diatomics(), H2_Reset(), and H2_RT_tau_inc().

◆ shortlabel

string diatomics::shortlabel

private

Referenced by diatomics(), and init().

◆ Solomon_dissoc_rate_g

double diatomics::Solomon_dissoc_rate_g

Referenced by diatomics(), H2_LevelPops(), H2_PunchDo(), H2_Solomon_rate(), and H2_zero_pops_too_low().

◆ Solomon_dissoc_rate_s

double diatomics::Solomon_dissoc_rate_s

Referenced by diatomics(), H2_LevelPops(), H2_PunchDo(), H2_Solomon_rate(), and H2_zero_pops_too_low().

◆ Solomon_elec_decay_g

double diatomics::Solomon_elec_decay_g

these are decay rates from electronic levels into g and s

Referenced by H2_PunchDo(), and H2_Solomon_rate().

◆ Solomon_elec_decay_s

double diatomics::Solomon_elec_decay_s

Referenced by H2_PunchDo(), and H2_Solomon_rate().

◆ sp

molecule* diatomics::sp

Referenced by H2_LevelPops(), H2_PunchDo(), H2_ReadEnergies(), H2_X_sink_and_source(), and init().

◆ sp_star

molecule* diatomics::sp_star

Referenced by H2_PunchDo(), H2_X_sink_and_source(), and init().

◆ spon_diss_tot

double diatomics::spon_diss_tot

Referenced by diatomics(), H2_X_sink_and_source(), and SolveExcitedElectronicLevels().

◆ stat_levn

vector<double> diatomics::stat_levn

private

Referenced by H2_Level_low_matrix().

◆ states

◆ TeUsedBoltz

double diatomics::TeUsedBoltz

Referenced by diatomics(), H2_Reset(), and mole_H2_LTE().

◆ TeUsedColl

double diatomics::TeUsedColl

Referenced by diatomics(), H2_LevelPops(), and H2_Reset().

◆ trans

◆ xMeanNoise

double diatomics::xMeanNoise

std and mean for the noise, log normal distribution

Referenced by init(), and ParseDatabaseH2().

◆ xSTDNoise

double diatomics::xSTDNoise

Referenced by init(), and ParseDatabaseH2().

The documentation for this class was generated from the following files:

Public Member Functions

Public Attributes

Private Attributes

Constructor & Destructor Documentation

◆ diatomics()

Member Function Documentation

◆ Abund()

◆ CalcPhotoionizationRate()

◆ Cont_Diss_Heat_Rate()

◆ GetDissociationRate()

◆ GetExcitedElecDensity()

◆ GetHeatRate()

◆ GetIndices()

◆ getLine()

◆ GetXColden()

◆ gs_rate()

◆ H2_Accel()

◆ H2_Calc_Average_Rates()

◆ H2_Colden()

◆ H2_CollidRateEvalAll()

◆ H2_CollidRateEvalOne()

◆ H2_CollidRateRead()

◆ H2_ContPoint()

◆ H2_Cooling()

◆ H2_DR()

◆ H2_InterEnergy()

◆ H2_itrzn()

◆ H2_Level_low_matrix()

◆ H2_LevelPops()

◆ H2_LinesAdd()

◆ H2_LineZero()

◆ H2_ParseSave()

◆ H2_Prt_column_density()

◆ H2_Prt_Zone()

◆ H2_PrtDepartCoef()

◆ H2_Punch_line_data()

◆ H2_PunchDo()

◆ H2_PunchLineStuff()

◆ H2_RadPress()

◆ H2_Read_hminus_distribution()

◆ H2_Read_LTE_cooling_per_H2()

◆ H2_ReadDissocEnergies()

◆ H2_ReadDissprob()

◆ H2_ReadEnergies() [1/2]

◆ H2_ReadEnergies() [2/2]

◆ H2_ReadTransprob()

◆ H2_Reset()

◆ H2_RT_diffuse()

◆ H2_RT_OTS()

◆ H2_RT_tau_inc()

◆ H2_RT_tau_reset()

◆ H2_RTMake()

◆ H2_Solomon_rate()

◆ H2_X_coll_rate_evaluate()

◆ H2_X_sink_and_source()

◆ H2_zero_pops_too_low()

◆ init()

◆ interpolate_LTE_Cooling()

◆ LTE_Cooling_per_H2()

◆ Mol_Photo_Diss_Rates()

◆ MolDissocOpacity()

◆ mole_H2_form()

◆ mole_H2_LTE()

◆ OpacityCreate()

◆ Read_Mol_Diss_cross_sections()

◆ set_numLevelsMatrix()

◆ SolveExcitedElectronicLevels()

◆ SolveSomeGroundElectronicLevels()

Member Data Documentation

◆ AulDest

◆ AulEscp

◆ AulPump

◆ Average_A

◆ Average_collH2_deexcit

◆ Average_collH2_dissoc_g

◆ Average_collH2_dissoc_s

◆ Average_collH2_excit

◆ Average_collH_deexcit

◆ Average_collH_dissoc_g

◆ Average_collH_dissoc_s