diff --git a/docs/examples/10_mesh2d_global_grid.ipynb b/docs/examples/10_mesh2d_global_grid.ipynb
index 9a23dfba..fcacb5cb 100644
--- a/docs/examples/10_mesh2d_global_grid.ipynb
+++ b/docs/examples/10_mesh2d_global_grid.ipynb
@@ -35,7 +35,7 @@
    "outputs": [],
    "source": [
     "mk = meshkernel.MeshKernel(projection=meshkernel.ProjectionType.SPHERICAL)\n",
-    "mk.mesh2d_make_global(num_longitude_nodes=192, num_latitude_nodes=100)\n",
+    "mk.mesh2d_make_global(num_longitude_nodes=192)\n",
     "mesh2d = mk.mesh2d_get()"
    ]
   },
diff --git a/meshkernel/meshkernel.py b/meshkernel/meshkernel.py
index 7d7be7a8..0c2a53f0 100644
--- a/meshkernel/meshkernel.py
+++ b/meshkernel/meshkernel.py
@@ -643,21 +643,20 @@ def mesh2d_delete_hanging_edges(self) -> None:
             self.lib.mkernel_mesh2d_delete_hanging_edges, self._meshkernelid
         )
 
-    def mesh2d_make_global(
-        self, num_longitude_nodes: int, num_latitude_nodes: int
-    ) -> None:
-        """Compute the global mesh with a given number of points along the longitude and latitude directions.
+    def mesh2d_make_global(self, num_longitude_nodes: int) -> None:
+        """Compute the global mesh with a given number of points along the longitude direction.
+        The number of points in the latitude direction will be automatically calculated, such
+        that the resulting cells are as square as possible and the resulting mesh extends up
+        to the poles.
 
         Args:
         num_longitude_nodes (int): The number of points along the longitude.
-        num_latitude_nodes (int):  The number of points along the latitude (half hemisphere)
         """
 
         self._execute_function(
             self.lib.mkernel_mesh2d_make_global,
             self._meshkernelid,
             c_int(num_longitude_nodes),
-            c_int(num_latitude_nodes),
         )
 
     def mesh2d_make_triangular_mesh_from_polygon(
diff --git a/tests/test_mesh2d_basics.py b/tests/test_mesh2d_basics.py
index fadcdb0b..cbf00943 100644
--- a/tests/test_mesh2d_basics.py
+++ b/tests/test_mesh2d_basics.py
@@ -683,9 +683,8 @@ def test_mesh2d_make_global():
     mk = MeshKernel(ProjectionType.SPHERICAL)
 
     num_longitude_nodes = 19
-    num_latitude_nodes = 25
 
-    mk.mesh2d_make_global(num_longitude_nodes, num_latitude_nodes)
+    mk.mesh2d_make_global(num_longitude_nodes)
     mesh2d = mk.mesh2d_get()
 
     assert mesh2d.edge_x.size == 1225
@@ -698,10 +697,9 @@ def test_mesh2d_make_global_with_cartesian_coordinates_should_throw():
     mk = MeshKernel(ProjectionType.CARTESIAN)
 
     num_longitude_nodes = 19
-    num_latitude_nodes = 25
 
     with pytest.raises(MeshKernelError):
-        mk.mesh2d_make_global(num_longitude_nodes, num_latitude_nodes)
+        mk.mesh2d_make_global(num_longitude_nodes)
 
 
 def test_mesh2d_make_triangular_mesh_from_polygon():