replace manual netcdf4 engine logic with _open_dataset_with_fallback utility

Author: rajeeja

uxarray/core/api.py

@@ -111,12 +111,9 @@ def open_grid(
 
     else:
         # Attempt to use Xarray directly for remaining input types
-        # Force netCDF4 engine for .nc files to avoid scipy reader issues
-        if isinstance(grid_filename_or_obj, (str, os.PathLike)) and str(
-            grid_filename_or_obj
-        ).endswith(".nc"):
-            kwargs.setdefault("engine", "netcdf4")
-        grid_ds = xr.open_dataset(grid_filename_or_obj, chunks=grid_chunks, **kwargs)
+        grid_ds = _open_dataset_with_fallback(
+            grid_filename_or_obj, chunks=grid_chunks, **kwargs
+        )
         grid = Grid.from_dataset(grid_ds, use_dual=use_dual)
 
     # Return the grid (and chunks, if requested) in a consistent manner.