rt_escprob.cpp File Reference

#include "cddefines.h"
#include "physconst.h"
#include "dense.h"
#include "conv.h"
#include "rfield.h"
#include "opacity.h"
#include "lines_service.h"
#include "taulines.h"
#include "doppvel.h"
#include "pressure.h"
#include "wind.h"
#include "rt.h"

Include dependency graph for rt_escprob.cpp:

Go to the source code of this file.

Defines

#define SCALE 2.

Functions

static double escmase (double tau)

static void RTesc_lya_1side (double taume, double beta, float *esc, float *dest, long ipLine)

double esc_PRD_1side (double tau, double a)

double esc_CRDwing_1side (double tau, double a)

double RTesc_lya (double *esin, double *dest, double abund, long int nelem)

double esc_PRD (double tau, double tau_out, double damp)

double esc_CRDwing (double tau_in, double tau_out, double damp)

double esc_CRDcore (double tau_in, double tau_out)

double esca0k2 (double taume)

static void FindNeg (void)

static double cone2 (double t)

double escConE2 (double x)

static double conrec (double x)

double esccon (double tau, double hnukt)

double RT_DestProb (double abund, double crsec, long int ipanu, double widl, double escp, int nCore)

double RT_LineWidth (EmLine *t)

double RT_DestHummer (double beta)

Variables

static double chnukt_ContTkt

static double chnukt_ctau

Define Documentation

#define SCALE 2.

Definition at line 15 of file rt_escprob.cpp.

Referenced by esc_CRDcore(), esc_CRDwing(), and esc_PRD().

Function Documentation

static double cone2 ( double t ) [static]

Definition at line 544 of file rt_escprob.cpp.

References DEBUG_ENTRY, and DEBUG_EXIT.

Referenced by escConE2().

static double conrec ( double x ) [static]

Definition at line 591 of file rt_escprob.cpp.

References chnukt_ContTkt, DEBUG_ENTRY, and DEBUG_EXIT.

Referenced by esccon().

double esc_CRDcore	(	double	tau,
		double	tout
	)

esc_CRDcore escape probability CRD with no wings, for subordinate lines

Parameters:

	tau
	tout

Definition at line 329 of file rt_escprob.cpp.

References ASSERT, DEBUG_ENTRY, DEBUG_EXIT, esca0k2(), t_rt::fracin, t_opac::lgTauOutOn, MIN2, opac, rt, SCALE, t_rt::wayin, and t_rt::wayout.

Referenced by RT_line_static().

Here is the call graph for this function:

double esc_CRDwing	(	double	tau,
		double	tout,
		double	damp
	)

esc_CRDwing escape probability CRD with wings, for subordinate lines

Parameters:

	tau
	tout
	damp

Definition at line 273 of file rt_escprob.cpp.

References ASSERT, DEBUG_ENTRY, DEBUG_EXIT, esc_CRDwing_1side(), t_rt::fracin, t_opac::lgTauOutOn, MIN2, opac, rt, SCALE, t_rt::wayin, and t_rt::wayout.

Referenced by RT_line_static().

Here is the call graph for this function:

double esc_CRDwing_1side	(	double	tau,
		double	a
	)

esc_CRDwing_1side fundamental escape probability radiative transfer routine, for complete redistribution

Definition at line 99 of file rt_escprob.cpp.

References DEBUG_ENTRY, DEBUG_EXIT, esca0k2(), and SQRTPI.

Referenced by atom_oi_calc(), DrvEscP(), esc_CRDwing(), RT_continuum_shield_fcn(), and RT_line_static().

Here is the call graph for this function:

double esc_PRD	(	double	tau,
		double	tout,
		double	damp
	)

esc_PRD escape probability radiative transfer for incomplete redistribution

Parameters:

	tau
	tout
	damp

Definition at line 228 of file rt_escprob.cpp.

