centroid.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	centroid3 (nel, length, nft, ixs, xgrid, elem_centroid, face_centroid)
subroutine	centroid3t (nel, length, nft, ixs, xgrid, elem_centroid, face_centroid)
subroutine	centroid2 (nel, length, nft, ixq, xgrid, elem_centroid, face_centroid)
subroutine	centroid2t (nel, length, nft, ixtg, xgrid, elem_centroid, face_centroid)

Function/Subroutine Documentation

centroid2()

subroutine centroid2	(	integer, intent(in)	nel,
		integer, intent(in)	length,
		integer, intent(in)	nft,
		integer, dimension(nixq, *), intent(in)	ixq,
		dimension(3, *), intent(in)	xgrid,
		dimension(3, nel), intent(out)	elem_centroid,
		dimension(3, 6, nel), intent(out)	face_centroid )

Definition at line 195 of file centroid.F.

      use element_mod , only : nixq
C-----------------------------------------------
C     D e s c r i p t i o n
C-----------------------------------------------
C     Computes the centroids of the elements
C-----------------------------------------------
C     I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C     C o m m o n   B l o c k s
C-----------------------------------------------
!     NIXS
C-----------------------------------------------
C     D u m m y   A r g u m e n t s
C-----------------------------------------------      
      INTEGER, INTENT(IN) :: NEL, NFT, IXQ(NIXQ, *), LENGTH
      my_real, INTENT(IN) :: xgrid(3, *)
      my_real, INTENT(OUT) :: elem_centroid(3, nel), face_centroid(3, 6, nel)
C-----------------------------------------------
C     L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I, II, NODE1, NODE2, NODE3, NODE4
      my_real
     .     x1(3), x2(3), x3(3), x4(3)
      DO ii = 1, nel
         i = ii + nft
C     Nodes of the element locally stored in NODE* for memory access
         node1 = ixq(2, i)
         node2 = ixq(3, i)
         node3 = ixq(4, i)
         node4 = ixq(5, i)
C     Node coordinates
         x1(1:3) = xgrid(1:3, node1)
         x2(1:3) = xgrid(1:3, node2)
         x3(1:3) = xgrid(1:3, node3)
         x4(1:3) = xgrid(1:3, node4)
C     Centroid
         elem_centroid(1:3, ii) = fourth * (x1(1:3) + x2(1:3) + x3(1:3) + x4(1:3))
C     Face centroid
C     Face 1 
         face_centroid(1:3, 1, ii) = half * (x1(1:3) + x2(1:3))
C     Face 2 
         face_centroid(1:3, 2, ii) = half * (x2(1:3) + x3(1:3))
C     Face 3 
         face_centroid(1:3, 3, ii) = half * (x3(1:3) + x4(1:3))
C     Face 4 
         face_centroid(1:3, 4, ii) = half * (x1(1:3) + x4(1:3))
 
      ENDDO

centroid2t()

subroutine centroid2t	(	integer, intent(in)	nel,
		integer, intent(in)	length,
		integer, intent(in)	nft,
		integer, dimension(nixtg, *), intent(in)	ixtg,
		dimension(3, *), intent(in)	xgrid,
		dimension(3, nel), intent(out)	elem_centroid,
		dimension(3, 6, nel), intent(out)	face_centroid )

Definition at line 256 of file centroid.F.

      use element_mod , only : nixtg
C-----------------------------------------------
C     D e s c r i p t i o n
C-----------------------------------------------
C     Computes the centroids of the elements
C-----------------------------------------------
C     I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C     C o m m o n   B l o c k s
C-----------------------------------------------
!     NIXS
C-----------------------------------------------
C     D u m m y   A r g u m e n t s
C-----------------------------------------------      
      INTEGER, INTENT(IN) :: NEL, NFT, IXTG(NIXTG, *), LENGTH
      my_real, INTENT(IN) :: xgrid(3, *)
      my_real, INTENT(OUT) :: elem_centroid(3, nel), face_centroid(3, 6, nel)
C-----------------------------------------------
C     L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I, II, NODE1, NODE2, NODE3
      my_real
     .     x1(3), x2(3), x3(3)
      DO ii = 1, nel
         i = ii + nft
