/home66/gary/public_html/cloudy/c10_branch/source/helike_einsta.cpp File Reference

#include "cddefines.h"
#include "physconst.h"
#include "taulines.h"
#include "dense.h"
#include "trace.h"
#include "hydro_bauman.h"
#include "iso.h"
#include "helike.h"
#include "helike_einsta.h"
#include "hydroeinsta.h"

Include dependency graph for helike_einsta.cpp:

Go to the source code of this file.

Functions
STATIC double	ritoa (long li, long lf, long nelem, double k, double RI2)
STATIC double	ForbiddenAuls (long ipHi, long ipLo, long nelem)
STATIC double	Jint (double theta)
STATIC double	AngerJ (double vv, double zz)
double	scqdri (double nstar, long int l, double npstar, long int lp, double iz)
double	he_1trans (long nelem, double Enerwn, double Eff_nupper, long lHi, long sHi, long jHi, double Eff_nlower, long lLo, long sLo, long jLo)
void	DoFSMixing (long nelem, long ipLoSing, long ipHiSing)
realnum	helike_transprob (long nelem, long ipHi, long ipLo)
void	HelikeTransProbSetup (void)
Variables
static double ***	TransProbs
static double	vJint
static double	zJint

Function Documentation

STATIC double AngerJ	(	double	vv,
		double	zz
	)

Definition at line 58 of file helike_einsta.cpp.

References DEBUG_ENTRY, Jint(), PI, qg32(), vJint, and zJint.

Referenced by scqdri().

Here is the call graph for this function:

void DoFSMixing	(	long	nelem,
		long	ipLoSing,
		long	ipHiSing
	)

Every bit of this routine is based upon the singlet-triplet mixing formalism given in >>refer He FSM Drake, G. W. F. 1996, in Atomic, Molecular, & Optical Physics Handbook, >>refercon ed. G. W. F. Drake (New York: AIP Press). That formalism mixes the levels themselves, but since this code is not J-resolved, we simulate that by mixing only the transition probabilities. We find results comparable to those calculated in the fully J-resolved model spearheaded by Rob Bauman, and described in >>refer He FSM Bauman, R. P., Porter, R. L., Ferland, G. J., & MacAdam, K. B. 2005, ApJ, accepted

Parameters:

	nelem
	ipLoSing
	ipHiSing

Definition at line 810 of file helike_einsta.cpp.

References ASSERT, DEBUG_ENTRY, ioQQQ, ipHE_LIKE, iso, PI, POW2, t_iso::QuantumNumbers2Index, StatesElemNEW, TRANS_PROB_CONST, and Transitions.

Referenced by iso_create().

STATIC double ForbiddenAuls	(	long	ipHi,
		long	ipLo,
		long	nelem
	)

Todo:: 2 find a transition probability for this 2^3P0 - 2^3P1 transition. It will require a bit of trickery to insert into the rate matrix, because of the fact that the lower level has a higher index. See discussion "Energy order within 2 3P" near the top of helike.c

Todo:: 2 find a transition probability for this 2^3P1 - 2^3P2 transition. It will require a bit of trickery to insert into the rate matrix, because of the fact that the lower level has a higher index. See discussion "Energy order within 2 3P" near the top of helike.c

Definition at line 221 of file helike_einsta.cpp.

References ASSERT, DEBUG_ENTRY, fixit(), ipARGON, ipHe1s1S, ipHe2p1P, ipHe2p3P0, ipHe2p3P1, ipHe2p3P2, ipHe2s1S, ipHe2s3S, ipHe3s3S, ipHE_LIKE, ipHELIUM, ipNEON, IPRAD, iso, iso_put_error(), L_, N_, pow(), S_, t_iso::SmallA, and TotalInsanity().

Referenced by he_1trans().

Here is the call graph for this function:

double he_1trans	(	long	nelem,
		double	Enerwn,
		double	Eff_nupper,
		long	lHi,
		long	sHi,
		long	jHi,
		double	Eff_nlower,
		long	lLo,
		long	sLo,
		long	jLo
	)

compute energy diffference in wn and Aul for given line return is 0 for success, 1 for failure

Parameters:

	nelem	charge on the C scale, 1 is helium
	Enerwn	energy difference in wavenumber
	Eff_nupper	upper quantum numbers
	Eff_nlower	lower quantum numbers
	lHi
	sHi
	jHi
	lLo
	sLo
	jLo