GEANT321/fluka/qnrg.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:19:58  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.44  by  S.Giani
  11 *-- Author :
  12 *$ CREATE QNRG.FOR
  13 *COPY QNRG
  14 *                                                                      *
  15 *=== qnrg =============================================================*
  16 *                                                                      *
  17       FUNCTION QNRG ( A1, Z1, A2, Z2 )
  18
  19 #include "geant321/dblprc.inc"
  20 #include "geant321/dimpar.inc"
  21 #include "geant321/iounit.inc"
  22 C--------------------------------------------------------------------
  23 C SUBNAME = QNRG
  24 C--------------------------------------------------------------------
  25 #include "geant321/eva0.inc"
  26 C--------------------------------------------------------------------
  27       KA1  = NINT (A1)
  28       KZ1  = NINT (Z1)
  29       KA2  = NINT (A2)
  30       KZ2  = NINT (Z2)
  31       N1= KA1 - KZ1
  32       N2= KA2 - KZ2
  33       IF (N1 .LE. 0) STOP 105
  34       IF (N2 .LE. 0) STOP 106
  35       IZZ1 = NINT ( WAPS (KA1,1) )
  36       IZZ2 = NINT ( WAPS (KA2,1) )
  37       IF ( KZ1 .LT. IZZ1 .OR. KZ1 .GT. IZZ1 + 9 ) THEN
  38          ENRG1 = ENRG ( A1, Z1 )
  39       ELSE
  40          IZ1   = KZ1 - IZZ1 + 2
  41          ENRG1 = WAPS ( KA1, IZ1 )
  42          IF ( ENRG1 .EQ. 0.D0 .AND. ( KA1 .NE. 12 .OR. KZ1 .NE. 6) )
  43      &        THEN
  44             ENRG1 = ENRG ( A1, Z1 )
  45          END IF
  46       END IF
  47       IF ( KZ2 .LT. IZZ2 .OR. KZ2 .GT. IZZ2 + 9 ) THEN
  48          ENRG2 = ENRG ( A2, Z2 )
  49       ELSE
  50          IZ2   = KZ2 - IZZ2 + 2
  51          ENRG2 = WAPS ( KA2, IZ2 )
  52          IF ( ENRG2 .EQ. 0.D0 .AND. ( KA2 .NE. 12 .OR. KZ2 .NE. 6) )
  53      &        THEN
  54             ENRG2 = ENRG ( A2, Z2 )
  55          END IF
  56       END IF
  57       QNRG = ENRG1 - ENRG2
  58       RETURN
  59       END