diff --git a/src/FEM_mech.f90 b/src/FEM_mech.f90
index e10071eb8..0fdaf6283 100644
--- a/src/FEM_mech.f90
+++ b/src/FEM_mech.f90
@@ -140,7 +140,7 @@ subroutine FEM_mech_init(fieldBC)
  nc = dimPlex
  call PetscQuadratureSetData(mechQuad,dimPlex,nc,nQuadrature,qPointsP,qWeightsP,ierr)
  CHKERRQ(ierr)
- call PetscFECreateDefault(PETSC_COMM_SELF,mech_mesh,dimPlex,nc,PETSC_TRUE,prefix, &
+ call PetscFECreateDefault(mech_mesh,dimPlex,nc,PETSC_TRUE,prefix, &
                            integrationOrder,mechFE,ierr); CHKERRQ(ierr)
  call PetscFESetQuadrature(mechFE,mechQuad,ierr); CHKERRQ(ierr)
  call PetscFEGetDimension(mechFE,nBasis,ierr); CHKERRQ(ierr)