diff --git a/src/delete_atoms.cpp b/src/delete_atoms.cpp
index a5b15f5f52a848894b35467639f772779c9d831e..4c7829662c6d85bd25df0bfe4bac8f5439f1d058 100644
--- a/src/delete_atoms.cpp
+++ b/src/delete_atoms.cpp
@@ -109,11 +109,14 @@ void DeleteAtoms::command(int narg, char **arg)
   // reset atom tags to be contiguous
   // set all atom IDs to 0, call tag_extend()
 
-  if (atom->molecular == 0 && compress_flag) {
-    tagint *tag = atom->tag;
-    int nlocal = atom->nlocal;
-    for (int i = 0; i < nlocal; i++) tag[i] = 0;
-    atom->tag_extend();
+  if (compress_flag) {
+    if (atom->molecular == 0) {
+      tagint *tag = atom->tag;
+      int nlocal = atom->nlocal;
+      for (int i = 0; i < nlocal; i++) tag[i] = 0;
+      atom->tag_extend();
+    } else if (comm->me == 0)
+        error->warning(FLERR,"Ignoring 'compress yes' for molecular system");
   }
 
   // reset atom->natoms and also topology counts