diff --git a/demos/reanalysis-forced.ipynb b/demos/reanalysis-forced.ipynb
index a4d161f2..b9b3dc04 100644
--- a/demos/reanalysis-forced.ipynb
+++ b/demos/reanalysis-forced.ipynb
@@ -191,7 +191,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "expt.get_glorys_rectangular(\n",
+    "expt.get_glorys(\n",
     "    raw_boundaries_path=glorys_path\n",
     ")"
    ]
diff --git a/regional_mom6/regional_mom6.py b/regional_mom6/regional_mom6.py
index fc666732..9cde4552 100644
--- a/regional_mom6/regional_mom6.py
+++ b/regional_mom6/regional_mom6.py
@@ -1486,7 +1486,7 @@ def setup_initial_condition(
 
         return
 
-    def get_glorys_rectangular(self, raw_boundaries_path):
+    def get_glorys(self, raw_boundaries_path):
         """
         This function is a wrapper for `get_glorys_data`, calling this function once for each of the rectangular boundary segments and the initial condition. For more complex boundary shapes, call `get_glorys_data` directly for each of your boundaries that aren't parallel to lines of constant latitude or longitude. For example, for an angled Northern boundary that spans multiple latitudes, you'll need to download a wider rectangle containing the entire boundary.