#include "module.h"

Include dependency graph for rt.h:

This graph shows which files directly or indirectly include this file:

Data Structures
struct	t_rt

Macros
#define	DEST0 SMALLFLOAT

Typedefs
typedef void(*	linefunc )(const TransitionProxy &t, bool lgShield_this_zone, realnum pestrk, realnum DopplerWidth)

Enumerations
enum	{ LINE_CONT_SHIELD_PESC = 1, LINE_CONT_SHIELD_FEDERMAN, LINE_CONT_SHIELD_FEDERMAN_BUG, LINE_CONT_SHIELD_FERLAND, LINE_CONT_SHIELD_RODGERS, LINE_CONT_SHIELD_INTEGRAL }

Functions
void	RT_line_one_escape (const TransitionProxy &t, bool lgShield_this_zone, realnum pestrk, realnum DopplerWidth)

void	RT_line_one_fine (const TransitionProxy &t, bool lgShield_this_zone, realnum pestrk, realnum DopplerWidth)

void	RT_fine_clear ()

void	RT_line_all (linefunc line_one)

void	RT_line_all_escape (realnum *error)

double	RT_line_driving (void)

double	RT_continuum_shield_fcn (const TransitionProxy &t, bool lgShieldThisZone, double dTau)

void	RT_diffuse (void)

void	RT_continuum (void)

void	RT_OTS (void)

void	RT_OTS_AddLine (double ots, long int ip)

void	RT_OTS_Update (double *SumOTS)

void	RT_OTS_Zero (void)

void	RT_OTS_ChkSum (long int ipPnt)

void	RT_line_one_tauinc (const TransitionProxy &t, long int mas_species, long int mas_ion, long int mas_hi, long int mas_lo, realnum DopplerWidth)

void	RT_tau_init (void)

void	RT_line_one_tau_reset (const TransitionProxy &t)

void	RT_tau_reset (void)

void	RT_tau_inc (void)

void	RT_OTS_PrtRate (double weak, int chFlag)

double	RT_recom_effic (long int ip)

void	RT_stark (void)

double	RT_EscLVG (double tau, double sigma)

double	DrvContPump (double tau, double damp)

void	DrivePump (double tau)

void	prt_trans_opc_debug (const char *LineGroup, const TransitionProxy &t)

Variables
t_rt	rt

Macro Definition Documentation

#define DEST0 SMALLFLOAT

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

Definition at line 154 of file rt.h.

Referenced by RT_OTS(), and RTesc_lya_1side().

Typedef Documentation

typedef void(* linefunc)(const TransitionProxy &t, bool lgShield_this_zone, realnum pestrk, realnum DopplerWidth)

Definition at line 35 of file rt.h.

Enumeration Type Documentation

anonymous enum

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

Enumerator
LINE_CONT_SHIELD_PESC
LINE_CONT_SHIELD_FEDERMAN
LINE_CONT_SHIELD_FEDERMAN_BUG
LINE_CONT_SHIELD_FERLAND
LINE_CONT_SHIELD_RODGERS
LINE_CONT_SHIELD_INTEGRAL

Definition at line 221 of file rt.h.

Function Documentation

void DrivePump ( double tau )

TestPump – test continuum pumping routines

Definition at line 455 of file rt_continuum_shield_fcn.cpp.

References conpmp(), conpmp_qg32(), conpmp_romb(), exp10(), fprintf(), growth_romb(), growthRodgers(), ioQQQ, open_data(), safe_div(), shieldFederman(), shieldRodgers(), and integrate::Romberg< T >::update().

Referenced by ParseDrive().

Here is the call graph for this function:

double DrvContPump	(	double	tau,
		double	damp
	)

DrvContPump local continuum pumping rate radiative transfer for all lines

Parameters

*t
DopplerWidth

Definition at line 449 of file rt_continuum_shield_fcn.cpp.

References conpmp(), and DEBUG_ENTRY.

Referenced by ParseDrive().

Here is the call graph for this function:

void prt_trans_opc_debug	(	const char *	LineGroup,
		const TransitionProxy &	t
	)

