GEANT321/fluka/energy.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1999/05/18 15:55:15  fca
   6 * AliRoot sources
   7 *
   8 * Revision 1.1.1.1  1995/10/24 10:19:55  cernlib
   9 * Geant
  10 *
  11 *
  12 #include "geant321/pilot.h"
  13 #if defined(CERNLIB_OLDNAME)
  14 *CMZ :  3.21/02 29/03/94  15.41.42  by  S.Giani
  15 *-- Author :
  16 *=== energy ===========================================================*
  17 *                                                                      *
  18       FUNCTION ENERGY(A,Z)
  19
  20 #include "geant321/dblprc.inc"
  21 #include "geant321/dimpar.inc"
  22 #include "geant321/iounit.inc"
  23 C--------------------------------------------------------------------
  24 C SUBNAME = ENERGY
  25 C--------------------------------------------------------------------
  26 #include "geant321/eva0.inc"
  27 C--------------------------------------------------------------------
  28       KA0 = NINT ( A )
  29       KZ0 = NINT ( Z )
  30       N   = KA0 - KZ0
  31       IF ( N .LE. 0 ) THEN
  32          IF ( KA0 .NE. 1 ) THEN
  33             IF ( N .LT. 0 ) THEN
  34                WRITE ( LUNOUT, * )
  35      &         ' Stopped in energy: mass number =< atomic number !!',
  36      &           KA0, KZ0
  37             END IF
  38          ELSE
  39             ENERGY = WAPS ( 1, 3 )
  40             RETURN
  41          END IF
  42       END IF
  43       IZZ = NINT ( WAPS ( KA0, 1 ) )
  44       IF ( KZ0 .LT. IZZ .OR. KZ0 .GT. IZZ + 9 ) THEN
  45          ENERGY = ENRG ( A, Z )
  46          RETURN
  47       ELSE
  48          IZ0 = KZ0 - IZZ + 2
  49          ENERGY = WAPS ( KA0, IZ0 )
  50          IF ( ENERGY .EQ. 0.D+00 .AND. (KA0 .NE. 12 .OR. KZ0 .NE. 6) )
  51      &        ENERGY = ENRG ( A, Z )
  52          RETURN
  53       END IF
  54 *     RETURN
  55       END
  56 #else
  57       SUBROUTINE ENERGY_DUMMY
  58       END
  59 #endif