diff --git a/tests/fixtures/sample_rules.py b/tests/fixtures/sample_rules.py
index 3bd3cb2..3b30cf3 100644
--- a/tests/fixtures/sample_rules.py
+++ b/tests/fixtures/sample_rules.py
@@ -4,6 +4,27 @@
 from pymorize.rule import Rule
 
 
+@pytest.fixture
+def fesom_2p6_esmtools_temp_rule(fesom_2p6_pimesh_esm_tools_data):
+    return Rule.from_dict(
+        {
+            "name": "temp",
+            "experiment_id": "piControl",
+            "output_directory": "./output",
+            "source_id": "FESOM",
+            "variant_label": "r1i1p1f1",
+            "inputs": [
+                {
+                    "path": fesom_2p6_pimesh_esm_tools_data / "outdata/fesom",
+                    "pattern": "temp.fesom..*.nc",
+                },
+            ],
+            "cmor_variable": "thetao",
+            "model_variable": "temp",
+        }
+    )
+
+
 @pytest.fixture
 def pi_uxarray_temp_rule(pi_uxarray_data):
     return Rule.from_dict(
diff --git a/tests/unit/test_gather_inputs.py b/tests/unit/test_gather_inputs.py
index 652be24..9266227 100644
--- a/tests/unit/test_gather_inputs.py
+++ b/tests/unit/test_gather_inputs.py
@@ -7,3 +7,9 @@ def test_load_mfdataset_pi_uxarray(pi_uxarray_temp_rule):
     data = load_mfdataset(None, pi_uxarray_temp_rule)
     # Check if load worked correctly and we got back a Dataset
     assert isinstance(data, xr.Dataset)
+
+
+def test_load_mfdataset_fesom_2p6_esmtools(fesom_2p6_esmtools_temp_rule):
+    data = load_mfdataset(None, fesom_2p6_esmtools_temp_rule)
+    # Check if load worked correctly and we got back a Dataset
+    assert isinstance(data, xr.Dataset)