Definition at line 24 of file rt_tau_inc.cpp.

References t_elementnames::chElementSym, TransitionProxy::chLabel(), dense, elementnames, TransitionProxy::Emis(), fprintf(), TransitionProxy::Hi(), ioQQQ, EmissionProxy::TauCon(), EmissionProxy::TauIn(), EmissionProxy::TauTot(), and t_dense::xIonDense.

Here is the call graph for this function:

void RT_continuum ( void )

RT_continuum attenuation of diffuse and beamed continua

Definition at line 70 of file rt_continuum.cpp.

References t_mesh::anu(), ASSERT, cdEXIT, t_rfield::chContLabel, t_rfield::chLineLabel, cmshft(), column(), t_rfield::ConEmitLocal, t_rfield::ConEmitOut, t_rfield::ConEmitReflec, t_rfield::ConInterOut, t_rfield::ConRefIncid, t_rfield::convoc, cosmology, t_hextra::CrsSecNeutron, DEBUG_ENTRY, dense, t_radius::depth, t_rfield::DiffuseEscape, t_geometry::DirectionalCosin, t_radius::drad, t_radius::drad_x_fillfac, t_radius::dRadSign, t_radius::dVolOutwrd, e2(), t_opac::E2TauAbsFace, t_opac::E2TauAbsOut, t_rfield::EnergyDiffCont, t_rfield::EnergyIncidCont, EXIT_FAILURE, t_opac::ExpmTau, t_opac::ExpZone, t_iso_sp::fb, findspecieslocal(), fixit, t_rfield::flux, t_rfield::flux_beam_const, t_rfield::flux_beam_time, t_rfield::flux_isotropic, t_rfield::flux_total_incident, fprintf(), t_dense::gas_phase, geometry, get_ptr(), hextra, ioQQQ, ipH_LIKE, ipHYDROGEN, t_rfield::ipPlasma, iso_sp, iteration, t_trace::lgConBug, t_cosmology::lgDo, t_hextra::lgNeutrnHeatOn, t_opac::lgOpacNeg, t_opac::lgScatON, t_geometry::lgSphere, t_trace::lgTrace, MAX2, t_rfield::nflux, nzone, t_rfield::OccNumbContEmitOut, t_rfield::OccNumbDiffCont, t_rfield::OccNumbIncidCont, opac, t_opac::opacity_abs, t_opac::opacity_sct, t_rfield::otscon, t_rfield::otslin, t_rfield::outlin, t_rfield::outlin_noplot, pnegopc(), POW2, pow2(), t_radius::r1r0sq, radius, t_radius::Radius, t_rfield::reflin, rfield, t_radius::rinner, SDIV(), t_rfield::setTrimming(), sexp(), ShowMe(), SMALLFLOAT, t_rfield::SummedCon, t_rfield::SummedDif, t_rfield::SummedOcc, t_opac::TauAbsFace, t_opac::TauAbsGeo, t_opac::TauAbsTotal, t_opac::TauScatFace, t_opac::TauScatGeo, t_opac::TauTotalGeo, t_opac::tmn, TotalInsanity(), t_hextra::totneu, trace, vexp(), t_rfield::vexp_arg, and t_dense::xIonDense.

Referenced by cloudy().

Here is the call graph for this function:

double RT_continuum_shield_fcn	(	const TransitionProxy &	t,
		bool	lgShieldThisZone,
		double	dTau
	)

rt_continuum_shield_fcn computing continuum shielding due to single line

Parameters

t

Definition at line 328 of file rt_continuum_shield_fcn.cpp.

References avg_shield(), chLineLbl(), DEBUG_ENTRY, TransitionProxy::Emis(), fprintf(), ioQQQ, TransitionProxy::ipCont(), nzone, RT_continuum_shield_fcn_point(), and EmissionProxy::TauCon().

Referenced by RT_line_pumping().

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 should use ConEmitLocal for all continua, but not followed by rfield.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 add fegrain to outward beams, but within main formalism by including grains in all x-ray processes

Definition at line 35 of file rt_diffuse.cpp.

