cloudy  trunk
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Macros | Functions | Variables
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.

Macros

#define SCALE   2.
 

Functions

STATIC double escmase (double tau)
 
STATIC void RTesc_lya_1side (double taume, double beta, realnum *esc, realnum *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 (const transition *t)
 
double RT_DestHummer (double beta)
 

Variables

static double chnukt_ContTkt
 
static double chnukt_ctau
 

Macro Definition 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)

Definition at line 528 of file rt_escprob.cpp.

References DEBUG_ENTRY.

Referenced by escConE2().

STATIC double conrec ( double  x)

Definition at line 571 of file rt_escprob.cpp.

References chnukt_ContTkt, and DEBUG_ENTRY.

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 315 of file rt_escprob.cpp.

References ASSERT, DEBUG_ENTRY, esca0k2(), t_rt::fracin, iteration, MIN2, 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 261 of file rt_escprob.cpp.

References ASSERT, DEBUG_ENTRY, esc_CRDwing_1side(), t_rt::fracin, iteration, MIN2, 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 95 of file rt_escprob.cpp.

References DEBUG_ENTRY, esca0k2(), and SQRTPI.

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

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 218 of file rt_escprob.cpp.

References ASSERT, DEBUG_ENTRY, esc_PRD_1side(), t_rt::fracin, iteration, 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 ASSERT, DEBUG_ENTRY, escmase(), MIN2, and rt.

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

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 368 of file rt_escprob.cpp.

References DEBUG_ENTRY, escmase(), and SQRTPI.

Referenced by DrvEscP(), esc_CRDcore(), esc_CRDwing_1side(), RT_continuum_shield_fcn(), 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 582 of file rt_escprob.cpp.

References chnukt_ContTkt, chnukt_ctau, conrec(), DEBUG_ENTRY, 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 516 of file rt_escprob.cpp.

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

Referenced by esccon().

Here is the call graph for this function:

STATIC double escmase ( double  tau)

Definition at line 481 of file rt_escprob.cpp.

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

Referenced by esc_PRD_1side(), and esca0k2().

Here is the call graph for this function:

STATIC void FindNeg ( void  )

Definition at line 426 of file rt_escprob.cpp.

References atmolEmis, C12O16Rotate, C13O16Rotate, DEBUG_ENTRY, DumpLine(), HFLines, t_transition::Hi, linesAdded2, nCORotate, t_quantumState::nelem, nHFLines, NISO, nLevel1, nWindLine, TauLine2, and TauLines.

Referenced by escmase().

Here is the call graph for this function:

double RT_DestHummer ( double  beta)
Parameters
betabeta is ratio of continuum to mean line opacity,
Returns
dest prob = beta F(beta)

Definition at line 981 of file rt_escprob.cpp.

References ASSERT, and DEBUG_ENTRY.

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
abundabundance of species
crsecits line absorption cross section
ipanupointer to energy within continuum array, to get background opacity, this is on the f not c scale
widlline width
escpescape probability
nCoretype 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 722 of file rt_escprob.cpp.

References t_opac::albedo, t_rfield::anu, ASSERT, cdEXIT, conv, DEBUG_ENTRY, fp_equal(), 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, POW2, rfield, SQRTPI, and Transitions.

Referenced by RT_line_static().

Here is the call graph for this function:

double RT_LineWidth ( const transition 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 that this checks on size not sign of windv

Definition at line 866 of file rt_escprob.cpp.

References ASSERT, t_emission::damp, DEBUG_ENTRY, dense, t_DoppVel::doppler, DoppVel, t_dense::eden, t_transition::Emis, esc_PRD_1side(), t_transition::Hi, iteration, t_pressure::lgPradDen, MAX2, MIN2, t_quantumState::nelem, opac, t_emission::Pdest, t_emission::Pelec_esc, PI, pressure, t_emission::TauIn, t_opac::taumin, t_emission::TauTot, wind, t_wind::windv, and t_wind::windv0.

Referenced by FeIIRadPress(), PressureRadiationLine(), and PresTotCurrent().

Here is the call graph for this function:

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
ipZ0 for H,

Definition at line 122 of file rt_escprob.cpp.

References ASSERT, DEBUG_ENTRY, t_DoppVel::doppler, DoppVel, t_rt::fracin, ipH1s, ipH2p, ipH_LIKE, LIMELM, MAX2, MIN2, opac, t_opac::opacity_abs, rt, RTesc_lya_1side(), SQRTPI, Transitions, 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,
realnum esc,
realnum dest,
long  ipLine 
)

Definition at line 624 of file rt_escprob.cpp.

References t_opac::albedo, DEBUG_ENTRY, 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().


Generated for cloudy by doxygen 1.8.3.1