diff --git a/src/mesh/mesh_mech_FEM.f90 b/src/mesh/mesh_mech_FEM.f90
index 60c7fd942..404238453 100644
--- a/src/mesh/mesh_mech_FEM.f90
+++ b/src/mesh/mesh_mech_FEM.f90
@@ -357,14 +357,10 @@ subroutine FEM_mech_formResidual(dm_local,xx_local,f_local,dummy,ierr)
 ! evaluate field derivatives
   do cell = cellStart, cellEnd-1                                                                    !< loop over all elements
     
-    ! BEGIN DEBUG
     call PetscSectionGetNumFields(section,numFields,ierr)
     CHKERRQ(ierr)
-    write(6,*) 'numFields', numFields
     call DMPlexGetDepth(dm_local,depth,ierr)
     CHKERRQ(ierr)
-    write(6,*) 'depth', depth
-    ! END DEBUG
     
     call DMPlexVecGetClosure(dm_local,section,x_local,cell,x_scal,ierr)                             !< get Dofs belonging to element
     CHKERRQ(ierr)