Remove netcdf_engine from Field and FieldSet

Author: VeckoTheGecko

parcels/field.py

@@ -226,7 +226,6 @@ def __init__(
 
         self._dimensions = kwargs.pop("dimensions", None)
         self._dataFiles = kwargs.pop("dataFiles", None)
-        self._netcdf_engine = kwargs.pop("netcdf_engine", "netcdf4")
         self._creation_log = kwargs.pop("creation_log", "")
 
         # data_full_zdim is the vertical dimension of the complete field data, ignoring the indices.
@@ -276,10 +275,6 @@ def interp_method(self, value):
     def gridindexingtype(self):
         return self._gridindexingtype
 
-    @property
-    def netcdf_engine(self):
-        return self._netcdf_engine
-
     @classmethod
     def _get_dim_filenames(cls, filenames, dim):
         if isinstance(filenames, str) or not isinstance(filenames, collections.abc.Iterable):
@@ -295,11 +290,11 @@ def _get_dim_filenames(cls, filenames, dim):
             return filenames
 
     @staticmethod
-    def _collect_timeslices(data_filenames, dimensions, indices, netcdf_engine):
+    def _collect_timeslices(data_filenames, dimensions, indices):
         timeslices = []
         dataFiles = []
         for fname in data_filenames:
-            with NetcdfFileBuffer(fname, dimensions, indices, netcdf_engine=netcdf_engine) as filebuffer:
+            with NetcdfFileBuffer(fname, dimensions, indices) as filebuffer:
                 ftime = filebuffer.time
                 timeslices.append(ftime)
                 dataFiles.append([fname] * len(ftime))
@@ -385,7 +380,6 @@ def from_netcdf(
                 raise NotImplementedError("Vertically adaptive meshes not implemented for from_netcdf()")
             depth_filename = depth_filename[0]
 
-        netcdf_engine = kwargs.pop("netcdf_engine", "netcdf4")
         gridindexingtype = kwargs.get("gridindexingtype", "nemo")
 
         indices = {}  # TODO Nick: Cleanup
@@ -440,9 +434,7 @@ def from_netcdf(
             # Concatenate time variable to determine overall dimension
             # across multiple files
             if "time" in dimensions:
-                time, time_origin, timeslices, dataFiles = cls._collect_timeslices(
-                    data_filenames, dimensions, indices, netcdf_engine
-                )
+                time, time_origin, timeslices, dataFiles = cls._collect_timeslices(data_filenames, dimensions, indices)
                 grid = Grid.create_grid(lon, lat, depth, time, time_origin=time_origin, mesh=mesh)
                 kwargs["dataFiles"] = dataFiles
             else:  # e.g. for the CROCO CS_w field, see https://github.com/OceanParcels/Parcels/issues/1831
@@ -451,7 +443,7 @@ def from_netcdf(
         elif grid is not None and ("dataFiles" not in kwargs or kwargs["dataFiles"] is None):
             # ==== means: the field has a shared grid, but may have different data files, so we need to collect the
             # ==== correct file time series again.
-            _, _, _, dataFiles = cls._collect_timeslices(data_filenames, dimensions, indices, netcdf_engine)
+            _, _, _, dataFiles = cls._collect_timeslices(data_filenames, dimensions, indices)
             kwargs["dataFiles"] = dataFiles
 
         if "time" in indices:
@@ -486,7 +478,6 @@ def from_netcdf(
             allow_time_extrapolation = False if "time" in dimensions else True
 
         kwargs["dimensions"] = dimensions.copy()
-        kwargs["netcdf_engine"] = netcdf_engine
 
         return cls(
             variable,

parcels/fieldfilebuffer.py

@@ -19,7 +19,6 @@ def __init__(
         interp_method: InterpMethodOption = "linear",
         data_full_zdim=None,
         gridindexingtype="nemo",
-        netcdf_engine="netcdf4",
     ):
         self.filename: PathLike | list[PathLike] = filename
         self.dimensions = dimensions  # Dict with dimension keys for file data
@@ -28,10 +27,9 @@ def __init__(
         self.interp_method = interp_method
         self.gridindexingtype = gridindexingtype
         self.data_full_zdim = data_full_zdim
-        self.netcdf_engine = netcdf_engine
 
     def __enter__(self):
-        self.dataset = open_xarray_dataset(self.filename, self.netcdf_engine)
+        self.dataset = open_xarray_dataset(self.filename)
         return self
 
     def __exit__(self, type, value, traceback):
@@ -159,12 +157,12 @@ def time_access(self):
         return time
 
 
-def open_xarray_dataset(filename: Path | str, netcdf_engine: str) -> xr.Dataset:
+def open_xarray_dataset(filename: Path | str) -> xr.Dataset:
     try:
         # Unfortunately we need to do if-else here, cause the lock-parameter is either False or a Lock-object
         # (which we would rather want to have being auto-managed).
         # If 'lock' is not specified, the Lock-object is auto-created and managed by xarray internally.
-        ds = xr.open_mfdataset(filename, decode_cf=True, engine=netcdf_engine)
+        ds = xr.open_mfdataset(filename, decode_cf=True)
         ds["decoded"] = True
     except:
         warnings.warn(  # TODO: Is this warning necessary? What cases does this except block get triggered - is it to do with the bare except???
@@ -174,7 +172,7 @@ def open_xarray_dataset(filename: Path | str, netcdf_engine: str) -> xr.Dataset:
             stacklevel=2,
         )
 
-        ds = xr.open_mfdataset(filename, decode_cf=False, engine=netcdf_engine)
+        ds = xr.open_mfdataset(filename, decode_cf=False)
         ds["decoded"] = False
     return ds
 

parcels/fieldset.py

@@ -356,9 +356,6 @@ def from_netcdf(
         gridindexingtype : str
             The type of gridindexing. Either 'nemo' (default), 'mitgcm', 'mom5', 'pop', or 'croco' are supported.
             See also the Grid indexing documentation on oceanparcels.org
-        netcdf_engine :
-            engine to use for netcdf reading in xarray. Default is 'netcdf',
-            but in cases where this doesn't work, setting netcdf_engine='scipy' could help. Accepted options are the same as the ``engine`` parameter in ``xarray.open_dataset()``.
         **kwargs :
             Keyword arguments passed to the :class:`parcels.Field` constructor.