rt.h File Reference

Go to the source code of this file.

Data Structures

struct  t_rt

Defines

#define DEST0   SMALLFLOAT
#define LINE_CONT_SHIELD_PESC   1
#define LINE_CONT_SHIELD_FEDERMAN   2
#define LINE_CONT_SHIELD_FERLAND   3

Functions

void RT_line_all (bool lgDoEsc, bool lgUpdateFineOpac)
void RT_radiative_acceleration (void)
double RT_line_driving (void)
void RT_line_one (EmLine *t, bool lgDoEsc, bool lgUpdateFineOpac, bool lgShield_this_zone)
float RT_continuum_shield_fcn (EmLine *t)
void RT_diffuse (void)
void RT_OTS (void)
void RT_OTS_AddLine (double ots, long int ip)
void RT_OTS_Update (double *SumOTS, double BigFrac)
void RT_OTS_Zero (void)
void RT_OTS_ChkSum (long int ipPnt)
void RT_line_one_tauinc (EmLine *t, long int mas_species, long int mas_ion, long int mas_hi, long int mas_lo)
void RT_tau_init (void)
void RT_line_one_tau_reset (EmLine *t, double f1)
void RT_tau_reset (void)
void RT_tau_inc (void)
void RT_OTS_PrtRate (double weak, int chFlag)
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 ipZ)
double esc_CRDwing (double tau, double tout, double damp)
double esc_CRDcore (double tau, double tout)
double esc_PRD (double tau, double tout, double damp)
double esca0k2 (double taume)
double esccon (double tau, double hnukt)
double escConE2 (double x)
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)
double RT_recom_effic (long int ip)
void RT_stark (void)

Variables

EXTERN struct t_rt rt


Define Documentation

#define DEST0   SMALLFLOAT

DEST0 is the smallest destruction probability to return in high metallicity models

Definition at line 235 of file rt.h.

Referenced by RT_OTS(), and RTesc_lya_1side().

#define LINE_CONT_SHIELD_FEDERMAN   2

Definition at line 296 of file rt.h.

Referenced by ParseSet(), RT_continuum_shield_fcn(), and Zerologic().

#define LINE_CONT_SHIELD_FERLAND   3

Definition at line 297 of file rt.h.

Referenced by ParseSet(), and RT_continuum_shield_fcn().

#define LINE_CONT_SHIELD_PESC   1

these are all possible values of rt.nLineContShield, first is default, these are set with set continuum shielding

Definition at line 295 of file rt.h.

Referenced by ParseSet(), and RT_continuum_shield_fcn().


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 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:

float RT_continuum_shield_fcn ( EmLine t  ) 

rt_continuum_shield_fcn computing continuum shielding due to single line

Parameters:
t 

Definition at line 16 of file rt_continuum_shield_fcn.cpp.

References ASSERT, conpmp(), EmLine::damp, DEBUG_ENTRY, DEBUG_EXIT, esc_CRDwing_1side(), esc_PRD_1side(), esca0k2(), ipCRD, ipCRDW, ipLY_A, ipPRD, EmLine::iRedisFun, LINE_CONT_SHIELD_FEDERMAN, LINE_CONT_SHIELD_FERLAND, LINE_CONT_SHIELD_PESC, MIN2, t_rt::nLineContShield, POW2, rt, SDIV(), sexp(), EmLine::TauCon, and TotalInsanity().

Referenced by RT_line_static().

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:

void RT_diffuse ( void   ) 

RT_diffuse fill in DiffCont array with diffuse emission for this zone

Todo:
2 This doesn't really seem to be the expression above!!!

Todo:
2 Steve Kraemer says we're lacking radio recombination continua. Is this why? --RLP

Todo:
2 should use ConEmitLocal for all continua, but not followed by DiffuseEscape - put that at the end. Once continua all bundled this way, it will be easy to save them as a function of depth and then do exact rt

Todo:
2 Replace this constant with the appropriate macro, if any

