From 2d366f8aca49526e0ad8b4ebb0ebba5d299c601c Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Sat, 18 May 2019 23:20:22 +0200 Subject: [PATCH] no clue where or when this was used --- src/mesh_marc.f90 | 27 --------------------------- 1 file changed, 27 deletions(-) diff --git a/src/mesh_marc.f90 b/src/mesh_marc.f90 index 2b605f257..78cedae21 100644 --- a/src/mesh_marc.f90 +++ b/src/mesh_marc.f90 @@ -215,7 +215,6 @@ integer, dimension(:,:), allocatable, private :: & mesh_build_cellnodes, & mesh_build_ipVolumes, & mesh_build_ipCoordinates, & - mesh_cellCenterCoordinates, & mesh_FEasCP @@ -1259,32 +1258,6 @@ subroutine mesh_build_ipCoordinates end subroutine mesh_build_ipCoordinates -!-------------------------------------------------------------------------------------------------- -!> @brief Calculates cell center coordinates. -!-------------------------------------------------------------------------------------------------- -pure function mesh_cellCenterCoordinates(ip,el) - - - integer, intent(in) :: el, & !< element number - ip !< integration point number - real(pReal), dimension(3) :: mesh_cellCenterCoordinates !< x,y,z coordinates of the cell center of the requested IP cell - integer :: t,g,c,n - - t = mesh_element(2,el) ! get element type - g = theMesh%elem%geomType - c = theMesh%elem%cellType - mesh_cellCenterCoordinates = 0.0_pReal - do n = 1,theMesh%elem%nCellnodesPerCell - mesh_cellCenterCoordinates = mesh_cellCenterCoordinates + mesh_cellnode(1:3,mesh_cell(n,ip,el)) - enddo - mesh_cellCenterCoordinates = mesh_cellCenterCoordinates / real(theMesh%elem%nCellnodesPerCell,pReal) - - end function mesh_cellCenterCoordinates - - - - - !-------------------------------------------------------------------------------------------------- !> @brief calculation of IP interface areas, allocate globals '_ipArea', and '_ipAreaNormal' !--------------------------------------------------------------------------------------------------