References t_mesh::anu(), t_mesh::anu2(), ASSERT, atmdat_outer_shell(), EmissionProxy::Aul(), t_radius::BeamInOut, t_radius::BeamOutOut, t_gaunt::brems_rt(), t_gaunt::brems_sum_ions(), CalcTwoPhotonEmission(), t_rfield::ConEmitLocal, t_rfield::ConInterOut, t_rfield::ConSourceFcnLocal, t_rfield::ContBoltz, t_rfield::convoc, CoolHeavy, dBaseSpecies, dBaseTrans, DEBUG_ENTRY, t_brems_den::den_Hep, t_brems_den::den_Hepp, t_brems_den::den_Hm, t_brems_den::den_Hp, t_brems_den::den_ion, dense, diatoms, t_rfield::DiffuseEscape, t_rfield::DiffuseLineEmission, DumpLine(), t_radius::dVolOutwrd, t_dense::eden, t_rfield::eeBremsDif, emergent_line(), TransitionProxy::Emis(), TransitionList::Emis(), TransitionProxy::EnergyErg(), t_opac::ExpmTau, t_rfield::flux, t_rfield::flux_beam_const, t_rfield::flux_beam_time, t_rfield::flux_isotropic, fprintf(), GrainVar::GrainEmission, GrainMakeDiffuse(), gv, Heavy, HFLines, hmi, Singleton< t_gaunt >::Inst(), t_dense::IonHigh, t_dense::IonLow, ioQQQ, t_Heavy::ipBalHeavy, TransitionProxy::ipCont(), t_opac::ipElement, ipH1s, ipH_LIKE, ipHE_LIKE, t_Heavy::ipHeavy, t_hmi::iphmin, t_opac::iphmop, ipHYDROGEN, t_Heavy::ipLyHeavy, t_rfield::ipMaxBolt, t_rfield::ipPlasma, ipSatelliteLines, iso_ctrl, iso_sp, lgAbort, t_isoCTRL::lgDielRecom, GrainVar::lgDustOn(), t_dense::lgElmtOn, t_CoolHeavy::lgFreeOn, GrainVar::lgGrainPhysicsOn, t_isoCTRL::lgInd2nu_On, t_rfield::lgInducProcess, t_rfield::lgKillOutCont, t_rfield::lgKillOutLine, t_trace::lgTrace, LIMELM, MAX2, MIN2, t_rfield::nflux, t_rfield::nflux_with_check, NISO, nSpecies, t_Heavy::nsShells, t_iso_sp::numLevels_local, 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, TransitionProxy::outline(), EmissionProxy::Pesc(), phots(), phycon, pow2(), PrtTwoPhotonEmissCoef(), radius, t_Heavy::RadRecCon, t_rfield::reflin, rfield, RT_iso_integrate_RRC(), safe_div(), SatelliteLines, SEXP_LIMIT, TransitionList::size(), SMALLFLOAT, t_phycon::sqrte, t_iso_sp::st, TauLine2, t_phycon::te, t_opac::tmn, t_rfield::TotDiff2Pht, trace, t_iso_sp::trans(), t_iso_sp::TwoNu, t_mesh::widflx(), EmissionProxy::xIntensity(), t_dense::xIonDense, and t_Heavy::xLyaHeavy.

Referenced by AGN_Hemis(), and cloudy().

Here is the call graph for this function:

double RT_EscLVG	(	double	tau,
		double	sigma
	)

Definition at line 1115 of file rt_escprob.cpp.

Referenced by RT_line_escape().

void RT_fine_clear ( )

Definition at line 402 of file rt_line_all.cpp.

References cosmology, DEBUG_ENTRY, t_rfield::fine_opac_velocity_width, t_rfield::fine_opac_zone, t_rfield::ipFineConVelShift, t_cosmology::lgDo, t_rfield::nfine, POW2, t_cosmology::redshift_current, t_cosmology::redshift_start, rfield, wind, Wind::windv, and Wind::windv0.

Referenced by RT_tau_inc().

void RT_line_all ( linefunc line_one )

MakeRT drive static or wind metal line radiative transfer,

