del bounds attr from coord in drop_attributes() (#733)

* del bounds attr from coord in drop_attributes()

* add TODO statement

src/preprocessor.py

@@ -849,9 +849,14 @@ def drop_attributes(self, xr_ds: xr.Dataset) -> xr.Dataset:
                      'average_T1',
                      'height',
                      'date']
+        # TODO: find a suitable answer to conflicts in xarray merging (i.e. nctoolkit)
         for att in drop_atts:
             if xr_ds.get(att, None) is not None:
                 xr_ds = xr_ds.drop_vars(att)
+                for coord in xr_ds.coords:
+                    if 'bounds' in xr_ds[coord].attrs:
+                        if xr_ds[coord].attrs['bounds'] == att:
+                            del xr_ds[coord].attrs['bounds']
         return xr_ds
 
     def check_multichunk(self, group_df: pd.DataFrame, case_dr, log) -> pd.DataFrame:

src/xr_parser.py

@@ -1072,7 +1072,7 @@ def reconcile_coord_bounds(self, our_coord, ds, ds_coord_name):
         expectations based on the model's convention (*our_var*), for the bounds
         on the dimension coordinate *our_coord*.
         """
-        if len(ds.cf.bounds) > 0:
+        if len(ds.cf.bounds.get(ds_coord_name, [])) > 0:
             bounds = ds.cf.get_bounds(ds_coord_name)
         elif hasattr(ds[ds_coord_name], 'attrs'):
             if ds[ds_coord_name].attrs.get('bounds', None):