DPMJET/flukapro/(BREANG)

   1 *$ CREATE BREANG.ADD
   2 *COPY BREANG
   3 *
   4 *=== breang ===========================================================*
   5 *
   6 *----------------------------------------------------------------------*
   7 *                                                                      *
   8 *     Created on  22  march 1991   by    Alfredo Ferrari & Paola Sala  *
   9 *                                                   Infn - Milan       *
  10 *                                                                      *
  11 *     Last change on 22-mar-91     by    Alfredo Ferrari               *
  12 *                                                                      *
  13 *     Included in the following routines:                              *
  14 *                                                                      *
  15 *                      BDANDI                                          *
  16 *                      BREMNW                                          *
  17 *                      BREMS                                           *
  18 *                                                                      *
  19 *         Mxphtb = number of intervals on which phi(y) is tabulated    *
  20 *         Yphimn = minimum value of y for which phi is taublated       *
  21 *         Yphirt = ratio of y between two tabulated points             *
  22 *         Yphimx = maximum value of y for which phi is taublated       *
  23 *                  (=Yphimn*Yphirt**(Mxphtb-1))                        *
  24 *         Phiytb = tabulated values of phi                             *
  25 *         Alymin = natural logarithm of Yphimn                         *
  26 *         Alyrat = natural logarithm of Yphirt                         *
  27 *         Alytra = natural logarithm of the max. y allowed for Xsitra  *
  28 *         Phia00 = used for asymptotic behaviour                       *
  29 *         Phib00 = used for asymptotic behaviour                       *
  30 *         Phic0  = used for asymptotic behaviour                       *
  31 *         Phid0  = used for asymptotic behaviour                       *
  32 *        ( for y > yphimx, Phi = Phia00 + Phib00 / y^2                 *
  33 *          for y < yphimn, Phi = log(y) + Phic0 + Phid0 x y^2 )        *
  34 *         Accrit = parameter used for the Migdal polarization effect,  *
  35 *                  given by Nel x Lambda_compt^2 x r0 / pi             *
  36 *         Zbrem  = approximate "effective" Z for bremsstrahlung        *
  37 *         Fcoul  = Coulomb factor fc(Z)                                *
  38 *         Gmolie = factor entering Moliere's expansion of Thomas-Fermi *
  39 *                  form factor ( = Z^1/3/121 )                         *
  40 *         Algmol = natural logarithm of Gmolie                         *
  41 *                                                                      *
  42 *----------------------------------------------------------------------*
  43 *
  44 *  Tpifsc = 2 x pi x fine structure constant
  45       PARAMETER ( TPIFSC = 4.5850621648360624 D-02 )
  46       PARAMETER ( MXPHTB = 100 )
  47       PARAMETER ( YPHIMN = 1.0  D-01 )
  48       PARAMETER ( YPHIRT = 1.07 D+00 )
  49 *
  50       COMMON / BREANG / PHIYTB (MXPHTB), YPHIMX, ALYMIN, ALYRAT,
  51      &                  ALYTRA, PHIA00, PHIB00, PHIC0, PHID0,
  52      &                  ACCRIT (MXXMDE), ZBREM  (MXXMDE),
  53      &                  FCOUL  (MXXMDE), GMOLIE (MXXMDE),
  54      &                  ALGMOL (MXXMDE)
  55