C LIRANG    SOURCE    CHAT      05/01/13    01:20:21     5004
      SUBROUTINE LIRANG(XTOT,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
*--------------------------------------------------------------------*
*                                                                    *
*     Operateur DYNE : algorithme de Fu - de Vogelaere               *
*     ________________________________________________               *
*                                                                    *
*     Repere l'indice NRG du tableau XABSCI representant la liste de *
*     reels des abscisses de la loi de comportement elasto-plastique.*
*     Pour une valeur XTOT de l'abscisse, cet indice sera tel que :  *
*              XABSCI(I,NRG)<= XTOT <= XABSCI(I,NRG+1)               *
*                                                                    *
*     Parametres:                                                    *
*                                                                    *
* e   XTOT    Valeur donnee de l'abscisse                            *
* e   XABSCI  Tableau contenant les abscisses de la loi plastique    *
*             pour la liaison point-point-plastique                  *
* e   I       numero de la liaison.                                  *
*             pour la liaison point-point-plastique                  *
* s   NRG     Rang du point repere                                   *
*                                                                    *                                                                          *
*     Auteur, date de creation:                                      *
*                                                                    *
*     Nicolas WECXSTEEN modifications 04/96 point-point-..-plastique *
*                                                                    *
*--------------------------------------------------------------------*
*
      REAL*8 XABSCI(NLIAB,*),XORDON(NLIAB,*)
*


      NRG = 1
      XVAL = XABSCI(I,2)
 10   IF (XTOT.GE.XVAL) THEN
        NRG = NRG + 1
        if (nrg.eq.NIP) then
          XPENTE = (XORDON(I,NRG)-XORDON(I,NRG-1))/
     &               (XABSCI(I,NRG)-XABSCI(I,NRG-1))
          go to 20
        endif

        XVAL = XABSCI(I,NRG+1)
        GOTO 10
      endif

        XPENTE = (XORDON(I,NRG+1)-XORDON(I,NRG))/
     &               (XABSCI(I,NRG+1)-XABSCI(I,NRG))

 20    continue

      END