diff --git a/BathymetryReco/Project.toml b/BathymetryReco/Project.toml
index a89cfac4a7b1c12a0f04a29a8a518c972b14d2b2..ab82890321c8c07a2823fabde14a8d6144b56b99 100644
--- a/BathymetryReco/Project.toml
+++ b/BathymetryReco/Project.toml
@@ -9,6 +9,7 @@ HDF5 = "f67ccb44-e63f-5c2f-98bd-6dc0ccc4ba2f"
 ProgressBars = "49802e3a-d2f1-5c88-81d8-b72133a6f568"
 PyCall = "438e738f-606a-5dbb-bf0a-cddfbfd45ab0"
 TOML = "fa267f1f-6049-4f14-aa54-33bafae1ed76"
+Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [compat]
 Distributions = "0.25.117"
@@ -16,3 +17,4 @@ HDF5 = "0.17.2"
 ProgressBars = "1.5.1"
 PyCall = "1.96.4"
 TOML = "1.0.3"
+Test = "1.11.0"
diff --git a/BathymetryReco/test/runtests.jl b/BathymetryReco/test/runtests.jl
new file mode 100644
index 0000000000000000000000000000000000000000..04b16cfd8d728a6b66fd4a6d0879c41786170ae6
--- /dev/null
+++ b/BathymetryReco/test/runtests.jl
@@ -0,0 +1,10 @@
+using Test
+using BathymetryReco
+
+@testset "Test bathymetry reconstruction" begin
+    # Load the observation data
+    obs_data, b = load_observation("../../data/toy_measurement/waterchannel_exact_bathy")
+
+    @test size(obs_data.H) == (length(obs_data.t), length(obs_data.x))
+    @test obs_data.tstart == 32.0
+end
\ No newline at end of file