vinter_smooth.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

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

Function/Subroutine Documentation

vinter_smooth()

subroutine vinter_smooth	(		tf,
		integer, dimension(*)	iad,
		integer, dimension(*)	ipos,
		integer, dimension(*)	ilen,
		integer	nel0,
			x,
			dydx,
			y )

Definition at line 29 of file vinter_smooth.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 x(*),dydx(*),y(*),tf(2,*)
      INTEGER I,J1,J,ICONT,J2,J_FIRST,J_LAST
      my_real dydx1,dydx2,dydx3,x_first,x_last
C-----------------------------------------------
      j = 0
      icont = 1
      DO WHILE (icont == 1)
! 
       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
         ELSEIF (ipos(i) >= 1 .AND. x(i) < tf(1,j1-1)) THEN
           ipos(i) = ipos(i) - 1
           icont = 1
         ENDIF
       ENDDO
!
      ENDDO
!
!  smooth interpolation
C
       DO i=1,nel0
C
        j_first = ipos(i)+iad(i)
        j_last  = j_first + 1
        x_first = tf(1,j_first)
        x_last  = tf(1,j_last)
C
        IF (x(i) <= x_first) THEN
          y(i)  = tf(2,j_first)
        ELSEIF (x(i) >= x_last) THEN
          y(i)  = tf(2,j_last)
        ELSE
! within interval
          j1   =ipos(i)+iad(i)
          j2   = j1+1
          dydx(i)=(x(i)-tf(1,j1))/(tf(1,j2)-tf(1,j1))
!
          dydx1 = dydx(i)
          dydx2 = dydx1*dydx1
          dydx3 = dydx1*dydx2
!
          y(i)  = tf(2,j1) + (tf(2,j2)-tf(2,j1))*dydx3*
     .                       (10. - 15.*dydx1 + 6.*dydx2)
!
!!        Y(I)  = TF(2,J1) + (TF(2,J2)-TF(2,J1))*DYDX(I)**3*
!!     .                     (10. - 15.*DYDX(I) + 6.*DYDX(I)**2)
        ENDIF ! IF (X <= X_FIRST)
      ENDDO
!---
      RETURN