gpamela/gptrd/xtrinter.F

*
* $Id$
*
* $Log$
*
*CMZ :  3.00/00 09/04/2001  15.38.41  by  Marialuigia Ambriola
*-- Author :    Marialuigia Ambriola   05/04/2001
      REAL FUNCTION XTRINTER(S,Y,X,N)
****************************************
* makes a linear interpolation
* called by GPUSTEP
****************************************
      REAL X(N),Y(N),S,XX
      INTEGER N
      XTRINTER=0.
      XX=0.
c      PRINT*,'N,S',N,S
c      PRINT*,'Y=',Y
      DO I=2,N
         IF(S.EQ.Y(I-1)) THEN
            XX=X(I-1)
            GOTO 100
         ENDIF
         IF(S.LT.Y(I).AND.S.GT.Y(I-1)) THEN
            XX=X(I-1)+(X(I)-X(I-1))*(S-Y(I-1))/(Y(I)-Y(I-1))
            GOTO 100
         ENDIF
      ENDDO
 100  XTRINTER=XX
c      PRINT*,'EXTRACTED ENERGY',XX
      RETURN
      END
1	*
2	* $Id$
3	*
4	* $Log$
5	*
6	*CMZ : 3.00/00 09/04/2001 15.38.41 by Marialuigia Ambriola
7	*-- Author : Marialuigia Ambriola 05/04/2001
8	REAL FUNCTION XTRINTER(S,Y,X,N)
9	****************************************
10	* makes a linear interpolation
11	* called by GPUSTEP
12	****************************************
13	REAL X(N),Y(N),S,XX
14	INTEGER N
15	XTRINTER=0.
16	XX=0.
17	c PRINT*,'N,S',N,S
18	c PRINT*,'Y=',Y
19	DO I=2,N
20	IF(S.EQ.Y(I-1)) THEN
21	XX=X(I-1)
22	GOTO 100
23	ENDIF
24	IF(S.LT.Y(I).AND.S.GT.Y(I-1)) THEN
25	XX=X(I-1)+(X(I)-X(I-1))*(S-Y(I-1))/(Y(I)-Y(I-1))
26	GOTO 100
27	ENDIF
28	ENDDO
29	100 XTRINTER=XX
30	c PRINT*,'EXTRACTED ENERGY',XX
31	RETURN
32	END