ninterp.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	ninterp (ifunc, npc, pld, npoint, xv, yv)

Function/Subroutine Documentation

ninterp()

subroutine ninterp	(	integer	ifunc,
		integer, dimension(*)	npc,
			pld,
		integer	npoint,
			xv,
			yv )

Definition at line 31 of file ninterp.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IFUNC,NPC(*),NPOINT
C     REAL
      my_real
     .  pld(*),xv(*),yv(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real
     .   pente, y
      INTEGER RANG(MVSIZ)
      INTEGER NDEB, NREST, I,
     .        IPT, NPT, II
C-----------------------------------------------
      ndeb=0
      nrest=npoint
 100  CONTINUE
C---------------------------------
      DO i=1,min(mvsiz,nrest)
        rang(i)=0
      ENDDO
C---------------------------------
      npt = (npc(ifunc+1)-npc(ifunc))/2
      DO ipt=1,npt-1
        ii=npc(ifunc)+2*(ipt-1)
        DO i=1,min(mvsiz,nrest)
          IF (pld(ii)<=xv(ndeb+i)
     .    .AND.xv(ndeb+i)<=pld(ii+2)) rang(i)=ipt
        ENDDO
      ENDDO
      ii=npc(ifunc)+2*(npt-1)
      DO i=1,min(mvsiz,nrest)
          IF (xv(ndeb+i)>pld(ii)) rang(i)=npt
      ENDDO
C---------------------------------
      DO i=1,min(mvsiz,nrest)
        IF (rang(i)==0) THEN
          ii=npc(ifunc)
          pente=(pld(ii+3)-pld(ii+1))/(pld(ii+2)-pld(ii))
          yv(ndeb+i) =pld(ii+1)+pente*(xv(ndeb+i)-pld(ii))
        ELSEIF (rang(i)==npt) THEN
          ii=npc(ifunc)+2*(npt-1)
          pente=(pld(ii+1)-pld(ii-1))/(pld(ii)-pld(ii-2))
          yv(ndeb+i) =pld(ii+1)
     .       +max(-pld(ii+1),pente*(xv(ndeb+i)-pld(ii)))
        ELSE
          ii=npc(ifunc)+2*(rang(i)-1)
          pente=(pld(ii+3)-pld(ii+1))/(pld(ii+2)-pld(ii))
          yv(ndeb+i) =pld(ii+1)
     .       +pente*(xv(ndeb+i)-pld(ii))
        ENDIF
      ENDDO
C---------------------------------
C     Groupe suivant.
C---------------------------------
      IF (nrest-mvsiz>0) THEN
        nrest=nrest-mvsiz
        ndeb =ndeb +mvsiz
        GOTO 100
      END IF
C-----------------------------------------------
      RETURN