diff --git a/pyop2/petsc_base.py b/pyop2/petsc_base.py
index 3575d7cb3..7db96b9fa 100644
--- a/pyop2/petsc_base.py
+++ b/pyop2/petsc_base.py
@@ -229,9 +229,9 @@ def vecscatters(self):
         scatters = []
         for d in self:
             size = d.size * d.cdim
-            vscat = PETSc.Scatter().create(d.layout_vec, None, self.layout_vec,
-                                           PETSc.IS().createStride(size, offset, 1,
-                                                                   comm=d.comm))
+            vscat = PETSc.Scatter().createWithData(d.layout_vec, None, self.layout_vec,
+                                                   PETSc.IS().createStride(size, offset, 1,
+                                                                           comm=d.comm))
             offset += size
             scatters.append(vscat)
         return tuple(scatters)