223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

no clue where or when this was used

This commit is contained in:
Martin Diehl 2019-05-18 23:20:22 +02:00
parent 2ccfc9f6da
commit 2d366f8aca
1 changed files with 0 additions and 27 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
src/mesh_marc.f90
View File

@ -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'
!--------------------------------------------------------------------------------------------------