Cloudy
Spectral Synthesis Code for Astrophysics
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
t_rfield Struct Reference

#include <rfield.h>

Inheritance diagram for t_rfield:
Collaboration diagram for t_rfield:

Public Member Functions

const char * chName () const
 
void zero ()
 
void comment (t_warnings &)
 
 t_rfield ()
 
const realnumgetCoarseTransCoef ()
 
void setCoarseTransCoefPtr (size_t size)
 
void resetCoarseTransCoef ()
 
void setTrimming ()
 
- Public Member Functions inherited from module
 module ()
 
virtual ~module ()
 
- Public Member Functions inherited from t_mesh
void InitMesh (bool lgUnitCell)
 
void ValidateEdges () const
 
void CheckMesh () const
 
bool lgMeshSetUp () const
 
long ncells () const
 
double emm () const
 
double egamry () const
 
double getResolutionScaleFactor () const
 
void setResolutionScaleFactor (double fac)
 
string mesh_md5sum () const
 
const double * anuptr () const
 
double anu (size_t i) const
 
double anu2 (size_t i) const
 
double anu3 (size_t i) const
 
double anuln (size_t i) const
 
const double * anulog10ptr () const
 
double anulog10 (size_t i) const
 
double anusqrt (size_t i) const
 
double anumin (size_t i) const
 
double anumax (size_t i) const
 
double widflx (size_t i) const
 
size_t ipointC (double anu) const
 
size_t ipointF (double anu) const
 
size_t ithreshC (double threshold) const
 
size_t ithreshF (double threshold) const
 
bool isEnergyBound (Energy en) const
 
 t_mesh ()
 

Public Attributes

long int nflux
 
long int nflux_with_check
 
long int nPositive
 
double FluxFaint
 
vector< long > line_count
 
vector< realnumOccNumbContEmitOut
 
Spectrum flux
 
vector< realnumflux_isotropic
 
vector< realnumflux_beam_time
 
vector< realnumflux_beam_const
 
vector< realnumflux_accum
 
vector< realnumExtinguishFactor
 
realnum ExtinguishLeakage
 
realnum ExtinguishColumnDensity
 
realnum ExtinguishLowEnergyLimit
 
realnum ExtinguishConvertColDen2OptDepth
 
realnum ExtinguishEnergyPowerLow
 
bool lgMustBlockHIon
 
bool lgBlockHIon
 
bool lgDoLineTrans
 
bool lgOpacityReevaluate
 
bool lgCMB_set
 
bool lgIonizReevaluate
 
vector< realnumconvoc
 
vector< realnumOccNumbIncidCont
 
vector< realnumOccNumbDiffCont
 
vector_avx< double > ContBoltz
 
vector_avx< double > ContBoltzHelp1
 
vector_avx< double > ContBoltzHelp2
 
vector_avx< double > vexp_arg
 
vector_avx< double > ContBoltzAvg
 
multi_arr< realnum, 2 > ConEmitLocal
 
multi_arr< realnum, 2 > ConSourceFcnLocal
 
Spectrum ConEmitReflec
 
Spectrum ConEmitOut
 
vector< realnumConInterOut
 
Spectrum ConRefIncid
 
vector< double > SummedCon
 
vector< realnumSummedDif
 
vector< realnumSummedOcc
 
vector< realnumSummedDifSave
 
vector< realnumConOTS_local_photons
 
vector< realnumConOTS_local_OTS_rate
 
vector< realnumDiffuseEscape
 
vector< realnumTotDiff2Pht
 
vector< realnumotslin
 
vector< realnumotscon
 
multi_arr< realnum, 2 > otssav
 
Spectrum outlin
 
vector< realnumoutlin_noplot
 
vector< realnumDiffuseLineEmission
 
Spectrum reflin
 
Spectrum flux_total_incident
 
vector< realnumflux_beam_const_save
 
vector< realnumflux_time_beam_save
 
vector< realnumflux_isotropic_save
 
realnum time_continuum_scale
 
realnum DiffPumpOn
 
vector< string > chLineLabel
 
vector< string > chContLabel
 
char chDffTrns [4]
 
bool lgOutOnly
 
long int ipEnergyBremsThin
 
realnum EnergyBremsThin
 
long int ipMaxBolt
 
bool lgInducProcess
 
vector< double > comup
 
vector< double > comdn
 
long int ipB_filter
 
long int ipV_filter
 
long int ipG0_TH85_lo
 
long int ipG0_TH85_hi
 
long int ipG0_DB96_lo
 