Todo:
2 add fegrain to outward beams, but within main formalism by including grains in all x-ray processes

Definition at line 34 of file rt_diffuse.cpp.

References t_rfield::anu, t_rfield::anu2, t_iso::As2nu, ASSERT, atmdat_outer_shell(), EmLine::Aul, t_radius::BeamInIn, t_radius::BeamInOut, t_radius::BeamOutOut, C12O16Rotate, C13O16Rotate, t_iso::CaseBCheck, EmLine::ColOvTot, t_rfield::ConEmitLocal, t_rfield::ConInterOut, t_rfield::ConOTS_local_photons, t_rfield::ContBoltz, t_rfield::convoc, CoolHeavy, DEBUG_ENTRY, DEBUG_EXIT, dense, t_radius::drad, DumpLine(), t_radius::dVolOutwrd, t_dense::eden, EmisLines, EmLine::EnergyWN, t_opac::ExpmTau, FeII_RT_Out(), t_rfield::flux, EmLine::FracInwd, t_rfield::gff, t_gv::GrainEmission, GrainMakeDiffuse(), gv, H2_RT_diffuse(), Heavy, HFLines, hmi, HPLANCK, t_dense::IonHigh, t_dense::IonLow, ioQQQ, t_Heavy::ipBalHeavy, EmLine::ipCont, t_opac::ipElement, t_rfield::ipEnergyBremsThin, ipH1s, ipH_LIKE, ipHe2s3S, ipHE_LIKE, t_Heavy::ipHeavy, ipHELIUM, iphi, t_hmi::iphmin, t_opac::iphmop, t_opac::iphmra, ipHYDROGEN, t_iso::ipIsoLevNIonCon, t_Heavy::ipLyHeavy, t_rfield::ipMaxBolt, ipoint(), t_iso::ipOpac, t_rfield::ipPlasma, ipRecEsc, t_iso::ipSym2nu, t_iso::ipTwoPhoE, iso, lgAbort, t_gv::lgDustOn, t_CoolHeavy::lgFreeOn, t_gv::lgGrainPhysicsOn, t_iso::lgInd2nu_On, t_rfield::lgKillOutCont, t_rfield::lgKillOutLine, t_trace::lgTrace, LIMELM, MALLOC, MAX2, MILNE_CONST, MIN2, nCORotate, t_rfield::nflux, nHFLines, NISO, nLevel1, t_Heavy::nsShells, nu, t_iso::numLevels_local, t_rfield::nupper, nWindLine, nzone, t_rfield::OccNumbDiffCont, opac, t_opac::opacity_abs, t_opac::OpacStack, t_rfield::otscon, t_rfield::otslin, t_rfield::outlin, t_rfield::outlin_noplot, outline(), EmLine::Pesc, EmLine::phots, phycon, t_iso::Pop2Ion, EmLine::PopHi, POW2, t_radius::Radius, radius, t_iso::RadRec_caseB, t_iso::RadRecomb, t_rfield::reflin, rfield, t_radius::rinner, RYD_INF, SEXP_LIMIT, SMALLFLOAT, t_phycon::sqrte, t_iso::stat, t_iso::stat_ion, t_rfield::SummedOcc, TauLine2, TauLines, t_phycon::te, TE1RYD, t_phycon::te_ryd, t_opac::tmn, t_rfield::TotDiff2Pht, trace, TwoS, t_rfield::widflx, t_dense::xIonDense, t_iso::xIsoLevNIonRyd, and t_Heavy::xLyaHeavy.

Referenced by AGN_Hemis(), and cloudy().

Here is the call graph for this function:

void RT_line_all ( bool  lgDoEsc,
bool  lgUpdateFineOpac 
)

MakeRT drive static or wind metal line radiative transfer, if argument is true then do both escape and deEXTERN struction probs, if false only update deEXTERN struction probabilities. routine only called by cloudy(true) and ConvBase(false)