Definition at line 200 of file rt_line_all.cpp.

References atmdat, t_dense::AtomicWeight, conv, dBaseSpecies, dBaseTrans, DEBUG_ENTRY, dense, diatoms, TransitionList::Emis(), ExtraLymanLines, fixit, fnzone, fprintf(), GetDopplerWidth(), HFLines, t_dense::IonHigh, ioQQQ, ipExtraLymanLines, ipH_LIKE, ipHE_LIKE, ipHYDROGEN, ipSatelliteLines, iso_ctrl, iso_sp, t_atmdat::lgChiantiHybrid, t_isoCTRL::lgDielRecom, t_rfield::lgDoLineTrans, t_dense::lgElmtOn, t_conv::lgFirstSweepThisZone, t_dense::lgIonChiantiOn, t_dense::lgIonStoutOn, t_conv::lgLastSweepThisZone, t_atmdat::lgStoutHybrid, t_trace::lgTrace, LIMELM, t_dense::maxWN, NISO, t_isoCTRL::nLyman, t_conv::nPres2Ioniz, nSpecies, t_iso_sp::numLevels_local, nWindLine, rfield, SatelliteLines, TransitionList::size(), t_iso_sp::st, TauLine2, trace, UTALines, and t_dense::xIonDense.

Referenced by RT_line_all_escape(), and RT_tau_inc().

Here is the call graph for this function:

void RT_line_all_escape ( realnum * error )

Definition at line 21 of file rt_line_all.cpp.

References AllTransitions, t_dense::AtomicWeight, conv, DEBUG_ENTRY, dense, t_radius::drad_x_fillfac, t_hydro::dstfe2lya, DumpLine(), TransitionProxy::Emis(), fnzone, fprintf(), frac(), GetDopplerWidth(), hydro, t_dense::IonHigh, ioQQQ, TransitionProxy::ipCont(), ipH1s, ipH2p, ipH_LIKE, ipHYDROGEN, iso_ctrl, iso_sp, t_opac::lgCaseB_no_pdest, t_dense::lgElmtOn, t_conv::lgFirstSweepThisZone, t_conv::lgIonStageTrimed, t_conv::lgLastSweepThisZone, t_hydro::lgLymanPumping, lgTauGood(), LIMELM, max(), NISO, t_isoCTRL::nLyaLevel, t_iso_sp::numLevels_max, nzone, opac, EmissionProxy::opacity(), t_opac::opacity_abs, EmissionProxy::Pdest(), EmissionProxy::PopOpc(), EmissionProxy::pump(), radius, RT_line_all(), RT_line_one_escape(), RT_stark(), SDIV(), SMALLFLOAT, and t_iso_sp::trans().

Referenced by ConvBase(), EdenChange(), EdenError(), 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 returns line radiative acceleration cm s-2

Definition at line 15 of file rt_line_driving.cpp.

References atmdat, dBaseSpecies, dBaseTrans, DEBUG_ENTRY, dense, diatoms, TransitionProxy::Emis(), TransitionList::Emis(), TransitionProxy::EnergyErg(), ExtraLymanLines, fprintf(), HFLines, t_dense::IonHigh, ioQQQ, ipCARBON, ipExtraLymanLines, ipH_LIKE, ipHE_LIKE, ipSatelliteLines, t_atmdat::ipSpecIon, iso_ctrl, iso_sp, t_isoCTRL::lgDielRecom, LIMELM, NISO, t_isoCTRL::nLyman, nSpecies, t_iso_sp::numLevels_local, nWindLine, EmissionProxy::PopOpc(), EmissionProxy::pump(), SatelliteLines, TransitionList::size(), TauLine2, t_iso_sp::trans(), and UTALines.

Referenced by PresTotCurrent().

Here is the call graph for this function:

void RT_line_one_escape	(	const TransitionProxy &	t,
		bool	lgShield_this_zone,
		realnum	pestrk,
		realnum	DopplerWidth
	)

RT_line_one do line radiative transfer

Parameters

t	line structure
pestrk	Stark escape probability
DopplerWidth
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 383 of file rt_line_one.cpp.

