Télécharger inter3.eso

Retour à la liste

Numérotation des lignes :

  1. C INTER3 SOURCE CHAT 05/01/13 00:40:22 5004
  2. SUBROUTINE INTER3(KTE,KTF,T0,FT0)
  3. C
  4. C =========================================================
  5. C = INTERPOLATION LINEAIRE =
  6. C = KTE : POINTEUR SUR OBJET DE TYPE LISTREEL DE T =
  7. C = KTF : POINTEUR SUR OBJET DE TYPE LISTREEL DE F(T) =
  8. C = CREATION : 22/02/85 =
  9. C = PROGRAMMEUR : GUILBAUD =
  10. C =========================================================
  11. C
  12. IMPLICIT INTEGER(I-N)
  13. IMPLICIT REAL*8(A-H,O-Z)
  14. -INC SMLREEL
  15. -INC CCOPTIO
  16. LOGICAL EGALDP
  17. MLREEL=KTE
  18. SEGACT MLREEL
  19. LON=PROG(/1)
  20. DO 40 L2=2,LON
  21. T2=PROG(L2)
  22. IF(T0.LE.T2) GOTO 60
  23. 40 CONTINUE
  24. L2=LON
  25. 60 CONTINUE
  26. L1=L2-1
  27. T1=PROG(L1)
  28. *
  29. *--- 1.D-15 A LA PLACE DE 1.D-5
  30. *
  31. IF (EGALDP(T1,T2,1.D-15)) THEN
  32. DT=0.D0
  33. ELSE
  34. DT=(T0-T1)/(T2-T1)
  35. ENDIF
  36. SEGDES MLREEL
  37. MLREEL=KTF
  38. SEGACT MLREEL
  39. FT1=PROG(L1)
  40. FT2=PROG(L2)
  41. SEGDES MLREEL
  42. FT0=(FT2-FT1)*DT+FT1
  43. RETURN
  44. END
  45.  
  46.  

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