Parameters:
lgDoEsc this is true if we want to do both escape and destruction probs, and false if only destruction probabilities are needed
lgUpdateFineOpac flag saying whether to update fine opacities

Definition at line 22 of file rt_line_all.cpp.

References atmdat_2phot_rate(), atom_oi_calc(), atoms, atoms_fe2ovr(), EmLine::Aul, C12O16Rotate, C13O16Rotate, conv, DEBUG_ENTRY, DEBUG_EXIT, dense, t_DoppVel::doppler, DoppVel, t_radius::drad_x_fillfac_mean, t_hydro::dstfe2lya, DumpLine(), EmisLines, t_opac::ExpmTau, t_iso::ExtraLymanLines, FeII_RT_Make(), t_rfield::fine_opac_zone, fnzone, H2_RTMake(), HFLines, HLineTransOpacSet(), hydro, t_dense::IonHigh, EmLine::IonStg, ioQQQ, EmLine::ipCont, ipH1s, ipH2p, ipH_LIKE, ipHYDROGEN, ipLY_A, EmLine::iRedisFun, iso, t_opac::lgCaseB_no_pdest, t_rfield::lgDoLineTrans, t_rfield::lgFine_opac_update, t_conv::lgIonStageTrimed, t_trace::lgIsoTraceFull, t_hydro::lgLymanPumping, lgTauGood(), t_opac::lgTauOutOn, t_trace::lgTrace, LIMELM, MAX2, MIN2, Elevels::n, nCORotate, EmLine::nelem, t_rfield::nfine, nHFLines, NISO, nLevel1, t_iso::nLyaLevel, t_iso::nLyman, t_conv::nPres2Ioniz, t_iso::numLevels_local, t_iso::numLevels_max, nUTA, nWindLine, nzone, opac, EmLine::opacity, t_opac::opacity_abs, EmLine::Pdest, EmLine::Pesc, t_hydro::pestrk, t_atoms::pmph31, EmLine::PopLo, EmLine::PopOpc, EmLine::pump, t_iso::quant_desig, radius, rfield, RT_line_one(), RT_stark(), SDIV(), SMALLFLOAT, TauLine2, TauLines, trace, UTALines, and t_dense::xIonDense.

Referenced by cloudy(), ConvBase(), RT_tau_inc(), and RT_tau_init().

Here is the call graph for this function:

double RT_line_driving ( void   ) 

rt_line_driving radiative acceleration due to line absorption of incident continuum

Definition at line 13 of file rt_line_driving.cpp.

References C12O16Rotate, C13O16Rotate, DEBUG_ENTRY, DEBUG_EXIT, dense, EmisLines, EN1RYD, EmLine::EnergyErg, FeIIAccel(), H2_Accel(), HFLines, t_dense::IonHigh, ipH_LIKE, ipHE_LIKE, iso, LIMELM, nCORotate, nHFLines, NISO, nLevel1, t_iso::numLevels_local, nUTA, nWindLine, EmLine::PopOpc, EmLine::pump, TauLine2, TauLines, UTALines, and t_dense::xIonDense.

Referenced by radius_increment(), and RT_radiative_acceleration().

Here is the call graph for this function:

void RT_line_one ( EmLine t,
bool  lgDoEsc,
bool  lgUpdateFineOpac,
bool  lgShield_this_zone 
)

RT_line_one do line radiative transfer

Parameters:
t 
lgDoEsc the logical variable says whether to do the escape probs too (true) or just the deEXTERN struction probs (false)
lgUpdateFineOpac 
lgShield_this_zone this is option to not include line self shielding across this zone. this can cause pump to depend on zone thickness, and leads to unstable feedback in some models with the large H2 molecule, due to Solomon process depending on zone thickness and level populations.

Definition at line 34 of file rt_line_one.cpp.