C     Nodes of the element locally stored in NODE* for memory access
         node1 = ixtg(2, i)
         node2 = ixtg(3, i)
         node3 = ixtg(4, i)
C     Node coordinates
         x1(1:3) = xgrid(1:3, node1)
         x2(1:3) = xgrid(1:3, node2)
         x3(1:3) = xgrid(1:3, node3)
C     Centroid
         elem_centroid(1:3, ii) = third * (x1(1:3) + x2(1:3) + x3(1:3))
C     Face centroid
C     Face 1 
         face_centroid(1:3, 1, ii) = half * (x1(1:3) + x2(1:3))
C     Face 2 
         face_centroid(1:3, 2, ii) = half * (x2(1:3) + x3(1:3))
C     Face 3 
         face_centroid(1:3, 3, ii) = half * (x3(1:3) + x1(1:3))
 
      ENDDO

centroid3()

subroutine centroid3	(	integer, intent(in)	nel,
		integer, intent(in)	length,
		integer, intent(in)	nft,
		integer, dimension(nixs, *), intent(in)	ixs,
		dimension(3, *), intent(in)	xgrid,
		dimension(3, nel), intent(out)	elem_centroid,
		dimension(3, 6, nel), intent(out)	face_centroid )

Definition at line 31 of file centroid.F.

      use element_mod , only : nixs
C-----------------------------------------------
C     D e s c r i p t i o n
C-----------------------------------------------
C     Computes the centroids of the elements
C-----------------------------------------------
C     I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C     C o m m o n   B l o c k s
C-----------------------------------------------
!     NIXS
C-----------------------------------------------
C     D u m m y   A r g u m e n t s
C-----------------------------------------------      
      INTEGER, INTENT(IN) :: NEL, NFT, IXS(NIXS, *), LENGTH
      my_real, INTENT(IN) :: xgrid(3, *)
      my_real, INTENT(OUT) :: elem_centroid(3, nel), face_centroid(3, 6, nel)
C-----------------------------------------------
C     L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I, II, NODE1, NODE2, NODE3, NODE4, NODE5, NODE6, NODE7, NODE8
      my_real
     .     x1(3), x2(3), x3(3), x4(3), x5(3), x6(3), x7(3), x8(3)
      DO ii = 1, nel
         i = ii + nft
C     Nodes of the element locally stored in NODE* for memory access
         node1 = ixs(2, i)
         node2 = ixs(3, i)
         node3 = ixs(4, i)
         node4 = ixs(5, i)
         node5 = ixs(6, i)
         node6 = ixs(7, i)
         node7 = ixs(8, i)
         node8 = ixs(9, i)  
C     Node coordinates
         x1(1:3) = xgrid(1:3, node1)
         x2(1:3) = xgrid(1:3, node2)
         x3(1:3) = xgrid(1:3, node3)
         x4(1:3) = xgrid(1:3, node4)
         x5(1:3) = xgrid(1:3, node5)
         x6(1:3) = xgrid(1:3, node6)
         x7(1:3) = xgrid(1:3, node7)
         x8(1:3) = xgrid(1:3, node8)
C     Centroid
         elem_centroid(1:3, ii) = one_over_8 * (x1(1:3) + x2(1:3) + x3(1:3) + x4(1:3) + 
     .        x5(1:3) + x6(1:3) + x7(1:3) + x8(1:3))
C     Face centroid
C     Face 1 
         face_centroid(1:3, 1, ii) = fourth * (x1(1:3) + x2(1:3) + x3(1:3) + x4(1:3))
C     Face 2 
         face_centroid(1:3, 2, ii) = fourth * (x3(1:3) + x4(1:3) + x8(1:3) + x7(1:3))
C     Face 3 
         face_centroid(1:3, 3, ii) = fourth * (x5(1:3) + x6(1:3) + x7(1:3) + x8(1:3))
C     Face 4 
         face_centroid(1:3, 4, ii) = fourth * (x1(1:3) + x2(1:3) + x6(1:3) + x5(1:3))
