Esope source files

   1 : C INTERE    SOURCE    CHAT      05/01/13    00:40:37     5004           
   2 : C  INTERE    SOURCE   NECH   99/9/12
   3 :       SUBROUTINE INTERE(T,TE,A,DF)
   4 : C
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 : C
   8 : *--------------------------------------------------------------------*
   9 : *                                                                    *
  10 : *          INTERPOLATION DES LOIS D'EVOLUTION                        * 
  11 : *                                                                    *
  12 : *--------------------------------------------------------------------*
  13 : C
  14 :       DIMENSION TE(6), A(6)
  15 : C
  16 :       
  17 :       IF(T.LE.TE(1))THEN
  18 :        DF=1.D0
  19 :       ELSE
  20 :       IF(T.LT.TE(2))THEN
  21 :        AA=(A(2)-A(1))/(TE(2)-TE(1))
  22 :        AB=(TE(2)*A(1)-TE(1)*A(2))/(TE(2)-TE(1))
  23 :        DF=AA*T+AB
  24 :       ELSE
  25 :       IF(T.LT.TE(3))THEN
  26 :        AA=(A(3)-A(2))/(TE(3)-TE(2))
  27 :        AB=(TE(3)*A(2)-TE(2)*A(3))/(TE(3)-TE(2))
  28 :        DF=AA*T+AB
  29 :       ELSE
  30 :       IF(T.LT.TE(4))THEN
  31 :        AA=(A(4)-A(3))/(TE(4)-TE(3))
  32 :        AB=(TE(4)*A(3)-TE(3)*A(4))/(TE(4)-TE(3))
  33 :        DF=AA*T+AB
  34 :       ELSE
  35 :       IF(T.LT.TE(5))THEN
  36 :        AA=(A(5)-A(4))/(TE(5)-TE(4))
  37 :        AB=(TE(5)*A(4)-TE(4)*A(5))/(TE(5)-TE(4))
  38 :        DF=AA*T+AB
  39 :       ELSE
  40 :        AA=(A(6)-A(5))/(TE(6)-TE(5))
  41 :        AB=(TE(6)*A(5)-TE(5)*A(6))/(TE(6)-TE(5))
  42 :        DF=AA*T+AB
  43 :       ENDIF
  44 :       ENDIF
  45 :       ENDIF
  46 :       ENDIF
  47 :       ENDIF        
  48 : C
  49 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  50 : C                       
  51 :         RETURN
  52 :         END 
  53 :  
  54 :