References ASSERT, DEBUG_ENTRY, DEBUG_EXIT, esc_PRD_1side(), t_rt::fracin, t_opac::lgTauOutOn, opac, rt, SCALE, t_rt::wayin, and t_rt::wayout.

Referenced by RT_line_static().

Here is the call graph for this function:

double esc_PRD_1side	(	double	tau,
		double	a
	)

esc_PRD_1side fundamental escape probability radiative transfer routine for incomplete redistribution

Parameters:

	tau
	a

Definition at line 42 of file rt_escprob.cpp.

References DEBUG_ENTRY, DEBUG_EXIT, escmase(), MIN2, and rt.

Referenced by DrvEscP(), esc_PRD(), hmole_reactions(), lines_lv1_k_zn(), RT_continuum_shield_fcn(), and RT_line_static().

Here is the call graph for this function:

double esca0k2 ( double taume )

esca0k2 derive Hummer's K2 escape probability for Doppler core only

Parameters:

taume

Definition at line 384 of file rt_escprob.cpp.

References DEBUG_ENTRY, DEBUG_EXIT, escmase(), and SQRTPI.

Referenced by DrvEscP(), esc_CRDcore(), esc_CRDwing_1side(), RT_continuum_shield_fcn(), RT_line_static(), and th85rate_co().

Here is the call graph for this function:

double esccon	(	double	tau,
		double	hnukt
	)

escpcn continuum escape probability

Parameters:

	tau
	hnukt

Definition at line 604 of file rt_escprob.cpp.

References chnukt_ContTkt, chnukt_ctau, conrec(), DEBUG_ENTRY, DEBUG_EXIT, escConE2(), and qg32().

Referenced by RT_recom_effic().

Here is the call graph for this function:

double escConE2 ( double x )

conesc one of the forms of the continuum escape probability

Parameters:

x

Definition at line 530 of file rt_escprob.cpp.

References chnukt_ContTkt, chnukt_ctau, cone2(), DEBUG_ENTRY, and DEBUG_EXIT.

Referenced by esccon().

Here is the call graph for this function:

static double escmase ( double tau ) [static]

Definition at line 492 of file rt_escprob.cpp.

References ASSERT, cdEXIT(), DEBUG_ENTRY, DEBUG_EXIT, FindNeg(), ioQQQ, nzone, and ShowMe().

Referenced by esc_PRD_1side(), and esca0k2().

Here is the call graph for this function:

static void FindNeg ( void ) [static]

Definition at line 444 of file rt_escprob.cpp.

References C12O16Rotate, C13O16Rotate, DEBUG_ENTRY, DEBUG_EXIT, DumpLine(), HFLines, nCORotate, nHFLines, NISO, nLevel1, nWindLine, TauLine2, and TauLines.

Referenced by escmase().

Here is the call graph for this function:

double RT_DestHummer ( double beta )

Parameters:

beta

beta is ratio of continuum to mean line opacity,

Returns:: dest prob = beta F(beta)

Definition at line 1157 of file rt_escprob.cpp.

References ASSERT, DEBUG_ENTRY, and DEBUG_EXIT.

double RT_DestProb	(	double	abund,
		double	crsec,
		long int	ipanu,
		double	widl,
		double	escp,
		int	nCore
	)

RT_DestProb returns line destruction probability due to continuum opacity

Parameters:

	abund	abundance of species
	crsec	its line absorption cross section
	ipanu	pointer to energy within continuum array, to get background opacity, this is on the f not c scale
	widl	line width
	escp	escape probability
	nCore	type of redistribution function

Todo:: 2 this min is because there are no calculations that show what to do for beta beyound this value

Definition at line 752 of file rt_escprob.cpp.

References t_opac::albedo, t_rfield::anu, ASSERT, cdEXIT(), conv, DEBUG_ENTRY, DEBUG_EXIT, EmisLines, ioQQQ, ipDEST_INCOM, ipDEST_K2, ipDEST_SIMPL, ipH1s, ipH2p, ipH_LIKE, ipHYDROGEN, MIN2, t_rfield::nflux, t_conv::nTotalIoniz, nzone, opac, t_opac::opacity_abs, EmLine::PopOpc, POW2, rfield, and SQRTPI.

