diff --git a/tests/test_bathymetric_gradient_function.py b/tests/test_bathymetric_gradient_function.py
index 03319f7..0bb0d72 100644
--- a/tests/test_bathymetric_gradient_function.py
+++ b/tests/test_bathymetric_gradient_function.py
@@ -54,7 +54,7 @@ def test_bathymetric_gradient_function():
     EPSG = 4326  # EPSG:4326 or WGS84
     bbox = (-74.4, -73.4, 40.2, 41.2)
     extent = om.Region(extent=bbox, crs=EPSG)
-    # dem = om.DEM(fdem, crs=4326)
+    dem = om.DEM(fdem, crs=4326)
 
     min_edge_length = 0.0025  # minimum mesh size in domain in projection
     max_edge_length = 0.10  # maximum mesh size in domain in projection