#include <dense.h>
variables dealing with pressure across model
Definition at line 29 of file dense.h.
t_dense::t_dense | ( | ) | [inline] |
Definition at line 31 of file dense.h.
References AtomicWeight.
this is lower limit to abundance of element that will be include in the calculation, default is zero, set with command ELEMENT LIMIT OFF XXX command
Definition at line 130 of file dense.h.
Referenced by AbundancesSet(), InitDefaultsPreparse(), and ParseElement().
vector of atomic weights for all elements, set in zerologic
Definition at line 72 of file dense.h.
Referenced by AbundancesSet(), atom_oi_calc(), t_fe2ovr_la::atoms_fe2ovr(), collision_strength_VF01(), CoolCalc(), CoolDima(), CoolEvaluate(), CS_l_mixing_PS64(), FeII_RT_Make(), FeII_RT_TauInc(), FeIILyaPump(), FeIIPunchLineStuff(), FeIIRadPress(), GrainChrgTransferRates(), GrainCollHeating(), GrainElecEmis1(), GrainElecRecomb1(), H2_PunchLineStuff(), H2_RadPress(), H2_RT_tau_inc(), H2_RTMake(), hmole_reactions(), hmole_step(), hydro_vs_coll_str(), iso_collide(), iso_satellite_update(), iso_state_lifetime(), mie_read_form(), mie_read_opc(), newspecies(), pah1_fun(), pah2_fun(), PresTotCurrent(), PrtLinePres(), radius_increment(), rfield_opac_malloc(), ritoa(), RT_line_all(), RT_tau_inc(), RT_tau_init(), S62_Therm_ave_coll_str(), SaveLineStuff(), states_nelemfill(), and t_dense().
double t_dense::cdsqte |
cdsqte is eden/sqrte times 8.629e-6 - this multiplies the collision strength to produce the deexcitation rate coefficient, s-1 8.629e-6 is COLL_CONST
Definition at line 224 of file dense.h.
Referenced by atom_level2(), atom_level3(), atom_levelN(), atom_pop2(), atom_pop3(), atom_pop5(), AtomSeqBeryllium(), ConvRate2CS(), CoolCalc(), CoolCarb(), CoolEvaluate(), CoolNitr(), CoolOxyg(), CoolSili(), emit_frac(), H21_cm_pops(), LineConvRate2CS(), lines_lv1_k_zn(), lines_lv1_li_ne(), oi_level_pops(), and tfidle().
Definition at line 241 of file dense.h.
Referenced by AbundancesSet(), lgConvPres(), ParseFluc(), and PressureChange().
char t_dense::chDenseLaw[5] |
label describing the density law for current calculation 'DLW2' is dense_tabden interpolated table
Definition at line 149 of file dense.h.
Referenced by AgeCheck(), ConvFail(), ConvInitSolution(), hmole(), InitDefaultsPreparse(), iter_end_check(), lgConvPres(), ParseCommands(), ParseConstant(), ParseDLaw(), ParseDynaWind(), ParseFluc(), ParseGlobule(), ParseHDEN(), ParseTLaw(), PresTotCurrent(), prt_smooth_predictions(), PrtComment(), PrtFinal(), radius_first(), radius_next(), and ZoneStart().
Definition at line 242 of file dense.h.
Referenced by AbundancesSet(), lgConvPres(), ParseFluc(), and PressureChange().
Definition at line 230 of file dense.h.
Referenced by hmole(), lgConvPres(), and ParseHDEN().
double t_dense::density_low_limit |
lowest allowed density for any ion = if density falls below this then set to zero in ion_trim
Definition at line 191 of file dense.h.
Referenced by ContSetIntensity(), InitCoreload(), ion_trim(), iso_collide(), and RT_tau_inc().
double t_dense::DensityLaw[10] |
parameters set by the dlaw command, used by dense_fabden (maybe)
Definition at line 164 of file dense.h.
Referenced by dense_fabden(), dense_parametric_wind(), ParseDLaw(), and ParseTLaw().
parameters dealing with hydrogen density scaling as power of radius DensityPower is power
Definition at line 228 of file dense.h.
Referenced by lgConvPres(), ParseCommands(), and ParseHDEN().
double t_dense::eden |
electron density, units cm-3
Definition at line 181 of file dense.h.
Referenced by AGN_Hemis(), atmdat_DielSupres(), atom_level3(), cdEDEN_last(), CO_update_species_cache(), collision_strength_VF01(), ContNegative(), ContSetIntensity(), ConvBase(), ConvEdenIoniz(), ConvFail(), ConvInitSolution(), ConvPresTempEdenIoniz(), CoolCarb(), CoolChlo(), CoolEvaluate(), CoolIron(), CoolNick(), CoolNitr(), CoolSili(), CoolSulf(), CS_l_mixing_PS64(), CS_l_mixing_VF01(), CS_ThermAve_PR78(), da(), dBase_solve(), DoSatelliteLines(), DumpLine(), DynaPunchTimeDep(), eden_sum(), EdenError(), Fe2_cooling(), FeIILevelPops(), find_solution(), fndstr(), GetLineRec(), GrainChargeTemp(), GrainCollHeating(), GrainDrift(), GrainDrive(), GrainElecRecomb1(), GrainTemperature(), H2_LevelPops(), H2_X_coll_rate_evaluate(), highen(), hmole(), hmole_reactions(), hmole_step(), HydroLevel(), HyperfineCreate(), ion_recomb(), ion_recombAGN(), IonCalci(), IonPotas(), iso_collapsed_bnl_set(), iso_continuum_lower(), iso_cool(), iso_departure_coefficients(), iso_ionize_recombine(), iso_satellite_update(), iter_end_check(), IterRestart(), IterStart(), lgCheckMonitors(), lgConvEden(), lines(), lines_continuum(), lines_helium(), lines_hydro(), lines_lv1_k_zn(), lines_lv1_li_ne(), lines_lv1_na_ar(), map_do(), t_mean::MeanInc(), OpacityAddTotal(), PressureChange(), PresTotCurrent(), PrintRates(), PrtHeader(), PrtHydroTrace1(), PrtHydroTrace1a(), PrtZone(), radius_increment(), radius_next(), RT_diffuse(), RT_line_electron_scatter(), RT_LineWidth(), RT_OTS(), RT_stark(), RT_tau_inc(), RT_tau_init(), Save_Line_RT(), SaveDo(), SaveLineData(), SaveSpecial(), tauff(), tfidle(), zero(), and ZoneStart().
electron density, units cm-3, floating point version for compatibility in CO_mole
Definition at line 187 of file dense.h.
Referenced by CO_Init(), and CO_update_species_cache().
double t_dense::eden_from_metals |
fraction of electron density due to ions rather than molecules and grains
Definition at line 213 of file dense.h.
Referenced by eden_sum(), and radius_next().
extra electron density, set with eden command
Definition at line 200 of file dense.h.
Referenced by ContSetIntensity(), eden_sum(), ParseEden(), and zero().
double t_dense::EdenHCorr |
EdenHCorr is eden + hi * 1.7e-4, includes correction for H atom collisions, units cm-3
Definition at line 207 of file dense.h.
Referenced by ContSetIntensity(), CoolCarb(), HydroLevel(), ion_collis(), iso_collide(), iso_cool(), iso_ionize_recombine(), iso_level(), IterRestart(), lines_helium(), oi_othercs(), PrintRates(), PrtHydroTrace1a(), SaveDo(), SaveLineData(), and tfidle().
double t_dense::EdenMax |
max and min eden over this iteration
Definition at line 184 of file dense.h.
Referenced by IterStart(), PrtComment(), and radius_increment().
double t_dense::EdenMin |
Definition at line 184 of file dense.h.
Referenced by IterStart(), and radius_increment().
EdenSet electron density set with set eden command
Definition at line 197 of file dense.h.
Referenced by ContSetIntensity(), eden_sum(), ParseSet(), and zero().
double t_dense::edensqte |
double t_dense::EdenTrue |
this is the true eden as set in eden_sum, we will try to converge eden to this
Definition at line 210 of file dense.h.
Referenced by ContSetIntensity(), ConvBase(), ConvEdenIoniz(), ConvFail(), CoolEvaluate(), eden_sum(), EdenError(), iter_end_check(), IterRestart(), lgConvEden(), lines_grains(), radius_increment(), RT_OTS(), and SaveDo().
Definition at line 172 of file dense.h.
Referenced by dense_tabden(), ParseDLaw(), and ParseTLaw().
Definition at line 243 of file dense.h.
Referenced by AbundancesSet(), InitDefaultsPreparse(), lgConvPres(), ParseFluc(), and PressureChange().
parameters for the density fluctuations command
Definition at line 240 of file dense.h.
Referenced by InitDefaultsPreparse(), lgConvPres(), ParseFluc(), PressureChange(), radius_first(), and radius_next().
frad is log radius in cm, fhden is log hden
Definition at line 171 of file dense.h.
Referenced by dense_tabden(), ParseDLaw(), and ParseTLaw().
dense.gas_phase is the total gas phase abundances, including anything within molecules, but not including grains
Definition at line 69 of file dense.h.
Referenced by AbundancesPrt(), AbundancesSet(), ChemImportance(), CO_drive(), CO_solve(), conorm(), ContRate(), ContSetIntensity(), ConvBase(), ConvFail(), ConvInitSolution(), ConvPresTempEdenIoniz(), ConvTempEdenIoniz(), CoolEvaluate(), DynaIonize(), DynaIterEnd(), DynaPresChngFactor(), DynaPrtZone(), DynaPunchTimeDep(), DynaSave(), DynaSaveLast(), DynaStartZone(), Fe2_cooling(), fill_array(), get_total_abundance_ions(), GrainChargeTemp(), GrainDrive(), GrainsInit(), GrainTemperature(), GrainUpdateRadius1(), GrnStdDpth(), GrnVryDpth(), H2_LevelPops(), H2_PunchDo(), t_dynamics::Heat(), highen(), hmole(), hmole_step(), ion_trim(), ion_wrapper(), IonAlumi(), IonBeryl(), IonBoron(), IonCarbo(), IonHelium(), IonHydro(), IonIron(), IonMagne(), IonNeon(), IonNitro(), IonOxyge(), IonSodiu(), iso_level(), iso_solve(), iter_end_check(), IterRestart(), IterStart(), lgConvPres(), lgHomogeneousSource(), lgIonizConverg(), lgMolecAver(), lines_general(), map_do(), t_mean::MeanInc(), OpacityAddTotal(), ParseCommands(), ParseDLaw(), ParseFluc(), ParseGlobule(), ParseHDEN(), pltopc(), PressureChange(), PresTotCurrent(), PrintRates(), PrtComment(), PrtFinal(), PrtZone(), radius_first(), radius_increment(), radius_next(), RT_continuum(), RT_OTS_PrtRate(), RT_tau_init(), save_opacity(), SaveDo(), SaveLineData(), SaveSpecial(), and zero().
sum of all H in CO network
Definition at line 79 of file dense.h.
Referenced by CO_solve(), DynaIonize(), IterRestart(), IterStart(), and zero().
this is scale factor that multiplies the correction factor for neutral hydrogen collisions, def 1, changed with set command
Definition at line 102 of file dense.h.
Referenced by CoolCarb(), ParseSet(), tfidle(), and zero().
long int t_dense::IonHigh[LIMELM+1] |
Definition at line 111 of file dense.h.
Referenced by ChargTranPun(), CO_solve(), ContSetIntensity(), ConvBase(), CoolEvaluate(), DynaIonize(), fill_array(), find_solution(), get_total_abundance_ions(), GrainChargeTemp(), ion_collis(), ion_photo(), ion_recomb(), ion_recombAGN(), ion_solver(), ion_trim(), IonIron(), iso_level(), iso_solve(), iso_update_rates(), IterRestart(), IterStart(), lgTrivialSolution(), lines_general(), lines_hydro(), map_do(), ParseElement(), PresTotCurrent(), PrintRates(), PrtLinePres(), PrtZone(), t_ionbal::RateIonizTot(), RT_diffuse(), RT_line_all(), RT_line_driving(), RT_OTS(), RT_tau_inc(), store_new_densities(), and tfidle().
long int t_dense::IonLow[LIMELM+1] |
indices for lowest stage of ionization of the elements on C scale, lowest is 0 for atom, -1 if element turned off, the first stage of ionization with positive abundance is [IonLow] where 0 is the atom, the highest stage of ionization with positive abundance is [IonHigh], NB NB so loops should be ion=IonLow, ion<=IonHigh
Definition at line 110 of file dense.h.
Referenced by ChargTranPun(), CO_drive(), CO_solve(), ContSetIntensity(), ConvBase(), CoolEvaluate(), DynaIonize(), fill_array(), find_solution(), get_total_abundance_ions(), GrainChargeTemp(), ion_collis(), ion_photo(), ion_recomb(), ion_recombAGN(), ion_solver(), ion_trim(), IonIron(), IonMagne(), IonOxyge(), iso_level(), iso_solve(), iso_update_rates(), IterRestart(), IterStart(), lines_general(), map_do(), PresTotCurrent(), RT_diffuse(), and store_new_densities().
bool t_dense::lgAsChoose[LIMELM][LIMELM] |
options on set atomic data command
Definition at line 167 of file dense.h.
Referenced by InitDefaultsPreparse().
bool t_dense::lgCSChoose[LIMELM][LIMELM] |
bool t_dense::lgDenFlucOn |
set true when density fluctuations are turned on
Definition at line 233 of file dense.h.
Referenced by AbundancesSet(), InitDefaultsPreparse(), ParseFluc(), PressureChange(), PrtZone(), and radius_next().
set false when fluctuations are over col den rather than radius, set with column options on fluctuations command
Definition at line 237 of file dense.h.
Referenced by InitDefaultsPreparse(), lgConvPres(), ParseFluc(), and PressureChange().
Definition at line 153 of file dense.h.
Referenced by ConvInitSolution(), InitDefaultsPreparse(), and ParseConstant().
bool t_dense::lgDLWDepth |
lg is true if depth, false if radius to be used
Definition at line 178 of file dense.h.
Referenced by dense_tabden(), ParseDLaw(), and ParseTLaw().
bool t_dense::lgEdenBad |
flags set when bad electron density is detected
Definition at line 216 of file dense.h.
Referenced by ConvInitSolution(), and PrtComment().
bool t_dense::lgElmtOn[LIMELM] |
array of logical variables saying whether an element is enable (true) or disabled (false). It is set totally true in zero and is set false with the "element off" command. In SetAbundances if can be reset so that an element that was disabled on the first model in a core load is not later enabled
Definition at line 137 of file dense.h.
Referenced by AbundancesPrt(), AbundancesSet(), atmdat_readin(), Badnell_rec_init(), ChemImportance(), CO_drive(), CO_solve(), ContCreatePointers(), ContRate(), ContSetIntensity(), ConvBase(), ConvIterCheck(), CoolEvaluate(), DoSatelliteLines(), DynaCreateArrays(), DynaIonize(), DynaNewStep(), DynaSaveLast(), DynaStartZone(), GetStandardHeLines(), GrainChargeTemp(), GrainChrgTransferRates(), GrainCollHeating(), GrainElecEmis1(), GrainElecRecomb1(), GrainRateDr(), HeCollidSetup(), HighestIonStage(), InitCoreload(), InitCoreloadPostparse(), InitSimPostparse(), ion_wrapper(), IonAlumi(), IonArgon(), IonBeryl(), IonBoron(), IonCalci(), IonCarbo(), IonChlor(), IonChrom(), IonCobal(), IonCoppe(), IonFluor(), IonHelium(), IonIron(), IonLithi(), IonMagne(), IonManga(), IonNeon(), IonNicke(), IonNitro(), IonOxyge(), IonPhosi(), IonPotas(), IonScand(), IonSilic(), IonSodiu(), IonSulph(), IonTitan(), IonVanad(), IonZinc(), iso_allocate(), iso_assign_quantum_numbers(), iso_cool(), iso_create(), iso_recomb_malloc(), iso_recomb_setup(), iso_satellite(), iso_satellite_update(), iso_solve(), iso_update_rates(), iso_zero(), iter_end_check(), IterRestart(), IterStart(), lgCheckMonitors(), lgIonizConverg(), lines_continuum(), lines_helium(), lines_hydro(), lines_lv1_k_zn(), map_do(), OpacityAddTotal(), OpacityCreateAll(), OpacityValenceRescale(), ParseCompile(), ParseElement(), PresTotCurrent(), prt_smooth_predictions(), PrtAllTau(), PrtColumns(), PrtComment(), PrtFinal(), PrtMeanIon(), prtmet(), PrtZone(), radius_increment(), radius_next(), RT_diffuse(), RT_line_all(), RT_OTS(), RT_stark(), RT_tau_init(), RT_tau_reset(), SanityCheckBegin(), SaveDo(), SaveLineData(), SaveLineStuff(), state_get_put(), and UpdateRecomZ0().
bool t_dense::lgIonChiantiOn[LIMELM][LIMELM+1] |
which ions have chianti enabled?
Definition at line 125 of file dense.h.
Referenced by CoolEvaluate(), InitCoreload(), and states_nelemfill().
Definition at line 156 of file dense.h.
Referenced by ConvInitSolution(), InitDefaultsPreparse(), and ParseConstant().
bool t_dense::lgSetIoniz[LIMELM] |
will we solve for ionization (false) or specify it with element ionization cmnd true
Definition at line 140 of file dense.h.
Referenced by ContSetIntensity(), ConvBase(), InitCoreload(), ion_trim(), IonHelium(), IonHydro(), iso_level(), iso_solve(), lgTrivialSolution(), ParseElement(), and PresTotCurrent().
long int t_dense::nvals |
number of values in above table
Definition at line 175 of file dense.h.
Referenced by dense_tabden(), ParseDLaw(), and ParseTLaw().
long int t_dense::nzEdenBad |
zone where bad electron density was detected
Definition at line 194 of file dense.h.
Referenced by ConvInitSolution(), PrtComment(), and radius_increment().
total number of particles per cubic centimeter
Definition at line 89 of file dense.h.
Referenced by ConvInitSolution(), IterStart(), PresTotCurrent(), radius_first(), radius_increment(), and SaveDo().
Definition at line 161 of file dense.h.
Referenced by ConvInitSolution(), InitDefaultsPreparse(), and ParseConstant().
Definition at line 159 of file dense.h.
Referenced by ConvInitSolution(), InitDefaultsPreparse(), and ParseConstant().
Definition at line 229 of file dense.h.
Referenced by lgConvPres(), and ParseHDEN().
dense.SetIoniz the ionization fractions that are set when lgSetIoniz set true, gas phase abundance is this times total aabundance Ionization fraction for [nelem][ion]
Definition at line 145 of file dense.h.
Referenced by ContSetIntensity(), ConvBase(), IonHelium(), IonHydro(), iso_solve(), lgTrivialSolution(), and ParseElement().
double t_dense::SqrtEden |
square root of electron density, set in tfidle
Definition at line 203 of file dense.h.
Referenced by GrainDrift(), and tfidle().
mean AMU per particle
Definition at line 92 of file dense.h.
Referenced by dense_parametric_wind(), PresTotCurrent(), PrtFinal(), radius_first(), and radius_increment().
dense.xIonDense[nelem][i] is density of ith ionization stage (cm^-3), [nelem][0] is atom, [][1]) the first ion nelem = 0 for H, 1 for he, etc
with this scheme, OIII is dense.xIonDense[ipOXYGEN][2] HII is dense.xIonDense[ipHYDROGEN][1] at the moment 12CO, 13CO, and H2 are xIonDense[LIMELM][2], [L+1][3] and [L+2][1]
Definition at line 122 of file dense.h.
Referenced by AGN_Hemis(), atom_level2(), atom_level3(), atom_oi_calc(), t_fe2ovr_la::atoms_fe2ovr(), AtomSeqBeryllium(), AtomSeqBoron(), ChargTranEval(), ChargTranSumHeat(), CO_drive(), CO_Init(), CO_solve(), ContSetIntensity(), ConvBase(), ConvInitSolution(), CoolAlum(), CoolArgo(), CoolCalc(), CoolCarb(), CoolChlo(), CoolChro(), CoolDima(), CoolEvaluate(), CoolIron(), CoolMagn(), CoolNeon(), CoolNick(), CoolNitr(), CoolOxyg(), CoolPhos(), CoolScan(), CoolSili(), CoolSodi(), CoolSulf(), CoolTita(), CoolVana(), dBase_solve(), DoSatelliteLines(), DynaIonize(), eden_sum(), Fe11Lev5(), Fe13Lev5(), Fe2_cooling(), Fe3Lev14(), Fe4Lev12(), Fe7Lev8(), FeII_RT_Out(), FeIILevelPops(), FeIIPunPop(), FeIISumBand(), fill_array(), get_total_abundance_ions(), GetLineRec(), GrainChargeTemp(), GrainCollHeating(), GrainDrift(), GrainElecEmis1(), GrainElecRecomb1(), GrainMakeDiffuse(), GrnStdDpth(), GrnVryDpth(), H2_Cooling(), H2_PunchDo(), H2_X_coll_rate_evaluate(), HeatSum(), highen(), HighestIonStage(), hmole(), hmole_step(), HydroLevel(), HydroRenorm(), HyperfineCreate(), ion_recomb(), ion_solver(), ion_trim(), ion_wrapper(), IonAlumi(), IonBeryl(), IonBoron(), IonCarbo(), IonHelium(), IonHydro(), IonIron(), IonMagne(), IonNeon(), IonNitro(), IonOxyge(), IonSodiu(), iso_charge_transfer_update(), iso_collide(), iso_continuum_lower(), iso_cool(), iso_departure_coefficients(), iso_ionize_recombine(), iso_level(), iso_satellite_update(), iso_solve(), iter_end_check(), IterRestart(), IterStart(), lgIonizConverg(), lgMolecAver(), lgOH_ChargeTransferDominant(), lgTrivialSolution(), ligbar(), lines(), lines_continuum(), lines_helium(), lines_hydro(), lines_lv1_k_zn(), lines_lv1_li_ne(), lines_lv1_na_ar(), lines_molecules(), MakeCS(), map_do(), t_mean::MeanInc(), mole_H2_form(), oi_cs(), oi_othercs(), OpacityAdd1Element(), OpacityAddTotal(), PressureChange(), PresTotCurrent(), PrintRates(), PrtZone(), radius_first(), radius_increment(), radius_next(), RT_continuum(), RT_diffuse(), RT_line_all(), RT_OTS(), RT_OTS_PrtRate(), RT_tau_inc(), SanityCheckBegin(), save_opacity(), SaveDo(), SaveHeat(), SaveLineData(), store_new_densities(), tfidle(), UpdateRecomZ0(), and zero().
xMassDensity grams per cc
Definition at line 82 of file dense.h.
Referenced by ConvInitSolution(), DynaIonize(), DynaPresChngFactor(), DynaPrtZone(), GravitationalPressure(), Magnetic_evaluate(), PressureChange(), PresTotCurrent(), PrtFinal(), radius_first(), radius_increment(), radius_next(), SaveDo(), and SaveSpecial().
WJH: fiducial value that corresponds to hden set in init file, this is used for setting the mass-flux in dynamic models
Definition at line 86 of file dense.h.
Referenced by DynaFlux(), PresTotCurrent(), and zero().
total mass in grams PER 4 pi rinner^2
Definition at line 98 of file dense.h.
Referenced by GravitationalPressure(), iter_end_check(), IterStart(), PrtFinal(), radius_increment(), and zero().
dense.xMolecules density of elements locked in molecules, this is included in gas_phase
Definition at line 76 of file dense.h.
Referenced by CO_solve(), ConvBase(), fill_array(), get_total_abundance_ions(), hmole(), hmole_step(), IterRestart(), IterStart(), lgHomogeneousSource(), PressureChange(), PrintRates(), and zero().
total number of nuclei, set in PressureTotal
Definition at line 95 of file dense.h.
Referenced by ConvBase(), dBase_solve(), iso_continuum_lower(), PresTotCurrent(), RT_tau_init(), and zero().