Esope source files

   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                                                    
  72 : *        CALL ERREUR(   )                                                       
  73 :          SEGDES MLREEL                                                    
  74 :          SEGDES MLREE1                                                    
  75 :          RETURN                                                                 
  76 :       ENDIF                                                                     
  77 :       JG = LON1                                                                 
  78 :       SEGINI MLREE2                                                             
  79 :       IPT3 = MLREE2                                                             
  80 : *                                                                               
  81 : *     INTERPOLATION LINEAIRE                                                    
  82 : *                                                                               
  83 :       DO 10 I = 1,LON1                                                          
  84 :          X1 = PROG(I)                                                           
  85 :          X2 = MLREE1.PROG(I)                                                    
  86 :          MLREE2.PROG(I) = ( ALPHA * (X2 - X1) ) + X1                            
  87 : 10       CONTINUE                                                               
  88 : *     END DO                                                                    
  89 : *                                                                               
  90 :       SEGDES MLREEL                                                       
  91 :       SEGDES MLREE1                                                       
  92 :       SEGDES MLREE2                                                             
  93 : *                                                                               
  94 :       END                                                                       
  95 :