long int ipG0_DB96_hi
 
long int ipG0_spec_lo
 
long int ipG0_spec_hi
 
long int ip1000A
 
double extin_mag_B_point
 
double extin_mag_V_point
 
double extin_mag_B_extended
 
double extin_mag_V_extended
 
double opac_mag_B_point
 
double opac_mag_V_point
 
double opac_mag_B_extended
 
double opac_mag_V_extended
 
vector< realnumcsigh
 
vector< realnumcsigc
 
vector< realnumeeBremsDif
 
double comtot
 
double cmheat
 
double cmcool
 
double cinrat
 
bool lgComptonOn
 
bool lgComUndr
 
double totpow [LIMSPC]
 
double slope [LIMSPC]
 
double cutoff [LIMSPC][3]
 
double spfac [LIMSPC]
 
bool lgTimeVary [LIMSPC]
 
bool lgBeamed [LIMSPC]
 
realnum OpticalDepthScaleFactor [LIMSPC]
 
Illuminate::IlluminationType Illumination [LIMSPC]
 
long int nShape
 
long int ipSpec
 
vector< EnergytNu [LIMSPC]
 
vector< realnumtslop [LIMSPC]
 
vector< realnumtFluxLog [LIMSPC]
 
double TableRadius [LIMSPC]
 
bool lgSphericalDilution [LIMSPC]
 
long ncont [LIMSPC]
 
double RSFCheck [LIMSPC]
 
double range [LIMSPC][2]
 
char chSpNorm [LIMSPC][5]
 
char chRSpec [LIMSPC][5]
 
char chSpType [LIMSPC][6]
 
char chCumuType [5]
 
realnum qhtot
 
realnum qhe
 
realnum qheii
 
realnum qbal
 
realnum qrad
 
realnum qtot
 
realnum uh
 
realnum uheii
 
bool lgUSphON
 
realnum rstrom
 
bool lgHabing
 
long int fine_opac_nelem
 
long int fine_opac_nresolv
 
realnum fine_opac_velocity_width
 
bool trans_coef_total_stale
 
vector< long > ipnt_coarse_2_fine
 
realnum fine_ener_lo
 
realnum fine_ener_hi
 
long nfine
 
double fine_resol
 
vector< realnumfine_opac_zone
 
vector< realnumfine_opt_depth
 
vector< realnumfine_anu
 
unordered_map< long, vector
< int > > 
fine_lstack
 
long int ipFineConVelShift
 
bool lgOpacityFine
 
bool lgSaveOpacityFine
 
bool lgLyaOTS
 
bool lgHeIIOTS
 
bool lgKillOutLine
 
bool lgKillOutCont
 
bool lgKillOTSLine
 
bool lgPlasNu
 
realnum plsfrq
 
realnum plsfrqmax
 
long int nZonePlsFrqEval
 
long int ipPlasma
 
long int ipPlasmax
 
bool lgMMok
 
bool lgHPhtOK
 
bool lgXRayOK
 
bool lgGamrOK
 
realnum EnerGammaRay
 
long int ipEnerGammaRay
 
bool lgHionRad
 
realnum occmax
 
realnum occmnu
 
realnum tbrmax
 
realnum tbrmnu
 
realnum tbr4nu
 
realnum occ1nu
 
bool lgOcc1Hi
 
realnum EnergyIncidCont
 
realnum EnergyDiffCont
 

Private Attributes

vector< realnumtrans_coef_total
 

Constructor & Destructor Documentation

t_rfield::t_rfield ( )
inline

Member Function Documentation

const char* t_rfield::chName ( ) const
inlinevirtual

Implements module.

void t_rfield::comment ( t_warnings )
inlinevirtual

Implements module.

const realnum * t_rfield::getCoarseTransCoef ( )

References fine_opt_depth, ipnt_coarse_2_fine, t_opac::lgScatON, nflux, opac, sexp(), trans_coef_total, and trans_coef_total_stale.

Referenced by cdSPEC2(), SaveDo(), and sum_radiation().

Here is the call graph for this function:

void t_rfield::resetCoarseTransCoef ( )
inline
void t_rfield::setCoarseTransCoefPtr ( size_t  size)
inline

References trans_coef_total.

Referenced by rfield_opac_alloc().

void t_rfield::setTrimming ( )

References t_mesh::anu(), DEBUG_ENTRY, FluxFaint, max(), nflux, nPositive, SummedCon, and t_mesh::widflx().

Referenced by RT_continuum(), RT_OTS_Update(), RT_OTS_Zero(), and ZoneEnd().

