diff --git a/examples/ZeldovichPancake_3D/plotSolution.py b/examples/ZeldovichPancake_3D/plotSolution.py
index dd1c6243e4637b25da187e93c89490a655303153..163e61eb9d3738c39912c8c42c0b6c6bb9990cac 100644
--- a/examples/ZeldovichPancake_3D/plotSolution.py
+++ b/examples/ZeldovichPancake_3D/plotSolution.py
@@ -86,10 +86,8 @@ phi = sim["/PartType0/Potential"][:]
 
 x -= 0.5 * boxSize
 
-if snap > 40:
-    filename_g = "snapshot_%03d.hdf5"%(snap+1)
-else:
-    filename_g = "snapshot_%03d.hdf5"%(snap)
+# Check for Gadget solution
+filename_g = "snapshot_%03d.hdf5"%(snap)
 if os.path.exists(filename_g):
     sim_g = h5py.File(filename_g, "r")
     x_g = sim_g["/PartType0/Coordinates"][:,0]
diff --git a/examples/ZeldovichPancake_3D/run.sh b/examples/ZeldovichPancake_3D/run.sh
index 196cd1f41a6dfe5f75fbd0dee4d5d418185dc287..9b6b8166ac0d084898b96e2de6b0fc1ef378aecd 100755
--- a/examples/ZeldovichPancake_3D/run.sh
+++ b/examples/ZeldovichPancake_3D/run.sh
@@ -8,7 +8,10 @@ then
 fi
 
 # Run SWIFT
-../swift -s -c -G -t 4 zeldovichPancake.yml 2>&1 | tee output.log
+../swift -a -s -c -G -t 8 zeldovichPancake.yml 2>&1 | tee output.log
 
 # Plot the result
-python plotSolution.py 2
+for i in {0..119}
+do 
+    python plotSolution.py $i
+done
diff --git a/examples/ZeldovichPancake_3D/zeldovichPancake.yml b/examples/ZeldovichPancake_3D/zeldovichPancake.yml
index d0bbe55222447ab5967c6c46b956fc1d6251ede4..03445aa666c4a13f400b2458f6a2fff471248ace 100644
--- a/examples/ZeldovichPancake_3D/zeldovichPancake.yml
+++ b/examples/ZeldovichPancake_3D/zeldovichPancake.yml
@@ -33,7 +33,8 @@ InitialConditions:
 
 Scheduler:
   max_top_level_cells: 8
-  cell_split_size:           400
+  cell_split_size:     50
+  tasks_per_cell:      125
   
 Cosmology:
   Omega_m: 1.
@@ -46,7 +47,7 @@ Cosmology:
 Gravity:
   mesh_side_length:   16
   eta: 0.025
-  theta: 0.85
+  theta: 0.3
   r_cut_max: 5.
-  comoving_softening: 0.0001
-  max_physical_softening: 0.0001
+  comoving_softening: 0.001
+  max_physical_softening: 0.001