interp.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	interp (tf, tt, npoint, f, tg)

Function/Subroutine Documentation

interp()

subroutine interp	(		tf,
			tt,
		integer	npoint,
			f,
			tg )

Definition at line 34 of file interp.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NPOINT
      my_real tt, f, tg
      my_real tf(2,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, M
      my_real dt1, dt2, d1, d2
C-----------------------------------------------
      IF (tt-tf(1,npoint)>zero) GO TO 50
      DO 20 i=2,npoint
      m=i
      IF (tt-tf(1,m)<=zero) GO TO 30
   20 CONTINUE
   50 f=zero
      tg=zero
      RETURN
   30 dt1=tt-tf(1,m-1)
      dt2=-tt+tf(1,m)
      d1=tf(1,m)-tf(1,m-1)
      d2=tf(2,m)-tf(2,m-1)
      tg=zero
      IF(d1/=zero)tg=d2/d1
      !TEST POUR PRECISION DU CALCUL
      IF(dt1<=dt2)THEN
       f =tf(2,m-1)+dt1*tg
      ELSE
       f =tf(2,m)-dt2*tg
      ENDIF
      RETURN