diff --git a/examples/plot.py b/examples/plot.py
index 09a07bd4b002e19918fc1ba930b04394ec6d508b..8f6bfed588e69c0a548722be3c5a932ddee1ca08 100644
--- a/examples/plot.py
+++ b/examples/plot.py
@@ -39,17 +39,17 @@ print "Plotting..."
 # Read Quickshed accelerations
 data=loadtxt("particle_dump.dat")
 id = data[:,0]
-accx_e=data[:,4]
-accy_e=data[:,5]
-accz_e=data[:,6]
+accx_e=data[:,5]
+accy_e=data[:,6]
+accz_e=data[:,7]
 
-accx_bh=data[:,7]
-accy_bh=data[:,8]
-accz_bh=data[:,9]
+accx_bh=data[:,8]
+accy_bh=data[:,9]
+accz_bh=data[:,10]
 
-accx_new=data[:,10]
-accy_new=data[:,11]
-accz_new=data[:,12]
+accx_new=data[:,11]
+accy_new=data[:,12]
+accz_new=data[:,13]
 
 # Sort accelerations
 rank = argsort(id)
diff --git a/examples/test_bh_sorted.c b/examples/test_bh_sorted.c
index 90c5e9ef05449dd42cebfdf6156f2fc6230c3f55..74b5db024642cef41ea826b9c751541e245de9ad 100644
--- a/examples/test_bh_sorted.c
+++ b/examples/test_bh_sorted.c
@@ -1900,10 +1900,10 @@ void test_bh(int N, int nr_threads, int runs, char *fileName) {
 
   /* Dump the particles to a file */
   file = fopen("particle_dump.dat", "w");
-  fprintf(file, "# x y z a_exact.x   a_exact.y    a_exact.z    a_legacy.x    "
+  fprintf(file, "# m x y z a_exact.x   a_exact.y    a_exact.z    a_legacy.x    "
                 "a_legacy.y    a_legacy.z    a_new.x     a_new.y    a_new.z\n");
   for (k = 0; k < N; ++k)
-    fprintf(file, "%d %e %e %e %e %e %e %e %e %e %e %e %e\n", parts[k].id,
+    fprintf(file, "%d %e %e %e %e %e %e %e %e %e %e %e %e %e\n", parts[k].id, parts[k].mass,
             parts[k].x[0], parts[k].x[1], parts[k].x[2], parts[k].a_exact[0],
             parts[k].a_exact[1], parts[k].a_exact[2], parts[k].a_legacy[0],
             parts[k].a_legacy[1], parts[k].a_legacy[2], parts[k].a[0],