/home66/gary/public_html/cloudy/c13_branch/source/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 (void)
double	RT_line_driving (void)
void	RT_line_one (const TransitionProxy &t, bool lgShield_this_zone, realnum pestrk, realnum DopplerWidth)
double	RT_continuum_shield_fcn (const TransitionProxy &t)
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	esc_PRD_1side (double tau, double a)
double	esc_CRDwing_1side (double tau, double a)
double	RTesc_lya (double esin, double dest, double abund, const TransitionProxy &t, realnum DopplerWidth)
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	RT_DestProb (double abund, double crsec, long int ipanu, double widl, double escp, int nCore)
double	RT_LineWidth (const TransitionProxy &t, realnum DopplerWidth)
double	RT_DestHummer (double beta)
double	RT_recom_effic (long int ip)
void	RT_stark (void)
void	RT_iso_integrate_RRC (const long ipISO, const long nelem, const bool lgUpdateContinuum)
Variables
t_rt	rt

Define Documentation

#define DEST0 SMALLFLOAT

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

Definition at line 227 of file rt.h.

Referenced by RT_OTS(), and RTesc_lya_1side().

#define LINE_CONT_SHIELD_FEDERMAN 2

Definition at line 291 of file rt.h.

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

#define LINE_CONT_SHIELD_FERLAND 3

Definition at line 292 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 290 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 445 of file rt_escprob.cpp.

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

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

References DEBUG_ENTRY, esca0k2(), ioQQQ, POW2, SQRTPI, Integrator< Integrand, Method >::sum(), and trapezium().

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

References ASSERT, DEBUG_ENTRY, esc_PRD_1side(), t_rt::fracin, iteration, rt, tau_from_here(), 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 97 of file rt_escprob.cpp.

References ASSERT, DEBUG_ENTRY, escmase(), MIN2, pow(), and rt.

Referenced by DrvEscP(), esc_PRD(), mole_h_reactions(), 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 490 of file rt_escprob.cpp.

References DEBUG_ENTRY, escmase(), and SQRTPI.

Referenced by DrvEscP(), esc_CRDcore(), esc_CRDwing_1side(), and RT_continuum_shield_fcn().

Here is the call graph for this function:

double esccon	(	double	tau,
		double	hnukt
	)

escpcn continuum escape probability

Parameters:

	tau
	hnukt

Definition at line 628 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:

void RT_continuum ( void )

RT_continuum attenuation of diffuse and beamed continua

Definition at line 67 of file rt_continuum.cpp.

References t_rfield::anu, ASSERT, cdEXIT, t_rfield::chContLabel, t_rfield::chLineLabel, cmshft(), t_colden::colden, colden, t_rfield::ConEmitLocal, t_rfield::ConEmitOut, t_rfield::ConEmitReflec, t_rfield::ConInterOut, t_rfield::ConRefIncid, t_rfield::convoc, 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, EN1RYD, t_rfield::EnergyDiffCont, t_rfield::EnergyIncidCont, EXIT_FAILURE, t_opac::ExpmTau, t_opac::ExpZone, t_iso_sp::fb, t_rfield::fine_opac_zone, t_rfield::fine_opt_depth, t_rfield::flux, t_rfield::flux_beam_const, t_rfield::flux_beam_time, t_rfield::flux_isotropic, t_rfield::flux_total_incident, t_dense::gas_phase, geometry, hextra, ioQQQ, ipCOL_H0, ipH_LIKE, ipHYDROGEN, t_rfield::ipPlasma, iso_sp, iteration, t_trace::lgConBug, t_hextra::lgNeutrnHeatOn, t_rfield::lgOpacityFine, t_opac::lgOpacNeg, t_opac::lgScatON, t_geometry::lgSphere, t_trace::lgTrace, MAX2, t_rfield::nfine, 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, t_radius::r1r0sq, t_radius::Radius, radius, t_rfield::reflin, rfield, t_radius::rinner, SDIV(), sexp(), ShowMe(), SMALLFLOAT, t_rfield::SummedCon, t_rfield::SummedDif, t_opac::TauAbsFace, t_opac::TauAbsGeo, t_opac::TauAbsTotal, t_opac::TauScatFace, t_opac::TauScatGeo, t_opac::TauTotalGeo, TestCode(), t_opac::tmn, TotalInsanity(), t_hextra::totneu, trace, t_rfield::trans_coef_total_stale, and t_dense::xIonDense.

Referenced by cloudy().

Here is the call graph for this function:

double RT_continuum_shield_fcn ( const TransitionProxy & t )

rt_continuum_shield_fcn computing continuum shielding due to single line

Parameters:

t

Definition at line 33 of file rt_continuum_shield_fcn.cpp.

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

Referenced by RT_line_pumping().

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