Here is the call graph for this function:

Member Data Documentation

vector<string> t_rfield::chContLabel

string identifying the first continuum edge that occurred at this energy

Referenced by AGN_Hemis(), ContCreatePointers(), ContNegative(), ContSetIntensity(), GammaPrt(), ipContEnergy(), rfield_opac_alloc(), RT_continuum(), RT_OTS_ChkSum(), RT_OTS_PrtRate(), save_opacity(), and SaveDo().

char t_rfield::chCumuType[5]

Referenced by IterEnd(), ParseSet(), PrintSpectrum(), and zero().

char t_rfield::chDffTrns[4]

method for transferring diffuse continuum either 'ots' or 'oux' where x is n for which

Referenced by InitDefaultsPreparse(), ParseCosmology(), ParseDiffuse(), and RT_recom_effic().

vector<string> t_rfield::chLineLabel
char t_rfield::chSpNorm[LIMSPC][5]

chSpNorm says how spectrum was normalized - intensity or luminosity case chRSpec says whether per unit area or tot 4 pi

Referenced by conorm(), ParseAbsMag(), ParseBackgrd(), ParseBlackbody(), ParseCMB(), ParseEnergy(), ParseF_nu(), ParseIntensity(), ParseIonPar(), ParseLuminosity(), ParsePhi(), ParseQH(), ParseRatio(), and ParseTable().

double t_rfield::cinrat
double t_rfield::cmcool
double t_rfield::cmheat
vector<double> t_rfield::comdn

Referenced by highen(), rfield_opac_alloc(), and SaveDo().

double t_rfield::comtot
vector<double> t_rfield::comup

saves for the upward and downward Compton shifts

Referenced by highen(), rfield_opac_alloc(), and SaveDo().

multi_arr<realnum,2> t_rfield::ConEmitLocal

ConEmitLocal is local diffuse continuum, units photons cm-3 s-1 cell-1, for each zone, evaluated in RT_diffuse

Referenced by AGN_Hemis(), ContCreateMesh(), lines_continuum(), RT_continuum(), RT_diffuse(), RT_iso_integrate_RRC(), and SaveDo().

Spectrum t_rfield::ConEmitOut

outward diffuse emission continuum (not the interactive one), this is incremented in radinc because of interplay between absorption and emission - get the outward bremsstrahlung right photons cell-1 cm-2 s-1

Referenced by ContSetIntensity(), flxCell(), IterEnd(), IterRestart(), IterStart(), lines_continuum(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), SaveDo(), and sum_radiation().

vector<realnum> t_rfield::ConOTS_local_OTS_rate

the local photoionization rate corresponding to above photons

Referenced by GammaPrt(), IterRestart(), IterStart(), rfield_opac_alloc(), rfield_opac_zero(), RT_OTS_ChkSum(), RT_OTS_Update(), and RT_OTS_Zero().

vector<realnum> t_rfield::ConOTS_local_photons

this will control array of locally destroyed continuum photons, zeroed and evaluated in RT_diffuse, currently only two-photon

Referenced by IterRestart(), IterStart(), rfield_opac_alloc(), rfield_opac_zero(), RT_OTS_Update(), and RT_OTS_Zero().

Spectrum t_rfield::ConRefIncid

ConRefIncid is reflected portion of incident continuum

Referenced by ContSetIntensity(), flxCell(), IterEnd(), IterRestart(), lines_continuum(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), SaveDo(), and sum_radiation().

multi_arr<realnum,2> t_rfield::ConSourceFcnLocal

the local source function - diffuse emission, photons cell-1 cm-2 s-1

Referenced by ContCreateMesh(), RT_diffuse(), and SaveDo().

vector_avx<double> t_rfield::ContBoltz
vector_avx<double> t_rfield::ContBoltzAvg

array of Boltzmann factors for the continuum energy grid and current temperature properly averaged over the width of the cell

Referenced by ContSetIntensity(), GammaBn(), rfield_opac_alloc(), rfield_opac_zero(), and tfidle().

vector_avx<double> t_rfield::ContBoltzHelp1

helper arrays for calculating Boltzmann factors averaged over the width of the cell

Referenced by ContSetIntensity(), iso_rad_rec_cooling_discrete(), rfield_opac_alloc(), rfield_opac_zero(), RT_iso_integrate_RRC(), and tfidle().

vector<realnum> t_rfield::convoc
vector<realnum> t_rfield::csigc
vector<realnum> t_rfield::csigh

