diff --git a/src/mesh/FEM_mech.f90 b/src/mesh/FEM_mech.f90
index 03b9758b0..9a490ba6b 100644
--- a/src/mesh/FEM_mech.f90
+++ b/src/mesh/FEM_mech.f90
@@ -190,6 +190,7 @@ subroutine FEM_mech_init(fieldBC)
 
 #endif
   CHKERRQ(ierr)
+  call DMSetSection(mech_mesh,section,ierr); CHKERRQ(ierr)
   do faceSet = 1, numBC
     call ISDestroy(pbcPoints(faceSet),ierr); CHKERRQ(ierr)
   enddo