Referenced by RT_line_static().

Here is the call graph for this function:

double RT_LineWidth ( EmLine * t )

RT_LineWidth compute line width (cm/sec), using optical depth array information

Parameters:

t

Todo:: 1 dynamics; this test assumes that neg vel are subsonic, so that sobolev length would overestimate the optical depth, since ion is at most present over computed slab, and possibly more.

Todo:: 1 rewrite so tht this checks on size not sign of windv

Definition at line 1020 of file rt_escprob.cpp.

References ASSERT, EmLine::damp, DEBUG_ENTRY, DEBUG_EXIT, dense, t_DoppVel::doppler, DoppVel, t_dense::eden, t_pressure::lgPradDen, t_opac::lgTauOutOn, MAX2, MIN2, EmLine::nelem, opac, EmLine::Pdest, EmLine::Pelec_esc, EmLine::Pesc, PI, pressure, EmLine::TauIn, t_opac::taumin, EmLine::TauTot, wind, t_wind::windv, and t_wind::windv0.

Referenced by FeIIRadPress(), H2_RadPress(), PresTotCurrent(), and PrtLinePres().

double RTesc_lya	(	double *	esin,
		double *	dest,
		double	abund,
		long int	ipZ
	)

RTesc_lya escape prob for hydrogen atom Lya, using Hummer and Kunasz results

Parameters:

	*esin
	*dest
	abund
	ipZ	0 for H,

Definition at line 128 of file rt_escprob.cpp.

References ASSERT, DEBUG_ENTRY, DEBUG_EXIT, t_DoppVel::doppler, DoppVel, EmisLines, t_rt::fracin, EmLine::FracInwd, EmLine::ipCont, ipH1s, ipH2p, ipH_LIKE, LIMELM, MAX2, MIN2, opac, EmLine::opacity, t_opac::opacity_abs, EmLine::Pdest, EmLine::Pesc, rt, RTesc_lya_1side(), SQRTPI, t_rt::wayin, and t_rt::wayout.

Referenced by RT_line_static().

Here is the call graph for this function:

static void RTesc_lya_1side	(	double	taume,
		double	beta,
		float *	esc,
		float *	dest,
		long	ipLine
	)			`[static]`

Definition at line 652 of file rt_escprob.cpp.

References t_opac::albedo, DEBUG_ENTRY, DEBUG_EXIT, DEST0, ioQQQ, MAX2, MIN2, opac, POW2, and SQRTPI.

Referenced by RTesc_lya().

Variable Documentation

double chnukt_ContTkt [static]

Definition at line 30 of file rt_escprob.cpp.

Referenced by conrec(), esccon(), and escConE2().

double chnukt_ctau [static]

Definition at line 30 of file rt_escprob.cpp.

Referenced by esccon(), and escConE2().


Defines
#define	SCALE 2.
Functions
static double	escmase (double tau)
static void	RTesc_lya_1side (double taume, double beta, float esc, float dest, long ipLine)
double	esc_PRD_1side (double tau, double a)
double	esc_CRDwing_1side (double tau, double a)
double	RTesc_lya (double esin, double dest, double abund, long int nelem)
double	esc_PRD (double tau, double tau_out, double damp)
double	esc_CRDwing (double tau_in, double tau_out, double damp)
double	esc_CRDcore (double tau_in, double tau_out)
double	esca0k2 (double taume)
static void	FindNeg (void)
static double	cone2 (double t)
double	escConE2 (double x)
static double	conrec (double x)
double	esccon (double tau, double hnukt)
double	RT_DestProb (double abund, double crsec, long int ipanu, double widl, double escp, int nCore)
double	RT_LineWidth (EmLine *t)
double	RT_DestHummer (double beta)
Variables
static double	chnukt_ContTkt
static double	chnukt_ctau