GEANT321/fluka/energy.F

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