diff --git a/examples/IsothermalPotential/GravityOnly/README b/examples/IsothermalPotential/GravityOnly/README
new file mode 100644
index 0000000000000000000000000000000000000000..90fb1872aa2301cab133b8a20a8cd8de724d4553
--- /dev/null
+++ b/examples/IsothermalPotential/GravityOnly/README
@@ -0,0 +1,6 @@
+;
+; this probelm generates a set of gravity particles in an isothermal
+; potential and follows their orbits. Tests verify consdevation of
+; energy and angular momentum
+;
+;
diff --git a/examples/IsothermalPotential/GravityOnly/externalGravity.yml b/examples/IsothermalPotential/GravityOnly/externalGravity.yml
index 150b06dd0f3fbcc0f3ef60be24dbb0554bd63e4d..0de99779f07591a5b71be11b75bc56ec741ddaed 100644
--- a/examples/IsothermalPotential/GravityOnly/externalGravity.yml
+++ b/examples/IsothermalPotential/GravityOnly/externalGravity.yml
@@ -20,55 +20,28 @@ Statistics:
 # Parameters governing the snapshots
 Snapshots:
   basename:            Isothermal # Common part of the name of output files
-  time_first:          0.        # Time of the first output (in internal units)
-  delta_time:          0.02      # Time difference between consecutive outputs (in internal units)
-
-# Parameters for the task scheduling
-Scheduler:
-  nr_threads:       8        # The number of threads per MPI rank to use.
-  nr_queues:        0        # The number of task queues to use. Use 0  to let the system decide.
-  cell_max_size:    8000000  # Maximal number of interactions per task (this is the default value).
-  cell_sub_size:    8000000  # Maximal number of interactions per sub-task  (this is the default value).
-  cell_split_size:  400      # Maximal number of particles per cell (this is the default value).
-
+  time_first:          0.         # Time of the first output (in internal units)
+  delta_time:          0.02       # Time difference between consecutive outputs (in internal units)
 
 # Parameters for the hydrodynamics scheme
 SPH:
   resolution_eta:        1.2349   # Target smoothing length in units of the mean inter-particle separation (1.2349 == 48Ngbs with the cubic spline kernel).
   delta_neighbours:      1.       # The tolerance for the targetted number of neighbours.
   CFL_condition:         0.1      # Courant-Friedrich-Levy condition for time integration.
-  max_ghost_iterations:  30       # Maximal number of iterations allowed to converge towards the smoothing length.
   max_smoothing_length:  40.      # Maximal smoothing length allowed (in internal units).
 
 # Parameters related to the initial conditions
 InitialConditions:
   file_name:  Isothermal.hdf5       # The file to read
-  h_scaling:  1.                    # A scaling factor to apply to all smoothing lengths in the ICs.
   shift_x:    100.                  # A shift to apply to all particles read from the ICs (in internal units).
   shift_y:    100.
   shift_z:    100.
-
-# Parameters govering domain decomposition
-DomainDecomposition:
-  initial_type:       m     # The initial strategy ("g", "m", "w", or "v"). See documentation for details.
-  initial_grid_x:    10     # Grid size if the 'g' strategy is chosen.
-  initial_grid_y:    10
-  initial_grid_z:    10
-  repartition_type:   b     # The re-decomposition strategy ("n", "b", "v", "e" or "x"). See documentation for details.
  
-
 # External potential parameters
-PointMass:
-	position_x:      50.     # location of external point mass in internal units
-	position_y:      50.
-	position_z:      50.	
-	mass:            1e10     # mass of external point mass in internal units
-	timestep_mult:   0.03     # controls time step
-
 IsothermalPotential:
-	position_x:      100.     # location of centre of isothermal potential in internal units
-	position_y:      100.
-	position_z:      100.	
-	vrot:            200.     # rotation speed of isothermal potential in internal units
-	timestep_mult:   0.03     # controls time step
+  position_x:      100.     # location of centre of isothermal potential in internal units
+  position_y:      100.
+  position_z:      100.	
+  vrot:            200.     # rotation speed of isothermal potential in internal units
+  timestep_mult:   0.03     # controls time step
 
diff --git a/examples/IsothermalPotential/GravityOnly/run.sh b/examples/IsothermalPotential/GravityOnly/run.sh
index 89e7c10ce19ce0576127551f63eea5c9e4e955ff..caa616f56ed18ec84209829b183ea08d89b71348 100755
--- a/examples/IsothermalPotential/GravityOnly/run.sh
+++ b/examples/IsothermalPotential/GravityOnly/run.sh
@@ -3,8 +3,8 @@
 # Generate the initial conditions if they are not present.
 if [ ! -e Isothermal.hdf5 ]
 then
-    echo "Generating initial conditions for the point mass potential box example..."
-    python makeIC.py 1000
+    echo "Generating initial conditions for the isothermal potential box example..."
+    python makeIC.py 1000 1
 fi
 
-../swift -g -t 2 externalGravity.yml
+../../swift -g -t 2 externalGravity.yml