DPMJET/flukapro/(NUCSFT)

   1 *$ CREATE NUCSFT.ADD
   2 *COPY NUCSFT
   3 *
   4 *=== nucsft ===========================================================*
   5 *
   6 *----------------------------------------------------------------------*
   7 *                                                                      *
   8 *     NUClear Symmetrized Fermi/woods-saxon density and potential      *
   9 *     tabulations                                                      *
  10 *                                                                      *
  11 *     Last change   on   20-jun-94   by   Alfredo Ferrari, INFN-Milan  *
  12 *                                                                      *
  13 *          Vparvp, cparvp, aparvp (bparvp auxiliary parameter          *
  14 *          for Cparvp calculations):                                   *
  15 *                                                                      *
  16 *          V0 (r)    = Vparvp * sinh (cparvp/aparvp)                   *
  17 *                    / ( cosh (cparvp/aparvp) + cosh (r/aparvp) )      *
  18 *          expvca    = exp ( -cparvp / aparvp )                        *
  19 *          expvra    = exp ( rcursf / aparvp )                         *
  20 *          bnvpsf    = "average" binding energy used when generating   *
  21 *                      the potential tabulations                       *
  22 *          v0cesf    = well depth at r=0 for the "average" potential   *
  23 *          evsfmx    = Max. Fermi energy for the ip nucleon potential  *
  24 *                      evsfmx = |v0cesf| - bnvpsf                      *
  25 *          pvsfmx    = Max. Fermi momentum for proton/neutron          *
  26 *          pfsksf    = Fermi momentum for proton/neutron computed at   *
  27 *                      radiu1 defined ad the (mrsfbn-5)th bin limit    *
  28 *          vcoumx    = maximum height of the Coulomb barrier           *
  29 *                      (for Z=Zvpcou), modified during tracking to     *
  30 *                      get the real barrier                            *
  31 *          vcousv    = maximum height of the Coulomb barrier           *
  32 *                      (for Z=Zvpcou)                                  *
  33 *          v0ccou    = Coulomb potential at R=Radtot                   *
  34 *          zvpcou    = Z for which the barrier has been computed       *
  35 *          zbourd    = Z reduction factor to be used for bound nucleons*
  36 *                      when computing Coulomb effects                  *
  37 *          ravpsf    = radius at which the "average" nuclear potential *
  38 *                      stops                                           *
  39 *          ravcou    = radius at which the Coulomb potential stops     *
  40 *          ravpot    = radius at which the nuclear potential stops     *
  41 *          dradsf    = radial increment used from radiu0 up to radtot  *
  42 *          dravsf    = radial increment used from radtot up to ravpot  *
  43 *          ivsfmn    = index of the innermost potential shell which    *
  44 *                      can be reached by the current particle accord-  *
  45 *                      ing to its angular momentum                     *
  46 *          ivsfmx    = index of the outermost potential shell which    *
  47 *                      can be reached by the current particle accord-  *
  48 *                      ing to its energy and angular momentum          *
  49 *          ivcomx    = bin index of the maximum height of the Coulomb  *
  50 *                      barrier                                         *
  51 *          isfint    = current tabulation index                        *
  52 *          ibsfmn    = index of the radial bin of the the minimum      *
  53 *                      approach point                                  *
  54 *          ibsfmx    = index of the radial bin of the the maximum      *
  55 *                      approach point                                  *
  56 *                                                                      *
  57 *     The actual potential is given by (Preeqmvax version with no      *
  58 *     energy dependence of the potential):                             *
  59 *     for r < radtot (note V_nuc < 0):                                 *
  60 *       neutral particles :                                            *
  61 *          V_true(r) = V_nuc(r,ip) + bnvpsf (ip) - bn_true (ip)        *
  62 *       charged particles :                                            *
  63 *          V_true(r) = V_nuc(r,ip) + bnvpsf (ip) - bn_true (ip)        *
  64 *                    + Zproj * V_coul(r)                               *
  65 *     for ravpsf > r > radtot (note V_nuc < 0):                        *
  66 *          V_true(r) = V_nuc(r,ip) / bnvpsf (ip) x bn_true (ip)        *
  67 *                    + Znow / Zvpcou x Zproj x V_coul(r)               *
  68 *     for r > ravpsf :                                                 *
  69 *          V_true(r) = Znow / Zvpcou x Zproj x V_coul(r)               *
  70 *     Please note that Ekfcen(ip)=Evsfmx(ip), for ip=1,2               *
  71 *     For the actual potential with the Preeqmvax version with the     *
  72 *     energy dependence of the potential, see the relevant routines    *
  73 *                                                                      *
  74 *----------------------------------------------------------------------*
  75 *
  76       COMMON / NUCSFT / VPARVP (2), CPARVP (2), APARVP (2), BPARVP (2),
  77      &                  EXPVCA (2), EXPVRA (2), RAVPSF (2), BNVPSF (2),
  78      &                  V0CESF (2), EVSFMX (2), VSFHLP (2), PVSFMX (2),
  79      &                  PFSKSF (2), RAVCOU, RAVPOT, VCOUMX, V0CCOU,
  80      &                  ZVPCOU, ZBOURD, VPSFA0, VPSFA1, VPSFAE, DRADSF,
  81      &                  DRAVSF, VCOUSV,
  82      &                  IVSFMN, IVSFMX, IVCOMX, ISFINT, IBSFMN, IBSFMX
  83