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