diff --git a/CHANGELOG.md b/CHANGELOG.md
index e6b3089..0a3adef 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -12,10 +12,11 @@ This project adheres to [Semantic Versioning](http://semver.org/).
   - Added Vertical TEC Instrument
   - Added documentation
 - Bug Fix
-  - Updated madrigal methods to simplify compound data types and enable
-    creation of netCDF4 files using `Instrument.to_netcdf4()`.
+  - Updated Madrigal methods to simplify compound data types and enable
+    creation of netCDF4 files using `Instrument.to_netcdf4()`
   - Updated load for multiple files in pandas format
   - Fixed remote listing routine to return filenames instead of experiments
+  - Fixed bug introduced by change in xarray requiring engine kwarg
 
 ## [0.0.3] - 2020-06-15
 - pypi compatibility
diff --git a/pysatMadrigal/instruments/methods/general.py b/pysatMadrigal/instruments/methods/general.py
index a09f006..17c115f 100644
--- a/pysatMadrigal/instruments/methods/general.py
+++ b/pysatMadrigal/instruments/methods/general.py
@@ -91,9 +91,10 @@ def load(fnames, tag=None, inst_id=None, xarray_coords=[], file_type='hdf5'):
     if file_type == "netCDF4":
         # Xarray natively opens netCDF data into a Dataset
         if len(fnames) == 1:
-            file_data = xr.open_dataset(fnames[0])
+            file_data = xr.open_dataset(fnames[0], engine="netcdf4")
         else:
-            file_data = xr.open_mfdataset(fnames, combine='by_coords')
+            file_data = xr.open_mfdataset(fnames, combine='by_coords',
+                                          engine="netcdf4")
 
         # Currently not saving file header data, as all metadata is at
         # the data variable level