/home66/gary/public_html/cloudy/c08_branch/source/mole_h_step.cpp File Reference

#include "cddefines.h"
#include "physconst.h"
#include "iso.h"
#include "atmdat.h"
#include "grainvar.h"
#include "ionbal.h"
#include "dense.h"
#include "secondaries.h"
#include "mole.h"
#include "mole_co_priv.h"
#include "opacity.h"
#include "rfield.h"
#include "thermal.h"
#include "timesc.h"
#include "trace.h"
#include "phycon.h"
#include "doppvel.h"
#include "thirdparty.h"
#include "gammas.h"
#include "h2.h"
#include "dynamics.h"
#include "conv.h"
#include "radius.h"
#include "hextra.h"
#include "hmi.h"

Include dependency graph for mole_h_step.cpp:

Go to the source code of this file.

Data Structures
struct	Hmole_rate_s
Defines
#define	ABSLIM   1e-12
#define	INTSZ(a)   (sizeof(a)/sizeof(int))
#define	CATALYST   true
#define	USE_MCCALL
#define	FACTOR   2.25
#define	OLD_FRAC   0.0
#define	CO_ON   true
#define	MOLMIN   1
#define	N_H_MAT   (N_H_MOLEC-MOLMIN)
#define	MAT(a, I_, J_)   (*((a)+(I_)*(N_H_MAT)+(J_)))
Typedefs
typedef Hmole_rate_s	reaction
Functions
reaction *	newreaction (int rindex, int *in, int nin, int *out, int nout, int *rate, int nrate)
void	hmole_step (int *nFixup, double *error)

---

Define Documentation

#define ABSLIM   1e-12

Definition at line 50 of file mole_h_step.cpp.

#define CATALYST   true

#define CO_ON   true

#define FACTOR   2.25

#define INTSZ

(

a

)

(sizeof(a)/sizeof(int))

Definition at line 53 of file mole_h_step.cpp.

#define MAT	(	a,
		I_,
		J_		)	(*((a)+(I_)*(N_H_MAT)+(J_)))

#define MOLMIN   1

#define N_H_MAT   (N_H_MOLEC-MOLMIN)

#define OLD_FRAC   0.0

#define USE_MCCALL

---

Typedef Documentation

typedef struct Hmole_rate_s reaction

Definition at line 64 of file mole_h_step.cpp.

---

Function Documentation

void hmole_step	(	int *	nFixup,
		double *	error
	)

Take one Newton step of the chemical network

Parameters:

	*nFixup
	*error

HeatH2Dexc_used is heating due to collisional deexcitation of vib-excited H2 actually used

these are derivative wrt temp for collisional processes within X

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

Todo:

1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

2 process is net source term for H(n=3) states, must be added in

Todo:

2 process is net ionization term for H(n=3) states

Todo:

1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

0 update photoelectric opacity for H2 to include real cross sections and energies. this is not a higher priority because when H2 is formed there can be very little ionizing radiation. this process must be trivial compared with the Solomon process following reference gives cross section for all energies >>refer H2 photo cs Yan, M., Sadeghpour, H.R., & Dalgarno, A., 1998, ApJ, 496, 1044 Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914

Todo:

0 must include heating, Compton ionization

Todo:

1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? reverse of this reaction i not in detailed balance,why? GS

Todo:

2 this rate drives numerical instability in such models as secondary1 and 2.in

Todo:

1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

2 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

2 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

2 must add process H2+ + H- => H2 + H, Dalgarno&Lepp 87

Todo:

2 put in H2+ + gamma => H + H+

Todo:

1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

2 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

2 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

2 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

2 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS

Todo:

0 had to comment following test out - NA change to hmole caused massive prints

Definition at line 153 of file mole_h_step.cpp.

References amat, DEBUG_ENTRY, dense, hmi, t_hmi::Hmolec, ipHYDROGEN, ipiv, t_hmi::lgNoH2Mole, MALLOC, MAXREACTANTS, N_H_MOLEC, Hmole_rate_s::rk, and t_dense::xMolecules.

Referenced by hmole().

reaction* newreaction	(	int	rindex,
		int *	in,
		int	nin,
		int *	out,
		int	nout,
		int *	rate,
		int	nrate
	)

Definition at line 67 of file mole_h_step.cpp.

References ASSERT, Hmole_rate_s::index, MALLOC, MAXPRODUCTS, MAXREACTANTS, Hmole_rate_s::nproducts, and Hmole_rate_s::nrates.

---