Added zarr support for input

tzis/tzis.py

@@ -152,7 +152,10 @@ class Tzis:
             if self.verbose:               
                 print("Calculating chunk size")
 
-            testds=xarray.open_dataset(mf[0])
+            if xarray_kwargs is not None:
+                testds=xarray.open_dataset(mf[0], **xarray_kwargs)
+            else:
+                testds=xarray.open_dataset(mf[0])
 
             if not varname:
                 self.varname = self._get_varname(testds)
@@ -195,7 +198,7 @@ class Tzis:
                 for var_id in mf_dset.variables:
                     mf_dset[var_id].unify_chunks()
 
-            conflict_attrs = self.get_conflict_attrs(mf, mf_dset)
+            conflict_attrs = self.get_conflict_attrs(mf, mf_dset, xarray_kwargs)
 
             self.init_verify_provenance(reset=True)
 
@@ -234,13 +237,17 @@ class Tzis:
             print("You did not provide anything to open.")
         return None
 
-    def get_conflict_attrs(self, mf, mf_dset):
+    def get_conflict_attrs(self, mf, mf_dset, xarray_kwargs):
         conflict_attrs = {}
         # try:
         maxdigits = len(str(len(mf)))
         digitstring = "{0:0" + str(maxdigits) + "}"
         for fileno, dset in enumerate(mf):
-            dset_attrs = xarray.open_dataset(dset).attrs
+            if xarray_kwargs is not None:
+                ds=xarray.open_dataset(dset,**xarray_kwargs)
+            else:
+                ds=xarray.open_dataset(dset)
+            dset_attrs = ds.attrs
             missing_attrs = {
                 k: v for k, v in dset_attrs.items() if k not in mf_dset.attrs
             }