#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, 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 |
| #define SCALE 2. |
Definition at line 15 of file rt_escprob.cpp.
Referenced by esc_CRDcore(), esc_CRDwing(), and esc_PRD().
| STATIC double cone2 | ( | double | t | ) |
| 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
| 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
| 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 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
| 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
| 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
| 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
| 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
| 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, DEBUG_ENTRY, DumpLine(), t_transition::Emis, t_transition::Hi, linesAdded2, t_quantumState::nelem, NISO, nLevel1, nWindLine, t_emission::TauIn, TauLine2, and TauLines.
Referenced by escmase().
Here is the call graph for this function:
| double RT_DestHummer | ( | double | beta | ) |
| beta | beta is ratio of continuum to mean line opacity, |
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
| 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 |
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
| t |
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 PressureRadiationLine().
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
| *esin | ||
| *dest | ||
| abund | ||
| ipZ | 0 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().
double chnukt_ContTkt [static] |
double chnukt_ctau [static] |
1.4.7