coefficients for fitting Tarter expressions for Compton heating and cooling, over full energy range of code

Referenced by ContCreateMesh(), ContSetIntensity(), highen(), PrtHeader(), and rfield_opac_alloc().

realnum t_rfield::DiffPumpOn

this is zero or one depending whether pumping by diffuse fields is turned on (1) or turned off (0)

Referenced by ParseDont(), and zero().

vector<realnum> t_rfield::DiffuseEscape

computed in rt_diffuse.cpp, escaping continuum emission added to beam in radius_increment

Referenced by IterRestart(), IterStart(), rfield_opac_alloc(), RT_continuum(), RT_diffuse(), and RT_iso_integrate_RRC().

vector<realnum> t_rfield::DiffuseLineEmission

local diffuse line emission, photons cm-3 s-1

Referenced by lines_continuum(), outline_base_bin(), rfield_opac_alloc(), RT_diffuse(), and SaveDo().

vector<realnum> t_rfield::eeBremsDif

ee brems emission, still needs to be multiplied with pow2(dense.eden)

Referenced by OpacityAddTotal(), rfield_opac_alloc(), and RT_diffuse().

realnum t_rfield::EnerGammaRay

lowest energy counted as gamma rays, 100 keV

Referenced by ContCreatePointers(), InitSimPostparse(), and ParseInterp().

realnum t_rfield::EnergyBremsThin
realnum t_rfield::EnergyDiffCont
realnum t_rfield::EnergyIncidCont

intensity [erg cm-2 s-1] in incident and diffuse continua in current zone

Referenced by lgCheckMonitors(), PrtZone(), and RT_continuum().

double t_rfield::extin_mag_B_extended

extinction in magnitudes at B and V filters for a resolved source, this does discount forward scattering, so is appropriate for a resolved source

Referenced by InitSimPostparse(), IterRestart(), IterStart(), and radius_increment().

double t_rfield::extin_mag_B_point

extinction in magnitudes at B and V filters for a point source, this does not discount forward scattering

Referenced by InitSimPostparse(), IterRestart(), IterStart(), PrtFinal(), and radius_increment().

realnum t_rfield::ExtinguishColumnDensity
realnum t_rfield::ExtinguishConvertColDen2OptDepth

the constant that multiplies the column density to get optical depth

Referenced by extin(), and t_rfield().

realnum t_rfield::ExtinguishEnergyPowerLow

the power on the energy

Referenced by extin(), and t_rfield().

vector<realnum> t_rfield::ExtinguishFactor

extinction factor set with extinguish command

Referenced by extin(), IncidentContinuumHere(), and rfield_opac_alloc().

realnum t_rfield::ExtinguishLeakage
realnum t_rfield::ExtinguishLowEnergyLimit
vector<realnum> t_rfield::fine_anu

energies at center of each bin for fine continuum

Referenced by ContCreateMesh(), ContCreatePointers(), print_emline_fine(), rfield_opac_alloc(), save_opacity(), and SaveDo().

realnum t_rfield::fine_ener_hi
realnum t_rfield::fine_ener_lo

low and high bounds of fine continuum - set by need to include all resonance lines

Referenced by ContCreateMesh(), ipFineCont(), ParseSet(), and rfield_opac_alloc().

unordered_map<long, vector<int> > t_rfield::fine_lstack

map of list of line stack indices (value) centered on fine opacity vector index (key)

Referenced by lincom(), LineStackCreate(), and SaveDo().

long int t_rfield::fine_opac_nelem

heaviest element to be considered - default is iron - the fine opacity array resolution depends on heaviest element and lowest temperature

Referenced by ParseSet(), rfield_opac_alloc(), and zero().

long int t_rfield::fine_opac_nresolv

number of resolution elements over width of this element, default is 4

Referenced by ParseSet(), rfield_opac_alloc(), and zero().

realnum t_rfield::fine_opac_velocity_width

number of cm/s of each cell in fine mesh

Referenced by PrtFinal(), rfield_opac_alloc(), RT_fine_clear(), and RT_line_fine_opacity().

vector<realnum> t_rfield::fine_opt_depth

total optical depth array for fine continuum

Referenced by getCoarseTransCoef(), rfield_opac_alloc(), rfield_opac_zero(), RT_tau_inc(), RT_tau_reset(), and SaveDo().

double t_rfield::fine_resol

the dimensionless resolution of the fine continuum - dE/E

Referenced by ipFineCont(), and rfield_opac_alloc().

vector<realnum> t_rfield::flux_accum