References ASSERT, cdEXIT(), conv, EmLine::damp, DEBUG_ENTRY, DEBUG_EXIT, t_radius::depth_x_fillfac, t_DoppVel::doppler, DoppVel, t_rfield::fine_opac_velocity_width, t_rfield::fine_opac_zone, ioQQQ, EmLine::ipFine, ipLY_A, EmLine::iRedisFun, t_rfield::lgDoLineTrans, t_rfield::lgOpacityFine, MAX2, MIN2, EmLine::nelem, t_rfield::nfine, t_conv::nTotalIoniz, nzone, EmLine::opacity, EmLine::PopLo, EmLine::PopOpc, POW2, radius, rfield, RT_line_static(), and SMALLFLOAT.

Referenced by FeII_RT_Make(), H2_RTMake(), and RT_line_all().

Here is the call graph for this function:

void RT_line_one_tau_reset ( EmLine t,
double  f1 
)

RT_line_one_tau_reset computes average of old and new optical depths for new scale at end of iter

Parameters:
t 
f1 

Definition at line 12 of file rt_line_one_tau_reset.cpp.

References chLineLbl(), EmLine::cool, DEBUG_ENTRY, DEBUG_EXIT, t_rt::DoubleTau, EmLine::FracInwd, geometry, EmLine::heat, ioQQQ, iteration, t_prt::lgPrtMaser, t_geometry::lgSphere, t_geometry::lgStatic, MAX2, MIN2, opac, EmLine::Pdest, EmLine::Pelec_esc, EmLine::Pesc, EmLine::phots, EmLine::PopHi, EmLine::PopLo, EmLine::PopOpc, prt, rt, EmLine::TauCon, EmLine::TauIn, t_opac::taumin, EmLine::TauTot, and EmLine::xIntensity.

Referenced by FeII_RT_tau_reset(), H2_RT_tau_reset(), and RT_tau_reset().

Here is the call graph for this function:

void RT_line_one_tauinc ( EmLine t,
long int  mas_species,
long int  mas_ion,
long int  mas_hi,
long int  mas_lo 
)

RT_line_one_tauinc increment optical depths for all heavy element lines, zone by zone

Parameters:
t 
mas_species 
mas_ion 
mas_hi 
mas_lo 

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.

Definition at line 13 of file rt_line_one_tauinc.cpp.

References t_wind::AccelTot, ASSERT, DEBUG_ENTRY, DEBUG_EXIT, t_radius::depth, t_DoppVel::doppler, DoppVel, t_radius::drad_x_fillfac, t_rt::dTauMase, t_geometry::FillFac, t_rfield::fine_opac_zone, geometry, EmLine::ipFine, t_rt::lgMaserCapHit, t_rfield::lgOpacityFine, t_rt::mas_hi, t_rt::mas_ion, t_rt::mas_lo, t_rt::mas_species, MAX2, MIN2, EmLine::nelem, t_rfield::nfine, EmLine::opacity, EmLine::PopOpc, radius, rfield, rt, SDIV(), EmLine::TauCon, EmLine::TauIn, EmLine::TauTot, wind, and t_wind::windv.

Referenced by FeII_RT_TauInc(), H2_RT_tau_inc(), and RT_tau_inc().

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

void RT_OTS ( void   ) 

void RT_OTS_AddLine ( double  ots,
long int  ip 
)

RT_OTS_AddLine add local destruction of lines to ots field

Parameters:
ots 
ip pointer on the f scale

Definition at line 415 of file rt_ots.cpp.

References ASSERT, DEBUG_ENTRY, DEBUG_EXIT, ioQQQ, iteration, t_rfield::nflux, nOTS1, nOTS2, nOTS_Line_type, opac, t_opac::opacity_abs, t_rfield::otslinNew, and rfield.

Referenced by atom_level3(), CO_OTS(), CO_PopsEmisCool(), FeII_OTS(), H2_RT_OTS(), and RT_OTS().

void RT_OTS_ChkSum ( long int  ipPnt  ) 

