Esope source files

   1 : C MAXINT    SOURCE    CHAT      05/01/13    01:36:56     5004           
   2 :        SUBROUTINE MAXINT(KET,IDIM,TEMPS,FTEMPS,IRET)
   3 : *
   4 : *
   5 : *
   6 : * INTERPOLATION LINEAIRE DU MODELE DE MAXWELL
   7 : * ----------------------------------------------------------
   8 : *
   9 : * entree
  10 : *        KET tableau de reels de dimension 2*IDIM
  11 : *        IDIM nombre de points definissant la courbe d'evolution
  12 : *        TEMPS valeur de l'abscisse a interpoler
  13 : *
  14 : * sortie
  15 : *        FTEMPS ordonnee interpolee correspondant a TEMPS
  16 : *        IRET retourne 
  17 : *                0 si il y a deux valeurs successives en abscisse
  18 : *                       consecutives
  19 : *                1 si l'interpolation est reussie
  20 : *
  21 : * -----------------------------------------------------------
  22 : *
  23 : *
  24 :       IMPLICIT INTEGER(I-N)
  25 :        IMPLICIT REAL*8(A-H,O-Z)
  26 :        REAL*8 KET(2,IDIM)
  27 : *
  28 :        DO 10 I=1,IDIM-1
  29 :          IF (KET(1,I+1).EQ.KET(1,I)) THEN
  30 :            IRET=0
  31 :            RETURN
  32 :          ENDIF
  33 :          IF (TEMPS.LT.KET(1,I+1)) THEN
  34 :           FTEMPS=(TEMPS-KET(1,I))/(KET(1,I+1)-KET(1,I))
  35 :      1         *(KET(2,I+1)-KET(2,I))+KET(2,I)
  36 :            IRET=1
  37 :            RETURN
  38 :          ENDIF
  39 : 10     CONTINUE
  40 :       END
  41 :  
  42 :