C     Face 5 
         face_centroid(1:3, 5, ii) = fourth * (x2(1:3) + x3(1:3) + x7(1:3) + x6(1:3))
C     Face 6 
         face_centroid(1:3, 6, ii) = fourth * (x1(1:3) + x4(1:3) + x8(1:3) + x5(1:3))
      ENDDO

centroid3t()

subroutine centroid3t	(	integer, intent(in)	nel,
		integer, intent(in)	length,
		integer, intent(in)	nft,
		integer, dimension(nixs, *), intent(in)	ixs,
		dimension(3, *), intent(in)	xgrid,
		dimension(3, nel), intent(out)	elem_centroid,
		dimension(3, 6, nel), intent(out)	face_centroid )

Definition at line 107 of file centroid.F.

      use element_mod , only : nixs
C-----------------------------------------------
C     D e s c r i p t i o n
C-----------------------------------------------
C     Computes the centroids of the elements
C-----------------------------------------------
C     I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C     C o m m o n   B l o c k s
C-----------------------------------------------
!     NIXS
C-----------------------------------------------
C     D u m m y   A r g u m e n t s
C-----------------------------------------------      
      INTEGER, INTENT(IN) :: NEL, NFT, IXS(NIXS, *), LENGTH
      my_real, INTENT(IN) :: xgrid(3, *)
      my_real, INTENT(OUT) :: elem_centroid(3, nel), face_centroid(3, 6, nel)
C-----------------------------------------------
C     L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I, II, NODE1, NODE2, NODE3, NODE4
      my_real
     .     x1(3), x2(3), x3(3), x4(3), vol
      my_real 
     .     a(3), b(3), c(3)
      DO ii = 1, nel
         i = ii + nft
C     Nodes of the element locally stored in NODE* for memory access
         node1 = ixs(2, i)
         node2 = ixs(4, i)
         node3 = ixs(7, i)
         node4 = ixs(6, i)
C     Node coordinates
         x1(1:3) = xgrid(1:3, node1)
         x2(1:3) = xgrid(1:3, node2)
         x3(1:3) = xgrid(1:3, node3)
         x4(1:3) = xgrid(1:3, node4)
         elem_centroid(1:3, ii) = zero
C     Centroid
         vol = zero
         a(1:3) = x1(1:3)
         b(1:3) = x2(1:3)
         c(1:3) = x4(1:3)
         CALL geom(a, b, c,elem_centroid(1,ii),elem_centroid(2,ii),elem_centroid(3,ii), vol)
 
         a(1:3) = x1(1:3)
         b(1:3) = x3(1:3)
         c(1:3) = x2(1:3)
         CALL geom(a, b, c,elem_centroid(1,ii),elem_centroid(2,ii),elem_centroid(3,ii), vol)
 
         a(1:3) = x1(1:3)
         b(1:3) = x4(1:3)
         c(1:3) = x3(1:3)
         CALL geom(a, b, c,elem_centroid(1,ii),elem_centroid(2,ii),elem_centroid(3,ii), vol)
 
         a(1:3) = x2(1:3)
         b(1:3) = x3(1:3)
         c(1:3) = x4(1:3)
         CALL geom(a, b, c,elem_centroid(1,ii),elem_centroid(2,ii),elem_centroid(3,ii), vol)
         
         vol = vol / 6.d0
         elem_centroid(1:3, ii) = elem_centroid(1:3, ii) / 12.d0 / vol
 
         !ELEM_CENTROID(1:3, II) = FOURTH * (X1(1:3) + X2(1:3) + X3(1:3) + X4(1:3))
C     Face centroid
C     Face 1 
         face_centroid(1:3, 5, ii) = third * (x1(1:3) + x2(1:3) + x3(1:3))
C     Face 2 
         face_centroid(1:3, 6, ii) = third * (x1(1:3) + x2(1:3) + x4(1:3))
C     Face 3 
         face_centroid(1:3, 2, ii) = third * (x2(1:3) + x3(1:3) + x4(1:3))
C     Face 4 
         face_centroid(1:3, 4, ii) = third * (x1(1:3) + x3(1:3) + x4(1:3))
      ENDDO