vinterdp.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	vinterdp (tf, iad, ipos, ilen, nel0, x, dydx, y)

Function/Subroutine Documentation

vinterdp()

subroutine vinterdp	(		tf,
		integer, dimension(*)	iad,
		integer, dimension(*)	ipos,
		integer, dimension(*)	ilen,
		integer	nel0,
		double precision, dimension(*)	x,
		double precision, dimension(*)	dydx,
		double precision, dimension(*)	y )

Definition at line 28 of file vinterdp.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C
      INTEGER ILEN(*),IPOS(*),IAD(*),NEL0
      my_real tf(2,*)
      INTEGER I,J1,J,ICONT,J2
      DOUBLE PRECISION :: X(*),DYDX(*),Y(*)
      DOUBLE PRECISION  TF2J2,TF2J1,TF1J2,TF1J1,
     .  XDP,YYDP, DYDXDP
C
C
      j = 0
 100  CONTINUE
C 
       j = j+1
       icont = 0
       DO i=1,nel0
         j1 = ipos(i)+iad(i)+1
         IF(j<=ilen(i)-1.AND.x(i)>tf(1,j1))THEN
           ipos(i)=ipos(i)+1
           icont = 1
         ENDIF
       ENDDO
C
      IF(icont==1)GOTO 100
C
      DO i=1,nel0
        j1   =ipos(i)+iad(i)
        j2   = j1+1
 
        tf2j2 = tf(2,j2)
        tf2j1 = tf(2,j1)
        tf1j2 = tf(1,j2)
        tf1j1 = tf(1,j1)
        xdp = x(i)
 
        dydxdp = (tf2j2-tf2j1)/(tf1j2-tf1j1)
        dydx(i) = dydxdp
 
        yydp = tf2j1 + dydxdp * (xdp - tf1j1)
        y(i) = yydp
 
!        DYDX(I)=(TF(2,J2)-TF(2,J1))/(TF(1,J2)-TF(1,J1))
!        Y(I)   = TF(2,J1) + DYDX(I)*(X(I)-TF(1,J1))
      ENDDO
C
      RETURN