References t_opac::albedo, t_rfield::anu, ASSERT, cdEXIT, conv, DEBUG_ENTRY, EXIT_FAILURE, fp_equal(), ioQQQ, ipDEST_INCOM, ipDEST_K2, ipDEST_SIMPL, ipH1s, ipH2p, ipH_LIKE, ipHYDROGEN, iso_sp, MIN2, t_rfield::nflux, t_conv::nTotalIoniz, nzone, opac, t_opac::opacity_abs, POW2, rfield, and SQRTPI.

Referenced by RT_line_escape().

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 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, ASSERT, atmdat_outer_shell(), EmissionProxy::Aul(), t_radius::BeamInOut, t_radius::BeamOutOut, CalcTwoPhotonEmission(), t_rfield::ConEmitLocal, t_rfield::ConInterOut, t_rfield::ConOTS_local_photons, t_rfield::ConSourceFcnLocal, t_rfield::ContBoltz, t_rfield::convoc, CoolHeavy, dBaseSpecies, dBaseTrans, DEBUG_ENTRY, dense, diatoms, t_rfield::DiffuseEscape, t_rfield::DiffuseLineEmission, DumpLine(), t_radius::dVolOutwrd, t_dense::eden, emergent_line(), TransitionList::Emis(), TransitionProxy::Emis(), EN1RYD, t_opac::ExpmTau, FeII_RT_Out(), t_rfield::flux, t_rfield::flux_beam_const, t_rfield::flux_beam_time, t_rfield::flux_isotropic, t_rfield::gff, GrainVar::GrainEmission, GrainMakeDiffuse(), gv, Heavy, HFLines, hmi, t_dense::IonHigh, t_dense::IonLow, ioQQQ, t_Heavy::ipBalHeavy, TransitionProxy::ipCont(), t_opac::ipElement, t_rfield::ipEnergyBremsThin, ipH1s, ipH_LIKE, ipHE_LIKE, t_Heavy::ipHeavy, t_hmi::iphmin, t_opac::iphmop, t_opac::iphmra, 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, t_mole_global::list, MAX2, MILNE_CONST, MIN2, mole, mole_global, t_rfield::nflux, nHFLines, NISO, nLevel1, nSpecies, t_Heavy::nsShells, t_mole_global::num_calc, t_iso_sp::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, TransitionProxy::outline(), EmissionProxy::Pesc(), EmissionProxy::phots(), phycon, POW2, PrtTwoPhotonEmissCoef(), radius, t_Heavy::RadRecCon, t_rfield::reflin, rfield, RT_iso_integrate_RRC(), safe_div(), SatelliteLines, SEXP_LIMIT, SMALLFLOAT, t_mole_local::species, t_phycon::sqrte, t_iso_sp::st, TauLine2, TauLines, t_phycon::te, TE1RYD, t_opac::tmn, t_rfield::TotDiff2Pht, trace, t_iso_sp::trans(), t_iso_sp::TwoNu, t_rfield::widflx, t_dense::xIonDense, and t_Heavy::xLyaHeavy.

Referenced by AGN_Hemis(), and cloudy().

Here is the call graph for this function:

void RT_iso_integrate_RRC	(	const long	ipISO,
		const long	nelem,
		const bool	lgUpdateContinuum
	)

Definition at line 556 of file rt_diffuse.cpp.

References t_rfield::anu, t_rfield::anu2, ASSERT, t_iso_sp::CaseBCheck, t_rfield::ConEmitLocal, conv, DEBUG_ENTRY, dense, t_rfield::DiffuseEscape, t_dense::eden, emergent_line(), EN1RYD, t_iso_sp::fb, fixit(), fp_equal(), t_dense::IonHigh, ipRecEsc, ipRecNetEsc, iso_ctrl, iso_sp, LIMELM, MAX2, MILNE_CONST, t_rfield::nflux, NISO, t_conv::nTotalIoniz, t_iso_sp::numLevels_local, nzone, opac, t_opac::OpacStack, phycon, t_iso_sp::RadRec_caseB, rfield, SEXP_LIMIT, t_phycon::sqrte, t_iso_sp::st, t_isoCTRL::stat_ion, t_phycon::te, t_phycon::te_ryd, t_rfield::widflx, and t_dense::xIonDense.

Referenced by iso_cool(), and RT_diffuse().

Here is the call graph for this function:

void RT_line_all ( void )

MakeRT drive static or wind metal line radiative transfer,

Definition at line 26 of file rt_line_all.cpp.

References t_dense::AtomicWeight, t_fe2ovr_la::atoms_fe2ovr(), conv, cosmology, dBaseSpecies, dBaseTrans, DEBUG_ENTRY, dense, diatoms, t_radius::drad_x_fillfac, t_hydro::dstfe2lya, DumpLine(), TransitionList::Emis(), TransitionProxy::Emis(), ex, ExtraLymanLines, t_rfield::fine_opac_velocity_width, t_rfield::fine_opac_zone, fixit(), fnzone, GetDopplerWidth(), hydro, Singleton< t_fe2ovr_la >::Inst(), t_dense::IonHigh, ioQQQ, TransitionProxy::ipCont(), ipExtraLymanLines, t_rfield::ipFineConVelShift, ipH1s, ipH2p, ipH_LIKE, ipHE_LIKE, ipHYDROGEN, ipSatelliteLines, ipTFe56, iso_ctrl, iso_sp, t_opac::lgCaseB_no_pdest, t_isoCTRL::lgDielRecom, t_cosmology::lgDo, t_rfield::lgDoLineTrans, t_dense::lgElmtOn, t_conv::lgFirstSweepThisZone, t_conv::lgIonStageTrimed, t_conv::lgLastSweepThisZone, t_hydro::lgLymanPumping, lgTauGood(), t_trace::lgTrace, LIMELM, t_rfield::nfine, NISO, nLevel1, t_isoCTRL::nLyaLevel, t_isoCTRL::nLyman, t_conv::nPres2Ioniz, nSpecies, t_iso_sp::numLevels_local, t_iso_sp::numLevels_max, nUTA, nzone, opac, EmissionProxy::opacity(), t_opac::opacity_abs, EmissionProxy::Pdest(), EmissionProxy::PopOpc(), POW2, EmissionProxy::pump(), radius, t_cosmology::redshift_current, t_cosmology::redshift_start, rfield, RT_line_one(), RT_stark(), SatelliteLines, SDIV(), SMALLFLOAT, SPEEDLIGHT, t_iso_sp::st, TauLines, trace, t_iso_sp::trans(), UTALines, wind, Wind::windv, Wind::windv0, and t_dense::xIonDense.

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

Definition at line 16 of file rt_line_driving.cpp.

References dBaseSpecies, dBaseTrans, DEBUG_ENTRY, dense, diatoms, TransitionList::Emis(), TransitionProxy::Emis(), TransitionProxy::EnergyErg(), ExtraLymanLines, FeIIAccel(), HFLines, t_dense::IonHigh, ipExtraLymanLines, ipH_LIKE, ipHE_LIKE, ipSatelliteLines, iso_ctrl, iso_sp, LIMELM, nHFLines, NISO, nLevel1, t_isoCTRL::nLyman, nSpecies, t_iso_sp::numLevels_local, nUTA, nWindLine, EmissionProxy::PopOpc(), EmissionProxy::pump(), SatelliteLines, TauLine2, TauLines, t_iso_sp::trans(), and UTALines.

Referenced by PresTotCurrent().

Here is the call graph for this function:

void RT_line_one	(	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 392 of file rt_line_one.cpp.

References ASSERT, cdEXIT, conv, EmissionProxy::damp(), EmissionProxy::dampXvel(), DEBUG_ENTRY, TransitionProxy::Emis(), EN1RYD, TransitionProxy::EnergyErg(), EXIT_SUCCESS, EmissionProxy::FracInwd(), ioQQQ, ipLY_A, EmissionProxy::iRedisFun(), t_rfield::lgDoLineTrans, t_conv::lgLastSweepThisZone, lgTauGood(), TransitionProxy::Lo(), MAX2, MIN2, nzone, EmissionProxy::Pdest(), EmissionProxy::Pelec_esc(), EmissionProxy::Pesc(), t_rfield::plsfrq, EmissionProxy::pump(), rfield, RT_line_electron_scatter(), RT_line_escape(), RT_line_fine_opacity(), RT_line_pumping(), and SMALLFLOAT.

Referenced by CoolDima(), CoolEvaluate(), FeII_RT_Make(), diatomics::H2_RTMake(), and RT_line_all().

Here is the call graph for this function:

void RT_line_one_tau_reset ( const TransitionProxy & t )

Definition at line 12 of file rt_line_one_tau_reset.cpp.

References chLineLbl(), TransitionProxy::Coll(), CollisionProxy::cool(), DEBUG_ENTRY, t_rt::DoubleTau, TransitionProxy::Emis(), EmissionProxy::FracInwd(), geometry, CollisionProxy::heat(), TransitionProxy::Hi(), ioQQQ, TransitionProxy::ipCont(), iteration, t_prt::lgPrtMaser, t_geometry::lgSphere, t_geometry::lgStatic, TransitionProxy::Lo(), MAX2, min(), MIN2, opac, EmissionProxy::Pdest(), EmissionProxy::Pelec_esc(), EmissionProxy::Pesc(), EmissionProxy::phots(), EmissionProxy::PopOpc(), prt, rt, EmissionProxy::TauCon(), EmissionProxy::TauIn(), t_opac::taumin, EmissionProxy::TauTot(), EmissionProxy::TauTrack(), and EmissionProxy::xIntensity().

Referenced by FeII_RT_tau_reset(), 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 16 of file rt_line_one_tauinc.cpp.

References Wind::AccelGravity, Wind::AccelTotalOutward, t_rfield::anu, 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, fixit(), geometry, GetHubbleFactor(), 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, t_rfield::nfine, EmissionProxy::opacity(), t_rfield::plsfrq, EmissionProxy::PopOpc(), radius, t_cosmology::redshift_current, rfield, rt, SDIV(), EmissionProxy::TauCon(), EmissionProxy::TauIn(), EmissionProxy::TauTot(), wind, and Wind::windv.

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

Here is the call graph for this function:

double RT_LineWidth	(	const TransitionProxy &	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 918 of file rt_escprob.cpp.

References ASSERT, EmissionProxy::damp(), EmissionProxy::dampXvel(), DEBUG_ENTRY, dense, t_dense::eden, TransitionProxy::Emis(), esc_PRD_1side(), iteration, Wind::lgBallistic(), t_pressure::lgPradDen, MAX2, MIN2, opac, EmissionProxy::Pdest(), EmissionProxy::Pelec_esc(), PI, pow(), pressure, EmissionProxy::TauIn(), t_opac::taumin, EmissionProxy::TauTot(), wind, and Wind::windv0.

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

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

References ASSERT, EmissionProxy::Aul(), dBaseSpecies, dBaseTrans, DEBUG_ENTRY, dense, DEST0, diatoms, t_dense::eden, t_dense::EdenTrue, TransitionList::Emis(), TransitionProxy::Emis(), t_opac::ExpmTau, t_iso_sp::fb, FeII_OTS(), findspecieslocal(), fnzone, 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, nLevel1, 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, TauLines, 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 402 of file rt_ots.cpp.

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

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

void RT_OTS_ChkSum ( long int ipPnt )

RT_OTS_ChkSum sanity check confirms summed continua reflect contents of individuals

Definition at line 609 of file rt_ots.cpp.

References t_rfield::anu, t_rfield::chContLabel, t_rfield::chLineLabel, t_rfield::ConInterOut, t_rfield::ConOTS_local_OTS_rate, DEBUG_ENTRY, t_rfield::flux, 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 690 of file rt_ots.cpp.

References t_rfield::anu, ASSERT, t_rfield::chContLabel, t_rfield::chLineLabel, conv, DEBUG_ENTRY, dense, fnzone, t_dense::gas_phase, ioQQQ, ipHELIUM, ipOXYGEN, t_conv::lgConvIoniz(), t_rfield::line_count, t_rfield::nupper, 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 488 of file rt_ots.cpp.

References t_rfield::anu2, ASSERT, t_rfield::ConInterOut, t_rfield::ConOTS_local_OTS_rate, t_rfield::ConOTS_local_photons, t_rfield::ContBoltz, t_rfield::convoc, DEBUG_ENTRY, fixit(), t_rfield::flux, t_rfield::flux_accum, t_rfield::ipEnergyBremsThin, t_rfield::ipPlasma, t_rfield::lgKillOTSLine, t_rfield::lgOutOnly, MAX2, MIN2, t_rfield::nflux, t_rfield::OccNumbBremsCont, opac, t_opac::opacity_abs, t_rfield::otscon, t_rfield::otslin, t_rfield::outlin, t_rfield::outlin_noplot, rfield, SDIV(), SMALLDOUBLE, 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 )

zero these things out, called in zero

Definition at line 577 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::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 11 of file rt_recom_effic.cpp.

References t_rfield::anu, ASSERT, cdEXIT, t_rfield::chDffTrns, t_rfield::ContBoltz, cosmology, DEBUG_ENTRY, t_opac::E2TauAbsFace, esccon(), EXIT_FAILURE, t_opac::ExpZone, t_rfield::fine_opac_zone, 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::nupper, opac, t_opac::opacity_abs, t_opac::otsmin, phycon, t_cosmology::redshift_current, rfield, sexp(), SMALLFLOAT, SPEEDLIGHT, t_opac::TauAbsGeo, t_opac::taumin, t_phycon::te, TE1RYD, t_phycon::te_ryd, 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 14 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, MAX2, MIN2, NISO, t_isoCTRL::nLyaLevel, t_iso_sp::numLevels_local, t_iso_sp::numLevels_max, nzone, phycon, pow(), rt, SDIV(), t_phycon::sqrte, t_iso_sp::st, strkar(), EmissionProxy::TauIn(), t_phycon::te01, t_phycon::te03, t_phycon::te10, and t_iso_sp::trans().