diff --git a/examples/SedovBlast/run.sh b/examples/SedovBlast/run.sh
new file mode 100755
index 0000000000000000000000000000000000000000..0fdec1e04d293aaa37c3a1e9d7bc6d5203e65ef9
--- /dev/null
+++ b/examples/SedovBlast/run.sh
@@ -0,0 +1,3 @@
+#!/bin/bash
+
+../swift -s sedov.yml
diff --git a/examples/SedovBlast/sedov.yml b/examples/SedovBlast/sedov.yml
new file mode 100644
index 0000000000000000000000000000000000000000..1fb0e9f42c7e50e4751b376ba8ab0a7f0a018bbd
--- /dev/null
+++ b/examples/SedovBlast/sedov.yml
@@ -0,0 +1,48 @@
+
+# Define the system of units to use internally.
+UnitSystem:
+ UnitMass_in_cgs: 1 # Grams
+ UnitLength_in_cgs: 1 # Centimeters
+ UnitVelocity_in_cgs: 1 # Centimeters per second
+ UnitCurrent_in_cgs: 1 # Amperes
+ UnitTemp_in_cgs: 1 # Kelvin
+
+# Parameters for the task scheduling
+Scheduler:
+ nr_threads: 16 # 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).
+
+# Parameters governing the time integration
+TimeIntegration:
+ time_begin: 0. # The starting time of the simulation (in internal units).
+ time_end: 1. # The end time of the simulation (in internal units).
+ dt_min: 1e-6 # The minimal time-step size of the simulation (in internal units).
+ dt_max: 1e-2 # The maximal time-step size of the simulation (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: 1. # Maximal smoothing length allowed (in internal units).
+
+# Parameters related to the initial conditions
+InitialConditions:
+ file_name: ./sedov.hdf5 # The file to read
+ h_scaling: 1. # A scaling factor to apply to all smoothing lengths in the ICs.
+ shift_x: 0. # A shift to apply to all particles read from the ICs (in internal units).
+ shift_y: 0.
+ shift_z: 0.
+
+# Parameters govering domain decomposition
+DomainDecomposition:
+ initial_type: m # The initial strategy ("g", "m", "w", or "v"). See documentation for details.
+ intitial_grid_x: 10 # Grid size if the 'g' strategy is chosen.
+ intitial_grid_y: 10
+ intitial_grid_z: 10
+ repartition_type: b # The re-decomposition strategy ("n", "b", "v", "e" or "x"). See documentation for details.
+
diff --git a/examples/SodShock/run.sh b/examples/SodShock/run.sh
new file mode 100755
index 0000000000000000000000000000000000000000..445367a02e37e141e7232e6b865c18f1d73c8b36
--- /dev/null
+++ b/examples/SodShock/run.sh
@@ -0,0 +1,3 @@
+#!/bin/bash
+
+../swift -s sodShock.yml
diff --git a/examples/SodShock/sodShock.yml b/examples/SodShock/sodShock.yml
new file mode 100644
index 0000000000000000000000000000000000000000..24e6f5d20c8d026f2c59960834607b5306a3665c
--- /dev/null
+++ b/examples/SodShock/sodShock.yml
@@ -0,0 +1,48 @@
+
+# Define the system of units to use internally.
+UnitSystem:
+ UnitMass_in_cgs: 1 # Grams
+ UnitLength_in_cgs: 1 # Centimeters
+ UnitVelocity_in_cgs: 1 # Centimeters per second
+ UnitCurrent_in_cgs: 1 # Amperes
+ UnitTemp_in_cgs: 1 # Kelvin
+
+# Parameters for the task scheduling
+Scheduler:
+ nr_threads: 16 # 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).
+
+# Parameters governing the time integration
+TimeIntegration:
+ time_begin: 0. # The starting time of the simulation (in internal units).
+ time_end: 1. # The end time of the simulation (in internal units).
+ dt_min: 1e-6 # The minimal time-step size of the simulation (in internal units).
+ dt_max: 1e-2 # The maximal time-step size of the simulation (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: 0.01 # Maximal smoothing length allowed (in internal units).
+
+# Parameters related to the initial conditions
+InitialConditions:
+ file_name: ./sodShock.hdf5 # The file to read
+ h_scaling: 1. # A scaling factor to apply to all smoothing lengths in the ICs.
+ shift_x: 0. # A shift to apply to all particles read from the ICs (in internal units).
+ shift_y: 0.
+ shift_z: 0.
+
+# Parameters govering domain decomposition
+DomainDecomposition:
+ initial_type: m # The initial strategy ("g", "m", "w", or "v"). See documentation for details.
+ intitial_grid_x: 10 # Grid size if the 'g' strategy is chosen.
+ intitial_grid_y: 10
+ intitial_grid_z: 10
+ repartition_type: b # The re-decomposition strategy ("n", "b", "v", "e" or "x"). See documentation for details.
+
diff --git a/examples/UniformBox/run.sh b/examples/UniformBox/run.sh
new file mode 100755
index 0000000000000000000000000000000000000000..d9a06997eb3c4d4ff480722883337571906c79bf
--- /dev/null
+++ b/examples/UniformBox/run.sh
@@ -0,0 +1,3 @@
+#!/bin/bash
+
+../swift -s uniformBox.yml
diff --git a/examples/UniformBox/uniformBox.yml b/examples/UniformBox/uniformBox.yml
new file mode 100644
index 0000000000000000000000000000000000000000..2c5588215e9a845391e7606bfbcb877772853537
--- /dev/null
+++ b/examples/UniformBox/uniformBox.yml
@@ -0,0 +1,48 @@
+
+# Define the system of units to use internally.
+UnitSystem:
+ UnitMass_in_cgs: 1 # Grams
+ UnitLength_in_cgs: 1 # Centimeters
+ UnitVelocity_in_cgs: 1 # Centimeters per second
+ UnitCurrent_in_cgs: 1 # Amperes
+ UnitTemp_in_cgs: 1 # Kelvin
+
+# Parameters for the task scheduling
+Scheduler:
+ nr_threads: 16 # 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).
+
+# Parameters governing the time integration
+TimeIntegration:
+ time_begin: 0. # The starting time of the simulation (in internal units).
+ time_end: 1. # The end time of the simulation (in internal units).
+ dt_min: 1e-6 # The minimal time-step size of the simulation (in internal units).
+ dt_max: 1e-2 # The maximal time-step size of the simulation (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: 0.1 # Maximal smoothing length allowed (in internal units).
+
+# Parameters related to the initial conditions
+InitialConditions:
+ file_name: ./uniformBox.hdf5 # The file to read
+ h_scaling: 1. # A scaling factor to apply to all smoothing lengths in the ICs.
+ shift_x: 0. # A shift to apply to all particles read from the ICs (in internal units).
+ shift_y: 0.
+ shift_z: 0.
+
+# Parameters govering domain decomposition
+DomainDecomposition:
+ initial_type: m # The initial strategy ("g", "m", "w", or "v"). See documentation for details.
+ intitial_grid_x: 10 # Grid size if the 'g' strategy is chosen.
+ intitial_grid_y: 10
+ intitial_grid_z: 10
+ repartition_type: b # The re-decomposition strategy ("n", "b", "v", "e" or "x"). See documentation for details.
+