Go to the source code of this file.
Variables | |
const double | E = 2.718281828459045235360287 |
const double | EULER = 0.577215664901532860606512090082 |
const double | PI = 3.141592653589793238462643 |
const double | PI2 = 6.283185307179586476925287 |
const double | PI4 = 12.56637061435917295385057 |
const double | PI8 = 25.13274122871834590770115 |
const double | SQRT2 = 1.414213562373095048801689 |
const double | SQRTPI = 1.772453850905516027298167 |
const double | SQRTPIBY2 = 1.253314137315500251207883 |
const double | LN_TWO = 0.6931471805599453094172321 |
const double | LN_TEN = 2.302585092994045684017991 |
const double | LOG10_E = 0.4342944819032518276511289 |
const double | OPTDEP2EXTIN = 1.085736204758129569127822 |
const double | RADIAN = 57.29577951308232087679815 |
const double | SOLAR_MASS = 1.9889e33 |
const double | SOLAR_LUMINOSITY = 3.846e33 |
const double | AU = 1.49597870691e13 |
const double | ATOMIC_MASS_UNIT = 1.66053886e-24 |
const double | ELECTRON_MASS = 9.1093826e-28 |
const double | PROTON_MASS = 1.67262171e-24 |
const double | BOLTZMANN = 1.3806505e-16 |
const double | SPEEDLIGHT = 2.99792458e10 |
const double | HPLANCK = 6.6260693e-27 |
const double | ELEM_CHARGE = 1.60217653e-19 |
const double | RYD_INF = 1.0973731568525e5 |
const double | HIONPOT = 0.999466508219 |
const double | PARSEC = AU*RADIAN*3600. |
const double | H_BAR = HPLANCK/(2.*PI) |
const double | ELEM_CHARGE_ESU = ELEM_CHARGE*SPEEDLIGHT/10. |
const double | ELECTRIC_CONST = 1.e11/(PI4*pow2(SPEEDLIGHT)) |
const double | HION_LTE_POP = pow2(HPLANCK)/(PI2*BOLTZMANN*ELECTRON_MASS) |
const double | SAHA = sqrt(pow3(HION_LTE_POP)) |
const double | ERG1CM = HPLANCK*SPEEDLIGHT |
const double | T1CM = HPLANCK*SPEEDLIGHT/BOLTZMANN |
const double | WAVNRYD = 1./RYD_INF |
const double | RYDLAM = 1.e8/RYD_INF |
const double | EN1RYD = HPLANCK*SPEEDLIGHT*RYD_INF |
const double | TE1RYD = HPLANCK*SPEEDLIGHT*RYD_INF/BOLTZMANN |
const double | EVDEGK = ELEM_CHARGE*1.e7/BOLTZMANN |
const double | EVRYD = HPLANCK*SPEEDLIGHT*RYD_INF/ELEM_CHARGE*1.e-7 |
const double | EN1EV = EN1RYD/EVRYD |
const double | FR1RYD = SPEEDLIGHT*RYD_INF |
const double | HNU3C2 = 2.*HPLANCK*SPEEDLIGHT*pow3(RYD_INF) |
const double | FR1RYDHYD = SPEEDLIGHT*RYD_INF*HIONPOT |
const double | HBAReV = H_BAR/EN1EV |
const double | RYDLAMHYD = RYDLAM/HIONPOT |
const double | STEFAN_BOLTZ = pow2(PI*pow2(BOLTZMANN))/(60.*pow3(H_BAR)*pow2(SPEEDLIGHT)) |
const double | FREQ_1EV = SPEEDLIGHT*RYD_INF/EVRYD |
const double | FINE_STRUCTURE = pow2(ELEM_CHARGE_ESU)/SPEEDLIGHT/H_BAR |
const double | FINE_STRUCTURE2 = pow2(FINE_STRUCTURE) |
const double | BOHR_RADIUS_CM = FINE_STRUCTURE/(PI4*RYD_INF) |
const double | TWO_PHOT_CONST = 9.*pow3(FINE_STRUCTURE2)*FR1RYD/2048. |
const double | COLL_CONST = SAHA*BOLTZMANN/HPLANCK |
const double | MILNE_CONST = SPEEDLIGHT*sqrt(pow3(FINE_STRUCTURE2)*pow3(TE1RYD)/PI) |
const double | TRANS_PROB_CONST = PI4*HPLANCK*FINE_STRUCTURE/ELECTRON_MASS |
* physical constants used by Cloudy, taken from * >>refer phys const Mohr P.J. & Taylor B.N., 2002 Codata, see Reviews of Modern Physics, 77, 1 (2005) * >>refercon or http://www.physics.nist.gov/constants * <BR><BR> * NB - these are all printed with the "print constants" command, * which is in parse_print.c, so any new constants * added here must also be added to the prt_constants routine * this is in the func_test case in the auto test suite
Definition in file physconst.h.
const double ATOMIC_MASS_UNIT = 1.66053886e-24 |
atomic mass unit, grams
Definition at line 86 of file physconst.h.
Referenced by collision_strength_VF01(), CS_l_mixing_PS64(), dense_fabden(), GrainCollHeating(), GrainsInit(), HeCollid(), HydroCollid(), lines_helium(), mie_write_opc(), newspecies(), pah1_fun(), PresTotCurrent(), prt_constants(), rfield_opac_malloc(), ritoa(), velset(), and zero().
const double AU = 1.49597870691e13 |
astronomical unit, nearly the length of the semimajor axis of the Earth's elliptical orbit around the sun
Definition at line 78 of file physconst.h.
Referenced by dense_fabden(), and prt_constants().
const double BOHR_RADIUS_CM = FINE_STRUCTURE/(PI4*RYD_INF) |
Bohr radius in cm, 5.29177249e-9
Definition at line 199 of file physconst.h.
Referenced by collision_strength_VF01(), hydro_vs_coll_str(), prt_constants(), and ritoa().
const double BOLTZMANN = 1.3806505e-16 |
this is the Boltzmann factor, ergs/K
Definition at line 95 of file physconst.h.
Referenced by AgeCheck(), atom_levelN(), CoolEvaluate(), Fe11Lev5(), Fe13Lev5(), Fe3Lev14(), Fe4Lev12(), GrainCollHeating(), GrainDrift(), GrainElecRecomb1(), GrainScreen(), GrainTemperature(), grn_abs(), HydroRecCool(), iso_cool(), IterStart(), lines_molecules(), noneq_offset(), PresTotCurrent(), prt_constants(), PrtZone(), PunchDo(), qheat_init(), radius_increment(), rfield_opac_malloc(), uderiv(), UpdatePot(), and velset().
const double COLL_CONST = SAHA*BOLTZMANN/HPLANCK |
this is the square of the value roughly equal to 8.629e-6 that appears in converting collision strengths to rates. The constant is h^2/((2PI*me)^3/2 * k^1/2).
Definition at line 206 of file physconst.h.
Referenced by CoolCarb(), CS_l_mixing_PS64(), FeIICollRatesBoltzmann(), H1cs123(), HeCollid(), HeCreate(), Hion_colldeexc_cs(), Hydcs123(), hydro_vs_deexcit(), HydroCollid(), HyperfineCreate(), lines_lv1_li_ne(), prt_constants(), Punch1LineData(), and tfidle().
const double E = 2.718281828459045235360287 |
the number e
Definition at line 20 of file physconst.h.
Referenced by atmdat_readin(), bhintegrand_log(), CO_create_react(), CoolCarb(), CoolChlo(), CoolIron(), CoolMang(), CoolNeon(), CoolOxyg(), Fe7Lev8(), ForbiddenAuls(), HeCSInterp(), HeLikeLevel(), HeRecom(), hmole_step(), and prt_constants().
const double ELECTRIC_CONST = 1.e11/(PI4*pow2(SPEEDLIGHT)) |
electric constant, in F/m, 8.854e-12
Definition at line 130 of file physconst.h.
Referenced by mie_read_opc(), and prt_constants().
const double ELECTRON_MASS = 9.1093826e-28 |
electron mass, grams
Definition at line 89 of file physconst.h.
Referenced by CO_Init(), ContCreatePointers(), CS_l_mixing_PS64(), GrainCollHeating(), GrainElecRecomb1(), HeCollid(), highen(), hv(), HydroCollid(), lines_helium(), prt_constants(), and ritoa().
const double ELEM_CHARGE = 1.60217653e-19 |
elementary charge, in C in SI units, to use this must convert to cgs
Definition at line 104 of file physconst.h.
Referenced by GrainScreen(), and prt_constants().
const double ELEM_CHARGE_ESU = ELEM_CHARGE*SPEEDLIGHT/10. |
elementary charge, in ESU, 4.8032e-10
Definition at line 127 of file physconst.h.
Referenced by collision_strength_VF01(), and prt_constants().
ergs per eV, 1.602176e-012
Definition at line 169 of file physconst.h.
Referenced by ChargTranSumHeat(), CS_l_mixing_S62(), GrainCollHeating(), H2_Cooling(), highen(), hmole_step(), and prt_constants().
const double EN1RYD = HPLANCK*SPEEDLIGHT*RYD_INF |
ergs per inf mass Ryd, 2.180e-11
Definition at line 156 of file physconst.h.
Referenced by AGN_Hemis(), collision_strength_VF01(), conorm(), ContCreatePointers(), ContSetIntensity(), CoolEvaluate(), GammaBn(), GammaBnPL(), GammaK(), GammaPL(), GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), GrainCollHeating(), GrainDrift(), GrainMakeDiffuse(), GrainTemperature(), HeatSum(), HeCollid(), HeCreate(), HeRecom(), highen(), hmole_reactions(), hv(), HydroCollid(), ion_collis(), iso_cool(), iso_create(), lindst(), lines(), lines_continuum(), lines_helium(), OscStr_f(), OscStr_f_log10(), outsum(), PresTotCurrent(), prt_constants(), PrtFinal(), PrtHeader(), PrtZone(), PunchDo(), PunchSpecial(), qheat(), qheat_init(), radius_increment(), RebinQHeatResults(), RT_line_driving(), RT_radiative_acceleration(), Spec_cont(), sumcon(), TryDoubleStep(), uderiv(), and UpdatePot().
const double ERG1CM = HPLANCK*SPEEDLIGHT |
number of ergs per wavenumber, 1.9864e-16
Definition at line 144 of file physconst.h.
Referenced by atmdat_readin(), ContCreatePointers(), Fe7Lev8(), FeIIPoint(), H2_Cooling(), H2_Create(), lines_helium(), and prt_constants().
const double EULER = 0.577215664901532860606512090082 |
the Euler constant (aka Euler-Mascheroni constant or gamma)
Definition at line 23 of file physconst.h.
Referenced by expn(), and prt_constants().
const double EVDEGK = ELEM_CHARGE*1.e7/BOLTZMANN |
Kelvins per eV, 1.1604e4
Definition at line 163 of file physconst.h.
Referenced by hmole_step(), Hydcs123(), ligbar(), prt_constants(), S62_Therm_ave_coll_str(), and tfidle().
const double EVRYD = HPLANCK*SPEEDLIGHT*RYD_INF/ELEM_CHARGE*1.e-7 |
eV per inf mass Ryd, 13.606
Definition at line 166 of file physconst.h.
Referenced by AnuUnit(), atmdat_dielrec_fe(), t_ADfA::coll_ion(), ContCreatePointers(), cross_section(), CS_l_mixing_S62(), CS_VS80(), GrainCharge(), GrainDrive(), H1cs123(), H2_LevelPops(), Hion_coll_ioniz_ratecoef(), Hion_colldeexc_cs(), Hydcs123(), hydro_vs_coll_str(), hydro_vs_deexcit(), hydro_vs_excit(), hydro_vs_ioniz(), t_yield::init_yield(), ipShells(), iso_create(), lines(), Opacity_iso_photo_cs(), OpacityCreate1Element(), OpacityCreateAll(), OpacityCreateReilMan(), pah1_fun(), PresTotCurrent(), prt_constants(), PrtZone(), punch_opacity(), PunchDo(), punchFITS_EnergyData(), radius_increment(), ReadTable(), S62_Therm_ave_coll_str(), Therm_ave_coll_str_int_VS80(), UpdatePot1(), and Yan_H2_CS().
const double FINE_STRUCTURE = pow2(ELEM_CHARGE_ESU)/SPEEDLIGHT/H_BAR |
the fine-structure constant a= 2pi e^2/hc 7.297 352 533 x 10-3
Definition at line 193 of file physconst.h.
Referenced by prt_constants().
const double FINE_STRUCTURE2 = pow2(FINE_STRUCTURE) |
the square of the fine-structure constant
Definition at line 196 of file physconst.h.
const double FR1RYD = SPEEDLIGHT*RYD_INF |
frequency of one Ryd for infinite mass nuclei, 3.289842e15
Definition at line 172 of file physconst.h.
Referenced by conorm(), ContSetIntensity(), GrainMakeDiffuse(), lgCompileAtmosphere(), lgValidModel(), ParseF_nu(), PlanckIntegral(), plankf(), prt_constants(), and rayleh().
const double FR1RYDHYD = SPEEDLIGHT*RYD_INF*HIONPOT |
frequency of ionization potential of H (not inf mass), 3.288087e15 - never used
Definition at line 178 of file physconst.h.
Referenced by prt_constants().
const double FREQ_1EV = SPEEDLIGHT*RYD_INF/EVRYD |
the frequency of one eV, 2.418e14
Definition at line 190 of file physconst.h.
Referenced by prt_constants().
h/2pi = 1.05457e-27
Definition at line 124 of file physconst.h.
Referenced by collision_strength_VF01(), prt_constants(), and ritoa().
H_BAR in eV sec, 6.582e-16
Definition at line 181 of file physconst.h.
Referenced by prt_constants().
const double HION_LTE_POP = pow2(HPLANCK)/(PI2*BOLTZMANN*ELECTRON_MASS) |
this is the factor that appears in front of Boltzmann factor to get LTE level populations for hydrogenic ions. It is given in the first parts of section 5 of part 2 of hazy, and is actually ( planck^2 / (2 pi m_e k ) )^3/2, but cannot evaluate powers here, must raise this to 3/2 when used, HION_LTE_POP = 5.556e-11
Definition at line 137 of file physconst.h.
Referenced by HeCollid(), HydroCollid(), lines_helium(), and prt_constants().
const double HIONPOT = 0.999466508219 |
ionization potential of real hydrogen atom, in inf mass ryd, based on Codata 98, uncertainty 24e-12, calculated by Peter van Hoof
Definition at line 111 of file physconst.h.
Referenced by ParseBackgrd(), ParseF_nu(), ParseRangeOption(), ParseTable(), prt_constants(), and zero().
const double HNU3C2 = 2.*HPLANCK*SPEEDLIGHT*pow3(RYD_INF) |
2 h FR1RYD^3 / c^2 for infinite mass nucleus, 0.5250
Definition at line 175 of file physconst.h.
Referenced by ContSetIntensity(), and prt_constants().
const double HPLANCK = 6.6260693e-27 |
Planck's constant
Definition at line 101 of file physconst.h.
Referenced by ContSetIntensity(), H_Einstein_A_lin(), H_Einstein_A_log10(), PlanckIntegral(), prt_constants(), and RT_diffuse().
const double LN_TEN = 2.302585092994045684017991 |
const double LN_TWO = 0.6931471805599453094172321 |
ln(2)
Definition at line 47 of file physconst.h.
Referenced by GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), prt_constants(), and qheat().
const double LOG10_E = 0.4342944819032518276511289 |
const double MILNE_CONST = SPEEDLIGHT*sqrt(pow3(FINE_STRUCTURE2)*pow3(TE1RYD)/PI) |
this is the square of the value roughly equal to 4.123e11 that appears in the integration of photoionization cross-sections to obtain recombination coefficients.
Definition at line 210 of file physconst.h.
Referenced by prt_constants(), radrecomb(), and RT_diffuse().
const double OPTDEP2EXTIN = 1.085736204758129569127822 |
factor that converts optical depth into extinction in mags, 2.5 log e
Definition at line 57 of file physconst.h.
Referenced by prt_constants(), and radius_increment().
parsec in cm, 3.085678e18
Definition at line 121 of file physconst.h.
Referenced by dense_fabden(), ParseCommands(), ParseRadius(), ParseStop(), and prt_constants().
const double PI = 3.141592653589793238462643 |
pi
Definition at line 26 of file physconst.h.
Referenced by AngerJ(), bhG(), bhG_mx(), collision_strength_VF01(), CompileAtmosphereCoStar(), conorm(), ContCreatePointers(), DoFSMixing(), FreeFreeGaunt(), GetProbDistr_HighLimit(), GrainCollHeating(), GrainElecRecomb1(), GrainScreen(), grn_abs(), HydroOscilStr(), Jint(), L_mix_integrand_VF01(), mie_auxiliary(), mie_cs(), mie_integrate(), pah1_fun(), ParseCommands(), prt_constants(), PrtFinal(), radius_increment(), RauchInitializeSub(), ritoa(), RT_LineWidth(), S62_Therm_ave_coll_str(), scqdri(), StarkCollTransProb_VF01(), Stognienko(), velset(), vib_evap(), and WaveFunction().
const double PI2 = 6.283185307179586476925287 |
2*pi
Definition at line 29 of file physconst.h.
Referenced by conorm(), ParseMagnet(), and prt_constants().
const double PI4 = 12.56637061435917295385057 |
4*pi
Definition at line 32 of file physconst.h.
Referenced by conorm(), ContCreatePointers(), ContSetIntensity(), FeIIPoint(), GrainMakeDiffuse(), H2_Create(), HeCreate(), iso_create(), Magnetic_evaluate(), mie_calc_ial(), mie_read_opc(), ParseBackgrd(), ParseTable(), PlanckIntegral(), and prt_constants().
const double PI8 = 25.13274122871834590770115 |
8*pi
Definition at line 35 of file physconst.h.
Referenced by CoolEvaluate(), Magnetic_evaluate(), MeanInc(), and prt_constants().
const double PROTON_MASS = 1.67262171e-24 |
proton mass, grams
Definition at line 92 of file physconst.h.
Referenced by collision_strength_VF01(), hv(), and prt_constants().
const double RADIAN = 57.29577951308232087679815 |
const double RYD_INF = 1.0973731568525e5 |
infinite mass rydberg constant, in cm^-1
Definition at line 107 of file physconst.h.
Referenced by AnuUnit(), defect(), he_1trans(), he_energy(), HeCreate(), hydro_vs_coll_str(), HydroEinstA(), lines_helium(), pah1_fun(), prt_constants(), and RT_diffuse().
Angstrom per infinite mass Ryd, 911.2671
Definition at line 153 of file physconst.h.
Referenced by AnuUnit(), AtlasCompile(), CompileAtmosphereCoStar(), ContBandsCreate(), ContCreatePointers(), CoStarCompile(), HeRecom(), HLineTransOpacSet(), lgCompileAtmosphere(), lines(), lines_continuum(), lines_helium(), ParseTable(), PntForLine(), prt_constants(), PrtComment(), PrtFinal(), PunchNewContinuum(), PunchSpecial(), read_hm05(), and ReadTable().
wavelength (A) of ionization potential of Hydrogen, 911.7535 - never used
Definition at line 184 of file physconst.h.
Referenced by prt_constants().
const double SAHA = sqrt(pow3(HION_LTE_POP)) |
SAHA is ( h^2/2/pi/m/k )^3/2, is correct constant for free electron SAHA = 4.14132e-16
Definition at line 141 of file physconst.h.
Referenced by H2_LevelPops(), H2_X_coll_rate_evaluate(), hmole_reactions(), and prt_constants().
const double SOLAR_LUMINOSITY = 3.846e33 |
solar luminosity erg s-1 >>refer phys const http://pdg.lbl.gov/2000/astrorpp.pdf
Definition at line 72 of file physconst.h.
Referenced by prt_constants(), and PrtHeader().
const double SOLAR_MASS = 1.9889e33 |
solar mass in grams >>refer phys const http://pdg.lbl.gov/2000/astrorpp.pdf
Definition at line 68 of file physconst.h.
Referenced by AbundancesSet(), lgConvPres(), prt_constants(), PrtComment(), and PrtFinal().
const double SPEEDLIGHT = 2.99792458e10 |
speed of light, cm/s
Definition at line 98 of file physconst.h.
Referenced by AnuUnit(), conorm(), ContCreatePointers(), ContSetIntensity(), GrainDrift(), GrainMakeDiffuse(), highen(), lgCompileAtmosphere(), ParseSet(), PlanckIntegral(), prt_constants(), PrtZone(), radius_increment(), rfield_opac_malloc(), ritoa(), RT_radiative_acceleration(), and zero().
const double SQRT2 = 1.414213562373095048801689 |
const double SQRTPI = 1.772453850905516027298167 |
sqrt(pi)
Definition at line 41 of file physconst.h.
Referenced by esc_CRDwing_1side(), esca0k2(), prt_constants(), RT_DestProb(), RTesc_lya(), RTesc_lya_1side(), and vfun().
const double SQRTPIBY2 = 1.253314137315500251207883 |
sqrt(pi/2)
Definition at line 44 of file physconst.h.
Referenced by bhG(), bhG_mx(), and prt_constants().
const double STEFAN_BOLTZ = pow2(PI*pow2(BOLTZMANN))/(60.*pow3(H_BAR)*pow2(SPEEDLIGHT)) |
Stefan-Boltzmann constant, 5.6704e-5
Definition at line 187 of file physconst.h.
Referenced by lgValidModel(), and prt_constants().
const double T1CM = HPLANCK*SPEEDLIGHT/BOLTZMANN |
degrees kelvin per unit wavenumber, 1.4388
Definition at line 147 of file physconst.h.
Referenced by atom_levelN(), atom_pop5(), ContCreatePointers(), CoolCarb(), Fe11Lev5(), Fe13Lev5(), Fe3Lev14(), Fe4Lev12(), FeIIPoint(), H21_cm_pops(), H2_Create(), H2_Level_low_matrix(), H2_PunchDo(), hmole_reactions(), lines_helium(), prt_constants(), and tfidle().
const double TE1RYD = HPLANCK*SPEEDLIGHT*RYD_INF/BOLTZMANN |
the temperature of 1 Rydberg te1ryd is h/k is temp of 1 Rydberg, 1.579e5
Definition at line 160 of file physconst.h.
Referenced by atmdat_dielrec_fe(), t_ADfA::coll_ion(), cont_gaunt_calc(), ContSetIntensity(), CS_l_mixing_VF01(), CS_VS80(), dgaunt(), ffun1(), FillGFF(), FreeFreeGaunt(), GrainChargeTemp(), GrainDrift(), GrainMakeDiffuse(), HeCreate(), Hion_coll_ioniz_ratecoef(), Hion_colldeexc_cs(), Hydcs123(), hydro_vs_deexcit(), hydro_vs_excit(), hydro_vs_ioniz(), iso_cool(), iso_create(), OpacityAddTotal(), ParseAgn(), ParseBlackbody(), ParsePowerlawContinuum(), PlanckIntegral(), prt_constants(), PrtHeader(), PunchGaunts(), radrecomb(), Recomb_Seaton59(), RT_diffuse(), RT_recom_effic(), SanityCheckBegin(), tfidle(), Therm_ave_coll_str_int_VS80(), and UpdatePot2().
const double TRANS_PROB_CONST = PI4*HPLANCK*FINE_STRUCTURE/ELECTRON_MASS |
This is the constant used in converting oscillator strengths to As. The formula is Aul = TRANS_PROB_CONST * f(u,l) * wavenumber^2. TRANS_PROB_CONST is 0.667025
Definition at line 214 of file physconst.h.
Referenced by CS_l_mixing_S62(), DoFSMixing(), eina(), FeIICreate(), GetGF(), he_1trans(), HydroEinstA(), lines_setup(), prt_constants(), and ritoa().
const double TWO_PHOT_CONST = 9.*pow3(FINE_STRUCTURE2)*FR1RYD/2048. |
the two photon constant as defined by Breit & Teller, as in equation 4 of Spitzer & Greenstein 51, 2.18313
Definition at line 202 of file physconst.h.
Referenced by prt_constants().
number of Ryd per wavenumber, 9.11267e-6
Definition at line 150 of file physconst.h.
Referenced by ColStrGBar(), ContCreatePointers(), CS_l_mixing_S62(), FeIIPoint(), H2_ContPoint(), H2_Create(), H2_LevelPops(), HeCreate(), HeLikeLevel(), HLineTransOpacSet(), HydroLevelPop(), iso_create(), lines_helium(), MakeCS(), mie_calc_ial(), mie_write_opc(), ParsePunch(), prt_constants(), pun1Line(), and tbl_fun().