Esope source files

   1 : C VALEUR    SOURCE    CHAT      05/01/13    04:00:32     5004           
   2 :       SUBROUTINE VALEUR(X,CURVE,NCURV,XLAC,Y,YPOI,KERRE)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)
   5 : C============================================================
   6 : C 29/06/95  D. C. Modif de la routine de Pierre
   7 : C  Routine qui lit la valeur d'une courbe et de sa derivee 
   8 : C  quand on lui donne l'abscisse
   9 : C
  10 : C  ENTREES:
  11 : C    NCURV : Nombre de points de la courbe
  12 : C    X     : Abscisse
  13 : C    CURVE : Courbe
  14 : C    XLAC  : Longueur caracteristique
  15 : C     
  16 : C SORTIES
  17 : C    Y     : Valeur de la fonction en X
  18 : C    YPOI : Valeur de la derivee en X
  19 : C=============================================================
  20 :       DIMENSION CURVE(2,NCURV)
  21 : C
  22 :       IF (X.LT.CURVE(1,1)) THEN
  23 :         YPOI = 0.D0 
  24 :         Y=1.D0
  25 :       ELSE
  26 :         XX = CURVE(1,2) + (1.D0/XLAC)*(X - CURVE(1,2))
  27 : C        XX = (1.D0/XLAC)*X 
  28 :         IF(XX.GE.CURVE(1,NCURV))THEN
  29 :           Y=CURVE(2,NCURV)
  30 :           YPOI = 0.D0
  31 :         ELSE
  32 :           DO IE1=3,NCURV
  33 :             IF (XX.LE.CURVE(1,IE1)) GOTO 1
  34 :           ENDDO
  35 :  1        YPOR=(CURVE(2,IE1)-CURVE(2,IE1-1))
  36 :      &        /(CURVE(1,IE1)-CURVE(1,IE1-1))
  37 :           Y=CURVE(2,IE1)+YPOR*(XX-CURVE(1,IE1))
  38 :           YPOI=(1.D0/XLAC)*YPOR
  39 :         ENDIF
  40 :       ENDIF
  41 : C
  42 :       RETURN
  43 :       END
  44 : 
  45 :  
  46 :