Cloudy
Spectral Synthesis Code for Astrophysics
|
#include <rfield.h>
Public Member Functions | |
const char * | chName () const |
void | zero () |
void | comment (t_warnings &) |
t_rfield () | |
const realnum * | getCoarseTransCoef () |
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 () | |
Private Attributes | |
vector< realnum > | trans_coef_total |
|
inline |
References ExtinguishConvertColDen2OptDepth, ExtinguishEnergyPowerLow, and nZonePlsFrqEval.
|
inlinevirtual |
Implements module.
|
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().
|
inline |
References nflux_with_check, trans_coef_total, and trans_coef_total_stale.
Referenced by ContCreateMesh(), IterRestart(), and RT_OTS_Zero().
|
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().
|
virtual |
Implements module.
References chCumuType, DEBUG_ENTRY, DiffPumpOn, fine_opac_nelem, fine_opac_nresolv, ipIRON, lgBlockHIon, lgDoLineTrans, lgHabing, lgHeIIOTS, lgIonizReevaluate, lgKillOTSLine, lgKillOutCont, lgKillOutLine, lgLyaOTS, lgMustBlockHIon, lgOpacityReevaluate, lgSaveOpacityFine, and time_continuum_scale.
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 |
string identifying the first line that occurred at this energy
Referenced by ContCreatePointers(), ContNegative(), ContSetIntensity(), FindStrongestLineLabels(), GammaPrt(), ipLineEnergy(), rfield_opac_alloc(), RT_continuum(), RT_OTS_ChkSum(), RT_OTS_PrtRate(), and SaveDo().
char t_rfield::chRSpec[LIMSPC][5] |
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().
char t_rfield::chSpType[LIMSPC][6] |
double t_rfield::cinrat |
Referenced by highen(), InitSimPostparse(), and PrtComment().
double t_rfield::cmcool |
Referenced by CoolEvaluate(), highen(), InitSimPostparse(), PrtHeader(), and PrtZone().
double t_rfield::cmheat |
Referenced by HeatSum(), highen(), lines(), lines_general(), PrtHeader(), and PrtZone().
vector<double> t_rfield::comdn |
Referenced by highen(), rfield_opac_alloc(), and SaveDo().
double t_rfield::comtot |
Referenced by InitSimPostparse(), IterStart(), lines(), and PrtComment().
vector<double> t_rfield::comup |
saves for the upward and downward Compton shifts
Referenced by highen(), rfield_opac_alloc(), and SaveDo().
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().
Spectrum t_rfield::ConEmitReflec |
reflected diffuse emission continuum
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::ConInterOut |
this is set in RT_diffuse and carries interactive continua
Referenced by ContNegative(), GammaBn(), GammaK(), GammaPrt(), GetDissociationRateCoeff(), diatomics::GetHeatRate(), GrainDrift(), diatomics::H2_LevelPops(), highen(), IterEnd(), IterRestart(), IterStart(), mole_h_reactions(), outsum(), PresTotCurrent(), PrtComment(), PrtFinal(), PrtHeader(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), RT_diffuse(), RT_OTS_ChkSum(), RT_OTS_Update(), RT_OTS_Zero(), and SaveDo().
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().
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 |
array of Boltzmann factors for the continuum energy grid and current temperature
Referenced by atom_level2(), ContSetIntensity(), eeBremsSpectrum(), H21_cm_pops(), OpacityAdd1SubshellInduc(), OpacityAddTotal(), t_gaunt::p_setup_brems(), plankf(), rfield_opac_alloc(), rfield_opac_zero(), RT_diffuse(), RT_recom_effic(), RT_tau_inc(), and tfidle().
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_avx<double> t_rfield::ContBoltzHelp2 |
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 |
convoc is the conversion factor from rfield to OccNumbIncidCont
Referenced by ContSetIntensity(), IterRestart(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), RT_diffuse(), RT_OTS_Update(), RT_OTS_Zero(), set_xIntensity(), ZoneEnd(), and ZoneStart().
vector<realnum> t_rfield::csigc |
Referenced by ContCreateMesh(), ContSetIntensity(), highen(), PrtHeader(), and rfield_opac_alloc().
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().
double t_rfield::cutoff[LIMSPC][3] |
Referenced by ffun1(), ParseAgn(), ParseBackgrd(), ParseBlackbody(), ParseBremsstrahlung(), ParseCMB(), ParseLaser(), and ParsePowerlawContinuum().
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 |
Referenced by InitSimPostparse(), IterRestart(), IterStart(), PrtComment(), PrtFinal(), and tauff().
realnum t_rfield::EnergyDiffCont |
Referenced by lgCheckMonitors(), PrtZone(), and RT_continuum().
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().
double t_rfield::extin_mag_V_extended |
Referenced by InitSimPostparse(), iter_end_check(), IterRestart(), IterStart(), matchGeneric(), PrtFinal(), radius_increment(), radius_next(), and SaveDo().
double t_rfield::extin_mag_V_point |
realnum t_rfield::ExtinguishColumnDensity |
Referenced by extin(), InitDefaultsPreparse(), and ParseExtinguish().
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 |
Referenced by extin(), InitDefaultsPreparse(), and ParseExtinguish().
realnum t_rfield::ExtinguishLowEnergyLimit |
Referenced by extin(), InitDefaultsPreparse(), and ParseExtinguish().
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 |
Referenced by ContCreateMesh(), ipFineCont(), ParseSet(), and rfield_opac_alloc().
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_opac_zone |
the fine continuum opacity array
Referenced by PressureRadiationLine(), rfield_opac_alloc(), rfield_opac_zero(), RT_fine_clear(), RT_line_escape(), RT_line_fine_opacity(), RT_line_one_tauinc(), RT_recom_effic(), RT_tau_inc(), RT_tau_reset(), save_opacity(), and SaveDo().
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().
Spectrum t_rfield::flux |
================================================================================= the following are the arrays containing the local radiation field flux is photons per cell N.B. width of cells vary with energy, given by widflx
Referenced by ContNegative(), ContRate(), ContSetIntensity(), CoolEvaluate(), extin(), flux_correct_isotropic(), GammaBn(), GammaK(), GammaPrt(), GetDissociationRateCoeff(), diatomics::GetHeatRate(), GrainCharge(), GrainDrift(), GrainElecEmis1(), GrainRateDr(), GrainTemperature(), diatomics::H2_LevelPops(), highen(), IncidentContinuumHere(), IterEnd(), IterRestart(), lines_continuum(), mole_h_reactions(), outsum(), PresTotCurrent(), PrtFinal(), PrtHeader(), qintr(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), RT_diffuse(), RT_OTS_ChkSum(), RT_OTS_Update(), RT_OTS_Zero(), SaveDo(), sum_radiation(), sumcon(), ZoneEnd(), and ZoneStart().
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 |
Referenced by ContSetIntensity(), extin(), flux_correct_isotropic(), IterRestart(), rfield_opac_alloc(), RT_continuum(), RT_diffuse(), ZoneEnd(), and ZoneStart().
vector<realnum> t_rfield::flux_beam_const_save |
Referenced by ContSetIntensity(), IterRestart(), rfield_opac_alloc(), and rfield_opac_zero().
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 |
Referenced by ContSetIntensity(), IterRestart(), rfield_opac_alloc(), and rfield_opac_zero().
vector<realnum> t_rfield::flux_time_beam_save |
Referenced by ContSetIntensity(), IterRestart(), rfield_opac_alloc(), and rfield_opac_zero().
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().
Illuminate::IlluminationType t_rfield::Illumination[LIMSPC] |
Referenced by ContSetIntensity(), InitDefaultsPreparse(), ParseBackgrd(), ParseBlackbody(), ParseCMB(), ParseIlluminate(), and ParseTable().
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 |
Referenced by ContCreatePointers(), and ContSetIntensity().
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 |
Referenced by ContCreatePointers(), ContSetIntensity(), and mole_h_reactions().
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 |
Referenced by ContCreatePointers(), and mole_h_reactions().
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 |
Referenced by ContCreatePointers(), and mole_h_reactions().
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 |
Referenced by conorm(), ContSetIntensity(), ffun(), ffun1(), and qintr().
long int t_rfield::ipV_filter |
Referenced by ContCreatePointers(), radius_increment(), and SaveDo().
bool t_rfield::lgBeamed[LIMSPC] |
Referenced by ContSetIntensity(), ffun(), InitDefaultsPreparse(), ParseBackgrd(), ParseCMB(), ParseIlluminate(), and ParseTable().
bool t_rfield::lgBlockHIon |
Referenced by ParseCommands(), ParseExtinguish(), and zero().
bool t_rfield::lgCMB_set |
this flag says that CMB has been set
Referenced by InitDefaultsPreparse(), ParseBackgrd(), ParseCMB(), and PrtComment().
bool t_rfield::lgComptonOn |
Referenced by cmshft(), highen(), InitDefaultsPreparse(), and ParseDont().
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 |
Referenced by InitDefaultsPreparse(), ParseInterp(), and PrtComment().
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 |
Referenced by InitDefaultsPreparse(), ParseInterp(), and PrtComment().
bool t_rfield::lgInducProcess |
turn off continuum pumping, set with 'no induced processes' command
Referenced by GammaBn(), diatomics::H2_LevelPops(), InitDefaultsPreparse(), iso_update_rates(), OpacityAdd1SubshellInduc(), ParseDont(), RT_diffuse(), RT_line_pumping(), RT_OTS_Update(), and setXtraRatesO1().
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 |
option to turn off fine opacities with no fine opacity command
Referenced by InitDefaultsPreparse(), ParseCommands(), ParseDont(), PressureRadiationLine(), RT_line_escape(), RT_line_fine_opacity(), RT_line_one_tauinc(), RT_recom_effic(), and RT_tau_inc().
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::lgTimeVary[LIMSPC] |
option to have continuum intensity be time dependent
Referenced by ContSetIntensity(), DynaIterStart(), ffun(), InitDefaultsPreparse(), ParseAbsMag(), ParseCMBOuter(), ParseEnergy(), ParseF_nu(), ParseIntensity(), ParseIonPar(), ParseLuminosity(), ParsePhi(), ParseQH(), ParseRatio(), 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 |
Referenced by InitDefaultsPreparse(), ParseInterp(), and PrtComment().
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] |
Referenced by conorm(), ffun1(), InitDefaultsPreparse(), ParseInterp(), ParseTable(), and ReadTable().
long t_rfield::nfine |
the number of fine continuum cells - will be several million
Referenced by ContCreateMesh(), ipFineCont(), PressureRadiationLine(), rfield_opac_alloc(), RT_line_escape(), RT_line_fine_opacity(), RT_line_one_tauinc(), RT_recom_effic(), RT_tau_inc(), save_opacity(), and SaveDo().
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::nflux_with_check |
number of frequency cells including unity test cell
Referenced by AtlasInterpolate(), diatomics::CalcPhotoionizationRate(), cdDrive(), ContCreateMesh(), ContCreatePointers(), ContSetIntensity(), CoolEvaluate(), CoStarInterpolate(), emergent_line(), extin(), ffun1(), GrainsInit(), GrainUpdateRadius2(), GridInterpolate(), HaardtMadauInterpolate(), InitSimPostparse(), InterpolateGridCoStar(), InterpolateModel(), InterpolateModelCoStar(), InterpolateRectGrid(), IterRestart(), IterStart(), KhaireSrianandInterpolate(), Kurucz79Interpolate(), lgReadAtmosphereTail(), lgValidModel(), mie_calc_ial(), mie_write_opc(), MihalasInterpolate(), OpacityAddTotal(), diatomics::OpacityCreate(), OpacityCreate1Element(), OpacityCreateAll(), OpacityZero(), OpacityZeroOld(), t_gaunt::p_gauntff_vec(), t_gaunt::p_setup_brems(), PlanckIntegral(), qheat(), qintr(), radius_next(), RauchInterpolateCOWD(), RauchInterpolateHCa(), RauchInterpolateHelium(), RauchInterpolateHNi(), RauchInterpolateHpHe(), RauchInterpolateHydr(), RauchInterpolatePG1159(), RebinAtmosphere(), resetCoarseTransCoef(), rfield_opac_alloc(), RT_diffuse(), RT_OTS_AddCont(), RT_OTS_PrtRate(), RT_recom_effic(), RT_tau_init(), RT_tau_reset(), save_opacity(), SaveDo(), tfidle(), TlustyInterpolate(), UpdatePot1(), WernerInterpolate(), WMBASICInterpolate(), and ZoneEnd().
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 |
outward emitted continuum
Referenced by iso_ionize_recombine(), IterRestart(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), and RT_line_pumping().
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().
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::otscon |
the local ots continuum rates
Referenced by ContNegative(), GammaPrt(), IterRestart(), IterStart(), OpacityAddTotal(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), RT_diffuse(), RT_OTS(), RT_OTS_AddCont(), RT_OTS_ChkSum(), RT_OTS_PrtRate(), RT_OTS_Update(), RT_OTS_Zero(), and SaveDo().
vector<realnum> t_rfield::otslin |
otsline and otscon - local ots fields for line and continua outlin outward line fields the local ots line rates
Referenced by ContNegative(), eden_sum(), GammaBn(), GammaK(), GammaPrt(), GrainTemperature(), ion_photo(), IterRestart(), IterStart(), OpacityAddTotal(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), RT_diffuse(), RT_OTS(), RT_OTS_AddLine(), RT_OTS_ChkSum(), RT_OTS_PrtRate(), RT_OTS_Update(), RT_OTS_Zero(), SaveDo(), and setXtraRatesCa2().
Referenced by IterRestart(), IterStart(), and rfield_opac_alloc().
Spectrum t_rfield::outlin |
outward directed line emission photons cm-2 s-1
Referenced by ContNegative(), flxCell(), GammaPrt(), GetDissociationRateCoeff(), diatomics::GetHeatRate(), GrainDrift(), diatomics::H2_LevelPops(), highen(), IterEnd(), IterRestart(), IterStart(), lines_continuum(), mole_h_reactions(), outline_base_bin(), outsum(), PresTotCurrent(), PrtFinal(), PrtHeader(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), RT_diffuse(), RT_OTS_ChkSum(), RT_OTS_Update(), RT_OTS_Zero(), SaveDo(), and sum_radiation().
vector<realnum> t_rfield::outlin_noplot |
Referenced by GammaPrt(), GetDissociationRateCoeff(), diatomics::GetHeatRate(), GrainDrift(), diatomics::H2_LevelPops(), highen(), IterEnd(), IterRestart(), IterStart(), mole_h_reactions(), outsum(), PresTotCurrent(), PrtFinal(), PrtHeader(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), RT_diffuse(), RT_OTS_ChkSum(), RT_OTS_Update(), RT_OTS_Zero(), SaveDo(), and sum_radiation().
realnum t_rfield::plsfrq |
plasma frequency for current position in slab
Referenced by ContSetIntensity(), EdenChange(), iso_level(), TransitionProxy::outline(), PressureRadiationLine(), RT_line_one_escape(), RT_line_one_fine(), RT_line_one_tauinc(), and tauff().
realnum t_rfield::plsfrqmax |
store highest energy plasma frequency encountered
Referenced by ContSetIntensity(), EdenChange(), and PrtComment().
realnum t_rfield::qbal |
Referenced by ContSetIntensity(), PrtHeader(), and radius_first().
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().
Spectrum t_rfield::reflin |
reflected line
Referenced by flxCell(), IterEnd(), IterRestart(), lines_continuum(), outline_base_bin(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), RT_diffuse(), SaveDo(), and sum_radiation().
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::slope[LIMSPC] |
Referenced by conorm(), ContSetIntensity(), ffun1(), ParseAgn(), ParseBackgrd(), ParseBlackbody(), ParseBremsstrahlung(), ParseCMB(), ParseLaser(), and ParsePowerlawContinuum().
vector<double> t_rfield::SummedCon |
these are energy-by-energy sums of various arrays, used to save time in evaluating rate integrals
Referenced by ContNegative(), GammaBn(), GammaK(), GammaPrt(), GrainElecEmis1(), GrainTemperature(), highen(), IterRestart(), qheat_init(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), RT_OTS_ChkSum(), RT_OTS_Update(), RT_OTS_Zero(), setTrimming(), ZoneEnd(), and ZoneStart().
vector<realnum> t_rfield::SummedDif |
Referenced by GrainTemperature(), IterRestart(), IterStart(), rfield_opac_alloc(), rfield_opac_zero(), RT_continuum(), RT_OTS_ChkSum(), RT_OTS_Update(), RT_OTS_Zero(), ZoneEnd(), and ZoneStart().
vector<realnum> t_rfield::SummedDifSave |
Referenced by IterRestart(), IterStart(), and rfield_opac_alloc().
vector<realnum> t_rfield::SummedOcc |
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().
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().
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().
double t_rfield::totpow[LIMSPC] |
|
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 |
Referenced by getCoarseTransCoef(), resetCoarseTransCoef(), and RT_tau_inc().
Referenced by conorm(), ffun1(), InitDefaultsPreparse(), InterpolateGridCoStar(), InterpolateRectGrid(), ParseInterp(), and ParseTable().
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().