diff --git a/testing/CMakeLists.txt b/testing/CMakeLists.txt
index 008e87a4511f33e0da62adf7447140982e8abbd2..c854ff564b1c1328256ee1bb581f30ca0cc9a004 100644
--- a/testing/CMakeLists.txt
+++ b/testing/CMakeLists.txt
@@ -13,7 +13,8 @@ endfunction()
 
 # dorie run
 #dorie_add_system_test(dorie ode.mini)
-dorie_add_system_test(dorie ode_higherorder.mini)
+#dorie_add_system_test(dorie ode_higherorder.mini)
+dorie_add_system_test(dorie ode_higherorder_subsampling.mini)
 # dorie_add_system_test(dorie ode_evaporation.mini)
 # dorie_add_system_test(dorie ode_linear_interpolator.mini)
 # dorie_add_system_test(dorie muphi.mini)
diff --git a/testing/ode_higherorder_subsampling.mini.in b/testing/ode_higherorder_subsampling.mini.in
new file mode 100644
index 0000000000000000000000000000000000000000..d368e3f6d861a513fec2ebcc63e6b3b1bae4ca98
--- /dev/null
+++ b/testing/ode_higherorder_subsampling.mini.in
@@ -0,0 +1,39 @@
+include ${CMAKE_BINARY_DIR}/doc/default_files/config.ini
+
+__name = ode_higherorder
+_test_command = run
+_asset_path = "${CMAKE_CURRENT_LIST_DIR}"
+_evaluation = ode
+
+output.fileName = ode_higherorder_subsampling | unique
+output.outputPath = ode_higherorder_subsampling | unique
+output.verbose = 0
+output.subsamplingLevel = 0, 1, 2 | expand
+
+time.end = 1E8
+time.maxTimestep = 1E7
+time.startTimestep = 1E4
+
+adaptivity.threshold = 1E-8
+adaptivity.useAdaptivity = false, true | expand
+adaptivity.maxLevel = 3
+
+grid.FEorder = 1, 2, 3 | expand
+grid.gridType = rectangular
+# higher FEorders seem to have a problem with asymmetric grids
+grid.cells = 1 50
+grid.initialLevel = 1
+grid.gridFile = "{_asset_path}/meshes/square.msh"
+
+boundary.file = "{_asset_path}/bcs/infiltration_2d.dat"
+
+parameters.arrayFile = "{_asset_path}/parfields/sand.h5"
+
+[_ode]
+flux = -3e-6
+head_abstol = 1E-5
+flux_abstol = 2E-8
+
+[NewtonParameters]
+AbsoluteLimit = 1E-10
+Reduction = 1E-10