#include "cdstd.h"
#include "cddefines.h"
#include "prt.h"
#include "version.h"
#include "hmi.h"
#include "conv.h"
#include "thermal.h"
#include "grainvar.h"
#include "ionbal.h"
#include "opacity.h"
#include "radius.h"
#include "atmdat.h"
#include "trace.h"
#include "deuterium.h"
#include "grains.h"
#include "mole_priv.h"
#include "gammas.h"
#include "mole.h"
#include "freebound.h"
#include "dense.h"
#include "ipoint.h"
#include "phycon.h"
#include "doppvel.h"
#include "rfield.h"
#include "secondaries.h"
#include "hextra.h"
#include "rt_escprob.h"
#include "iso.h"
#include "h2.h"

Include dependency graph for mole_reactions.cpp:

Macros
#define	FLTEQ(a, b) (fabs((a)-(b)) <= 1e-6*fabs((a)+(b)))

Enumerations
enum	{ UDFA =0 }

enum

enum	{ BUFSIZE =256 }

Functions
STATIC void	newreact (const char label[], const char fun[], double a, double b, double c)

STATIC long	parse_reaction (shared_ptr< mole_reaction > &rate, const char label[])

STATIC string	canonicalize_reaction_label (const char label[])

STATIC void	canonicalize_reaction (shared_ptr< mole_reaction > &rate)

STATIC bool	lgReactionTrivial (shared_ptr< mole_reaction > &rate)

STATIC void	register_reaction_vectors (shared_ptr< mole_reaction > rate)

STATIC void	mole_check_reverse_reactions (void)

STATIC double	mole_get_equilibrium_condition (const char buf[])

STATIC double	mole_get_equilibrium_condition (const mole_reaction *const rate)

STATIC double	mole_partition_function (const molecule *const sp)

STATIC void	mole_generate_isotopologue_reactions (string atom_old, string atom_new)

STATIC double	sticking_probability_H_func (double T_gas, double T_grain)

STATIC double	sticking_probability_H_HM79 (double T_gas, double T_grain)

STATIC void	plot_sparsity (void)

STATIC bool	lgReactBalance (const shared_ptr< mole_reaction > &rate)

STATIC void	read_data (const char file[], void(parse)(char s))

STATIC void	parse_base (char *s)

STATIC void	parse_udfa (char *s)

STATIC void	compare_udfa (const shared_ptr< mole_reaction > &rate)

double	frac_H2star_hminus (void)

double	hmirat (double te)

STATIC void	mole_h2_grain_form (void)

STATIC void	mole_h_reactions (void)

void	mole_create_react (void)

void	mole_cmp_num_in_out_reactions ()

STATIC char *	getcsvfield (char **s, char c)

void	mole_update_rks (void)

void	mole_rk_bigchange (void)

mole_reaction *	mole_findrate_s (const char buf[])

Variables
static realnum	albedo = 0.5

static map< formula_species, shared_ptr< udfa_reaction > >	udfatab

Macro Definition Documentation

#define FLTEQ	(	a,
		b
	)	(fabs((a)-(b)) <= 1e-6*fabs((a)+(b)))

Referenced by compare_udfa(), and parse_udfa().

Enumeration Type Documentation

anonymous enum

Enumerator
BUFSIZE

anonymous enum

set UDFA = 1 to make UDFA comparison and stop http://udfa.ajmarkwick.net/index.php

Enumerator
UDFA

Function Documentation

STATIC void canonicalize_reaction ( shared_ptr< mole_reaction > & rate )

References DEBUG_ENTRY, and t_mole_global::sort().

Referenced by canonicalize_reaction_label(), and newreact().

Here is the call graph for this function:

STATIC string canonicalize_reaction_label ( const char label[] )

References canonicalize_reaction(), DEBUG_ENTRY, and parse_reaction().

Referenced by mole_create_react(), mole_findrate_s(), and mole_generate_isotopologue_reactions().

Here is the call graph for this function:

STATIC void compare_udfa ( const shared_ptr< mole_reaction > & rate )

References CONFLICT, CORRECT, FLTEQ, MAXPRODUCTS, and MAXREACTANTS.

Referenced by newreact().

double frac_H2star_hminus ( void )

References h2, t_hmi::H2_forms_hminus, t_hmi::H2star_forms_hminus, hmi, diatomics::lgEnabled, diatomics::lgEvaluated, t_hmi::lgH2_Chemistry_BigH2, and SDIV().

Referenced by mole_h_rate_diagnostics().

Here is the call graph for this function:

STATIC char * getcsvfield	(	char **	s,
		char	c
	)

Referenced by parse_base(), and parse_udfa().

double hmirat ( double te )

hmirat computes radiative association rate for H-

Parameters

te

