facnor.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	facnor (x, d, ii, xnorm, cdg, invert)

Function/Subroutine Documentation

facnor()

subroutine facnor	(		x,
			d,
		integer, dimension(4)	ii,
			xnorm,
			cdg,
		integer	invert )

Definition at line 28 of file facnor.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
     .   x(3,*),d(3,*),xnorm(3,*),cdg(*)
      INTEGER II(4),INVERT
C     REAL
      my_real
     .   x0(4),y0(4),z0(4),nx,ny,nz,xc,yc,zc,s
      INTEGER I,I1,I2,I3,I4
C-----------------------------------------------
C       FACE NORMAL
C-----------------------------------------------
      DO 100 i = 1,4
            x0(i) = x(1,ii(i))-d(1,ii(i))
            y0(i) = x(2,ii(i))-d(2,ii(i))
            z0(i) = x(3,ii(i))-d(3,ii(i))
 100  CONTINUE
C
      nx = (y0(3) - y0(1)) * (z0(4) - z0(2)) - 
     .     (z0(3) - z0(1)) * (y0(4) - y0(2))
      ny = (z0(3) - z0(1)) * (x0(4) - x0(2)) - 
     .     (x0(3) - x0(1)) * (z0(4) - z0(2))
      nz = (x0(3) - x0(1)) * (y0(4) - y0(2)) -
     .     (y0(3) - y0(1)) * (x0(4) - x0(2))
C
      xc = .25*(x0(1) + x0(2) + x0(3) + x0(4)) 
      yc = .25*(y0(1) + y0(2) + y0(3) + y0(4)) 
      zc = .25*(z0(1) + z0(2) + z0(3) + z0(4)) 
      xc = xc - cdg(1)
      yc = yc - cdg(2)
      zc = zc - cdg(3)
      s  = nx * xc + ny * yc + nz * zc
      invert = 1
c      IF(S<ZERO)THEN
c             INVERT = -1
c             I1    = II(2)
c             I2    = II(1)
c             I3    = II(4)
c             I4    = II(3)
c             II(1) = I1
c             II(2) = I2
c             II(3) = I3
c             II(4) = I4
c      ENDIF
C
      DO 200 i = 1,4
            x0(i) = x(1,ii(i))
            y0(i) = x(2,ii(i))
            z0(i) = x(3,ii(i))
 200  CONTINUE
C
      nx = (y0(3) - y0(1)) * (z0(4) - z0(2)) - 
     .     (z0(3) - z0(1)) * (y0(4) - y0(2))
      ny = (z0(3) - z0(1)) * (x0(4) - x0(2)) - 
     .     (x0(3) - x0(1)) * (z0(4) - z0(2))
      nz = (x0(3) - x0(1)) * (y0(4) - y0(2)) -
     .     (y0(3) - y0(1)) * (x0(4) - x0(2))
C
      DO 300 i = 1,4
          xnorm(1,ii(i)) = xnorm(1,ii(i)) + nx
          xnorm(2,ii(i)) = xnorm(2,ii(i)) + ny
          xnorm(3,ii(i)) = xnorm(3,ii(i)) + nz
 300  CONTINUE
C
C-----------------------------------------------
      RETURN