References ASSERT, EmissionProxy::damp(), EmissionProxy::dampXvel(), DEBUG_ENTRY, TransitionProxy::Emis(), TransitionProxy::EnergyErg(), EmissionProxy::FracInwd(), ipLY_A, EmissionProxy::iRedisFun(), t_rfield::lgDoLineTrans, lgTauGood(), TransitionProxy::Lo(), EmissionProxy::Pdest(), EmissionProxy::Pelec_esc(), EmissionProxy::Pesc(), t_rfield::plsfrq, EmissionProxy::pump(), rfield, RT_line_escape(), RT_line_pumping(), and SMALLFLOAT.

Referenced by CoolDima(), CoolHyperfine(), and RT_line_all_escape().

Here is the call graph for this function:

void RT_line_one_fine	(	const TransitionProxy &	t,
		bool	lgShield_this_zone,
		realnum	pestrk,
		realnum	DopplerWidth
	)

Definition at line 449 of file rt_line_one.cpp.

References ASSERT, cdEXIT, EmissionProxy::damp(), EmissionProxy::dampXvel(), DEBUG_ENTRY, TransitionProxy::Emis(), TransitionProxy::EnergyErg(), EXIT_SUCCESS, fprintf(), ioQQQ, ipLY_A, EmissionProxy::iRedisFun(), t_rfield::lgDoLineTrans, TransitionProxy::Lo(), MAX2, MIN2, nzone, t_rfield::plsfrq, rfield, RT_line_fine_opacity(), and SMALLFLOAT.

Referenced by RT_tau_inc().

Here is the call graph for this function:

void RT_line_one_tau_reset ( const TransitionProxy & t )

Definition at line 13 of file rt_line_one_tau_reset.cpp.

References chLineLbl(), TransitionProxy::Coll(), CollisionProxy::cool(), cosmology, DEBUG_ENTRY, t_rt::DoubleTau, TransitionProxy::Emis(), fprintf(), EmissionProxy::FracInwd(), geometry, CollisionProxy::heat(), TransitionProxy::Hi(), ioQQQ, TransitionProxy::ipCont(), iteration, t_cosmology::lgDo, t_prt::lgPrtMaser, t_geometry::lgSphere, t_geometry::lgStatic, TransitionProxy::Lo(), MAX2, MIN2, EmissionProxy::mult_opac(), secant_track::next_val(), opac, EmissionProxy::Pdest(), EmissionProxy::Pelec_esc(), EmissionProxy::Pesc(), EmissionProxy::PopOpc(), secant_track::PREV_ITER, prt, rt, EmissionProxy::TauCon(), EmissionProxy::TauIn(), EmissionProxy::TauInSpecific(), t_opac::taumin, EmissionProxy::TauTot(), EmissionProxy::TauTrack(), EmissionProxy::xIntensity(), and EmissionProxy::xObsIntensity().

Referenced by diatomics::H2_RT_tau_reset(), and RT_tau_reset().

Here is the call graph for this function:

void RT_line_one_tauinc	(	const TransitionProxy &	t,
		long int	mas_species,
		long int	mas_ion,
		long int	mas_hi,
		long int	mas_lo,
		realnum	DopplerWidth
	)

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
DopplerWidth

Definition at line 14 of file rt_line_one_tauinc.cpp.

References Wind::AccelGravity, Wind::AccelTotalOutward, t_mesh::anu(), TransitionProxy::chLabel(), cosmology, DEBUG_ENTRY, t_radius::depth, t_radius::drad_x_fillfac, t_rt::dTauMase, Wind::dvdr, TransitionProxy::Emis(), t_geometry::FillFac, t_rfield::fine_opac_zone, fprintf(), geometry, GetHubbleFactor(), ioQQQ, TransitionProxy::ipCont(), EmissionProxy::ipFine(), t_rfield::ipFineConVelShift, Wind::lgBallistic(), t_cosmology::lgDo, t_rt::lgMaserCapHit, t_rfield::lgOpacityFine, t_rt::mas_hi, t_rt::mas_ion, t_rt::mas_lo, t_rt::mas_species, MIN2, EmissionProxy::mult_opac(), t_rfield::nfine, nzone, EmissionProxy::opacity(), t_rfield::plsfrq, EmissionProxy::PopOpc(), radius, t_cosmology::redshift_current, rfield, rt, SDIV(), EmissionProxy::TauCon(), EmissionProxy::TauIn(), EmissionProxy::TauInSpecific(), EmissionProxy::TauTot(), EmissionProxy::VoigtLineCen(), wind, and Wind::windv.

