diff --git a/compareECLFiles.cmake b/compareECLFiles.cmake
index 45d68acca..c5985c19d 100755
--- a/compareECLFiles.cmake
+++ b/compareECLFiles.cmake
@@ -283,6 +283,22 @@ add_test_compareECLFiles(CASENAME fetkovich_2d
                          REL_TOL ${rel_tol}
                          DIR aquifer-fetkovich)
 
+add_test_compareECLFiles(CASENAME numerical_aquifer_3d_2aqu
+                         FILENAME 3D_2AQU_NUM
+                         SIMULATOR flow
+                         ABS_TOL ${abs_tol}
+                         REL_TOL ${rel_tol}
+                         DIR aquifer-num
+                         TEST_ARGS --enable-tuning=true)
+
+add_test_compareECLFiles(CASENAME numerical_aquifer_3d_1aqu
+                         FILENAME 3D_1AQU_3CELLS
+                         SIMULATOR flow
+                         ABS_TOL ${abs_tol}
+                         REL_TOL ${rel_tol}
+                         DIR aquifer-num
+                         TEST_ARGS --enable-tuning=true)
+
 add_test_compareECLFiles(CASENAME spe3
                          FILENAME SPE3CASE1
                          SIMULATOR flow
@@ -991,4 +1007,20 @@ if(MPI_FOUND)
                                        REL_TOL ${rel_tol_parallel}
                                        DIR parallel_fieldprops
                                        TEST_ARGS --linear-solver-reduction=1e-7 --tolerance-cnv=5e-6 --tolerance-mb=1e-6)
+
+  add_test_compare_parallel_simulation(CASENAME numerical_aquifer_3d_2aqu
+                                       FILENAME 3D_2AQU_NUM
+                                       SIMULATOR flow
+                                       ABS_TOL ${abs_tol_parallel}
+                                       REL_TOL ${rel_tol_parallel}
+                                       DIR aquifer-num
+                                       TEST_ARGS --enable-tuning=true)
+
+  add_test_compare_parallel_simulation(CASENAME numerical_aquifer_3d_1aqu
+                                       FILENAME 3D_1AQU_3CELLS
+                                       SIMULATOR flow
+                                       ABS_TOL ${abs_tol_parallel}
+                                       REL_TOL ${rel_tol_parallel}
+                                       DIR aquifer-num
+                                       TEST_ARGS --enable-tuning=true)
 endif()
diff --git a/tests/update_reference_data.sh b/tests/update_reference_data.sh
index 1ac622474..9329c9676 100755
--- a/tests/update_reference_data.sh
+++ b/tests/update_reference_data.sh
@@ -63,6 +63,8 @@ tests[spe1_brine]="flow spe1_brine SPE1CASE1_BRINE"
 tests[spe1_water]="flow_onephase spe1 SPE1CASE1_WATER"
 tests[ctaquifer_2d_oilwater]="flow aquifer-oilwater 2D_OW_CTAQUIFER"
 tests[fetkovich_2d]="flow aquifer-fetkovich 2D_FETKOVICHAQUIFER"
+tests[numerical_aquifer_3d_1aqu]="flow aquifer-num 3D_1AQU_3CELLS"
+tests[numerical_aquifer_3d_2aqu]="flow aquifer-num 3D_2AQU_NUM"
 tests[msw_2d_h]="flow msw_2d_h 2D_H__"
 tests[msw_3d_hfa]="flow msw_3d_hfa 3D_MSW"
 tests[polymer_oilwater]="flow polymer_oilwater 2D_OILWATER_POLYMER"