diff --git a/earthdiagnostics/ocean/moc.py b/earthdiagnostics/ocean/moc.py
index b6f845064a3d25970a6efbeb8a2302e81d99d899..f5931fc8b622bb7a0d9c40007cd8ecfe18ef91f3 100644
--- a/earthdiagnostics/ocean/moc.py
+++ b/earthdiagnostics/ocean/moc.py
@@ -91,6 +91,7 @@ class Moc(Diagnostic):
         if not basins:
             Log.error("Basins not recognized")
             return ()
+        basins.sort()
 
         job_list = list()
         for (
@@ -128,9 +129,9 @@ class Moc(Diagnostic):
         e3v = mesh.get_k_length(cell_point="V")
 
         masks = {}
-        self.basins.sort()
+        self.basins
         for basin in self.basins:
-            if basin == "Global":
+            if basin == Basins().Global:
                 global_mask = mesh.get_landsea_mask(cell_point="V")
                 global_mask[..., 0] = 0.0
                 global_mask[..., -1] = 0.0