Esope source files

   1 : C VARIFO    SOURCE    CHAT      05/01/13    04:01:57     5004           
   2 :       SUBROUTINE VARIFO(XX,YY,NB,X1,Y1,KERRE)                                   
   3 : C----------------------------------------------------------------               
   4 : C  ENTREES:                                                                     
   5 : C     XX YY NB X1                                                               
   6 : C  SORTIE:                                                                      
   7 : C     Y1 KERRE                                                                  
   8 : C----------------------------------------------------------------               
   9 :       IMPLICIT INTEGER(I-N)
  10 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  11 :       DIMENSION XX(*),YY(*)                                                     
  12 :       KERRE=0                                                                   
  13 : *
  14 : *   ON VERIFIE QUE LA LISTE DES ABSCISSES EST ORDONNEE
  15 : *
  16 :       DO 1 I=2,NB
  17 :       IF(XX(I).LE.XX(I-1)) GO TO 2
  18 :   1   CONTINUE
  19 :       JORDO=1
  20 :       GO TO 10
  21 : *
  22 :   2   CONTINUE
  23 :       DO 3 I=2,NB
  24 :       IF(XX(I).GE.XX(I-1)) GO TO 4
  25 :   3   CONTINUE
  26 :       JORDO=-1
  27 :       GO TO 10
  28 : *
  29 :   4   CONTINUE
  30 :       KERRE=211                                                                 
  31 :       RETURN                                                                    
  32 : *
  33 : *   ON INTERPOLE LINEAIREMENT
  34 : *
  35 :   10  CONTINUE
  36 :       IF (JORDO.EQ.1) THEN
  37 :         DO 11 I=2,NB
  38 :         IP=I
  39 :         IF(XX(I).GT.X1) GO TO 13
  40 :   11    CONTINUE
  41 : *
  42 :       ELSE
  43 :         DO 12 I=2,NB
  44 :         IP=I
  45 :         IF(XX(I).LT.X1) GO TO 13
  46 :   12    CONTINUE
  47 :       ENDIF
  48 :   13  CONTINUE
  49 :       IP1=IP-1
  50 :       PENTE=(YY(IP)-YY(IP1))/(XX(IP)-XX(IP1))
  51 :       Y1=YY(IP1)+PENTE*(X1-XX(IP1))
  52 :       RETURN                                                                    
  53 :       END                                                                       
  54 :