1 : C INTLIN SOURCE CHAT 05/01/13 00:41:24 5004
2 : SUBROUTINE INTLIN(TEMPS,KTE,KFT,LON,L1,L2, FT0,IRET)
3 : IMPLICIT INTEGER(I-N)
4 : IMPLICIT REAL*8(A-H,O-Z)
5 : *--------------------------------------------------------------------*
6 : * *
7 : * INTERPOLATION LINEAIRE D'UNE VALEUR *
8 : * *
9 : * Param}tres: *
10 : * *
11 : * e TEMPS Abscisse @ interpoler *
12 : * e KTE LISTREEL des abscisses *
13 : * e KFT LISTREEL des ordonn{es *
14 : * e LON Longueur de la LISTREEL des abscisses *
15 : * s FT0 Valeur interpol{e *
16 : * s IRET IRET = 1, succ{s *
17 : * s IRET IRET = 0, {chec *
18 : * *
19 : *--------------------------------------------------------------------*
20 : -INC SMLREEL
21 : FT0 = 0.D0
22 : MLREE1 = KTE
23 : MLREE2 = KFT
24 : T1 = MLREE1.PROG(L1)
25 : TL = MLREE1.PROG(LON)
26 : IF (TEMPS.LT.T1 .OR. TEMPS.GT.TL) THEN
27 : * temps en dehors de la liste
28 : IRET = 0
29 : RETURN
30 : ENDIF
31 : DO 10 L = L2,LON
32 : T2 = MLREE1.PROG(L)
33 : L2 = L
34 : IF (TEMPS.LT.T2) GOTO 20
35 : T1 = T2
36 : L1 = L2
37 : 10 CONTINUE
38 : * end do
39 : 20 FT1 = MLREE2.PROG(L1)
40 : FT2 = MLREE2.PROG(L2)
41 : DT = (TEMPS - T1) / (T2 - T1)
42 : FT0 = (FT2 - FT1) * DT + FT1
43 : IRET = 1
44 : *
45 : END
46 :
© Cast3M 2003 - All rights reserved.
Disclaimer