revert to open_dataset from open_mfdataset

to avoid TypeError: The dtype and signature arguments to ufuncs only select the general DType and not details such as the byte order or time unit. You can avoid this error by using the scalar types np.float64 or the dtype string notation. xarray 2022.3.0 was okay with how I handled xarray DataArrays with units (pint quantities), but xarray 2023.5.0 is not okay.
dtcenter · Aug 24, 2023 · 97a2c8e · 97a2c8e
1 parent e2d5693
commit 97a2c8e
Show file tree

Hide file tree

Showing 3 changed files with 9 additions and 9 deletions.
diff --git a/metplotpy/contributed/fv3_physics_tend/cross_section_vert.py b/metplotpy/contributed/fv3_physics_tend/cross_section_vert.py
@@ -103,18 +103,18 @@ def main():
 
     # Read lat/lon from gfile
     logging.debug(f"read lat/lon from {gfile}")
-    gds = xarray.open_mfdataset(gfile)
+    gds = xarray.open_dataset(gfile)
     lont = gds[fv3["lon_name"]]
     latt = gds[fv3["lat_name"]]
 
     # Open input file
     logging.debug("open %s", ifile)
-    fv3ds = xarray.open_mfdataset(ifile)
+    fv3ds = xarray.open_dataset(ifile)
 
     if subtract:
         logging.info("subtracting %s", subtract)
         with xarray.set_options(keep_attrs=True):
-            fv3ds -= xarray.open_mfdataset(subtract)
+            fv3ds -= xarray.open_dataset(subtract)
 
     datetimeindex = fv3ds.indexes['time']
     if hasattr(datetimeindex, "to_datetimeindex"):

diff --git a/metplotpy/contributed/fv3_physics_tend/planview_fv3.py b/metplotpy/contributed/fv3_physics_tend/planview_fv3.py
@@ -110,19 +110,19 @@ def main():
 
     # Read lat/lon/area from gfile
     logging.debug(f"read lat/lon/area from {gfile}")
-    gds = xarray.open_mfdataset(gfile)
+    gds = xarray.open_dataset(gfile)
     lont = gds[fv3["lon_name"]]
     latt = gds[fv3["lat_name"]]
     area = gds["area"]
 
     # Open input file
     logging.debug("open %s", ifile)
-    fv3ds = xarray.open_mfdataset(ifile)
+    fv3ds = xarray.open_dataset(ifile)
 
     if subtract:
         logging.info("subtracting %s", subtract)
         with xarray.set_options(keep_attrs=True):
-            fv3ds -= xarray.open_mfdataset(subtract)
+            fv3ds -= xarray.open_dataset(subtract)
 
     datetimeindex = fv3ds.indexes['time']
     if hasattr(datetimeindex, "to_datetimeindex"):

diff --git a/metplotpy/contributed/fv3_physics_tend/vert_profile_fv3.py b/metplotpy/contributed/fv3_physics_tend/vert_profile_fv3.py
@@ -97,19 +97,19 @@ def main():
 
     # Read lat/lon/area from gfile
     logging.debug(f"read lat/lon/area from {gfile}")
-    gds = xarray.open_mfdataset(gfile)
+    gds = xarray.open_dataset(gfile)
     lont = gds[fv3["lon_name"]]
     latt = gds[fv3["lat_name"]]
     area = gds["area"]
 
     # Open input file
     logging.debug("open %s", ifile)
-    fv3ds = xarray.open_mfdataset(ifile)
+    fv3ds = xarray.open_dataset(ifile)
 
     if subtract:
         logging.info("subtracting %s", subtract)
         with xarray.set_options(keep_attrs=True):
-            fv3ds -= xarray.open_mfdataset(subtract)
+            fv3ds -= xarray.open_dataset(subtract)
 
     datetimeindex = fv3ds.indexes['time']
     if hasattr(datetimeindex, "to_datetimeindex"):