Télécharger inter3.eso

Retour à la liste

Numérotation des lignes :

inter3
  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.  
  16. -INC PPARAM
  17. -INC CCOPTIO
  18. LOGICAL EGALDP
  19. MLREEL=KTE
  20. SEGACT MLREEL
  21. LON=PROG(/1)
  22. DO 40 L2=2,LON
  23. T2=PROG(L2)
  24. IF(T0.LE.T2) GOTO 60
  25. 40 CONTINUE
  26. L2=LON
  27. 60 CONTINUE
  28. L1=L2-1
  29. T1=PROG(L1)
  30. *
  31. *--- 1.D-15 A LA PLACE DE 1.D-5
  32. *
  33. IF (EGALDP(T1,T2,1.D-15)) THEN
  34. DT=0.D0
  35. ELSE
  36. DT=(T0-T1)/(T2-T1)
  37. ENDIF
  38. SEGDES MLREEL
  39. MLREEL=KTF
  40. SEGACT MLREEL
  41. FT1=PROG(L1)
  42. FT2=PROG(L2)
  43. SEGDES MLREEL
  44. FT0=(FT2-FT1)*DT+FT1
  45. RETURN
  46. END
  47.  
  48.  

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