the accumulated flux, sum from this energy to infinity

Referenced by diatomics::H2_X_sink_and_source(), rfield_opac_alloc(), rfield_opac_zero(), and RT_OTS_Update().

vector<realnum> t_rfield::flux_beam_const_save
vector<realnum> t_rfield::flux_beam_time

this is the variable and constant parts of the above

Referenced by ContSetIntensity(), extin(), flux_correct_isotropic(), IterRestart(), rfield_opac_alloc(), RT_continuum(), RT_diffuse(), ZoneEnd(), and ZoneStart().

vector<realnum> t_rfield::flux_isotropic

this is the isotropic part of the constant continuum

Referenced by ContSetIntensity(), extin(), IterRestart(), rfield_opac_alloc(), RT_continuum(), RT_diffuse(), set_xIntensity(), ZoneEnd(), and ZoneStart().

vector<realnum> t_rfield::flux_isotropic_save
vector<realnum> t_rfield::flux_time_beam_save
Spectrum t_rfield::flux_total_incident

save incident continuum for later iterations

Referenced by ContSetIntensity(), flxCell(), IterEnd(), IterRestart(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), and SaveDo().

double t_rfield::FluxFaint

faintest high energy flux to consider, set with set flxfnt command

Referenced by InitDefaultsPreparse(), ParseSet(), setTrimming(), and ZoneEnd().

long int t_rfield::ip1000A

this is the wavelength where Bertoldi & Draine estimate the Habing field

Referenced by ContCreatePointers(), and mole_h_reactions().

long int t_rfield::ipB_filter

array indices for centers of B and V filters

Referenced by ContCreatePointers(), and radius_increment().

long int t_rfield::ipEnerGammaRay
long int t_rfield::ipEnergyBremsThin

ipEnergyBremsThin is index for lowest energy thin to ff abs and plasma frequency EnergyBremsThin is energy there, Ryd

Referenced by t_gaunt::brems_cool(), ContRate(), CoolEvaluate(), GrainRateDr(), InitSimPostparse(), IterRestart(), IterStart(), RT_OTS_ChkSum(), and tauff().

long int t_rfield::ipFineConVelShift

shift in fine continuum rest frame scale due to velocity gradient, rest frame is velocity of first zone, positive means that first zone is blue shifted relative to current zone. This is a decelerating flow. Evaluated in RT_line_all

Referenced by PressureRadiationLine(), rfield_opac_alloc(), RT_fine_clear(), RT_line_escape(), RT_line_fine_opacity(), RT_line_one_tauinc(), and SaveDo().

long int t_rfield::ipG0_DB96_hi
long int t_rfield::ipG0_DB96_lo

these are the lower and upper bounds for the G0 radiation field used by Tielens & Hollenbach in their PDR work

Referenced by ContCreatePointers(), ContSetIntensity(), and mole_h_reactions().

long int t_rfield::ipG0_spec_hi
long int t_rfield::ipG0_spec_lo

these are the lower and upper bounds for the special G0 radiation field

Referenced by ContCreatePointers(), and mole_h_reactions().

long int t_rfield::ipG0_TH85_hi
long int t_rfield::ipG0_TH85_lo

these are the lower and upper bounds for the G0 radiation field used by Tielens & Hollenbach in their PDR work

Referenced by ContCreatePointers(), and mole_h_reactions().

long int t_rfield::ipMaxBolt

index of highest cell with positive Boltzmann factor

Referenced by t_gaunt::brems_cool(), t_gaunt::brems_rt(), CoolEvaluate(), t_gaunt::p_setup_brems(), RT_diffuse(), SaveDo(), and tfidle().

vector<long> t_rfield::ipnt_coarse_2_fine

and array indices for lower and upper bounds of each coarse continuum mapped onto the fine continuum 0 (false) if fine continuum does not extend to these cells

Referenced by ContCreateMesh(), getCoarseTransCoef(), rfield_opac_alloc(), and RT_recom_effic().

long int t_rfield::ipPlasma

pointer to current plasma frequency

Referenced by ContRate(), ContSetIntensity(), EdenChange(), RT_continuum(), RT_diffuse(), and RT_OTS_Update().

long int t_rfield::ipPlasmax

pointer to largest plasma freq encountered

Referenced by ContSetIntensity().

long int t_rfield::ipSpec
long int t_rfield::ipV_filter
bool t_rfield::lgBlockHIon
bool t_rfield::lgCMB_set

this flag says that CMB has been set

Referenced by InitDefaultsPreparse(), ParseBackgrd(), ParseCMB(), and PrtComment().

