*$ CREATE BREANG.ADD
*COPY BREANG
*
*=== breang ===========================================================*
*
*----------------------------------------------------------------------*
*                                                                      *
*     Copyright (C) 1991-2005      by    Alfredo Ferrari & Paola Sala  *
*     All Rights Reserved.                                             *
*                                                                      *
*                                                                      *
*     Created on  22  march 1991   by    Alfredo Ferrari & Paola Sala  *
*                                                   Infn - Milan       *
*                                                                      *
*     Last change on 22-mar-91     by    Alfredo Ferrari               *
*                                                                      *
*     Included in the following routines:                              *
*                                                                      *
*                      BDANDI                                          *
*                      BREMNW                                          *
*                      BREMS                                           *
*                                                                      *
*         Mxphtb = number of intervals on which phi(y) is tabulated    *
*         Yphimn = minimum value of y for which phi is taublated       *
*         Yphirt = ratio of y between two tabulated points             *
*         Yphimx = maximum value of y for which phi is taublated       *
*                  (=Yphimn*Yphirt**(Mxphtb-1))                        *
*         Phiytb = tabulated values of phi                             *
*         Alymin = natural logarithm of Yphimn                         *
*         Alyrat = natural logarithm of Yphirt                         *
*         Alytra = natural logarithm of the max. y allowed for Xsitra  *
*         Phia00 = used for asymptotic behaviour                       *
*         Phib00 = used for asymptotic behaviour                       *
*         Phic0  = used for asymptotic behaviour                       *
*         Phid0  = used for asymptotic behaviour                       *
*        ( for y > yphimx, Phi = Phia00 + Phib00 / y^2                 *
*          for y < yphimn, Phi = log(y) + Phic0 + Phid0 x y^2 )        *
*         Accrit = parameter used for the Migdal polarization effect,  *
*                  given by Nel x Lambda_compt^2 x r0 / pi             *
*         Zbrem  = approximate "effective" Z for bremsstrahlung        *
*         Fcoul  = Coulomb factor fc(Z)                                *
*         Gmolie = factor entering Moliere's expansion of Thomas-Fermi *
*                  form factor ( = Z^1/3/121 )                         *
*         Algmol = natural logarithm of Gmolie                         *
*                                                                      *
*----------------------------------------------------------------------*
*
*  Tpifsc = 2 x pi x fine structure constant
      PARAMETER ( TPIFSC = 4.5850621648360624 D-02 )
      PARAMETER ( MXPHTB = 100 )
      PARAMETER ( YPHIMN = 1.0  D-01 )
      PARAMETER ( YPHIRT = 1.07 D+00 )
*
      COMMON / BREANG / PHIYTB (MXPHTB), YPHIMX, ALYMIN, ALYRAT,
     &                  ALYTRA, PHIA00, PHIB00, PHIC0, PHID0,
     &                  ACCRIT (MXXMDE), ZBREM  (MXXMDE),
     &                  FCOUL  (MXXMDE), GMOLIE (MXXMDE),
     &                  ALGMOL (MXXMDE)
      SAVE / BREANG /