DPMJET/flukapro/(SLNUCM)

   1 *$ CREATE SLNUCM.ADD
   2 *COPY SLNUCM
   3 *
   4 *=== slnucm ===========================================================*
   5 *
   6 *----------------------------------------------------------------------*
   7 *                                                                      *
   8 *     SoLar ( and supernova) NeUtrino CoMmon:                          *
   9 *                                                                      *
  10 *     Created on 29 february 1996  by    Alfredo Ferrari & Paola Sala  *
  11 *                                                   Infn - Milan       *
  12 *                                                                      *
  13 *     Last change on 03-feb-01     by    Alfredo Ferrari               *
  14 *                                                                      *
  15 *                                                                      *
  16 *     Slnasg(j,i) = absorption cross section for jth level and ith     *
  17 *                   energy                                             *
  18 *          Qugtfr = reaction Q for absorption on 40-Ar (40-Ar -> 40-K) *
  19 *       Spsngt(i) = differential spectrum x cross section for solar    *
  20 *                   for Gamow-Teller vue absorption on 40-Ar           *
  21 *       Cusngt(i) = cumulative spectrum x cross section for solar neu- *
  22 *                   trino for Gamow-Teller vue absorption on 40-Ar     *
  23 *          Sngtcn = total interaction rate per year per kTon x f(40-Ar)*
  24 *                   for Gamow-Teller vue absorption on 40-Ar           *
  25 *          Ppcogt = pp    interaction rate per year per kTon x f(40-Ar)*
  26 *                   for Gamow-Teller vue absorption on 40-Ar           *
  27 *          B8cogt = 8-B   interaction rate per year per kTon x f(40-Ar)*
  28 *                   for Gamow-Teller vue absorption on 40-Ar           *
  29 *          Hepcog = hep   interaction rate per year per kTon x f(40-Ar)*
  30 *                   for Gamow-Teller vue absorption on 40-Ar           *
  31 *          Be71cg = 7-Be (1st) inter. rate per year per kTon x f(40-Ar)*
  32 *                   for Gamow-Teller vue absorption on 40-Ar           *
  33 *          Be72cg = 7-Be (2nd) inter. rate per year per kTon x f(40-Ar)*
  34 *                   for Gamow-Teller vue absorption on 40-Ar           *
  35 *          Pepcog = pep   interaction rate per year per kTon x f(40-Ar)*
  36 *                   for Gamow-Teller vue absorption on 40-Ar           *
  37 *       Spsnfr(i) = differential spectrum x cross section for solar    *
  38 *                   for Fermi vue absorption on 40-Ar                  *
  39 *       Cusnfr(i) = cumulative spectrum x cross section for solar neu- *
  40 *                   trino for Fermi vue absorption on 40-Ar            *
  41 *          Snfrcn = total interaction rate per year per kTon x f(40-Ar)*
  42 *                   for Fermi vue absorption on 40-Ar                  *
  43 *          Ppcofr = pp    interaction rate per year per kTon x f(40-Ar)*
  44 *                   for Fermi vue absorption on 40-Ar                  *
  45 *          B8cofr = 8-B   interaction rate per year per kTon x f(40-Ar)*
  46 *                   for Fermi vue absorption on 40-Ar                  *
  47 *          Hepcof = hep   interaction rate per year per kTon x f(40-Ar)*
  48 *                   for Fermi vue absorption on 40-Ar                  *
  49 *          Be71cf = 7-Be (1st) inter. rate per year per kTon x f(40-Ar)*
  50 *                   for Fermi vue absorption on 40-Ar                  *
  51 *          Be72cf = 7-Be (2nd) inter. rate per year per kTon x f(40-Ar)*
  52 *                   for Fermi vue absorption on 40-Ar                  *
  53 *          Pepcof = pep   interaction rate per year per kTon x f(40-Ar)*
  54 *                   for Fermi vue absorption on 40-Ar                  *
  55 *       Spsntt(i) = differential spectrum x cross section for solar    *
  56 *                   neutrino nue-electron elastic interactions         *
  57 *       Cusntt(i) = cumulative spectrum x cross section for solar neu- *
  58 *                   trino nue-electron elastic interactions            *
  59 *          Snttcn = total interaction rate per year per kTon x <Z/A>   *
  60 *                   for elastic vue-electron interactions              *
  61 *          Ppcont = pp    interaction rate per year per kTon x <Z/A>   *
  62 *                   for elastic vue-electron interactions              *
  63 *          B8cont = 8-B   interaction rate per year per kTon x <Z/A>   *
  64 *                   for elastic vue-electron interactions              *
  65 *          Hepcon = hep   interaction rate per year per kTon x <Z/A>   *
  66 *                   for elastic vue-electron interactions              *
  67 *          Be71cn = 7-Be (1st) inter. rate per year per kTon x <Z/A>   *
  68 *                   for elastic vue-electron interactions              *
  69 *          Be72cn = 7-Be (2nd) inter. rate per year per kTon x <Z/A>   *
  70 *                   for elastic vue-electron interactions              *
  71 *          Pepcon = pep   interaction rate per year per kTon x <Z/A>   *
  72 *                   for elastic vue-electron interactions              *
  73 *       Spsnot(i) = differential spectrum x cross section for oscilla- *
  74 *                   ted numu-electron elastic interactions             *
  75 *       Cusnot(i) = cumulative spectrum x cross section for oscillated *
  76 *                   numu-electron elastic interactions                 *
  77 *          Snotcn = total interaction rate per year per kTon x <Z/A>   *
  78 *                   for elastic vue-electron interactions              *
  79 *          Ppocon = pp    interaction rate per year per kTon x <Z/A>   *
  80 *                   for elastic oscillated vumu-electron interactions  *
  81 *          B8ocon = 8-B   interaction rate per year per kTon x <Z/A>   *
  82 *                   for elastic oscillated vumu-electron interactions  *
  83 *          Hepocn = hep   interaction rate per year per kTon x <Z/A>   *
  84 *                   for elastic oscillated vumu-electron interactions  *
  85 *          Be71oc = 7-Be (1st) inter. rate per year per kTon x <Z/A>   *
  86 *                   for elastic oscillated vumu-electron interactions  *
  87 *          Be72oc = 7-Be (2nd) inter. rate per year per kTon x <Z/A>   *
  88 *                   for elastic oscillated vumu-electron interactions  *
  89 *          Pepocn = pep   interaction rate per year per kTon x <Z/A>   *
  90 *                   for elastic oscillated vumu-electron interactions  *
  91 *          Tesnmn = secondary electron minimum kinetic energy (GeV)    *
  92 *       Elv40a(k) = kth energy level of 40-K (produced by vue absorp-  *
  93 *                   tion on 40-Ar)                                     *
  94 *          Inuosc = oscillation option flag                            *
  95 *       Ior40a(j) = index of the jth level accessible for vue absorp-  *
  96 *                   tion on 40-Ar (levels are of 40-K) in the 40-K     *
  97 *                   level scheme (see also Elv40a)                     *
  98 *                                                                      *
  99 *----------------------------------------------------------------------*
 100 *
 101 *  Number of 40-Ar levels for which absorption cross section data are
 102 *  given:
 103       PARAMETER ( NFGTAR =   21 )
 104 *  Number of (neutrino) energy bins for Gamow-Teller and Fermi
 105 *  absorption cross sections
 106       PARAMETER ( NFGTEN =  200 )
 107 *  Index of the Fermi transition for 40-Ar
 108       PARAMETER ( IFR40A =    8 )
 109 *  Initial (0 index) energy for Gamow-Teller and Fermi
 110 *  absorption cross sections
 111       PARAMETER ( E0FRGT = ZERZER )
 112 *  Energy step for Gamow-Teller and Fermi
 113 *  absorption cross sections
 114       PARAMETER ( DEFRGT = 0.0001  D+00 )
 115 *  Recorded levels of 40-K (produced by vue absoprtion on 40-Ar)
 116       PARAMETER ( NLV40A =   40 )
 117 *  Maximum number of branching to other levels of 40-K (produced
 118 *  by vue absoprtion on 40-Ar)
 119       PARAMETER ( NBR40A =   10 )
 120 *  Number of energy bins for elastic neutrino scattering
 121       PARAMETER ( NFSNTT = 1880 )
 122 *  Energy step for elastic neutrino scattering
 123       PARAMETER ( DEFLSN = 0.00001 D+00 )
 124 *  Energy of 1st 7-Be neutrino line
 125       PARAMETER ( ENBE71 = 0.000384D+00 )
 126 *  Energy of 2nd 7-Be neutrino line
 127       PARAMETER ( ENBE72 = 0.000862D+00 )
 128 *  Energy of pep neutrino line
 129       PARAMETER ( ENLPEP = 0.001442D+00 )
 130 *
 131       COMMON / SLNUCM / SLNASG (NFGTAR,NFGTEN), BRT40A (NBR40A,NLV40A),
 132      &                  SPSNGT (0:NFGTEN), CUSNGT (0:NFGTEN),   SNGTCN,
 133      &                  PPCOGT, B8COGT, HEPCOG, BE71CG, BE72CG, PEPCOG,
 134      &                  SPSNFR (0:NFGTEN), CUSNFR (0:NFGTEN),   SNFRCN,
 135      &                  PPCOFR, B8COFR, HEPCOF, BE71CF, BE72CF, PEPCOF,
 136      &                  SPSNTT (0:NFSNTT), CUSNTT (0:NFSNTT),   SNTTCN,
 137      &                  PPCONT, B8CONT, HEPCON, BE71CN, BE72CN, PEPCON,
 138      &                  SPSNOT (0:NFSNTT), CUSNOT (0:NFSNTT),   SNOTCN,
 139      &                  PPOCON, B8OCON, HEPOCN, BE71OC, BE72OC, PEPOCN,
 140      &                  ELV40A (NLV40A), QUGTFR, TESNMN, INUOSC,
 141      &                  IBG40A (NBR40A,NLV40A), IOR40A (NFGTAR)
 142