Referenced by diatomics::H2_RT_tau_inc(), and RT_tau_inc().

Here is the call graph for this function:

void RT_OTS ( void )

RT_OTS compute diffuse fields due to helium atom, metals, hydro done in HydroOTS

Todo:: 2 should also check whether IonLo is in bounds - in func_set_ion test He0 is set to zero, so this does not do anything. as NISO grows this will become larger waste of time

Definition at line 41 of file rt_ots.cpp.

References ASSERT, EmissionProxy::Aul(), dBaseSpecies, dBaseTrans, DEBUG_ENTRY, dense, DEST0, diatoms, t_dense::eden, t_dense::EdenTrue, TransitionProxy::Emis(), TransitionList::Emis(), t_opac::ExpmTau, t_iso_sp::fb, findspecieslocal(), fnzone, fprintf(), he, Heavy, hmi, t_hmi::HMinus_photo_rate, t_dense::IonHigh, ioQQQ, t_he::ip660, t_Heavy::ipBalHeavy, TransitionProxy::ipCont(), ipH1s, ipH2p, ipH_LIKE, ipHELIUM, ipHYDROGEN, t_Heavy::ipLyHeavy, ipRecEsc, ipRecRad, iso_sp, t_dense::lgElmtOn, t_rfield::lgHeIIOTS, t_rfield::lgLyaOTS, t_trace::lgTrace, LIMELM, MAX2, t_rfield::nflux, NISO, nOTS1, nOTS2, nOTS_Line_type, nSpecies, t_iso_sp::numLevels_local, nWindLine, nzone, opac, t_opac::opacity_abs, EmissionProxy::ots(), t_rfield::otscon, t_rfield::otslin, EmissionProxy::Pdest(), rfield, RT_OTS_AddCont(), RT_OTS_AddLine(), SMALLFLOAT, t_iso_sp::st, TauLine2, trace, t_iso_sp::trans(), t_dense::xIonDense, and t_Heavy::xLyaHeavy.

Referenced by ConvBase().

Here is the call graph for this function:

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 390 of file rt_ots.cpp.

References ASSERT, DEBUG_ENTRY, fprintf(), ioQQQ, iteration, t_rfield::nflux, nOTS1, nOTS2, nOTS_Line_type, opac, t_opac::opacity_abs, t_rfield::otslin, and rfield.

Referenced by diatomics::H2_RT_OTS(), and RT_OTS().

Here is the call graph for this function:

void RT_OTS_ChkSum ( long int ipPnt )

RT_OTS_ChkSum sanity check confirms summed continua reflect contents of individuals

Definition at line 606 of file rt_ots.cpp.

References t_mesh::anu(), t_rfield::chContLabel, t_rfield::chLineLabel, t_rfield::ConInterOut, t_rfield::ConOTS_local_OTS_rate, DEBUG_ENTRY, t_rfield::flux, fprintf(), ioQQQ, t_rfield::ipEnergyBremsThin, t_rfield::lgOutOnly, t_rfield::nflux, t_rfield::otscon, t_rfield::otslin, t_rfield::outlin, t_rfield::outlin_noplot, rfield, t_rfield::SummedCon, t_rfield::SummedDif, and TotalInsanity().

Referenced by ConvPresTempEdenIoniz().

Here is the call graph for this function:

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 687 of file rt_ots.cpp.

References t_mesh::anu(), ASSERT, t_rfield::chContLabel, t_rfield::chLineLabel, conv, DEBUG_ENTRY, dense, fnzone, fprintf(), t_dense::gas_phase, ioQQQ, ipHELIUM, ipOXYGEN, t_conv::lgConvIoniz(), t_rfield::line_count, t_rfield::nflux_with_check, opac, t_opac::opacity_abs, t_rfield::otscon, t_rfield::otslin, rfield, TorF(), and t_dense::xIonDense.

