diff --git a/server/map/__pycache__/get_relation.cpython-311.pyc b/server/map/__pycache__/get_relation.cpython-311.pyc
index fb8ec21176d8e459a50cb57a1e64c0b8458ef3d8..61a991605c9616d12954c7a67b412e78d41c1477 100644
Binary files a/server/map/__pycache__/get_relation.cpython-311.pyc and b/server/map/__pycache__/get_relation.cpython-311.pyc differ
diff --git a/server/map/get_relation.py b/server/map/get_relation.py
index 7b411154a8030113cb9d0867e45659d316c4fba2..0ca4075bf6f793c10b8e017498986b78b2f1989c 100644
--- a/server/map/get_relation.py
+++ b/server/map/get_relation.py
@@ -33,7 +33,10 @@ def get_relation(self, body_of_water: str):
         #hrz_lines = grid_lines[0]  # Horizontal coordinates
         #vrt_lines = grid_lines[1]  # Vertical coordinates
 
-        hrz_lines = [60.7530, 60.5931]  # Horizontal coordinates
+        hrz_lines = grid_lines[1]  # Horizontal coordinates
+        for line in hrz_lines:
+            print("line: ", line)
+
         vrt_lines = [11.0413]  # Vertical coordinates
 
         # Cut polygon into horizontal sections
@@ -49,6 +52,8 @@ def get_relation(self, body_of_water: str):
 
         divided_map.append(polygon)
 
+        break
+
     tiles = gpd.GeoDataFrame(geometry=divided_map)
 
     # NB test plot
@@ -88,12 +93,12 @@ def cut_polygon_in_two(polygon: Polygon, divisor: float, horizontal: bool):
     for point in exterior_coords:
         point = Point(point)  # Convert coordinates to Shapely Point object
         if horizontal:  # Horizontal split
-            if point.y > divisor:
+            if point.y < divisor:
                 split_shape.append(point)
             else:
                 remaining_shape.append(point)
         else:  # Vertical split
-            if point.x > divisor:
+            if point.x < divisor:
                 split_shape.append(point)
             else:
                 remaining_shape.append(point)
@@ -125,8 +130,12 @@ def create_grid_coords(polygon: Polygon, cell_size: float):
     x_coords = np.arange(min_x, max_x, cell_size)
     y_coords = np.arange(min_y, max_y, cell_size)
 
-    if len(x_coords) == 0 or len(y_coords) == 0:
+    # Round the coordinates to 4 decimals
+    x_coords_rounded = np.around(x_coords, decimals=4)
+    y_coords_rounded = np.around(y_coords, decimals=4)
+
+    if len(x_coords_rounded) == 0 or len(y_coords_rounded) == 0:
         raise ValueError("No grid points generated")
 
     # Return tuple of list of x coordinates and list of y coordinates
-    return x_coords, y_coords
+    return x_coords_rounded, y_coords_rounded