bool t_rfield::lgComptonOn
bool t_rfield::lgComUndr

set true if Compton cooling underflows

Referenced by PrtComment(), and PrtHeader().

bool t_rfield::lgDoLineTrans

option to not do line transfer, set false with no line transfer command

Referenced by outline_base(), ParseCommands(), ParseDont(), PresTotCurrent(), RT_line_all(), RT_line_one_escape(), RT_line_one_fine(), and zero().

bool t_rfield::lgGamrOK
bool t_rfield::lgHabing

flag set if incident radiation field less than 10x the Habing ISM field

Referenced by ContSetIntensity(), PrtComment(), and zero().

bool t_rfield::lgHeIIOTS

flag saying whether to include HeII Lya and rec cont ots in the radiation field - - usually true, set false with no HeII ots command

Referenced by ParseDont(), RT_OTS(), and zero().

bool t_rfield::lgHionRad

lgHionRad set to .true. if no hydrogen ionizing radiation

Referenced by ContSetIntensity(), and PrtComment().

bool t_rfield::lgHPhtOK
bool t_rfield::lgInducProcess
bool t_rfield::lgIonizReevaluate

says whether to constantly reevaluate ionization, normally true, set false with no ionization reevaluate command

Referenced by ConvBase(), ParseCommands(), ParseDont(), and zero().

bool t_rfield::lgKillOTSLine

flag saying whether to kill ots lines

Referenced by ParseDont(), RT_OTS_Update(), and zero().

bool t_rfield::lgKillOutCont

flag saying whether to kill outward only continuum

Referenced by ParseDont(), RT_diffuse(), and zero().

bool t_rfield::lgKillOutLine

flag saying whether to kill outward only lines

Referenced by ParseDont(), RT_diffuse(), and zero().

bool t_rfield::lgLyaOTS

flag saying whether to include H1 Lya ots in the radiation field - usually true, set false with no Lya ots command

Referenced by ParseDont(), RT_OTS(), and zero().

bool t_rfield::lgMMok

these are series of flags that say whether different parts of the continuum where entered ok

Referenced by InitDefaultsPreparse(), ParseInterp(), and PrtComment().

bool t_rfield::lgMustBlockHIon

this is set true is one of incident continua is expected to have all H-ionizing radiation blocked. That is done separately with the extinguish command.

Referenced by ParseCommands(), ParseTable(), and zero().

bool t_rfield::lgOcc1Hi

flag saying that occupation number at 1 Ryd > 1

Referenced by ContSetIntensity(), and PrtComment().

bool t_rfield::lgOpacityFine
bool t_rfield::lgOpacityReevaluate

says whether to constantly reevaluate opacities, normally true, set false with no opacity reevaluate command

Referenced by ConvBase(), ParseCommands(), ParseDont(), and zero().

bool t_rfield::lgOutOnly

another flag, true if outward only - used to multiply the ConInterOut continuum for creating the interactive continuum - when false, no outward only, this is not added

Referenced by GammaBn(), GammaK(), GammaPrt(), InitDefaultsPreparse(), ParseDiffuse(), RT_OTS_ChkSum(), and RT_OTS_Update().

bool t_rfield::lgPlasNu

following deal with plasma frequency, which enters the continuum array for even moderate densities due to very low frequencies considered set true if plasma freq enters energy array

Referenced by ContSetIntensity(), EdenChange(), iso_level(), and PrtComment().

bool t_rfield::lgSaveOpacityFine

says that fine optical depths will be saved, so save them

Referenced by ParseSave(), and zero().

bool t_rfield::lgSphericalDilution[LIMSPC]

does TABLE READ SED need additional spherical dilution?

Referenced by InitDefaultsPreparse(), ParseCommands(), and ParseTable().

bool t_rfield::lgUSphON

lgUSphON flag set when we hit Stromgren radius in spherical geometry

Referenced by dense_fabden(), InitSimPostparse(), IterRestart(), optimize_func(), PrtFinal(), and radius_increment().

bool t_rfield::lgXRayOK
vector<long> t_rfield::line_count

used to keep track of number of lines per freq interval

Referenced by ipLineEnergy(), rfield_opac_alloc(), RT_OTS_PrtRate(), and SaveDo().

long t_rfield::ncont[LIMSPC]
long t_rfield::nfine
long int t_rfield::nflux

================================================================================= the following define the continuum energy scale and its limits nflux is number of continuum points needed to get to high energy end of continuum. this excludes the unity test cell.

