Sources Esope | Cast3M

varifo
C VARIFO    SOURCE    FANDEUR   16/09/21    21:15:16     9098           
 
      SUBROUTINE VARIFO (XX,YY,NB,JORDO,X1, Y1)
 
C-----------------------------------------------------------------------
C  ENTREES :
C     XX YY NB JORDO X1
C  SORTIE :
C     Y1
C-----------------------------------------------------------------------
C JORDO indique que le tableau de valeurs XX est ordonne de maniere
C strictement croissante (=+1) ou strictement decroissante (=-1).
C Les autres cas ne sont pas consideres ici (deja traites en erreur).
C On pourrait faire une recherche par dichotomie de IP, mais on va 
C supposer ici que la dimension des tableaux (NB) est "petite".
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
      DIMENSION XX(*),YY(*)
*
*   ON INTERPOLE LINEAIREMENT MEME EN DEHORS DES BORNES !!!
*
      IF (JORDO.EQ.+1) THEN
        DO I = 2, NB
          IP = I
          IF (XX(I).GT.X1) GOTO 10
        ENDDO
*
**    ELSE IF (JORDO.EQ.-1) THEN
      ELSE
        DO I = 2, NB
          IP = I
          IF (XX(I).LT.X1) GOTO 10
        ENDDO
      ENDIF
 
  10  CONTINUE
      IP1 = IP-1
      PENTE = (YY(IP)-YY(IP1)) / (XX(IP)-XX(IP1))
      Y1 = YY(IP1) + PENTE * (X1-XX(IP1))
 
      RETURN
      END