void RT_OTS_PrtRate ( double  weak,
int  chFlag 
)

RT_OTS_PrtRate print ots arrays, called by ionize

Parameters:
weak arg is weakest rate to print
chFlag flag, 'c' continuum, 'l' line, 'b' both

Definition at line 1055 of file rt_ots.cpp.

References t_rfield::anu, ASSERT, t_rfield::chContLabel, t_rfield::chLineLabel, conv, DEBUG_ENTRY, DEBUG_EXIT, fnzone, ioQQQ, t_conv::lgConvIoniz, t_rfield::line_count, t_rfield::nupper, opac, t_opac::opacity_abs, t_rfield::otscon, t_rfield::otslin, rfield, and TorF().

Referenced by ConvBase().

Here is the call graph for this function:

void RT_OTS_Update ( double *  SumOTS,
double  BigFrac 
)

RTOTSUpdata sum flux, otscon, otslin, ConInterOut, outlin, to form SummeDif, SummedCon SummedOcc, int * is pointer to energy in continuum where this change happened,

Parameters:
BigFrac is largest relative change in ots allowed, 0 for no cap
SumOTS 
Returns:
sum of ots rates

Definition at line 545 of file rt_ots.cpp.

References t_rfield::anu, t_rfield::anu2, ASSERT, t_rfield::ConInterOut, t_rfield::ConOTS_local_OTS_rate, t_rfield::ConOTS_local_photons, t_rfield::ContBoltz, conv, t_rfield::convoc, DEBUG_ENTRY, DEBUG_EXIT, t_rfield::flux, t_rfield::flux_accum, ioQQQ, t_rfield::ipEnergyBremsThin, t_rfield::ipPlasma, t_rfield::lgKillOTSLine, t_rfield::lgOutOnly, t_conv::lgSearch, MALLOC, MAX2, MIN2, t_rfield::nflux, t_conv::nPres2Ioniz, t_conv::nTotalIoniz, t_rfield::nupper, nzone, t_rfield::OccNumbBremsCont, opac, t_opac::opacity_abs, t_rfield::otscon, t_rfield::otsconNew, t_rfield::otslin, t_rfield::otslinNew, t_rfield::outlin, t_rfield::outlin_noplot, rfield, SDIV(), SMALLDOUBLE, SMALLFLOAT, t_rfield::SummedCon, t_rfield::SummedDif, and t_rfield::SummedOcc.

Referenced by ConvBase(), and IterRestart().

Here is the call graph for this function:

void RT_OTS_Zero ( void   ) 

void RT_radiative_acceleration ( void   ) 

rt_radiative_acceleration radiative accel due to incident continuum, calls rt_line_driving for lines

Definition at line 12 of file rt_radiative_acceleration.cpp.

References t_wind::AccelCont, t_wind::AccelLine, t_wind::AccelPres, t_wind::AccelTot, t_rfield::anu, t_rfield::ConInterOut, DEBUG_ENTRY, DEBUG_EXIT, dense, EN1RYD, t_rfield::flux, t_rfield::nflux, opac, t_opac::opacity_abs, t_opac::opacity_sct, t_rfield::outlin, t_rfield::outlin_noplot, rfield, RT_line_driving(), SPEEDLIGHT, wind, and t_dense::xMassDensity.

Referenced by DynaPresChngFactor(), lgConvPres(), and radius_first().

Here is the call graph for this function:

double RT_recom_effic ( long int  ip  ) 

rt_recom_effic generate escape probability function for continua,

Parameters:
ip 

Definition at line 11 of file rt_recom_effic.cpp.

References t_rfield::anu, ASSERT, cdEXIT(), t_rfield::chDffTrns, t_rfield::ContBoltz, DEBUG_ENTRY, DEBUG_EXIT, t_opac::E2TauAbsFace, esccon(), t_opac::ExpZone, ioQQQ, t_opac::lgTauOutOn, MAX2, MIN2, t_rfield::nflux, t_rfield::nupper, opac, t_opac::otsmin, phycon, rfield, SMALLFLOAT, t_opac::TauAbsGeo, t_phycon::te, TE1RYD, and t_opac::tmn.

