Sources Esope | Cast3M

valeur
C VALEUR    SOURCE    CHAT      05/01/13    04:00:32     5004
      SUBROUTINE VALEUR(X,CURVE,NCURV,XLAC,Y,YPOI,KERRE)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C============================================================
C 29/06/95  D. C. Modif de la routine de Pierre
C  Routine qui lit la valeur d'une courbe et de sa derivee
C  quand on lui donne l'abscisse
C
C  ENTREES:
C    NCURV : Nombre de points de la courbe
C    X     : Abscisse
C    CURVE : Courbe
C    XLAC  : Longueur caracteristique
C
C SORTIES
C    Y     : Valeur de la fonction en X
C    YPOI : Valeur de la derivee en X
C=============================================================
      DIMENSION CURVE(2,NCURV)
C
      IF (X.LT.CURVE(1,1)) THEN
        YPOI = 0.D0
        Y=1.D0
      ELSE
        XX = CURVE(1,2) + (1.D0/XLAC)*(X - CURVE(1,2))
C        XX = (1.D0/XLAC)*X
        IF(XX.GE.CURVE(1,NCURV))THEN
          Y=CURVE(2,NCURV)
          YPOI = 0.D0
        ELSE
          DO IE1=3,NCURV
            IF (XX.LE.CURVE(1,IE1)) GOTO 1
          ENDDO
 1        YPOR=(CURVE(2,IE1)-CURVE(2,IE1-1))
     &        /(CURVE(1,IE1)-CURVE(1,IE1-1))
          Y=CURVE(2,IE1)+YPOR*(XX-CURVE(1,IE1))
          YPOI=(1.D0/XLAC)*YPOR
        ENDIF
      ENDIF
C
      RETURN
      END