/home66/gary/public_html/cloudy/c10_branch/source/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.

Data Structures

class  my_Integrand_escConE2
class  my_Integrand_conrec

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, realnum DopplerWidth)
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)
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, realnum DopplerWidth)
double RT_DestHummer (double beta)

Variables

static const double SCALE = 2.

Function Documentation

double esc_CRDcore ( double  tau,
double  tout 
)

esc_CRDcore escape probability CRD with no wings, for subordinate lines

Parameters:
tau 
tout 

Definition at line 339 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_escape().

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 285 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_escape().

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

References DEBUG_ENTRY, esca0k2(), and SQRTPI.

Referenced by 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 240 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_escape().

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

References ASSERT, DEBUG_ENTRY, escmase(), MIN2, pow(), 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 392 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 527 of file rt_escprob.cpp.

References my_Integrand_escConE2::chnukt_ContTkt, my_Integrand_conrec::chnukt_ContTkt, my_Integrand_escConE2::chnukt_ctau, DEBUG_ENTRY, and Integrator< Integrand, Method >::sum().

Referenced by RT_recom_effic().

Here is the call graph for this function:

STATIC double escmase ( double  tau  ) 

Definition at line 495 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 450 of file rt_escprob.cpp.

References dBaseLines, DEBUG_ENTRY, DumpLine(), HFLines, transition::Hi, linesAdded2, 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:
beta beta is ratio of continuum to mean line opacity,
Returns:
dest prob = beta F(beta)

Definition at line 939 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:
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 673 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_escape().

Here is the call graph for this function:

double RT_LineWidth ( const transition t,
realnum  DopplerWidth 
)

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

Parameters:
t 
DopplerWidth 

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

References ASSERT, BIGFLOAT, t_emission::damp, t_emission::dampXvel, DEBUG_ENTRY, dense, t_dense::eden, transition::Emis, esc_PRD_1side(), iteration, Wind::lgBallistic(), t_pressure::lgPradDen, MAX2, MIN2, opac, t_emission::Pdest, t_emission::Pelec_esc, PI, pow(), pressure, t_emission::TauIn, t_opac::taumin, t_emission::TauTot, wind, and 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,
realnum  DopplerWidth 
)

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

Parameters:
*esin 
*dest 
abund 
ipZ 0 for H
DopplerWidth 

Definition at line 143 of file rt_escprob.cpp.

References ASSERT, DEBUG_ENTRY, 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_escape().

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

References t_opac::albedo, DEBUG_ENTRY, DEST0, ioQQQ, MAX2, MIN2, opac, pow(), POW2, and SQRTPI.

Referenced by RTesc_lya().

Here is the call graph for this function:


Variable Documentation

const double SCALE = 2. [static]

Definition at line 26 of file rt_escprob.cpp.

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

 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines

Generated on 15 Nov 2012 for cloudy by  doxygen 1.6.1