diff --git a/tests/interpolor_particles_harmonic_average.prm b/tests/interpolor_particles_harmonic_average.prm
new file mode 100644
index 00000000000..a2edd2fa67d
--- /dev/null
+++ b/tests/interpolor_particles_harmonic_average.prm
@@ -0,0 +1,21 @@
+# Test whether the interpolator reads the input correctly.
+include $ASPECT_SOURCE_DIR/tests/time_stepping_repeat_particles.prm
+
+set End time                               = 0.0
+
+subsection Particles
+  set Interpolation scheme = harmonic average
+  set List of particle properties = initial position
+end
+
+subsection Compositional fields
+  set Number of fields = 0
+  set Names of fields =
+  set Compositional field methods =
+  set Mapped particle properties  =
+end
+
+
+subsection Initial composition model
+  set Model name = unspecified
+end
diff --git a/tests/interpolor_particles_harmonic_average/screen-output b/tests/interpolor_particles_harmonic_average/screen-output
new file mode 100644
index 00000000000..70364d6d64c
--- /dev/null
+++ b/tests/interpolor_particles_harmonic_average/screen-output
@@ -0,0 +1,27 @@
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------
+
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------
+Number of active cells: 128 (on 4 levels)
+Number of degrees of freedom: 1,836 (1,122+153+561)
+
+*** Timestep 0:  t=0 seconds, dt=0 seconds
+   Solving temperature system... 0 iterations.
+   Advecting particles...  done.
+   Solving Stokes system... 11+0 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-interpolor_particles_harmonic_average/solution/solution-00000
+     Writing particle output:  output-interpolor_particles_harmonic_average/particles/particles-00000
+
+Termination requested by criterion: end time
+
+
++----------------------------------------------+------------+------------+
++----------------------------------+-----------+------------+------------+
++----------------------------------+-----------+------------+------------+
+
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------
diff --git a/tests/interpolor_particles_nearest_neighbor.prm b/tests/interpolor_particles_nearest_neighbor.prm
new file mode 100644
index 00000000000..3938cab7dd8
--- /dev/null
+++ b/tests/interpolor_particles_nearest_neighbor.prm
@@ -0,0 +1,8 @@
+# Test whether the interpolator reads the input correctly.
+include $ASPECT_SOURCE_DIR/tests/time_stepping_repeat_particles.prm
+
+set End time                               = 0.0
+
+subsection Particles
+  set Interpolation scheme = nearest neighbor
+end
diff --git a/tests/interpolor_particles_nearest_neighbor/screen-output b/tests/interpolor_particles_nearest_neighbor/screen-output
new file mode 100644
index 00000000000..ffae82c5471
--- /dev/null
+++ b/tests/interpolor_particles_nearest_neighbor/screen-output
@@ -0,0 +1,28 @@
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------
+
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------
+Number of active cells: 128 (on 4 levels)
+Number of degrees of freedom: 2,958 (1,122+153+561+561+561)
+
+*** Timestep 0:  t=0 seconds, dt=0 seconds
+   Solving temperature system... 0 iterations.
+   Advecting particles...  done.
+   Solving field system ... 0 iterations.
+   Solving Stokes system... 11+0 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-interpolor_particles_nearest_neighbor/solution/solution-00000
+     Writing particle output:  output-interpolor_particles_nearest_neighbor/particles/particles-00000
+
+Termination requested by criterion: end time
+
+
++----------------------------------------------+------------+------------+
++----------------------------------+-----------+------------+------------+
++----------------------------------+-----------+------------+------------+
+
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------
diff --git a/tests/interpolor_particles_quadratic_least_squares.prm b/tests/interpolor_particles_quadratic_least_squares.prm
new file mode 100644
index 00000000000..96d437503aa
--- /dev/null
+++ b/tests/interpolor_particles_quadratic_least_squares.prm
@@ -0,0 +1,8 @@
+# Test whether the interpolator reads the input correctly.
+include $ASPECT_SOURCE_DIR/tests/time_stepping_repeat_particles.prm
+
+set End time                               = 0.0
+
+subsection Particles
+  set Interpolation scheme = quadratic least squares
+end
diff --git a/tests/interpolor_particles_quadratic_least_squares/screen-output b/tests/interpolor_particles_quadratic_least_squares/screen-output
new file mode 100644
index 00000000000..fbb13508c82
--- /dev/null
+++ b/tests/interpolor_particles_quadratic_least_squares/screen-output
@@ -0,0 +1,28 @@
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------
+
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------
+Number of active cells: 128 (on 4 levels)
+Number of degrees of freedom: 2,958 (1,122+153+561+561+561)
+
+*** Timestep 0:  t=0 seconds, dt=0 seconds
+   Solving temperature system... 0 iterations.
+   Advecting particles...  done.
+   Solving field system ... 0 iterations.
+   Solving Stokes system... 11+0 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-interpolor_particles_quadratic_least_squares/solution/solution-00000
+     Writing particle output:  output-interpolor_particles_quadratic_least_squares/particles/particles-00000
+
+Termination requested by criterion: end time
+
+
++----------------------------------------------+------------+------------+
++----------------------------------+-----------+------------+------------+
++----------------------------------+-----------+------------+------------+
+
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------