assert that λ and φ are monotonically increasing (#127)

COSIMA · Mar 19, 2024 · 9b5600e · 9b5600e
1 parent 89e57ce
commit 9b5600e
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/regional_mom6/regional_mom6.py b/regional_mom6/regional_mom6.py
@@ -410,6 +410,9 @@ def rectangular_hgrid(λ, φ):
         xarray.Dataset: An FMS-compatible ``hgrid`` that includes all required attributes.
     """
 
+    assert np.all(np.diff(λ) > 0), "longitudes array λ must be monotonically increasing"
+    assert np.all(np.diff(φ) > 0), "latitudes array φ must be monotonically increasing"
+
     R = 6371e3  # mean radius of the Earth; https://en.wikipedia.org/wiki/Earth_radius
 
     # compute longitude spacing and ensure that longitudes are uniformly spaced