DPMJET/flukapro/(ATNUBF)

   1 *$ CREATE ATNUBF.ADD
   2 *COPY ATNUBF
   3 *
   4 *=== atnubf ===========================================================*
   5 *
   6 *----------------------------------------------------------------------*
   7 *                                                                      *
   8 *     ATmospheric NeUtrino BuFfer:                                     *
   9 *                                                                      *
  10 *     Created on   29 may 1996     by    Alfredo Ferrari & Paola Sala  *
  11 *                                                   Infn - Milan       *
  12 *                                                                      *
  13 *     Last change on 15-dec-99     by    Alfredo Ferrari               *
  14 *                                                                      *
  15 *     Description of variables:                                        *
  16 *                                                                      *
  17 *          Idatnu = neutrino id (Paprop numbering)                     *
  18 *          Lgatnu = neutrino generation                                *
  19 *          Enatnu = neutrino energy (GeV)                              *
  20 *          Diatnu = neutrino production height or distance (cm)        *
  21 *          Thatnu = neutrino direction polar   (theta) angle (rad)     *
  22 *          Phatnu = neutrino direction azimuthal (phi) angle (rad)     *
  23 *          Wtatnu = neutrino weight (such to be automatically norm-    *
  24 *                   alized to fluence per unit time and area)          *
  25 *          Ipatnu = parent cosmic ray id (Z + 100 x A)                 *
  26 *          Ifatnu = "father" hadron/muon Fluka id                      *
  27 *          Igatnu = "grandfather" hadron Fluka id                      *
  28 *          Ictfnu = i1 + j1 x 3 + i2 x 9 + j2 x 27 + i3 x 81 + j3 x 243*
  29 *                 + i4 x 729 + j4 x 2187 + i5 x 6561 + j5 x 19683      *
  30 *                   ik = cutoff flag for the k_th location, direct nu  *
  31 *                   jk = cutoff flag for the k_th location, mirror nu  *
  32 *                           = 0 <-> not yet checked                    *
  33 *                           = 1 <-> neutrino not cutoffed              *
  34 *                           = 2 <-> neutrino cutoffed                  *
  35 *          Pmatnu = parent cosmic ray momentum (GeV/c/amu)             *
  36 *          Pfatnu = "father" hadron/muon momentum (GeV/c)              *
  37 *          Pgatnu = "grandfather" hadron momentum (GeV/c)              *
  38 *          Xpatnu = parent cosmic ray 1st interaction x coord. (cm)    *
  39 *          Ypatnu = parent cosmic ray 1st interaction y coord. (cm)    *
  40 *          Zpatnu = parent cosmic ray 1st interaction z coord. (cm)    *
  41 *          Tpatnu = parent cosmic ray direction polar     angle (rad)  *
  42 *          Ppatnu = parent cosmic ray direction azimuthal angle (rad)  *
  43 *          Wpatnu = accumulated primary weight at the previous buffer  *
  44 *                   flush                                              *
  45 *          Npatnu = accumulated primary number at the previous buffer  *
  46 *                   flush                                              *
  47 *          Ncatnu = current pointer in the buffer                      *
  48 *          Lbatnu = logical flag for atmospheric neutrino buffering    *
  49 *          Lunatn = logical unit for the atmospheric neutrino file     *
  50 *                                                                      *
  51 *   ALL VARIABLES ARE MEANT IN THE FRAME WHERE THE NEUTRINO POSITION   *
  52 *   IS ALONG (0,0,1), THAT IS Z IS THE LOCAL ZENITH AXIS, X IS POIN-   *
  53 *   TING NORTH, AND Y IS POINTING WEST                                 *
  54 *                                                                      *
  55 *----------------------------------------------------------------------*
  56 *
  57       PARAMETER ( MXATNU = 2000 )
  58       LOGICAL LBATNU
  59 *
  60       COMMON / ATNUBF / ENATNU (MXATNU), DIATNU (MXATNU),
  61      &                  THATNU (MXATNU), PHATNU (MXATNU),
  62      &                  WTATNU (MXATNU), PMATNU (MXATNU),
  63      &                  PFATNU (MXATNU), PGATNU (MXATNU),
  64      &                  XPATNU (MXATNU), YPATNU (MXATNU),
  65      &                  ZPATNU (MXATNU), TPATNU (MXATNU),
  66      &                  PPATNU (MXATNU), WPATNU,
  67      &                  IDATNU (MXATNU), LGATNU (MXATNU),
  68      &                  IPATNU (MXATNU), IFATNU (MXATNU),
  69      &                  IGATNU (MXATNU), ICTFNU (MXATNU),
  70      &                  NPATNU, NCATNU, LBATNU, LUNATN