From e3f2ad34b23c7225c7e5977493c19a0cfb6e3fca Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Sun, 23 Sep 2018 16:32:13 +0200 Subject: [PATCH] not needed --- src/meshFEM.f90 | 4 +--- 1 file changed, 1 insertion(+), 3 deletions(-) diff --git a/src/meshFEM.f90 b/src/meshFEM.f90 index 8553fb2a6..22568c99d 100644 --- a/src/meshFEM.f90 +++ b/src/meshFEM.f90 @@ -25,7 +25,6 @@ use PETScis mesh_NcpElems, & !< total number of CP elements in mesh mesh_NcpElemsGlobal, & mesh_Nnodes, & !< total number of nodes in mesh - mesh_maxNnodes, & !< max number of nodes in any CP element mesh_maxNips, & !< max number of IPs in any CP element mesh_maxNipNeighbors @@ -219,12 +218,11 @@ subroutine mesh_init(ip,el) CHKERRQ(ierr) FE_Nips(FE_geomtype(1_pInt)) = FEM_Zoo_nQuadrature(dimPlex,integrationOrder) - mesh_maxNnodes = FE_Nnodes(1_pInt) mesh_maxNips = FE_Nips(1_pInt) call mesh_FEM_build_ipCoordinates(dimPlex,FEM_Zoo_QuadraturePoints(dimPlex,integrationOrder)%p) call mesh_FEM_build_ipVolumes(dimPlex) - allocate (mesh_element (4_pInt+mesh_maxNnodes,mesh_NcpElems)); mesh_element = 0_pInt + allocate (mesh_element (4_pInt+FE_nodes(1_pInt),mesh_NcpElems)); mesh_element = 0_pInt do j = 1, mesh_NcpElems mesh_element( 1,j) = j mesh_element( 2,j) = mesh_elemType ! elem type