Referenced by ConvBase().

Here is the call graph for this function:

void RT_OTS_Update ( double * SumOTS )

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

Parameters

SumOTS

Returns: sum of ots rates

Definition at line 476 of file rt_ots.cpp.

References ASSERT, CalcTwoPhotonEmission(), t_rfield::ConInterOut, t_rfield::ConOTS_local_OTS_rate, t_rfield::ConOTS_local_photons, t_rfield::convoc, cosmology, DEBUG_ENTRY, dense, t_opac::ExpmTau, t_rfield::flux, t_rfield::flux_accum, t_dense::IonHigh, ipH_LIKE, t_rfield::ipPlasma, iso_ctrl, iso_sp, t_cosmology::lgDo, t_isoCTRL::lgInd2nu_On, t_rfield::lgInducProcess, t_rfield::lgKillOTSLine, t_rfield::lgOutOnly, LIMELM, MAX2, t_rfield::nflux, t_rfield::nflux_with_check, NISO, opac, t_opac::opacity_abs, t_rfield::otscon, t_rfield::otslin, t_rfield::outlin, t_rfield::outlin_noplot, rfield, t_rfield::setTrimming(), SMALLDOUBLE, t_rfield::SummedCon, t_rfield::SummedDif, t_rfield::SummedOcc, and t_iso_sp::TwoNu.

Referenced by ConvBase(), EdenError(), and IterRestart().

Here is the call graph for this function:

void RT_OTS_Zero ( void )

zero these things out, called in zero

Definition at line 571 of file rt_ots.cpp.

References t_rfield::ConInterOut, t_rfield::ConOTS_local_OTS_rate, t_rfield::ConOTS_local_photons, t_rfield::convoc, DEBUG_ENTRY, t_rfield::flux, t_rfield::nflux, t_rfield::otscon, t_rfield::otslin, t_rfield::outlin, t_rfield::outlin_noplot, t_rfield::resetCoarseTransCoef(), rfield, t_rfield::setTrimming(), t_rfield::SummedCon, t_rfield::SummedDif, and t_rfield::SummedOcc.

Referenced by ContSetIntensity().

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 12 of file rt_recom_effic.cpp.

References t_mesh::anu(), ASSERT, cdEXIT, t_rfield::chDffTrns, t_rfield::ContBoltz, cosmology, DEBUG_ENTRY, dsexp(), t_opac::E2TauAbsFace, esccon(), EXIT_FAILURE, t_opac::ExpZone, t_rfield::fine_opac_zone, fprintf(), GetHubbleFactor(), ioQQQ, t_rfield::ipnt_coarse_2_fine, iteration, t_cosmology::lgDo, t_rfield::lgOpacityFine, MAX2, MIN2, t_rfield::nfine, t_rfield::nflux, t_rfield::nflux_with_check, opac, t_opac::opacity_abs, t_opac::otsmin, phycon, t_cosmology::redshift_current, rfield, SMALLFLOAT, t_opac::TauAbsGeo, t_opac::taumin, t_phycon::te, and t_opac::tmn.

Referenced by iso_radiative_recomb().

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 EmissionProxy::Aul(), DEBUG_ENTRY, dense, t_dense::eden, TransitionProxy::Emis(), t_iso_sp::ex, TransitionProxy::ipCont(), ipH2s, ipH_LIKE, iso_ctrl, iso_sp, t_dense::lgElmtOn, t_rt::lgStarkON, LIMELM, max(), MAX2, MIN2, NISO, t_isoCTRL::nLyaLevel, t_iso_sp::numLevels_local, t_iso_sp::numLevels_max, nzone, phycon, powpq(), rt, t_phycon::sqrte, t_iso_sp::st, strkar(), EmissionProxy::TauIn(), t_phycon::te01, t_phycon::te03, t_phycon::te10, and t_iso_sp::trans().