Referenced by AGN_Hemis(), t_gaunt::brems_cool(), t_gaunt::brems_opac(), t_gaunt::brems_rt(), cdDrive(), cdSPEC2(), cmshft(), ContCreateMesh(), ContNegative(), ContRate(), ContSetIntensity(), ConvInitSolution(), CoolEvaluate(), eeBremsSpectrum(), ffun1(), FindStrongestLineLabels(), flux_correct_isotropic(), GammaBn(), GammaK(), GammaPrt(), getCoarseTransCoef(), GetDissociationRateCoeff(), diatomics::GetHeatRate(), GrainCharge(), GrainMakeDiffuse(), GrainUpdateRadius2(), GridRetrieveXSPECData(), highen(), IncidentContinuumHere(), ion_photo(), ion_trim(), ion_trim_init(), iso_photo(), iso_rad_rec_cooling_discrete(), IterEnd(), lincom(), OpacityAdd1Element(), OpacityAdd1Subshell(), OpacityAdd1SubshellInduc(), OpacityAddTotal(), OpacityCreateAll(), outsum(), t_gaunt::p_setup_brems(), pnegopc(), PresTotCurrent(), process_output(), PrtComment(), PrtFinal(), PrtHeader(), punchFITS_GenericData(), punchFITS_GenericHeader(), radius_first(), RT_continuum(), RT_DestProb(), RT_diffuse(), RT_iso_integrate_RRC(), RT_OTS(), RT_OTS_AddCont(), RT_OTS_AddLine(), RT_OTS_ChkSum(), RT_OTS_Update(), RT_OTS_Zero(), RT_recom_effic(), SanityCheckBegin(), save_opacity(), SaveDo(), saveFITSfile(), setTrimming(), sum_radiation(), sumcon(), tauff(), UpdatePot(), UpdatePot1(), ZoneEnd(), and ZoneStart().

long int t_rfield::nPositive

number of cells to include highest continuum cell with non-zero photon flux, loops over continuum should go to <nPositive to include positive cells

Referenced by cdDrive(), ContCreateMesh(), ContSetIntensity(), ConvInitSolution(), GammaBn(), GammaK(), GrainDrift(), GrainElecEmis1(), GrainsInit(), GrainTemperature(), highen(), PE_init(), qheat_init(), setTrimming(), UpdatePot(), UpdatePot1(), and ZoneEnd().

long int t_rfield::nShape

nShape is SED shape index number, this must equal the number of field intensities that are specified ipSpec is radiation field source number

Referenced by conorm(), ContSetIntensity(), DynaIterStart(), ffun(), ffun1(), InterpolateGridCoStar(), InterpolateRectGrid(), ParseAgn(), ParseBackgrd(), ParseBlackbody(), ParseBremsstrahlung(), ParseCMB(), ParseCommands(), ParseIlluminate(), ParseInterp(), ParseLaser(), ParsePowerlawContinuum(), ParseTable(), and ReadTable().

long int t_rfield::nZonePlsFrqEval

Referenced by EdenChange(), and t_rfield().

realnum t_rfield::occ1nu

Referenced by ContSetIntensity(), and PrtComment().

realnum t_rfield::occmax

these store photon occupation numbers at various energies in the continuum

Referenced by ContSetIntensity(), and PrtComment().

realnum t_rfield::occmnu

Referenced by ContSetIntensity(), and PrtComment().

vector<realnum> t_rfield::OccNumbContEmitOut
vector<realnum> t_rfield::OccNumbDiffCont

OccNumbDiffCont is the continuum occupation number, for local diffuse continuum

Referenced by IterRestart(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), and RT_diffuse().

vector<realnum> t_rfield::OccNumbIncidCont

OccNumbIncidCont is the continuum occupation number for the attenuated incident ONLY

Referenced by ContSetIntensity(), highen(), iso_ionize_recombine(), IterRestart(), PrtHeader(), rfield_opac_alloc(), RT_continuum(), RT_line_pumping(), and SaveDo().

double t_rfield::opac_mag_B_extended

Referenced by radius_increment().

double t_rfield::opac_mag_B_point

these are total opacities at these wavelengths, used to stop at exact Av

Referenced by radius_increment().

double t_rfield::opac_mag_V_extended

Referenced by radius_increment(), and radius_next().

double t_rfield::opac_mag_V_point

Referenced by radius_increment(), and radius_next().

realnum t_rfield::OpticalDepthScaleFactor[LIMSPC]

