Sources Esope | Cast3M

maxint
C MAXINT    SOURCE    CHAT      05/01/13    01:36:56     5004
       SUBROUTINE MAXINT(KET,IDIM,TEMPS,FTEMPS,IRET)
*
*
*
* INTERPOLATION LINEAIRE DU MODELE DE MAXWELL
* ----------------------------------------------------------
*
* entree
*        KET tableau de reels de dimension 2*IDIM
*        IDIM nombre de points definissant la courbe d'evolution
*        TEMPS valeur de l'abscisse a interpoler
*
* sortie
*        FTEMPS ordonnee interpolee correspondant a TEMPS
*        IRET retourne
*                0 si il y a deux valeurs successives en abscisse
*                       consecutives
*                1 si l'interpolation est reussie
*
* -----------------------------------------------------------
*
*
      IMPLICIT INTEGER(I-N)
       IMPLICIT REAL*8(A-H,O-Z)
       REAL*8 KET(2,IDIM)
*
       DO 10 I=1,IDIM-1
         IF (KET(1,I+1).EQ.KET(1,I)) THEN
           IRET=0
           RETURN
         ENDIF
         IF (TEMPS.LT.KET(1,I+1)) THEN
          FTEMPS=(TEMPS-KET(1,I))/(KET(1,I+1)-KET(1,I))
     1         *(KET(2,I+1)-KET(2,I))+KET(2,I)
           IRET=1
           RETURN
         ENDIF
10     CONTINUE
      END