/*
 * DO NOT MANUALLY UPDATE THIS FILE.
 * If needed: make cpp_headers
 *
 * fortran2c.py - automatic conversion
 *       dblprc.inc -> dblprc.h
 */

#ifndef __DBLPRC__
#define __DBLPRC__

#include "types.h"

#ifdef __cplusplus
	extern "C" {
#endif

//                                                                      *
//=== dblprc ===========================================================*
//                                                                      *
//----------------------------------------------------------------------*
//                                                                      *
//     Copyright (C) 2003-2019:  CERN & INFN                            *
//     All Rights Reserved.                                             *
//                                                                      *
//      Dblprc: included in any routine, machine, mathematical and      *
//              physical constants plus global declarations             *
//                                                                      *
//  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!   *
//  !!!! O N   M A C H I N E S   W H E R E   T H E   D O U B L E !!!!   *
//  !!!! P R E C I S I O N   I S   N O T   R E Q U I R E D  R E -!!!!   *
//  !!!! M O V E   T H E   D O U B L E   P R E C I S I O N       !!!!   *
//  !!!! S T A T E M E N T,  S E T   K A L G N M = 1   A N D     !!!!   *
//  !!!! C H A N G E   A L L   N U M E R I C A L   C O N S -     !!!!   *
//  !!!! T A N T S   T O   S I N G L E   P R E C I S I O N       !!!!   *
//  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!   *
//                                                                      *
//         Kalgnm = real address alignment, 2 for double precision,     *
//                  1 for single precision                              *
//         Kalch8 = character*8 address alignment wrt the precision     *
//                  defined by Kalgnm (mostly 1 in all situations)      *
//         I2algn = integer*2 address alignment wrt the normal integer  *
//                  precision (mostly 2, 4 for 64 bit integers)         *
//         Anglgb = this parameter should be set equal to the machine   *
//                  "zero" with respect to unit                         *
//         Anglsq = this parameter should be set equal to the square    *
//                  of Anglgb                                           *
//         Axcssv = this parameter should be set equal to the number    *
//                  for which unity is negligible for the machine       *
//                  accuracy                                            *
//         Andrfl = "underflow" of the machine for floating point       *
//                  operation                                           *
//         Avrflw = "overflow"  of the machine for floating point       *
//                  operation                                           *
//         Ainfnt = code "infinite"                                     *
//         Azrzrz = code "zero"                                         *
//         Einfnt = natural logarithm of the code "infinite"            *
//         Ezrzrz = natural logarithm of the code "zero"                *
//         Excssv = natural logarithm of the code number for which      *
//                  unit is negligible                                  *
//         Englgb = natural logarithm of the code "zero" with respect   *
//                  to unit                                             *
//         Onemns = 1- of the machine, it is 1 - 2 x Anglgb             *
//         Onepls = 1+ of the machine, it is 1 + 2 x Anglgb             *
//         Csnnrm = maximum tolerable error on cosine normalization,    *
//                  u**2+v**2+w**2: assuming a typical anglgb relative  *
//                  error on each component we would get 2xanglgb: use  *
//                  4xanglgb to avoid too many normalizations           *
//         Dmxtrn = "infinite" distance for transport (cm)              *
//         Rhflmn = minimal density for Fluka (g/cm^3)                  *
//                                                                      *
//   "Global" declarations:                                             *
//         Lfluka = set to true for a real (full) Fluka run             *
//         Lgbias = set to true for a fully biased run                  *
//         Lgbana = set to true for a fully analogue run                *
//         Lflgeo = set to true when using the standard Fluka geometry  *
//         Loflts = set to true for special off-line testing of speci-  *
//                  fic routines                                        *
//         Lusrin = set to true if the user dependent initialization    *
//                  routine Usrini has been called at least onec        *
//         Lnmgeo = set to true for a name-based geometry input         *
//         Lnminp = set to true for a name-based Fluka input            *
//         Lfrfmt = set to true for a free-format based Fluka input     *
//         Ldpmcr = set to true for allowing dump cores inside Flabrt   *
//----------------------------------------------------------------------*
//                                                                      *
const int    KALGNM = 2;
const int    KALCH8 = 1;
const int    I2ALGN = 2;
const double ANGLGB = 5.0E-16;
const double ANGLSQ = 2.5E-31;
const double AXCSSV = 0.2E+16;
const double ANDRFL = 1.0E-38;
const double AVRFLW = 1.0E+38;
const double AINFNT = 1.0E+30;
const double AZRZRZ = 1.0E-30;
const double EINFNT = +69.07755278982137E+00;
const double EZRZRZ = -69.07755278982137E+00;
const double EXCSSV = +35.23192357547063E+00;
const double ENGLGB = -35.23192357547063E+00;
const double ONEMNS = 0.999999999999999E+00;
const double ONEPLS = 1.000000000000001E+00;
const double CSNNRM = 2.0E-15;
const double DMXTRN = 1.0E+15*1.000001378E+00;
//      PARAMETER ( DMXTRN = 1.0D+15 )
//      PARAMETER ( RHFLMN = 1.0D-13 )
//
//======================================================================*
//======================================================================*
//=========                                                   ==========*
//=========    M A T H E M A T I C A L   C O N S T A N T S    ==========*
//=========                                                   ==========*
//======================================================================*
//======================================================================*
//                                                                      *
//   Numerical constants (long integer):                                *
//                                                                      *
//         I8zero = "0" for Integer*8                                   *
//                                                                      *
//   Numerical constants (single precision):                            *
//                                                                      *
//         Zersng = 0                                                   *
//                                                                      *
//   Numerical constants (double precision):                            *
//                                                                      *
//         Zerzer = 0                                                   *
//         Oneone = 1                                                   *
//         Twotwo = 2                                                   *
//         Thrthr = 3                                                   *
//         Foufou = 4                                                   *
//         Fivfiv = 5                                                   *
//         Sixsix = 6                                                   *
//         Sevsev = 7                                                   *
//         Eigeig = 8                                                   *
//         Aninen = 9                                                   *
//         Tenten = 10                                                  *
//         Eleven = 11                                                  *
//         Twelve = 12                                                  *
//         Fiften = 15                                                  *
//         Sixten = 16                                                  *
//         Hlfhlf = 1/2                                                 *
//         Onethi = 1/3                                                 *
//         Onefou = 1/4                                                 *
//         Onefiv = 1/5                                                 *
//         Onesix = 1/6                                                 *
//         Onesev = 1/7                                                 *
//         Oneeig = 1/8                                                 *
//         Twothi = 2/3                                                 *
//         Thrfou = 3/4                                                 *
//         Thrtwo = 3/2                                                 *
//         Fouthr = 4/3                                                 *
//         Pipipi = Circumference / diameter                            *
//         Twopip = 2 x Pipipi                                          *
//         Pip5o2 = 5/2 x Pipipi                                        *
//         Pipisq = Pipipi x Pipipi                                     *
//         Pihalf = 1/2 x Pipipi                                        *
//         R3tovl = 4/3 x Pipipi                                        *
//         Erfa00 = Erf (oo) = 1/2 x square root of pi                  *
//         Sqrtpi = square root of pi                                   *
//         Sqtwpi = square root of 2xpi                                 *
//         Eulero = Eulero's constant                                   *
//         Eulexp = exp ( Eulero )                                      *
//         E1m2eu = exp ( 1 - 2 eulero )                                *
//         Eneper = "e", base of natural logarithm                      *
//         Sqrent = square root of "e"                                  *
//         Sqrtwo = square root of  2                                   *
//         Sqrthr = square root of  3                                   *
//         Sqrfiv = square root of  5                                   *
//         Sqrsix = square root of  6                                   *
//         Sqrsev = square root of  7                                   *
//         Sqrt12 = square root of 12                                   *
//         S2fwhm = 2 x square root of 2 x logarithm of 2               *
//         Twolog = logarithm of 2                                      *
//         Two2o3 = 2^2/3                                               *
//         Tenlog = logarithm of 10                                     *
//         Atnfou = Atan[4]                                             *
//         Atnsix = Atan[6]                                             *
//         Alpemd = 2 / exp(Eulero) / Sqrt[e]                           *
//                                                                      *
//----------------------------------------------------------------------*
//
const double RHFLMN = 1.0E-10;
const long long I8ZERO = 0;
const float ZERSNG = 0.E+00;
const double ZERZER = 0.E+00;
const double ONEONE = 1.E+00;
const double TWOTWO = 2.E+00;
const double THRTHR = 3.E+00;
const double FOUFOU = 4.E+00;
const double FIVFIV = 5.E+00;
const double SIXSIX = 6.E+00;
const double SEVSEV = 7.E+00;
const double EIGEIG = 8.E+00;
const double ANINEN = 9.E+00;
const double TENTEN = 10.E+00;
const double ELEVEN = 11.E+00;
const double TWELVE = 12.E+00;
const double FIFTEN = 15.E+00;
const double SIXTEN = 16.E+00;
const double HLFHLF = 0.5E+00;
const double ONETHI = ONEONE/THRTHR;
const double ONEFOU = ONEONE/FOUFOU;
const double ONEFIV = ONEONE/FIVFIV;
const double ONESIX = ONEONE/SIXSIX;
const double ONESEV = ONEONE/SEVSEV;
const double ONEEIG = ONEONE/EIGEIG;
const double TWOTHI = TWOTWO/THRTHR;
const double THRFOU = THRTHR/FOUFOU;
const double THRTWO = THRTHR/TWOTWO;
const double FOUTHR = FOUFOU/THRTHR;
const double PIPIPI = 3.141592653589793238462643383279E+00;
const double TWOPIP = 6.283185307179586476925286766559E+00;
const double PIP5O2 = 7.853981633974483096156608458199E+00;
const double PIPISQ = 9.869604401089358618834490999876E+00;
const double PIHALF = 1.570796326794896619231321691640E+00;
const double R3TOVL = FOUFOU*PIPIPI/THRTHR;
const double ERFA00 = 0.886226925452758013649083741671E+00;
const double SQRTPI = 1.772453850905516027298167483341E+00;
const double SQTWPI = 2.506628274631000502415765284811E+00;
const double EULERO = 0.577215664901532860606512E+00;
const double EULEXP = 1.781072417990197985236504E+00;
const double EULLOG = -0.5495393129816448223376619E+00;
const double E1M2EU = 0.8569023337737540831433017E+00;
const double ENEPER = 2.718281828459045235360287471353E+00;
const double SQRENT = 1.648721270700128146848650787814E+00;
const double SQRTWO = 1.414213562373095048801688724210E+00;
const double SQRTHR = 1.732050807568877293527446341506E+00;
const double SQRFIV = 2.236067977499789696409173668731E+00;
const double SQRSIX = 2.449489742783178098197284074706E+00;
const double SQRSEV = 2.645751311064590590501615753639E+00;
const double SQRT12 = 3.464101615137754587054892683012E+00;
const double S2FWHM = 2.354820045030949382023138652919E+00;
const double TWOLOG = 0.693147180559945309417232121458E+00;
const double TWO2O3 = 1.587401051968199474751705639272E+00;
const double TENLOG = 2.302585092994045684017991454684E+00;
const double ATNFOU = 1.3258176636680326E+00;
const double ATNSIX = 1.4056476493802699E+00;
//
//======================================================================*
//======================================================================*
//=========                                                   ==========*
//=========       P H Y S I C A L   C O N S T A N T S         ==========*
//=========                                                   ==========*
//======================================================================*
//======================================================================*
//                                                                      *
//   Primary constants:                                                 *
//                                                                      *
//         Clight = speed of light in cm s-1                            *
//         Avogad = Avogadro number                                     *
//         Boltzm = k Boltzmann constant (J K-1)                        *
//         Amelgr = electron mass (g)                                   *
//         Plckbr = reduced Planck constant (erg s)                     *
//         Elccgs = elementary charge (CGS unit)                        *
//         Elcmks = elementary charge (MKS unit)                        *
//         Pvcmu0 = vacuum permeability constant (MKS unit)             *
//         Epsil0 = vacuum dielectric constant (MKS unit)               *
//         Amugrm = Atomic mass unit (g)                                *
//         Ammumu = Muon    mass (amu)                                  *
//         Amprmu = Proton  mass (amu)                                  *
//         Amnemu = Neutron mass (amu)                                  *
//                                                                      *
//   Derived constants:                                                 *
//                                                                      *
//         Alamb0 = Compton wavelength = 2 pi r0 / fsc                  *
//         Alpfsc = Fine structure constant  = e^2/(hbar c) (CGS units) *
//         Note:                                                        *
//          e[MKS]^2/ (4 pi epsilon0 hbar[MKS] c[MKS]) = alpha = 1/137  *
//          e[CGS]^2/ (hbar[CGS] c[CGS]) = alpha = 1/137                *
//          with c[MKs] = 10^-2 c[CGS], hbar[MKS] = 10^-7 hbar[CGS]     *
//         Amelct = Electron mass (GeV) = 10^-16Amelgr Clight^2 / Elcmks*
//         Amugev = Atomic mass unit (GeV) = 10^-16Amugrm Clight^2      *
//                                           / Elcmks                   *
//         Ammuon = Muon    mass (GeV) = Ammumu * Amugev                *
//         Amprtn = Proton  mass (GeV) = Amprmu * Amugev                *
//         Amntrn = Neutron mass (GeV) = Amnemu * Amugev                *
//         Amdeut = Deuteron mass (GeV)                                 *
//         Bnddtn = Deuteron binding energy (GeV)                       *
//         Amtrtn = Triton   mass (GeV) (derived from the CODATA value  *
//                  in kg, and the amu value in kg->GeV)                *
//         Am3hel = 3-He     mass (GeV) (derived from the CODATA value  *
//                  in kg, and the amu value in kg->GeV)                *
//         Amalph = Alpha    mass (GeV) (derived from the excess mass   *
//                  and an (approximate) atomic binding not a really    *
//                  measured constant)                                  *
//         Ampiop = Pion+  mass (GeV)                                   *
//         Ampiom = Pion-  mass (GeV)                                   *
//         Ampio0 = Pion0  mass (GeV)                                   *
//         Cougfm = e^2 (GeV fm) = Elccgs^2 / Elcmks * 10^-7 * 10^-9    *
//                * 10^13 (10^..=erg cm->joule cm->GeV cm->GeV fm       *
//                  it is equal to 0.00144 GeV fm                       *
//         Fscto2 = (Fine structure constant)^2                         *
//         Fscto3 = (Fine structure constant)^3                         *
//         Fscto4 = (Fine structure constant)^4                         *
//         Plabrc = Reduced Planck constant times the light velocity    *
//                  expressed in GeV fm                                 *
//         Rclsel = Classical electron radius (cm) = e^2 / (m_e c^2)    *
//         Bltzmn = k Boltzmann constant in GeV K-1                     *
//         A0bohr = Bohr radius, hbar^2 / ( m_e e^2) (fm) = Plabrc**2   *
//                / Amelct / Cougfm, or equivalently,                   *
//                Plabrc / Alpfsc / Amelct                              *
//         Gfohb3 = Fermi constant, G_f/(hbar c)^3, in GeV^-2           *
//         Gfermi = Fermi constant in GeV fm^3                          *
//         Sin2tw = sin^2 theta_Weinberg                                *
//         Prmgnm = proton  magnetic moment (magneton)                  *
//         Anmgnm = neutron magnetic moment (magneton)                  *
//         S0Thms = sigma_0 Thomson, 8/3 pi r_e^2 (mb)                  *
//         Plfrcs = plasma frequency constant = 28.816 eV               *
//                = Sqrt[ 4 pi alpha hbar^3 / m_e] with hbar in [eV cm] *
//                  and m_e in [eV]                                     *
//                                                                      *
//   Astronomical constants:                                            *
//                                                                      *
//         Rearth = Earth equatorial radius (cm)                        *
//         Auastu = Astronomical Unit       (cm)                        *
//                                                                      *
//   Conversion constants:                                              *
//                                                                      *
//         GeVMeV = from GeV to MeV                                     *
//         eMVGeV = from MeV to GeV                                     *
//         GeV2eV = from GeV to  eV                                     *
//         eV2GeV = from eV  to GeV                                     *
//         alGVMV = from GeV to MeV, log                                *
//         Raddeg = from radians to degrees                             *
//         Degrad = from degrees to radians                             *
//         GeVOmg = from (photon) energy [GeV] in 2pi x frequency [s^-1]*
//         Cmq2mb = from square centimetres to millibarns               *
//         Fmb2ba = from millibarns to barns                            *
//         Bar2mb = from barns to millibarns                            *
//         Fmb2fs = from millibarns to fermi^2                          *
//         Fms2mb = from fermi^2 to millibarns                          *
//         Ba2cmq = from barns to square centimetres                    *
//         Cmq2ba = from square centimetres to barns                    *
//         Fmi2cm = from fermi to cm                                    *
//         Cm2fmi = from cm to fermi                                    *
//                                                                      *
//   Useful constants:                                                  *
//                                                                      *
//         Fertho = constant to be used in the Fermi-Thomas approxima-  *
//                  ted expression for atomic binding energies          *
//         Expebn = exponent to be used in the Fermi-Thomas approxima-  *
//                  ted expression for atomic binding energies          *
//                    B_atomic (Z) = Fertho x Z^Expebn (GeV)            *
//         Bex1hy = 1-H atomic binding energies (GeV)                   *
//         Bex2hy = 2-H atomic binding energy   (GeV) as obtained from  *
//                  the 2-H (nuclear) mass, the 2-H atomic excess mass  *
//                  and the amu and electron masses. Consistent within  *
//                  the p and 2-H mass uncertainties with Bex1hy        *
//         Bex3hy = 3-H atomic binding energy   (GeV) as obtained from  *
//                  the 3-H (nuclear) mass, the 3-H atomic excess mass  *
//                  and the amu and electron masses. Consistent within  *
//         Bex3he = 3-He atomic binding energy  (GeV) as obtained from  *
//                  the 3-He (nuclear) mass, the 3-He atomic excess mass*
//                  and the amu and electron masses. Consistent within  *
//                  the alpha and 3-He mass uncertainties with Bex4he   *
//         Bex4he = 4-He atomic binding energy  (GeV) as obtained from  *
//                  the 4-He (nuclear) mass, the 4-He atomic excess mass*
//                  and the amu and electron masses. Consistent within  *
//                  the alpha and 3-He mass uncertainties with Bex4he   *
//         Bexc12 = Fermi-Thomas approximated expression for 12-C ato-  *
//                  mic binding energies (GeV)                          *
//         Bhfzxx = Hartree-Fock atomic binding energies for Z=xx, from *
//                  Adndt18 243 (1976) and Density-Functional Theory of *
//                  Atoms and Molecules, R.G.Parr and W.Yang            *
//                  Please note that the original values are in units   *
//                  of 2 x Bex1hy                                       *
//         Amunmu = difference between the atomic and nuclear mass units*
//         Amuc12 = "Nuclear" mass unit = 1/12 M_nucl (12-C),           *
//                  M_nucl (12-C) = M_atom (12-C) - 6 m_e + B_atom(12-C)*
//         Ex1hmv = 1-H (atomic) excess mass in MeV,                    *
//                  [(Amprtn - Amugev) + Amelct - Bex1hy] x GeVMeV      *
//         Exnemv = Neutron excess mass in MeV,                         *
//                  [Amntrn - Amugev] x GeVMeV                          *
//                                                                      *
//----------------------------------------------------------------------*
//
//  Clight is exact:
const double ALPEMD = TWOTWO/EULEXP/SQRENT;
const double CLIGHT = 2.99792458E+10;
const double AVOGAD = 6.0221367E+23;
const double BOLTZM = 1.380658E-23;
const double AMELGR = 9.1093897E-28;
const double PLCKBR = 1.05457266E-27;
const double ELCCGS = 4.8032068E-10;
const double ELCMKS = 1.60217733E-19;
const double AMUGRM = 1.6605402E-24;
const double AMMUMU = 0.113428913E+00;
const double AMPRMU = 1.007276470E+00;
const double AMNEMU = 1.008664904E+00;
//     PARAMETER ( ALPFSC = 1.D+00 / 137.035989561D+00 )
//  Pvcmu0 is exact:
const double EPSIL0 = 8.854187817E-12;
//  Epsilon0 is exact since c has been fixed:
//     PARAMETER ( FSCTO2 = ALPFSC * ALPFSC )
//     PARAMETER ( FSCTO3 = FSCTO2 * ALPFSC )
//     PARAMETER ( FSCTO4 = FSCTO3 * ALPFSC )
//    It is important to set the electron mass exactly with the same
//    rounding as in the mass tables, so use the explicit expression
//     PARAMETER ( AMELCT = 1.D-16 * AMELGR * CLIGHT * CLIGHT / ELCMKS )
//    It is important to set the amu mass exactly with the same
//    rounding as in the mass tables, so use the explicit expression
//     PARAMETER ( AMUGEV = 1.D-16 * AMUGRM * CLIGHT * CLIGHT / ELCMKS )
//    It is important to set the muon,proton,neutron masses exactly with
//    the same rounding as in the mass tables, so use the explicit
//    expression
//     PARAMETER ( AMMUON = AMMUMU * AMUGEV )
//     PARAMETER ( AMPRTN = AMPRMU * AMUGEV )
//     PARAMETER ( AMNTRN = AMNEMU * AMUGEV )
//     PARAMETER ( RCLSEL = ELCCGS * ELCCGS / CLIGHT / CLIGHT / AMELGR )
//     PARAMETER ( BLTZMN = BOLTZM / ELCMKS * 1.D-09 )
const double PVCMU0 = TWOTWO*TWOPIP*1.E-07;
const double ALPFSC = 7.2973530791728595E-03;
const double FSCTO2 = 5.3251361962113614E-05;
const double FSCTO3 = 3.8859399018437826E-07;
const double FSCTO4 = 2.8357075508200407E-09;
const double PLABRC = 0.197327053E+00;
const double AMELCT = 0.51099906E-03;
const double AMUGEV = 0.93149432E+00;
const double AMMUON = 0.105658389E+00;
const double AMPRTN = 0.93827231E+00;
const double AMNTRN = 0.93956563E+00;
const double AMDEUT = 1.87561339E+00;
const double BNDDTN = 2.250000000000000E-03;
const double AMTRTN = 2.808921005E+00;
//      PARAMETER ( AMALPH = 3.727379239        D+00 )
const double AM3HEL = 2.808391482E+00;
const double AMALPH = 3.72738025692891E+00;
const double AMPIOP = 1.3957018E-01;
const double AMPIOM = AMPIOP;
const double AMPIO0 = 1.349766E-01;
//      PARAMETER ( RCLSEL = 2.8179403267       D-13 )
const double COUGFM = ELCCGS*ELCCGS/ELCMKS*1.E-07*1.E+13*1.E-09;
const double RCLSEL = 2.8179409183694872E-13;
//      PARAMETER ( BLTZMN = 8.6173324          D-14 )
const double ALAMB0 = TWOTWO*PIPIPI*RCLSEL/ALPFSC;
const double BLTZMN = 8.617385E-14;
//      PARAMETER ( GFOHB3 = 1.1663787          D-05 )
const double A0BOHR = PLABRC/ALPFSC/AMELCT;
const double GFOHB3 = 1.16639E-05;
//      PARAMETER ( SIN2TW = 0.23126            D+00 )
//      PARAMETER ( SIN2TW = 0.23155            D+00 )
//      PARAMETER ( PRMGNM = 2.792847356        D+00 )
//      PARAMETER ( ANMGNM =-1.9130427          D+00 )
const double GFERMI = GFOHB3*PLABRC*PLABRC*PLABRC;
const double SIN2TW = 0.2319E+00;
const double PRMGNM = 2.792847386E+00;
const double ANMGNM = -1.91304275E+00;
const double REARTH = 6.378140E+08;
const double AUASTU = 1.4959787066E+13;
const double GEVMEV = 1.0E+03;
const double EV2GEV = 1.0E-09;
const double GEV2EV = 1.0E+09;
const double EMVGEV = 1.0E-03;
const double CMQ2MB = 1.0E+27;
const double FMB2BA = 1.0E-03;
const double BAR2MB = 1.0E+03;
const double FMB2FS = 1.0E-01;
const double FMS2MB = 1.0E+01;
const double BA2CMQ = 1.0E-24;
const double CMQ2BA = 1.0E+24;
const double FMI2CM = 1.0E-13;
const double CM2FMI = 1.0E+13;
const double AKG2GR = 1.0E+03;
const double ALGVMV = 6.90775527898214E+00;
const double RADDEG = 180.E+00/PIPIPI;
const double DEGRAD = PIPIPI/180.E+00;
const double GEVGRA = ELCMKS*GEV2EV*AKG2GR;
const double GEVOMG = CLIGHT*1.E+13/PLABRC;
//  Plasma frequency constant, to be recomputed every time alpha, m_e,
//  hbar or N_av change, it requires a Sqrt ergo it cannot be
//  parametrized explicitly as a function of the above constants
//     PARAMETER ( PLFRCS = SQRT ( ALPFSC * FOUFOU * PIPIPI
//    &                   * ( PLABRC * FMI2CM )**3 / AMELCT * AVOGAD )
//    &                   * GEV2EV )
const double S0THMS = EIGEIG/THRTHR*PIPIPI*RCLSEL*RCLSEL*CMQ2MB;
//  Old Fermi-Thomas parametrization of atomic binding energies:
//     PARAMETER ( FERTHO = 15.73       D-09 )
//     PARAMETER ( EXPEBN = 7.D+00 / 3.D+00  )
//     PARAMETER ( BEXC12 = FERTHO * 65.41634134195703D+00 )
//  New Fermi-Thomas parametrization of atomic binding energies:
const double PLFRCS = 2.8815937226648276E+01;
const double FERTHO = 14.33E-09;
const double EXPEBN = 2.39E+00;
const double BEX1HY = 13.61E-09;
const double BEX2HY = 14.188E-09;
const double BEX3HY = 14.818E-09;
const double BEX3HE = 81.346E-09;
const double BEX4HE = 77.246E-09;
const double BHFZ02 = TWOTWO*BEX1HY*2.8615E+00;
const double BHFZ03 = TWOTWO*BEX1HY*7.4290E+00;
const double BHFZ04 = TWOTWO*BEX1HY*14.569E+00;
const double BHFZ05 = TWOTWO*BEX1HY*24.530E+00;
const double BEXC12 = FERTHO*72.40715579499394E+00;
const double BHFZ06 = TWOTWO*BEX1HY*37.645E+00;
const double BHFZ07 = TWOTWO*BEX1HY*54.303E+00;
const double BHFZ08 = TWOTWO*BEX1HY*74.801E+00;
const double BHFZ09 = TWOTWO*BEX1HY*99.471E+00;
const double BHFZ10 = TWOTWO*BEX1HY*128.55E+00;
const double BHFZ18 = TWOTWO*BEX1HY*526.82E+00;
const double BHFZ36 = TWOTWO*BEX1HY*2752.0E+00;
//      PARAMETER ( AMUNMU = HLFHLF * AMELCT - BHFZ06 / 12.D+00 )
const double BHFZ54 = TWOTWO*BEX1HY*7232.1E+00;
const double AMUNMU = HLFHLF*AMELCT-BEXC12/12.E+00;
//  p and n excess mass in MeV:
const double AMUC12 = AMUGEV-AMUNMU;
const double EX1HMV = ((AMPRTN-AMUGEV)+AMELCT-BEX1HY)*GEVMEV;
//  Old MeV units:
const double EXNEMV = (AMNTRN-AMUGEV)*GEVMEV;
//  Infinite and "zero" T1/2:
const double AMEMEV = GEVMEV*AMELCT;
const double T12INF = 1.E+30;
//
const double T12ZER = 1.E-15;
extern struct {
	logical	lfluka;
	logical	lgbias;
	logical	lgbana;
	logical	lflgeo;
	logical	loflts;
	logical	lusrin;
	logical	lusrgl;
	logical	lnmgeo;
	logical	lnminp;
	logical	lfrfmt;
	logical	ldmpcr;
	logical	lldflr;
	logical	ldpbdy;
	int	kflgeo;
	int	kfldnr;
} global_;

extern struct {
	char	crvrfl[8];
} globch_;

#ifdef __cplusplus
}
#endif
#endif