diff --git a/mvcm/main.py b/mvcm/main.py
index a93981663df116a3422da0e6b5f5102d1732f528..44412e025fc404fa6fa13218d2220963a92c4898 100644
--- a/mvcm/main.py
+++ b/mvcm/main.py
@@ -581,8 +581,8 @@ def main(
         sf = {k: {"dims": ("x", "y"), "data": scene_flags[k]} for k in scene_flags}
         old_debug_data = old_debug_data | sf
 
-    debug_output = xr.Dataset.from_dict(old_debug_data)
-    debug_output.to_netcdf("/ships19/hercules/pveglio/debug_data.nc", mode="w")
+    # debug_output = xr.Dataset.from_dict(old_debug_data)
+    # debug_output.to_netcdf("/ships19/hercules/pveglio/debug_data.nc", mode="w")
 
     all_bits = np.zeros((bit["01"].shape[0], bit["01"].shape[1], 22))
     for i in range(len(_bitlist)):
diff --git a/mvcm/write_output.py b/mvcm/write_output.py
index 05dab585a8fb090ec3720bedf17c0e06cddb4b3e..ee933255232aab709f44368f2f617cd789faf176 100644
--- a/mvcm/write_output.py
+++ b/mvcm/write_output.py
@@ -173,7 +173,7 @@ def save_output(data: xr.Dataset, attrs: dict, fname: str, compression: int, deb
     if debug is True:
         debug_data = {
             "bits": {
-                "dims": ("number_of_lines", "number_of_pixels", "bits"),
+                "dims": ("number_of_lines", "number_of_pixels", "bits_dimension"),
                 "data": data.bits.values,
             },
             "bands": {