diff --git a/ess/TraceWin.py b/ess/TraceWin.py
index 8d7bead7cd7c587ef6a0dff2ddec80db83785307..c4c9f0f5fe03723af5cf64d5dfd4cf3d2b3b1bfc 100644
--- a/ess/TraceWin.py
+++ b/ess/TraceWin.py
@@ -740,15 +740,21 @@ class density_file:
             self.Np=self._sum_merge(o,'Np')
             self.Nrun+=o.Nrun
 
-    def savetohdf(self,filename='Density.h5',group='TraceWin'):
+    def savetohdf(self,filename='Density.h5',group='TraceWin', force=False):
         '''
         Saves data to HDF5
         '''
-        import h5py
+        import h5py,sys
 
         fout = h5py.File(filename,'a')
         if group in fout:
-            del fout[group]
+            if force:
+                del fout[group]
+            else:
+                if sys.flags.debug:
+                    print("Group {} already exist in {}".format(group,filename))
+                return
+
         group = fout.create_group(group)
 
         # header attributes..