t_dense Struct Reference

#include <dense.h>

Collaboration diagram for t_dense:

Public Member Functions
	t_dense ()
Data Fields
realnum	gas_phase [LIMELM]
realnum	AtomicWeight [LIMELM]
realnum	xMolecules [LIMELM]
realnum	H_sum_in_CO
realnum	xMassDensity
realnum	xMassDensity0
realnum	pden
realnum	wmole
realnum	xNucleiTotal
realnum	xMassTotal
realnum	HCorrFac
long int	IonLow [LIMELM+1]
long int	IonHigh [LIMELM+1]
realnum	xIonDense [LIMELM+3][LIMELM+1]
bool	lgIonChiantiOn [LIMELM][LIMELM+1]
realnum	AbundanceLimit
bool	lgElmtOn [LIMELM]
bool	lgSetIoniz [LIMELM]
realnum	SetIoniz [LIMELM][LIMELM+1]
char	chDenseLaw [5]
bool	lgDenseInitConstant
bool	lgPressureVaryTime
double	PressureVaryTimeTimescale
double	PressureVaryTimeIndex
double	DensityLaw [10]
bool	lgAsChoose [LIMELM][LIMELM]
bool	lgCSChoose [LIMELM][LIMELM]
realnum	frad [LIMTABDLAW]
realnum	fhden [LIMTABDLAW]
long int	nvals
bool	lgDLWDepth
double	eden
double	EdenMax
double	EdenMin
realnum	eden_f
double	density_low_limit
long int	nzEdenBad
realnum	EdenSet
realnum	EdenExtra
double	SqrtEden
double	EdenHCorr
double	EdenTrue
double	eden_from_metals
bool	lgEdenBad
double	edensqte
double	cdsqte
realnum	DensityPower
realnum	rscale
realnum	den0
bool	lgDenFlucOn
bool	lgDenFlucRadius
realnum	flong
realnum	cfirst
realnum	csecnd
realnum	flcPhase

Detailed Description

variables dealing with pressure across model

Definition at line 29 of file dense.h.

Constructor & Destructor Documentation

t_dense::t_dense ( ) [inline]

Definition at line 31 of file dense.h.

References AtomicWeight.

Field Documentation

realnum t_dense::AbundanceLimit

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().

realnum t_dense::AtomicWeight[LIMELM]

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().

realnum t_dense::cfirst

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().

realnum t_dense::csecnd

Definition at line 242 of file dense.h.

Referenced by AbundancesSet(), lgConvPres(), ParseFluc(), and PressureChange().

realnum t_dense::den0

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().

realnum t_dense::DensityPower

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().

realnum t_dense::eden_f

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().

realnum t_dense::EdenExtra

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().

realnum t_dense::EdenSet

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

edensqte is eden/sqrte

Definition at line 219 of file dense.h.

Referenced by tfidle().

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().

realnum t_dense::fhden[LIMTABDLAW]

Definition at line 172 of file dense.h.

Referenced by dense_tabden(), ParseDLaw(), and ParseTLaw().

realnum t_dense::flcPhase

Definition at line 243 of file dense.h.

Referenced by AbundancesSet(), InitDefaultsPreparse(), lgConvPres(), ParseFluc(), and PressureChange().

realnum t_dense::flong

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().

realnum t_dense::frad[LIMTABDLAW]

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().

realnum t_dense::gas_phase[LIMELM]

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().

realnum t_dense::H_sum_in_CO

sum of all H in CO network

Definition at line 79 of file dense.h.

Referenced by CO_solve(), DynaIonize(), IterRestart(), IterStart(), and zero().

realnum t_dense::HCorrFac

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]

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]

Definition at line 168 of file dense.h.

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().

bool t_dense::lgDenFlucRadius

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().

bool t_dense::lgDenseInitConstant

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().

bool t_dense::lgPressureVaryTime

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().

realnum t_dense::pden

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().

double t_dense::PressureVaryTimeIndex

Definition at line 161 of file dense.h.

Referenced by ConvInitSolution(), InitDefaultsPreparse(), and ParseConstant().

double t_dense::PressureVaryTimeTimescale

Definition at line 159 of file dense.h.

Referenced by ConvInitSolution(), InitDefaultsPreparse(), and ParseConstant().

realnum t_dense::rscale

Definition at line 229 of file dense.h.

Referenced by lgConvPres(), and ParseHDEN().

realnum t_dense::SetIoniz[LIMELM][LIMELM+1]

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().

realnum t_dense::wmole

mean AMU per particle

Definition at line 92 of file dense.h.

Referenced by dense_parametric_wind(), PresTotCurrent(), PrtFinal(), radius_first(), and radius_increment().

realnum t_dense::xIonDense[LIMELM+3][LIMELM+1]

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().