Update create_regular methods to remove cellsize divisibility check

cf/dimensioncoordinate.py

@@ -252,10 +252,7 @@ def create_regular(cls, args, units=None, standard_name=None, bounds=True):
         :Parameters:
 
             args: sequence of numbers
-                A sequence of three numeric values. The first two values in the
-                sequence represent the coordinate range (see the bounds
-                parameter for details), and the third value represents the
-                cellsize.
+                {{regular args}}
 
             bounds: `bool`, optional
                 If True (the default) then the given range represents
@@ -311,12 +308,6 @@ def create_regular(cls, args, units=None, standard_name=None, bounds=True):
                 f"negative cellsize ({cellsize})"
             )
 
-        if range_diff % cellsize != 0:
-            raise ValueError(
-                f"The range of the dimension ({range_diff}) must be "
-                f"divisible by the cellsize ({cellsize})"
-            )
-
         if standard_name is not None and not isinstance(standard_name, str):
             raise ValueError("standard_name must be either None or a string.")
 

cf/docstring/docstring.py

@@ -513,6 +513,18 @@
                 specified with or without the ``${...}`` syntax. For
                 instance, the following are equivalent: ``'base'`` and
                 ``'${base}'``.""",
+    # regular args
+    "{{regular args}}": 
+                """A sequence of three numeric values. The first two values in 
+                the sequence represent the coordinate range (see the bounds
+                parameter for details), and the third value represents the
+                cellsize.
+
+                .. note:: The cellsize does not have to explicitly divide into 
+                          the range of the given dimension. But as it follows 
+                          `numpy.arange` while creating the points, one should
+                          verify that that the number of grid points are 
+                          returned as expected.""",  
     # ----------------------------------------------------------------
     # Method description substitutions (4 levels of indentation)
     # ----------------------------------------------------------------

cf/domain.py

@@ -343,16 +343,10 @@ def create_regular(cls, x_args, y_args, bounds=True):
         :Parameters:
 
             x_args: sequence of numbers
-                A sequence of three numeric values. The first two values in the
-                sequence represent the coordinate range (see the bounds
-                parameter for details), and the third value represents the
-                cellsize.
+                {{regular args}}
 
             y_args: sequence of numbers
-                A sequence of three numeric values. The first two values in the
-                sequence represent the coordinate range (see the bounds
-                parameter for details), and the third value represents the
-                cellsize.
+                {{regular args}}
 
             bounds: `bool`, optional
                 If True (default), bounds will be created

cf/test/test_DimensionCoordinate.py

@@ -658,14 +658,6 @@ def test_DimensionCoordinate_create_regular(self):
                 standard_name="longitude",
             )
 
-        # Cellsize not divisible by range
-        with self.assertRaises(ValueError):
-            cf.DimensionCoordinate.create_regular(
-                (-180, 180, 50),
-                units="degrees_east",
-                standard_name="longitude",
-            )
-
         # Test decreasing case
         longitude_decreasing = cf.DimensionCoordinate.create_regular(
             (180, -180, -1), units="degrees_east", standard_name="longitude"

cf/test/test_Domain.py

@@ -322,10 +322,6 @@ def test_Domain_create_regular(self):
         self.assertTrue(np.allclose(longitude.array, x_points))
         self.assertTrue(np.allclose(latitude.array, y_points))
 
-        # Test dx and dy not divisors of the range
-        with self.assertRaises(ValueError):
-            cf.Domain.create_regular((-180, 180, 61), (-90, 90, 46))
-
         # Test if range difference in x_range is greater than 360
         with self.assertRaises(ValueError):
             cf.Domain.create_regular((-180, 190, 1), (-90, 90, 1))
@@ -366,6 +362,35 @@ def test_Domain_create_regular(self):
             np.allclose(latitude_no_bounds.array, y_points_no_bounds)
         )
 
+        # Test for the given specific domain
+        ymin, ymax, dy = 45.0, 90.0, 0.0083333
+        xmin, xmax, dx = 250.0, 360.0, 0.0083333
+
+        domain_specific = cf.Domain.create_regular(
+            (xmin, xmax, dx), (ymin, ymax, dy)
+        )
+        self.assertIsInstance(domain_specific, cf.Domain)
+
+        x_bounds_specific = np.arange(xmin, xmax + dx, dx)
+        y_bounds_specific = np.arange(ymin, ymax + dy, dy)
+
+        x_points_specific = (
+            x_bounds_specific[:-1] + x_bounds_specific[1:]
+        ) / 2
+        y_points_specific = (
+            y_bounds_specific[:-1] + y_bounds_specific[1:]
+        ) / 2
+
+        longitude_specific = domain_specific.construct("longitude")
+        latitude_specific = domain_specific.construct("latitude")
+
+        self.assertTrue(
+            np.allclose(longitude_specific.array - x_points_specific, 0)
+        )
+        self.assertTrue(
+            np.allclose(latitude_specific.array - y_points_specific, 0)
+        )
+
 
 if __name__ == "__main__":
     print("Run date:", datetime.datetime.now())