From 8ef3985da7d721df83fe0f4c28c6610ce4c2f5aa Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Tue, 21 Jan 2020 14:58:35 +0100 Subject: [PATCH] debug info --- src/mesh/FEM_mech.f90 | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/src/mesh/FEM_mech.f90 b/src/mesh/FEM_mech.f90 index 267d2f606..014ef1eba 100644 --- a/src/mesh/FEM_mech.f90 +++ b/src/mesh/FEM_mech.f90 @@ -131,11 +131,13 @@ subroutine FEM_mech_init(fieldBC) call DMPlexLabelComplete(mech_mesh,BCLabel,ierr); CHKERRQ(ierr) call DMGetSection(mech_mesh,section,ierr); CHKERRQ(ierr) allocate(pnumComp(1), source=dimPlex) - allocate(pnumDof(dimPlex+1), source = 0) + allocate(pnumDof(0:dimPlex), source = 0) do topologDim = 0, dimPlex call DMPlexGetDepthStratum(mech_mesh,topologDim,cellStart,cellEnd,ierr) CHKERRQ(ierr) - call PetscSectionGetDof(section,cellStart,pnumDof(topologDim+1),ierr) + call PetscSectionGetDof(section,cellStart,pnumDof(topologDim),ierr) + write(6,*) 'start',cellStart,'end',cellEnd + write(6,*) 'topologDim',topologDim,'numDOF',pNumDOF(topologDim) CHKERRQ(ierr) enddo numBC = 0