DPMJET/flukapro/(NCDNVP)

   1 *$ CREATE NCDNVP.ADD
   2 *COPY NCDNVP
   3 *
   4 *=== ncdnvp ===========================================================*
   5 *
   6 *----------------------------------------------------------------------*
   7 *                                                                      *
   8 *     Used only by PReequilibrium INitialization routines              *
   9 *                                                                      *
  10 *     Created on   22 june 1994    by    Alfredo Ferrari & Paola Sala  *
  11 *                                                   Infn - Milan       *
  12 *                                                                      *
  13 *     Last change on 02-oct-95     by    Alfredo Ferrari               *
  14 *                                                                      *
  15 *     Rhapsf (i,l) = a (skin) parameter for the density distribution,  *
  16 *                    or s-shell armonic oscillator constant (a_s) for  *
  17 *                    nuclei using shell model density distributions    *
  18 *                    already including (possible) folding distributions*
  19 *     Rhcpsf (i,l) = c (R1/2) parameter for the density distribution,  *
  20 *                    or p-shell armonic oscillator constant (a_p) for  *
  21 *                    nuclei using shell model density distributions    *
  22 *                    already including (possible) folding distributions*
  23 *     Rh0psf (i,l) = Rh0 (central density) parameter for the density   *
  24 *                    distribution, or p-shell armonic oscillator       *
  25 *                    constant (a_p) for nuclei using shell model dens- *
  26 *                    ity distributions NOT including (possible) fold-  *
  27 *                    ing distributions                                 *
  28 *     Rhcesf (i,l) = Rhocen, central density for the p/n density       *
  29 *                    distribution                                      *
  30 *     Rhrdsf (i,l) = radius at which the i,l density distribution drops*
  31 *                    at Rhfrsf (i,l) times the central value           *
  32 *     Rrmssf (i,l) = r-rms (<r^2>) parameter for the density distri-   *
  33 *                    bution                                            *
  34 *     Rreqsf (i,l) = equivalent sharp sphere radius for the density    *
  35 *                    distributions                                     *
  36 *                  ( i = 1 : proton, i = 2 : neutron                   *
  37 *                    l = 1 charge (matter) proton (neutron) distribu-  *
  38 *                        tion, integral = Z (N)                        *
  39 *                    l = 2 point/centre proton or neutron distribu-    *
  40 *                        tion, folded with folrms, integral = Z (N)    *
  41 *                    l = 3 point/centre proton or neutron distribu-    *
  42 *                        tion, folded with folrms, integral = Z (N)    *
  43 *                    l = 4 point/centre proton or neutron distribu-    *
  44 *                        tion, folded with folrms, integral = Z (N)    *
  45 *                    l = 5 point/centre proton (neutron) distribution, *
  46 *                        integral = Z (N)                              *
  47 *                    l = 6 proton (neutron) distribution, folded for   *
  48 *                        pion potential, integral = Z (N) )            *
  49 *       Sshlnc (i) = Number of nucleons of type i in the s shell (for  *
  50 *                    shell model like calculations, i = 1 : proton,    *
  51 *                    i = 2 : neutron )                                 *
  52 *       Pshlnc (i) = Number of nucleons of type i in the p shell (for  *
  53 *                    shell model like calculations, i = 1 : proton,    *
  54 *                    i = 2 : neutron )                                 *
  55 *     Vpapsf (i,k) = a (skin) parameter for the potential distribution *
  56 *                    or s-shell armonic oscillator constant (a_s) for  *
  57 *                    nuclei using shell model density like potential   *
  58 *                    distributions already including (possible) fold-  *
  59 *                    ing distributions                                 *
  60 *     Vpcpsf (i,k) = c (R1/2) parameter for the potential distribution *
  61 *                    or p-shell armonic oscillator constant (a_p) for  *
  62 *                    nuclei using shell model density like potential   *
  63 *                    distributions already including (possible) fold-  *
  64 *                    ing distributions                                 *
  65 *     Vp0psf (i,k) = Vp0 (central potential) parameter for the poten-  *
  66 *                    tial distribution distribution, or p-shell armo-  *
  67 *                    nic oscillator constant (a_p) for nuclei using    *
  68 *                    shell model density distributions NOT including   *
  69 *                    (possible) folding distributions                  *
  70 *     Vpdpsf (i,k) = Vp (central potential) depth for the potential    *
  71 *                    distribution                                      *
  72 *     Vpedsf (i,k) = Vp "edge" depth at the nucleus p/n radius         *
  73 *     Vprdsf (i,k) = Radius defining the potential "edge"              *
  74 *     Rcvpsf (i,k) = correlation r_0 for saturation correction for po- *
  75 *                    tential distributions                             *
  76 *     Vpbpsf (i,k) = b (skin-like) parameter for the potential distri- *
  77 *                    bution                                            *
  78 *     Vrmssf (i,k) = r-rms (<r^2>) parameter for the potential distri- *
  79 *                    bution                                            *
  80 *     Vreqsf (i,k) = sharp radius parameter for the potential distri-  *
  81 *                    bution                                            *
  82 *       Pcvpsf (k) = density exponent for saturation correction for po-*
  83 *                    tential distributions                             *
  84 *       Vpexsf (k) = volume integral of the contact like (Pauli) part  *
  85 *                    of the exchange interaction                       *
  86 *           Akvppp = relative strength of pp(nn) interaction           *
  87 *           Akvppn = relative strength of pn(np) interaction           *
  88 *           Bkvppp = relative strength of pp(nn) exchange interaction  *
  89 *           Bkvppn = relative strength of pn(np) exchange interaction  *
  90 *       Vbarnn (k) = volume integral of the NN interaction (GeV x fm^3)*
  91 *           Jnrhvp = present nucleon index (1=proton, 2=neutron)       *
  92 *           Jbouvp = present bound or unbound nucleon index (1=bound,  *
  93 *                    2=unbound)                                        *
  94 *           Jrhodn = present density index 1-6                         *
  95 *           Jrhofl = density index (1-6) to be (un)folded              *
  96 *           Lshmdu = flag for shell model like density calculations    *
  97 *                                                                      *
  98 *----------------------------------------------------------------------*
  99 *
 100 *  R_0 for the equivalent sharp sphere radius for the proton density
 101 *  (or charge) distribution
 102       PARAMETER ( R0PRCH = 1.128 D+00 )
 103 *  R_0 fundamental one according to Myers
 104       PARAMETER ( R0DRPM = 1.18  D+00 )
 105 *  Various parameters of the Droplet model *
 106       PARAMETER ( CKDRPM = 240.   D+00 )
 107       PARAMETER ( CJDRPM = 36.8   D+00 )
 108       PARAMETER ( CLDRPM = 100.   D+00 )
 109       PARAMETER ( CMDRPM = 0.     D+00 )
 110       PARAMETER ( CQDRPM = 17.    D+00 )
 111       PARAMETER ( CA1DRP = 15.96  D+00 )
 112       PARAMETER ( CA2DRP = 20.69  D+00 )
 113       PARAMETER ( CA3DRP = 0.     D+00 )
 114       PARAMETER ( CC1DRP = THRTHR / FIVFIV * COUGFM / R0DRPM * GEVMEV )
 115       PARAMETER ( CC2DRP = CC1DRP / 336.D+00 * ( ONEONE / CJDRPM
 116      &                   + 18.D+00 / CKDRPM ) )
 117       PARAMETER ( HC3DRP = CC1DRP * HLFHLF * FIVFIV / R0DRPM / R0DRPM )
 118 *  ( 3 / 2 pi )^2/3
 119       PARAMETER ( HC4DRP = 0.610887057710857 D+00 )
 120       PARAMETER ( CC4DRP = FIVFIV / FOUFOU * HC4DRP * CC1DRP )
 121       PARAMETER ( CC5DRP = CC1DRP * CC1DRP / CQDRPM / 64.D+00 )
 122       PARAMETER ( R3TOVL = FOUFOU * PIPIPI / THRTHR )
 123       PARAMETER ( R0SSHM = HLFHLF / PIPIPI / ERFA00 )
 124       PARAMETER ( R0PSHM = ONETHI / PIPIPI / ERFA00 )
 125       LOGICAL LSHMDU, LUNFRU
 126       COMMON / NCDNVP / RHCPSF (2,8), RHAPSF (2,8), RH0PSF (2,8),
 127      &                  RHCESF (2,8), RHRDSF (2,8), RHFRSF (2,8),
 128      &                  RRMSSF (2,8), RREQSF (2,8), VPCPSF (2,2),
 129      &                  VPAPSF (2,2), VP0PSF (2,2), VPDPSF (2,2),
 130      &                  VPEDSF (2,2), VPRDSF (2,2), VRMSSF (2,2),
 131      &                  VREQSF (2,2), VPBPSF (2,2), RCVPSF (2,2),
 132      &                  PCVPSF  (2) , VBARNN  (2) , VPEXSF  (2),
 133      &                  SSHLNC  (2) , PSHLNC  (2) , AKVPPP, AKVPPN,
 134      &                  BKVPPP, BKVPPN,
 135      &                  JRHDNF  (2) , JNRHVP, JBOUVP, JRHODN, JRHOFL,
 136      &                  LSHMDU, LUNFRU
 137