Esope source files

   1 : C INTE22    SOURCE    CHAT      05/01/13    00:39:53     5004           
   2 :       SUBROUTINE INTE22(PRT,PRF,T0,FT0,LON)                                     
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 : C                                                                               
   6 : C  ====================================================================         
   7 : C  =         INTERPOLATION LINEAIRE    RECHERCHE PAR DICHOTOMIE       =         
   8 : C  = CREATION    14/01/86                                             =         
   9 : C  = PROGRAMMEUR GUILBAUD                                             =         
  10 : C  ====================================================================         
  11 : C                                                                               
  12 : -INC CCOPTIO                                                                    
  13 :       DIMENSION PRT(1),PRF(2)                                                   
  14 :       DATA EPSI/1.D-10/                                                         
  15 :       N1=1                                                                      
  16 :       N2=LON                                                                    
  17 :     1 CONTINUE                                                                  
  18 :       N2N1=N2-N1                                                                
  19 :       IF(N2N1.EQ.1) GOTO 2                                                      
  20 :       NC=(N1+N2)/2                                                              
  21 :       IF(T0.LT.PRT(NC)) THEN                                                    
  22 :         N2=NC                                                                   
  23 :       ELSE                                                                      
  24 :         N1=NC                                                                   
  25 :       ENDIF                                                                     
  26 :       GOTO 1                                                                    
  27 :     2 CONTINUE                                                                  
  28 :       ALPHA=(T0-PRT(N1))/(PRT(N2)-PRT(N1))                                      
  29 :       FT0=ALPHA*(PRF(N2)-PRF(N1))+PRF(N1)                                       
  30 :       BETA=1.D0-ALPHA                                                           
  31 :       IF(BETA.LT.EPSI)FT0=PRF(N2)                                               
  32 :       RETURN                                                                    
  33 :       END                                                                       
  34 :