Télécharger inte44.eso

Retour à la liste

Numérotation des lignes :

  1. C INTE44 SOURCE CHAT 05/01/13 00:40:04 5004
  2. SUBROUTINE INTE44(IPT1,IPT2,ALPHA, IPT3)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. ************************************************************************
  6. *
  7. * I N T E 4 4
  8. * -----------
  9. *
  10. * FONCTION:
  11. * ---------
  12. *
  13. * INTERPOLATION ENTRE DEUX LISTREELS.
  14. *
  15. * MODULES UTILISES:
  16. * -----------------
  17. *
  18. -INC SMLREEL
  19. *
  20. * PARAMETRES: (E)=ENTREE (S)=SORTIE (+ = CONTENU DANS UN COMMUN)
  21. * -----------
  22. *
  23. * IPT1 (E) POINTEUR SUR LE PREMIER LISTREEL.
  24. * IPT2 (E) POINTEUR SUR LE DEUXIEME LISTREEL.
  25. * ALPHA (E) COEFFICIENT D'INTERPOLATION.
  26. * IPT3 (S) POINTEUR SUR LE LISTREEL RESULTANT.
  27. *
  28. *
  29. * MODE DE FONCTIONNEMENT:
  30. * -----------------------
  31. * X1 : LA VALEUR DE IPT1 A LA POSITION I.
  32. * X2 : LA VALEUR DE IPT2 A LA POSITION I.
  33. * X3 : LA VALEUR DE IPT3 A LA POSITION I.
  34. *
  35. * ON A :
  36. *
  37. * X3 = ( (1 - ALPHA) * X1 ) + ( ALPHA * X2 )
  38. *
  39. *
  40. * REMARQUES:
  41. * ----------
  42. *
  43. * LES DEUX LISTREELS DOIVENT ETRE DE MEME LONGUEUR.
  44. *
  45. *
  46. * AUTEUR, DATE DE CREATION:
  47. * -------------------------
  48. *
  49. * LIONEL VIVAN SEPTEMBRE 1988
  50. *
  51. * LANGAGE:
  52. * --------
  53. *
  54. * ESOPE + FORTRAN77
  55. *
  56. ************************************************************************
  57. *
  58. * ORDONNEES NUMERO 1
  59. *
  60. MLREEL=IPT1
  61. SEGACT MLREEL
  62. *
  63. * ORDONNEES NUMERO 2
  64. *
  65. MLREE1=IPT2
  66. SEGACT MLREE1
  67. *
  68. * ORDONNEES RECHERCHEES
  69. *
  70. LON1 = PROG(/1)
  71. LON2 = MLREE1.PROG(/1)
  72. IF (LON1.NE.LON2) THEN
  73. * CALL ERREUR( )
  74. SEGDES MLREEL
  75. SEGDES MLREE1
  76. RETURN
  77. ENDIF
  78. JG = LON1
  79. SEGINI MLREE2
  80. IPT3 = MLREE2
  81. *
  82. * INTERPOLATION LINEAIRE
  83. *
  84. DO 10 I = 1,LON1
  85. X1 = PROG(I)
  86. X2 = MLREE1.PROG(I)
  87. MLREE2.PROG(I) = ( ALPHA * (X2 - X1) ) + X1
  88. 10 CONTINUE
  89. * END DO
  90. *
  91. SEGDES MLREEL
  92. SEGDES MLREE1
  93. SEGDES MLREE2
  94. *
  95. END
  96.  
  97.  

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