diff --git a/particle.py b/particle.py
index cb4aea6977cb8bbd4d7e520d215f9716edf1a80c..b3ae6cfd60057428affe718c0dbe8c5d268bd036 100644
--- a/particle.py
+++ b/particle.py
@@ -14,8 +14,8 @@ position = psim.add_vector_property('position')
 velocity = psim.add_vector_property('velocity')
 force = psim.add_vector_property('force', vol=True)
 
-grid = psim.grid_3d(0.0, 4.0, 0.0, 4.0, 0.0, 4.0)
-psim.create_particle_lattice(grid, spacing=[1.0, 1.0, 1.0])
+grid = psim.grid_3d(0.0, 16.0, 0.0, 16.0, 0.0, 16.0)
+psim.create_particle_lattice(grid, spacing=[1.2, 1.2, 1.2])
 psim.create_cell_lists(grid, 2.8, 2.8)
 psim.periodic(2.8)
 
diff --git a/sim/particle_simulation.py b/sim/particle_simulation.py
index 0eae3a74acde241e8af24ce80e4ae207776b1f65..33ce6215b6f92aa9da68a1459c789f3ef402947c 100644
--- a/sim/particle_simulation.py
+++ b/sim/particle_simulation.py
@@ -26,7 +26,7 @@ class ParticleSimulation:
         self.properties = Properties(self)
         self.vars = Variables(self)
         self.arrays = Arrays(self)
-        self.particle_capacity = self.add_var('particle_capacity', Type_Int, 100)
+        self.particle_capacity = self.add_var('particle_capacity', Type_Int, 10000)
         self.nlocal = self.add_var('nlocal', Type_Int)
         self.nghost = self.add_var('nghost', Type_Int)
         self.nparticles = self.nlocal + self.nghost
@@ -167,8 +167,8 @@ class ParticleSimulation:
             Timestep(self, self.ntimesteps, [
                 (EnforcePBC(self.pbc).lower(), 20),
                 (SetupPBC(self.pbc).lower(), 20),
-                (CellListsBuild(self.cell_lists).lower(), 20),
                 UpdatePBC(self.pbc).lower(),
+                (CellListsBuild(self.cell_lists).lower(), 20),
                 PropertiesResetVolatile(self).lower(),
                 self.kernels.lower()
             ]).as_block()