cloudy  trunk
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Data Structures | Macros | Typedefs | Enumerations | Functions | Variables
rt.h File Reference
#include "module.h"
Include dependency graph for rt.h:
This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Data Structures

struct  t_rt




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


enum  {


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)


t_rt rt

Macro Definition Documentation


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


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


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


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

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

Referenced by RT_line_escape().

void RT_fine_clear ( )
void RT_line_all ( linefunc  line_one)
void RT_line_all_escape ( realnum error)
double RT_line_driving ( void  )
void RT_line_one_escape ( const TransitionProxy t,
bool  lgShield_this_zone,
realnum  pestrk,
realnum  DopplerWidth 

RT_line_one do line radiative transfer

tline structure
pestrkStark escape probability
lgShield_this_zonethis 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 381 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 
void RT_line_one_tau_reset ( const TransitionProxy t)
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_OTS ( void  )
void RT_OTS_AddLine ( double  ots,
long int  ip 

RT_OTS_AddLine add local destruction of lines to ots field

ippointer 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)
void RT_OTS_PrtRate ( double  weak,
int  chFlag 

RT_OTS_PrtRate print ots arrays, called by ionize

weakarg is weakest rate to print
chFlagflag, '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)
void RT_OTS_Zero ( void  )
double RT_recom_effic ( long int  ip)
void RT_stark ( void  )

rt_stark compute stark broadening escape probabilities using Puetter formalism

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 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, MAX2, MIN2, NISO, t_isoCTRL::nLyaLevel, t_iso_sp::numLevels_local, t_iso_sp::numLevels_max, nzone, phycon, pow(), powpq(), 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().

Referenced by RT_line_all_escape().

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 28 of file rt_tau_init.cpp.

References t_mesh::anu(), ASSERT, t_dense::AtomicWeight, EmissionProxy::Aul(), cdEXIT, t_StopCalc::colnut, conv, dBaseTrans, DEBUG_ENTRY, dense, diatoms, t_rt::DoubleTau, t_dense::eden, TransitionProxy::Emis(), EXIT_FAILURE, ExtraLymanLines, t_iso_sp::fb, fixit, fprintf(), EmissionProxy::FracInwd(), t_dense::gas_phase, geometry, GetDopplerWidth(), t_StopCalc::HColStop, HFLines, ioQQQ, ipExtraLymanLines, ipH1s, ipH2p, ipH2s, ipH3d, ipH3p, ipH3s, ipH_LIKE, ipHe1s1S, ipHe2p1P, ipHE_LIKE, ipHELIUM, ipHYDROGEN, t_trace::ipIsoTrace, ipoint(), ipSatelliteLines, t_StopCalc::iptnu, t_rt::ipxry, iso_ctrl, iso_sp, t_opac::lgCaseB, t_isoCTRL::lgDielRecom, t_dense::lgElmtOn, t_conv::lgFirstSweepThisZone, t_trace::lgIsoTraceFull, t_conv::lgOscilOTS, t_geometry::lgSphere, t_thermal::lgTemperatureConstant, t_trace::lgTrace, LIMELM, MAX2, MIN2, N_, t_iso_sp::n_HighestResolved_max, t_rfield::nflux_with_check, NISO, t_isoCTRL::nLyman, nSpecies, t_iso_sp::numLevels_local, t_iso_sp::numLevels_max, nWindLine, opac, EmissionProxy::opacity(), EmissionProxy::Ploss(), pow(), PrintEfmt, t_rfield::qhtot, t_iso_sp::QuantumNumbers2Index, rfield, rt, RT_line_all_escape(), S, SatelliteLines, ShowMe(), TransitionList::size(), t_isoCTRL::SmallA, SMALLFLOAT, StopCalc, SumDensities(), t_opac::TauAbsGeo, EmissionProxy::TauCon(), t_StopCalc::tauend, EmissionProxy::TauIn(), TAULIM, TauLine2, t_opac::taumin, t_StopCalc::taunu, t_opac::TauScatGeo, EmissionProxy::TauTot(), t_opac::TauTotalGeo, t_rt::tauxry, TempChange(), thermal, t_opac::tlamin, trace, t_iso_sp::trans(), UTALines, t_dense::xNucleiTotal, and TransitionProxy::Zero().

Referenced by cloudy().

Here is the call graph for this function:

void RT_tau_reset ( void  )

Variable Documentation

t_rt rt