diff --git a/metplotpy/contributed/fv3_physics_tend/cross_section_vert.py b/metplotpy/contributed/fv3_physics_tend/cross_section_vert.py
index 7f922056..2a4a0501 100644
--- 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
index 8cfec71e..4239fad9 100644
--- 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
index 7c1170b5..be20359f 100644
--- 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"):