Referenced by HeRecom(), and HydroRecom().

Here is the call graph for this function:

void RT_stark ( void   ) 

rt_stark compute stark broadening escape probabilities using Puetter formalism

Todo:
2 - Stark is disabled for now since Lya escape causes density dependent feedback on the radiative transfer. Would need to redo the escape probs every time the electron density is updated - see blr89.in for an example

Definition at line 12 of file rt_stark.cpp.

References EmLine::Aul, DEBUG_ENTRY, DEBUG_EXIT, dense, t_dense::eden, EmisLines, hydro, ipH1s, ipH2p, ipH2s, ipH_LIKE, ipHYDROGEN, iso, t_rt::lgStarkON, MAX2, MIN2, t_iso::numLevels_local, t_iso::numLevels_max, t_hydro::pestrk, phycon, rt, SDIV(), t_phycon::sqrte, t_hydro::strkar, t_phycon::te01, t_phycon::te03, and t_phycon::te10.

Referenced by RT_line_all().

Here is the call graph for this function:

void RT_tau_inc ( void   ) 

void RT_tau_init ( void   ) 

RT_tau_init set initial outward optical depths at start of first iteration

Definition at line 30 of file rt_tau_init.cpp.

References abscf(), abund, t_rfield::anu, ASSERT, EmLine::Aul, C12O16Rotate, C13O16Rotate, cdEXIT(), t_StopCalc::colnut, conv, DEBUG_ENTRY, DEBUG_EXIT, dense, t_DoppVel::doppler, DoppVel, t_rt::DoubleTau, t_dense::eden, EmisLines, EmLineZero(), t_iso::ExtraLymanLines, FeII_LineZero(), EmLine::FracInwd, t_dense::gas_phase, geometry, H2_LineZero(), t_StopCalc::HColStop, EmLine::heat, HFLines, HLineTransOpacSet(), ioQQQ, ipH1s, ipH2p, ipH2s, ipH_LIKE, ipHe1s1S, ipHe2p1P, ipHE_LIKE, ipHYDROGEN, t_iso::ipIsoLevNIonCon, t_trace::ipIsoTrace, ipoint(), t_StopCalc::iptnu, t_rt::ipxry, iso, t_opac::lgCaseB, t_dense::lgElmtOn, t_trace::lgHBug, t_trace::lgIsoTraceFull, t_conv::lgOscilOTS, t_geometry::lgSphere, t_trace::lgTrace, t_thermal::lgTSetOn, LIMELM, MIN2, nCORotate, nHFLines, NISO, nLevel1, t_iso::nLyman, t_iso::numLevels_max, t_rfield::nupper, nUTA, nWindLine, opac, EmLine::opacity, EmLine::Pdest, EmLine::Pelec_esc, EmLine::Pesc, phycon, PrintEfmt(), t_rfield::qhtot, t_iso::quant_desig, rfield, rt, RT_line_all(), Elevels::s, ShowMe(), t_abund::solar, t_iso::stat, StopCalc, t_rt::TAddHLya, t_opac::TauAbsGeo, EmLine::TauCon, TauDummy, t_StopCalc::tauend, EmLine::TauIn, TAULIM, TauLine2, TauLines, t_opac::taumin, t_StopCalc::taunu, t_opac::TauScatGeo, EmLine::TauTot, t_opac::TauTotalGeo, t_rt::tauxry, t_phycon::te, tfidle(), thermal, t_opac::tlamin, trace, and UTALines.

Referenced by cloudy().

Here is the call graph for this function:

void RT_tau_reset ( void   ) 

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:


Variable Documentation

EXTERN struct t_rt rt


Generated for cloudy by doxygen 1.5.9