i7lin3.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	i7lin3 (xi, yi, zi, xa, ya, za, xb, yb, zb, nx, ny, nz, p, alp)

Function/Subroutine Documentation

i7lin3()

subroutine i7lin3	(		xi,
			yi,
			zi,
			xa,
			ya,
			za,
			xb,
			yb,
			zb,
			nx,
			ny,
			nz,
			p,
			alp )

Definition at line 28 of file i7lin3.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-----------------------------------------------
C     REAL
      my_real
     .    xi,yi,zi,xa,ya,za,xb,yb,zb,nx,ny,nz,p,alp
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
C     REAL
      my_real
     .    xab,yab,zab,
     .    xp,yp,zp,norm
C
      xab = xb - xa
      yab = yb - ya
      zab = zb - za
C
      alp = (xi-xa)*xab + (yi-ya)*yab + (zi-za)*zab
      norm = max(em20, xab**2 + yab**2 + zab**2)
      alp = alp / norm
      alp = max(zero,alp)
      alp = min(one,alp)
C
      xp = xa + alp*xab 
      yp = ya + alp*yab 
      zp = za + alp*zab
C
      nx = xi -xp 
      ny = yi -yp 
      nz = zi -zp
C
      p = sqrt(nx**2 + ny**2 +nz**2)
      norm = max(em20,p)
C
      nx = nx / norm
      ny = ny / norm
      nz = nz / norm
C
      RETURN