GEANT321/fluka/eexlvl.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 *CMZ :  3.21/02 29/03/94  15.41.42  by  S.Giani
  11 *-- Author :
  12 *$ CREATE EEXLVL.FOR
  13 *COPY EEXLVL
  14 *
  15 *=== eexlvl ===========================================================*
  16 *
  17       SUBROUTINE EEXLVL ( JA, JZ, EEX1ST, EEX2ND, EEXCON )
  18
  19 #include "geant321/dblprc.inc"
  20 #include "geant321/dimpar.inc"
  21 #include "geant321/iounit.inc"
  22 *
  23 *----------------------------------------------------------------------*
  24 *                                                                      *
  25 *     Created on 06 december 1991  by    Alfredo Ferrari & Paola Sala  *
  26 *                                                   Infn - Milan       *
  27 *                                                                      *
  28 *     Last change on 28-apr-92     by    Alfredo Ferrari               *
  29 *                                                                      *
  30 *                                                                      *
  31 *----------------------------------------------------------------------*
  32 *
  33 #include "geant321/eva0.inc"
  34 *
  35       SAVE SQATAR, JAOLD
  36       DATA JAOLD / 0 /
  37 *
  38       IF ( JA .EQ. JZ ) THEN
  39          EEXCON = 0.D+00
  40          EEX1ST = 0.D+00
  41          EEX2ND = 0.D+00
  42          RETURN
  43       END IF
  44       IF ( JZ .EQ. 0 .OR. JA .EQ. JZ ) THEN
  45          EEXCON = 0.D+00
  46       ELSE
  47          EEXCON = 1.D-03 * ( CAM4 (JZ) + CAM5 (JA-JZ) )
  48       END IF
  49 *  **** Very tentative selection of Eex1st, Eex2nd, based on pairing
  50 *       energies according to delta = 12 MeV / A^1/2 ****
  51       IZODD  = 1 - MOD ( JZ, 2 )
  52       INODD  = 1 - MOD ( JA-JZ, 2 )
  53       IODD   = IZODD + INODD
  54 *  +-------------------------------------------------------------------*
  55 *  |  Even-even nucleus
  56       IF ( IODD .GE. 2 ) THEN
  57 *  |  +----------------------------------------------------------------*
  58 *  |  |
  59          IF ( JA .NE. JAOLD ) THEN
  60             JAOLD  = JA
  61             DJA = JA
  62             SQATAR = SQRT ( DJA )
  63          END IF
  64 *  |  |
  65 *  |  +----------------------------------------------------------------*
  66          EEX1ST = 12.0D-03 / SQATAR
  67          EEX2ND = 2.D+00 * EEX1ST
  68 *  |
  69 *  +-------------------------------------------------------------------*
  70 *  |  even-odd nucleus
  71       ELSE IF ( IODD .GT. 0 ) THEN
  72 *  |  +----------------------------------------------------------------*
  73 *  |  |
  74          IF ( JA .NE. JAOLD ) THEN
  75             JAOLD  = JA
  76             DJA    = JA
  77             SQATAR = SQRT ( DJA )
  78          END IF
  79 *  |  |
  80 *  |  +----------------------------------------------------------------*
  81          EEX1ST = 12.0D-03 / SQATAR
  82          EEX2ND = EEX1ST
  83 *  |
  84 *  +-------------------------------------------------------------------*
  85 *  |  odd-odd nucleus
  86       ELSE
  87          EEX2ND = 0.D+00
  88          EEX1ST = 0.D+00
  89       END IF
  90 *  |
  91 *  +-------------------------------------------------------------------*
  92       RETURN
  93 *=== End of subroutine eexlvl =========================================*
  94       END