geom.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	geom (a, b, c, center_x, center_y, center_z, vol)

Function/Subroutine Documentation

geom()

subroutine geom	(	dimension(3), intent(in)	a,
		dimension(3), intent(in)	b,
		dimension(3), intent(in)	c,
		intent(inout)	center_x,
		intent(inout)	center_y,
		intent(inout)	center_z,
		intent(inout)	vol )

Definition at line 29 of file geom.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, INTENT(IN) :: a(3), b(3), c(3)
      my_real, INTENT(INOUT) :: center_x,center_y,center_z, vol
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real :: norm(3), ab, bc, ca
C-----------------------------------------------
C   S o u r c e   L i n e s 
C-----------------------------------------------
 
      norm(1) = (b(2) - a(2)) * (c(3) - a(3)) - (b(3) - a(3)) * (c(2) - a(2))
      norm(2) = (b(3) - a(3)) * (c(1) - a(1)) - (b(1) - a(1)) * (c(3) - a(3))
      norm(3) = (b(1) - a(1)) * (c(2) - a(2)) - (b(2) - a(2)) * (c(1) - a(1))
      
      vol = vol + dot_product(a, norm)
      
      ab = dot_product(half * (a + b), half * (a + b))
      bc = dot_product(half * (b + c), half * (b + c))
      ca = dot_product(half * (c + a), half * (c + a))
 
      center_x = center_x + norm(1) * (ab + bc + ca)
      center_y = center_y + norm(2) * (ab + bc + ca)
      center_z = center_z + norm(3) * (ab + bc + ca)
 
C-----------------------------------------------