Remove FieldSet.dimrange

Author: VeckoTheGecko

parcels/fieldset.py

@@ -82,34 +82,6 @@ def time_interval(self):
             return None
         return functools.reduce(lambda x, y: x.intersection(y), time_intervals)
 
-    def dimrange(self, dim):
-        """Returns maximum value of a dimension (lon, lat, depth or time)
-        on 'left' side and minimum value on 'right' side for all grids
-        in a gridset. Useful for finding e.g. longitude range that
-        overlaps on all grids in a gridset.
-        """
-        maxleft, minright = (-np.inf, np.inf)
-        dim2ds = {
-            "depth": ["nz1", "nz"],
-            "lat": ["node_lat", "face_lat", "edge_lat"],
-            "lon": ["node_lon", "face_lon", "edge_lon"],
-            "time": ["time"],
-        }
-        for ds in self.datasets:
-            for field in ds.data_vars:
-                for d in dim2ds[dim]:  # check all possible dimensions
-                    if d in ds[field].dims:
-                        if dim == "depth":
-                            maxleft = max(maxleft, ds[field][d].min().data)
-                            minright = min(minright, ds[field][d].max().data)
-                        else:
-                            maxleft = max(maxleft, ds[field][d].data[0])
-                            minright = min(minright, ds[field][d].data[-1])
-        maxleft = 0 if maxleft == -np.inf else maxleft  # if all len(dim) == 1
-        minright = 0 if minright == np.inf else minright  # if all len(dim) == 1
-
-        return maxleft, minright
-
     def add_field(self, field: Field, name: str | None = None):
         """Add a :class:`parcels.field.Field` object to the FieldSet.