Télécharger intere.eso

Retour à la liste

Numérotation des lignes :

  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.  
  55.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales