Télécharger la plaquette Cast3M


fr en

Télécharger inte44.eso

Retour à la liste

Numérotation des lignes : 

inte44

  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