GEANT321/neutron/xsecnu.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:22:00  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.48  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE XSECNU(BUF,LEN,E,XSC,L1,L2)
  13 C       THIS ROUTINE DETERMINES A CROSS SECTION AT A GIVEN ENERGY
  14 C       FROM A CROSS SECTION VERSUS ENERGY TABLE
  15       DIMENSION BUF(*)
  16       SAVE
  17       IF(E.LT.BUF(L1))GO TO 40
  18       DO 10 J=1,LEN
  19          N=L1+2*(J-1)
  20          IF(E.LE.BUF(N))GO TO 20
  21    10 CONTINUE
  22       XSC=BUF(L2)
  23       RETURN
  24    20 IF(J.EQ.1)GO TO 30
  25       XSC=BUF(N-1)+(E-BUF(N-2))*(BUF(N+1)-BUF(N-1))/
  26      +(BUF(N)-BUF(N-2))
  27       RETURN
  28    30 XSC=BUF(N+1)
  29       RETURN
  30    40 XSC=0.0
  31       RETURN
  32       END