GEANT321/neutron/tbsplt.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:59  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 TBSPLT(A,E,NP,Y)
  13 C       THIS ROUTINE DETERMINES A CROSS SECTION AT A GIVEN
  14 C       ENERGY FROM A CROSS SECTION VERSUS ENERGY TABLE USING
  15 C       A TABLE SPLITTING METHOD
  16       DIMENSION A(1)
  17       SAVE
  18       IPP=1
  19       IF(E.LE.A(1))GO TO 40
  20       IF(E.GE.A(2*NP-1))GO TO 50
  21       INDXH=NP
  22       INDXL=0
  23    10 IF(INDXL+1.EQ.INDXH)GO TO 30
  24       J=(INDXH+INDXL)/2
  25       N=2*J-1
  26       IF(E.LE.A(N))GO TO 20
  27       INDXL=J
  28       GO TO 10
  29    20 INDXH=J
  30       GO TO 10
  31    30 N=2*INDXH-1
  32       Y=A(N-1)+(E-A(N-2))*(A(N+1)-A(N-1))/(A(N)-A(N-2))
  33       RETURN
  34    40 Y=A(IPP+1)
  35       RETURN
  36    50 Y=A(2*NP)
  37       RETURN
  38       END