idege8.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	idege8 (x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, amax, fac, it)

Function/Subroutine Documentation

idege8()

subroutine idege8	(		x1,
			x2,
			x3,
			x4,
			y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4,
			amax,
			fac,
		integer	it )

Definition at line 29 of file idege8.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-----------------------------------------------
      my_real
     .   x1, x2, x3, x4, y1, y2, y3, y4,
     .   z1, z2, z3, z4,amax,fac
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IDE,IT
      my_real
     .   x12,y12,z12,x23,y23,z23,x34,y34,z34,x41,y41,z41,a,
     .   rx,ry,rz,sx,sy,sz,nx,ny,nz
C-----------------------------------------------
       rx = x2+x3-x1-x4
       ry = y2+y3-y1-y4
       rz = z2+z3-z1-z4
       sx = x3+x4-x1-x2
       sy = y3+y4-y1-y2
       sz = z3+z4-z1-z2
       nx = ry * sz - rz * sy 
       ny = rz * sx - rx * sz 
       nz = rx * sy - ry * sx 
C--------A: (2*AREA)^2-----       
       a = nx*nx+ny*ny+nz*nz
       ide = 0
       x12 = x2 - x1
       y12 = y2 - y1
       z12 = z2 - z1
       IF (x12==zero.AND.y12==zero.AND.z12==zero) THEN
        ide = 1
          GOTO 10
       ENDIF
C       
       x23 = x2 - x3
       y23 = y2 - y3
       z23 = z2 - z3
       IF (x23==zero.AND.y23==zero.AND.z23==zero) THEN
        ide = 1
          GOTO 10
       ENDIF
C       
       x34 = x4 - x3
       y34 = y4 - y3
       z34 = z4 - z3
       IF (x34==zero.AND.y34==zero.AND.z34==zero) THEN
        ide = 1
          GOTO 10
       ENDIF
C       
       x41 = x4 - x1
       y41 = y4 - y1
       z41 = z4 - z1
       IF (x41==zero.AND.y41==zero.AND.z41==zero) THEN
        ide = 1
          GOTO 10
       ENDIF
C       
   10   CONTINUE
C     
       IF (a > amax )THEN
              it=ide
          amax = a
             END IF
      RETURN