tetra4v.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	tetra4v (det, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4)

Function/Subroutine Documentation

tetra4v()

subroutine tetra4v	(		det,
			x1,
			x2,
			x3,
			x4,
			y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4 )

Definition at line 28 of file tetra4v.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
     .   det,
     .   x1, x2, x3, x4, y1, y2, y3, y4,z1, z2, z3, z4
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
C                                                                     12
      my_real
     .   x41, y41, z41, x42, y42, z42, x43, y43, z43
      my_real
     .   rx, ry, rz, sx, sy, sz,tx, ty, tz,
     .   tsx, tsy,tsz
C-----------------------------------------------
       x43 = x4 - x3
       y43 = y4 - y3
       z43 = z4 - z3
       x41 = x4 - x1
       y41 = y4 - y1
       z41 = z4 - z1
       x42 = x4 - x2
       y42 = y4 - y2
       z42 = z4 - z2
C
       rx =  -x41
       ry =  -y41
       rz =  -z41
C
       sx =  -x42
       sy =  -y42
       sz =  -z42
C
       tx =  -x43
       ty =  -y43
       tz =  -z43
C
       tsx =  y43*z42 - y42*z43
       tsy =  z43*x42 - z42*x43
       tsz =  x43*y42 - x42*y43
C------V=r*(s^t)/6
       det = abs(rx*tsx + ry*tsy + rz*tsz)*one_over_6
C      
      RETURN