From f2d1343cdf5045b3c022d762a50b75dbd07c03aa Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Wed, 13 May 2020 12:58:39 +0200 Subject: [PATCH] need to select the whole mesh --- src/mesh/mesh_mech_FEM.f90 | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/src/mesh/mesh_mech_FEM.f90 b/src/mesh/mesh_mech_FEM.f90 index 8206f4174..0deca3e03 100644 --- a/src/mesh/mesh_mech_FEM.f90 +++ b/src/mesh/mesh_mech_FEM.f90 @@ -73,7 +73,7 @@ subroutine FEM_mech_init(fieldBC) PetscQuadrature :: mechQuad, functional PetscDS :: mechDS PetscDualSpace :: mechDualSpace - DMLabel, dimension(:),pointer :: pBCLabel + DMLabel, dimension(:),pointer :: nolabel=> NULL() DMLabel :: BCLabel PetscInt, dimension(:), pointer :: pNumComp, pNumDof, pBcField, pBcPoint @@ -182,8 +182,7 @@ subroutine FEM_mech_init(fieldBC) call DMPlexCreateSection(mech_mesh,dimPlex,1,pNumComp,pNumDof, & numBC,pBcField,pBcComps,pBcPoints,PETSC_NULL_IS,section,ierr) #else - allocate(pBClabel(1),source=BClabel) - call DMPlexCreateSection(mech_mesh,pBClabel,pNumComp,pNumDof, & + call DMPlexCreateSection(mech_mesh,nolabel,pNumComp,pNumDof, & numBC,pBcField,pBcComps,pBcPoints,PETSC_NULL_IS,section,ierr) #endif