Télécharger precn0.eso

Retour à la liste

Numérotation des lignes :

  1. C PRECN0 SOURCE CHAT 05/01/13 02:25:15 5004
  2. SUBROUTINE PRECN0(TEMPS,TE,FT,LON,FT0,IRET)
  3. C
  4. C INTERPOLATION LINEAIRE
  5. C IRET=POSITION DANS LE TABLEAU SI SUCCES , 0 SINON
  6. C
  7. IMPLICIT INTEGER(I-N)
  8. IMPLICIT REAL*8(A-H,O-Z)
  9. -INC CCOPTIO
  10. DIMENSION TE(*),FT(*)
  11. FT0=0.D0
  12. T2=TE(1)
  13. L2=1
  14. DT=TE(2)-T2
  15. IF(TEMPS.LT.T2) THEN
  16. IF(TEMPS.LT.T2-1.D-5*DT) GOTO 50
  17. FT0=FT(L2)
  18. IRET=L2
  19. RETURN
  20. ENDIF
  21. DO 40 L=2,LON
  22. L1=L2
  23. T2=TE(L)
  24. L2=L
  25. IF(T1.GT.T2) GOTO 55
  26. IF(TEMPS.LT.T2) GOTO 60
  27. 40 CONTINUE
  28. C TEMPS EN DEHORS DE LA TABLE
  29. DT=T2-T1
  30. IF(TEMPS.LT.T2+1.D-5*DT) THEN
  31. FT0=FT(L2)
  32. IRET=L2
  33. RETURN
  34. ENDIF
  35. C
  36. 50 IRET=-210
  37. RETURN
  38. C TEMPS NON CROISSANTS DANS LA TABLE
  39. 55 IRET=-211
  40. RETURN
  41. 60 DT=(TEMPS-T1)/(T2-T1)
  42. FT1=FT(L1)
  43. FT2=FT(L2)
  44. FT0=(FT2-FT1)*DT+FT1
  45. IRET=L1
  46. RETURN
  47. END
  48.  
  49.  
  50.  

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