1 : C INTER3 SOURCE CHAT 05/01/13 00:40:22 5004 2 : SUBROUTINE INTER3(KTE,KTF,T0,FT0) 3 : C 4 : C ========================================================= 5 : C = INTERPOLATION LINEAIRE = 6 : C = KTE : POINTEUR SUR OBJET DE TYPE LISTREEL DE T = 7 : C = KTF : POINTEUR SUR OBJET DE TYPE LISTREEL DE F(T) = 8 : C = CREATION : 22/02/85 = 9 : C = PROGRAMMEUR : GUILBAUD = 10 : C ========================================================= 11 : C 12 : IMPLICIT INTEGER(I-N) 13 : IMPLICIT REAL*8(A-H,O-Z) 14 : -INC SMLREEL 15 : -INC CCOPTIO 16 : LOGICAL EGALDP 17 : MLREEL=KTE 18 : SEGACT MLREEL 19 : LON=PROG(/1) 20 : DO 40 L2=2,LON 21 : T2=PROG(L2) 22 : IF(T0.LE.T2) GOTO 60 23 : 40 CONTINUE 24 : L2=LON 25 : 60 CONTINUE 26 : L1=L2-1 27 : T1=PROG(L1) 28 : * 29 : *--- 1.D-15 A LA PLACE DE 1.D-5 30 : * 31 : IF (EGALDP(T1,T2,1.D-15)) THEN 32 : DT=0.D0 33 : ELSE 34 : DT=(T0-T1)/(T2-T1) 35 : ENDIF 36 : SEGDES MLREEL 37 : MLREEL=KTF 38 : SEGACT MLREEL 39 : FT1=PROG(L1) 40 : FT2=PROG(L2) 41 : SEGDES MLREEL 42 : FT0=(FT2-FT1)*DT+FT1 43 : RETURN 44 : END 45 :
© Cast3M 2003 - All rights reserved.
Disclaimer