n2edge3l.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	n2edge3l (xxi, yyi, zzi, xxj, yyj, zzj, nx, ny, nz, xi, yi, zi, bbb)

Function/Subroutine Documentation

n2edge3l()

subroutine n2edge3l	(	intent(in)	xxi,
		intent(in)	yyi,
		intent(in)	zzi,
		intent(in)	xxj,
		intent(in)	yyj,
		intent(in)	zzj,
		intent(in)	nx,
		intent(in)	ny,
		intent(in)	nz,
		intent(in)	xi,
		intent(in)	yi,
		intent(in)	zi,
		intent(out)	bbb )

Definition at line 28 of file n2edge3l.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 , INTENT(IN)     :: 
     1                           xxi    ,yyi      ,zzi    ,
     2                           xxj    ,yyj      ,zzj    ,
     3                           nx     ,ny       ,nz     ,
     4                           xi     ,yi       ,zi     
      my_real , INTENT(OUT)    :: bbb
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER I
      my_real
     .     xmij(3),nij(3),di(3),aa,norm
C-----------------------------------------------
             xmij(1)=xxj-xxi
             xmij(2)=yyj-yyi
             xmij(3)=zzj-zzi
             nij(1)= ny*xmij(3) - nz*xmij(2)
             nij(2)= nz*xmij(1) - nx*xmij(3)
             nij(3)= nx*xmij(2) - ny*xmij(1)
             aa = nij(1)*nij(1)+nij(2)*nij(2)+nij(3)*nij(3)
             norm=max(em20,sqrt(aa))
             di(1)=half*(xxi+xxj)-xi
             di(2)=half*(yyi+yyj)-yi
             di(3)=half*(zzi+zzj)-zi
             bbb = (di(1)*nij(1)+di(2)*nij(2)+di(3)*nij(3))/norm
             IF (abs(bbb)<em3) then
               xmij(1)=xxj-xi
               xmij(2)=yyj-yi
               xmij(3)=zzj-zi
               aa = di(1)*xmij(1)+di(2)*xmij(2)+di(3)*xmij(3)
               IF (aa>abs(bbb)) bbb = aa
             END IF
C
      RETURN