References DEBUG_ENTRY, phycon, t_phycon::sqrte, t_phycon::te001, t_phycon::te003, t_phycon::te01, t_phycon::te03, t_phycon::te10, t_phycon::te20, and t_phycon::te70.

Referenced by mole_h_reactions().

STATIC bool lgReactBalance ( const shared_ptr< mole_reaction > & rate )

References fprintf(), and sign().

Referenced by newreact().

Here is the call graph for this function:

STATIC bool lgReactionTrivial ( shared_ptr< mole_reaction > & rate )

References DEBUG_ENTRY.

Referenced by newreact().

STATIC void mole_check_reverse_reactions ( void )

References DEBUG_ENTRY, exists(), fixit, fprintf(), ioQQQ, t_trace::lgTraceMole, NULL, mole_priv::reactab, and trace.

Referenced by mole_create_react().

Here is the call graph for this function:

void mole_cmp_num_in_out_reactions ( void )

References CHARS_SPECIES, DEBUG_ENTRY, fprintf(), ioQQQ, t_mole_global::list, mole_global, t_mole_global::num_total, and mole_priv::reactab.

Here is the call graph for this function:

void mole_create_react ( void )

mole_create_react build reaction structures

References canonicalize_reaction_label(), cdEXIT, DEBUG_ENTRY, dense, deut, EXIT_FAILURE, fixit, fprintf(), hmi, ipHELIUM, t_deuterium::lgElmtOn, t_dense::lgElmtOn, t_mole_global::lgFederman, t_hmi::lgLeiden_Keep_ipMH2s, t_mole_global::lgLeidenHack, t_mole_global::lgProtElim, mole, mole_check_reverse_reactions(), mole_generate_isotopologue_reactions(), mole_global, newreact(), nuclide_list, t_mole_global::offReactions, t_mole_local::old_zone, parse_base(), parse_udfa(), plot_sparsity(), pow(), mole_priv::reactab, t_mole_local::reaction_rks, read_data(), register_reaction_vectors(), and UDFA.

Referenced by InitSimPostparse().

Here is the call graph for this function:

mole_reaction* mole_findrate_s ( const char buf[] )

References canonicalize_reaction_label(), DEBUG_ENTRY, NULL, and mole_priv::reactab.

Referenced by t_mole_local::findrate(), t_mole_local::findrk(), mole_get_equilibrium_condition(), and t_mole_local::sink_rate().

Here is the call graph for this function:

STATIC void mole_generate_isotopologue_reactions	(	string	atom_old,
		string	atom_new
	)

References ASSERT, canonicalize_reaction_label(), create_isotopologues(), create_isotopologues_one_position(), DEBUG_ENTRY, findnuclide(), fixit, fprintf(), ioQQQ, newreact(), parse_species_label(), mole_priv::reactab, and TotalInsanity().

Referenced by mole_create_react().

Here is the call graph for this function:

STATIC double mole_get_equilibrium_condition ( const char buf[] )

References DEBUG_ENTRY, and mole_findrate_s().

Here is the call graph for this function:

STATIC double mole_get_equilibrium_condition ( const mole_reaction *const rate )

References BIGFLOAT, DEBUG_ENTRY, MIN2, mole_partition_function(), mole_reaction::nproducts, mole_reaction::nreactants, mole_reaction::products, mole_reaction::reactants, and SEXP_LIMIT.

Here is the call graph for this function:

STATIC void mole_h2_grain_form ( void )

Todo:: 2 - put in actual composition dependent Tad - this is only valid for bare surfaces - not ice - for ice Tad is 555K hmi.Tad is binding energy expressed as a temperature note that hmi.Tad is set to 800. in zero tau_nu the first equation in section 2.5 equation one paragraph before equation 2 at low grain temperatures all end in para, J=0

References ASSERT, t_dense::AtomicWeight, GrainVar::bin, t_hmi::chJura, DEBUG_ENTRY, dense, ENABLE_QUANTUM_HEATING, fixit, fnzone, fprintf(), t_dense::gas_phase, GetAveVelocity(), GrainDrive(), gv, h2, hmi, ioQQQ, ipHYDROGEN, GrainVar::lgDustOn(), GrainVar::lgGrainPhysicsOn, diatomics::lgH2_grain_deexcitation, MAT_CAR, MAT_CAR2, MAT_PAH, MAT_PAH2, MAT_SIC, MAT_SIL, MAT_SIL2, NQGRID, phycon, POW2, qheat(), diatomics::rate_grain_J1_to_J0, diatomics::rate_grain_op_conserve, t_hmi::rate_h2_form_grains_set, GrainVar::rate_h2_form_grains_used_total, t_hmi::ScaleJura, SDIV(), sexp(), t_phycon::sqrte, sticking_probability_H_func(), t_hmi::Tad, and t_phycon::te.