diff --git a/syconn/extraction/cs_processing_steps.py b/syconn/extraction/cs_processing_steps.py
index 61045d38e6e8c2ab0f9c5b503996505218fbee03..8bf8f51a282f5d547a30571309b30e7660bdb85d 100755
--- a/syconn/extraction/cs_processing_steps.py
+++ b/syconn/extraction/cs_processing_steps.py
@@ -694,7 +694,7 @@ def _combine_and_split_cs_thread(args):
     cs_ssv_id = base_id
 
     attr_dc = AttributeDict(base_dir + "/attr_dict.pkl", read_only=False)
-    mesh_dc = MeshStorage(base_dir + "/mesh.pkl", read_only=False)
+    mesh_dc = MeshStorage(base_dir + "/mesh.pkl", read_only=False, compress=True)
 
     # iterate over cell partners and their contact site IDs (each contact site is between two supervoxels
     # of the partner cells)
@@ -730,8 +730,7 @@ def _combine_and_split_cs_thread(args):
             csssv_attr_dc["mesh_area"] = cs_ssv.mesh_area
             csssv_attr_dc["bounding_box"] = cs_ssv.mesh_bb // scaling
             csssv_attr_dc["rep_coord"] = seghelp.calc_center_of_mass(mesh_cc[1].reshape((-1, 3))) // scaling
-            # TODO: PS remove
-            raise()
+            csssv_attr_dc["cs_ids"] = list(cs_ids)
             # create open3d mesh instance to compute volume
             # # TODO: add this as soon open3d >= 0.11 is supported (glibc error on cluster prevents upgrade)
             # tm = o3d.geometry.TriangleMesh