physconst.h

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2013 by Gary J. Ferland and
 * others.  For conditions of distribution and use see copyright notice in license.txt */
 
#ifndef PHYSCONST_H_
#define PHYSCONST_H_
 
/*#include "physconst.h"*/
 
/*********************************************************************
 * first come math constants                                         *
 *********************************************************************/
 
UNUSED const double EE = 2.718281828459045235360287;
 
UNUSED const double EULER = 0.577215664901532860606512090082;
 
UNUSED const double PI = 3.141592653589793238462643;
 
UNUSED const double PI2 = 6.283185307179586476925287;
 
UNUSED const double PI4 = 12.56637061435917295385057;
 
UNUSED const double PI8 = 25.13274122871834590770115;
 
UNUSED const double SQRT2 = 1.414213562373095048801689;
 
UNUSED const double SQRTPI = 1.772453850905516027298167;
 
UNUSED const double SQRTPIBY2 = 1.253314137315500251207883;
 
UNUSED const double LN_TWO = 0.6931471805599453094172321;
 
UNUSED const double LN_TEN = 2.302585092994045684017991;
 
UNUSED const double LOG10_E = 0.4342944819032518276511289;
 
UNUSED const double OPTDEP2EXTIN = 1.085736204758129569127822;
 
UNUSED const double RADIAN = 57.29577951308232087679815;
 
/*********************************************************************
 * astronomical constants go here                                    *
 *********************************************************************/
 
UNUSED const double SOLAR_MASS = 1.9884e33;
 
UNUSED const double SOLAR_LUMINOSITY = 3.8427e33;
 
/* >>refer      phys    const   http://pdg.lbl.gov/2010/reviews/rpp2010-rev-astrophysical-constants.pdf */
UNUSED const double AU = 1.49597870700e13;
 
/*********************************************************************
 * fundamental constants go next, eventually rest should be defined  *
 * in terms of these, these are Codata 2010 values.                  *
 *********************************************************************/
 
UNUSED const double ATOMIC_MASS_UNIT = 1.660538921e-24;
 
UNUSED const double ELECTRON_MASS = 9.10938291e-28;
 
UNUSED const double PROTON_MASS = 1.672621777e-24;
 
UNUSED const double BOLTZMANN = 1.3806488e-16;
 
UNUSED const double SPEEDLIGHT = 2.99792458e10;
 
UNUSED const double HPLANCK = 6.62606957e-27;
 
UNUSED const double AVOGADRO = 6.0221415e23;
 
UNUSED const double GRAV_CONST = 6.67384e-8;
 
UNUSED const double ELEM_CHARGE = 1.602176565e-19;
 
UNUSED const double RYD_INF = 1.0973731568539e5;
 
UNUSED const double HIONPOT = 0.999466508345;
 
/*********************************************************************
 * below here should be derived constants                            *
 *                                                                   *
 * NB - explicit values in comments are approximate                  *
 *      and are not maintained !                                     *
 *********************************************************************/
 
UNUSED const double AS1RAD = RADIAN*3600.;
 
UNUSED const double SQAS1SR = pow2(AS1RAD);
 
UNUSED const double SQAS_SKY = PI4*SQAS1SR;
 
UNUSED const double PARSEC = AU*AS1RAD;
 
UNUSED const double MEGAPARSEC = 1.e6*PARSEC;
 
UNUSED const double H_BAR = HPLANCK/(2.*PI);
 
UNUSED const double ELEM_CHARGE_ESU = ELEM_CHARGE*SPEEDLIGHT/10.;
 
UNUSED const double ELECTRIC_CONST = 1.e11/(PI4*pow2(SPEEDLIGHT));
 
UNUSED const double HION_LTE_POP = pow2(HPLANCK)/(PI2*BOLTZMANN*ELECTRON_MASS);
 
UNUSED const double SAHA = sqrt(pow3(HION_LTE_POP));
 
UNUSED const double ERG1CM = HPLANCK*SPEEDLIGHT;
 
UNUSED const double T1CM = HPLANCK*SPEEDLIGHT/BOLTZMANN;
 
UNUSED const double KJMOL1CM = ERG1CM*AVOGADRO/1e10;
 
UNUSED const double WAVNRYD = 1./RYD_INF;
 
UNUSED const double RYDLAM = 1.e8/RYD_INF;
 
UNUSED const double EN1RYD = HPLANCK*SPEEDLIGHT*RYD_INF;
 
UNUSED const double TE1RYD = HPLANCK*SPEEDLIGHT*RYD_INF/BOLTZMANN;
 
UNUSED const double EVDEGK = ELEM_CHARGE*1.e7/BOLTZMANN;
 
UNUSED const double EVRYD = HPLANCK*SPEEDLIGHT*RYD_INF/ELEM_CHARGE*1.e-7;
 
UNUSED const double EN1EV = EN1RYD/EVRYD;
 
UNUSED const double FR1RYD = SPEEDLIGHT*RYD_INF;
 
UNUSED const double HNU3C2 = 2.*HPLANCK*SPEEDLIGHT*pow3(RYD_INF);
 
UNUSED const double FR1RYDHYD = SPEEDLIGHT*RYD_INF*HIONPOT;
 
UNUSED const double HBAReV = H_BAR/EN1EV;  
 
UNUSED const double RYDLAMHYD = RYDLAM/HIONPOT;
 
UNUSED const double STEFAN_BOLTZ = pow2(PI*pow2(BOLTZMANN))/(60.*pow3(H_BAR)*pow2(SPEEDLIGHT));
 
UNUSED const double FREQ_1EV = SPEEDLIGHT*RYD_INF/EVRYD;
 
UNUSED const double FINE_STRUCTURE = pow2(ELEM_CHARGE_ESU)/SPEEDLIGHT/H_BAR;
 
UNUSED const double FINE_STRUCTURE2 = pow2(FINE_STRUCTURE);
 
UNUSED const double BOHR_RADIUS_CM = FINE_STRUCTURE/(PI4*RYD_INF);
 
UNUSED const double TWO_PHOT_CONST = 9.*pow3(FINE_STRUCTURE2)*FR1RYD/2048.;
 
UNUSED const double COLL_CONST = SAHA*BOLTZMANN/HPLANCK;
 
UNUSED const double MILNE_CONST = SPEEDLIGHT*sqrt(pow3(FINE_STRUCTURE2)*pow3(TE1RYD)/PI);
 
UNUSED const double TRANS_PROB_CONST = PI4*HPLANCK*FINE_STRUCTURE/ELECTRON_MASS;
 
#endif /* PHYSCONST_H_ */