1 / cos( illumination angle, angle measured from normal, default is angle=zero, normal illumination, DirectCos = 1

Referenced by ContSetIntensity(), InitDefaultsPreparse(), and ParseIlluminate().

vector<realnum> t_rfield::otslin
multi_arr<realnum,2> t_rfield::otssav
realnum t_rfield::plsfrqmax

store highest energy plasma frequency encountered

Referenced by ContSetIntensity(), EdenChange(), and PrtComment().

realnum t_rfield::qbal
realnum t_rfield::qhe

Referenced by ContSetIntensity(), and PrtHeader().

realnum t_rfield::qheii

Referenced by ContSetIntensity(), and PrtHeader().

realnum t_rfield::qhtot

these are total numbers of photons over various energy ranges

Referenced by ContSetIntensity(), lines_hydro(), PrtFinal(), PrtHeader(), radius_first(), and RT_tau_init().

realnum t_rfield::qrad

Referenced by ContSetIntensity().

realnum t_rfield::qtot

Referenced by ContSetIntensity(), and IterStart().

double t_rfield::range[LIMSPC][2]

energy range over which the intensity is integrated for normalizing each continuum source that contributes to the total source

Referenced by conorm(), InitDefaultsPreparse(), ParseAbsMag(), ParseBackgrd(), ParseBlackbody(), ParseCMB(), ParseEnergy(), ParseF_nu(), ParseIntensity(), ParseIonPar(), ParseLuminosity(), ParsePhi(), ParseQH(), ParseRangeOption(), ParseRatio(), and ParseTable().

double t_rfield::RSFCheck[LIMSPC]

used to store a check on the continuum mesh resolution scale factor this is used for the output of the SAVE TRANSMITTED CONTINUUM command

Referenced by conorm(), InitDefaultsPreparse(), and ReadTable().

realnum t_rfield::rstrom

the Stromgren radius var set to get u spherical

Referenced by dense_fabden(), optimize_func(), PrtFinal(), and radius_increment().

double t_rfield::spfac[LIMSPC]

Referenced by conorm(), and ffun().

vector<double> t_rfield::SummedCon
vector<realnum> t_rfield::SummedDifSave
double t_rfield::TableRadius[LIMSPC]

TABLE READ input is normalized on this radius

Referenced by InitDefaultsPreparse(), ParseCommands(), and ReadTable().

realnum t_rfield::tbr4nu

Referenced by ContSetIntensity(), and PrtComment().

realnum t_rfield::tbrmax

Referenced by ContSetIntensity(), and PrtComment().

realnum t_rfield::tbrmnu

Referenced by ContSetIntensity(), and PrtComment().

vector<realnum> t_rfield::tFluxLog[LIMSPC]

this is the log of f_nu of continuum point

Referenced by ffun1(), InitDefaultsPreparse(), ParseInterp(), ParseTable(), and ReadTable().

realnum t_rfield::time_continuum_scale

==1 for time steady, when continuum varies with time, is scale factor

Referenced by DynaIterEnd(), DynaIterStart(), DynaPunchTimeDep(), IterRestart(), and zero().

vector<Energy> t_rfield::tNu[LIMSPC]

these are used for interpolate and table commands, all have two indices, continuum and frequency tNuRyd is the linear energy Rydberg of continuum point must be explicit arrays again so that table commands will work before continuum defined.

Referenced by conorm(), ffun1(), InitDefaultsPreparse(), InterpolateGridCoStar(), InterpolateRectGrid(), ParseInterp(), ParseTable(), and ReadTable().

vector<realnum> t_rfield::TotDiff2Pht

saves total two photon continuum for debugging, set in RT_diffuse

Referenced by rfield_opac_alloc(), RT_diffuse(), and SaveDo().

vector<realnum> t_rfield::trans_coef_total
private

these are TOTAL transmission coefficients for fine opacity degraded to coarse continuum

Referenced by getCoarseTransCoef(), resetCoarseTransCoef(), and setCoarseTransCoefPtr().

bool t_rfield::trans_coef_total_stale
realnum t_rfield::uh

hydrogen ionization parameter

Referenced by ContSetIntensity(), ion_trim_init(), PrtComment(), PrtFinal(), PrtHeader(), and radius_first().

realnum t_rfield::uheii

helium ion ionization parameter

Referenced by ContSetIntensity(), PrtFinal(), and PrtHeader().

vector_avx<double> t_rfield::vexp_arg

helper array for vexp() calls

Referenced by rfield_opac_alloc(), RT_continuum(), and tfidle().


The documentation for this struct was generated from the following files: