diff --git a/src/lagrangian/intermediate/parcels/Templates/CollidingParcel/CollidingParcel.C b/src/lagrangian/intermediate/parcels/Templates/CollidingParcel/CollidingParcel.C
index f422f7ea372fc1c9904125e6b01e1fe2de8362e5..36f146d77af4e7892309ec000bc70512c76e925a 100644
--- a/src/lagrangian/intermediate/parcels/Templates/CollidingParcel/CollidingParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/CollidingParcel/CollidingParcel.C
@@ -80,6 +80,8 @@ bool Foam::CollidingParcel<ParcelType>::move
 
             p.angularMomentum() += 0.5*trackTime*p.torque();
 
+            td.keepParticle = true;
+
             break;
         }
 
diff --git a/src/lagrangian/intermediate/parcels/Templates/MPPICParcel/MPPICParcel.C b/src/lagrangian/intermediate/parcels/Templates/MPPICParcel/MPPICParcel.C
index 3844d22f673ebc0185168e10bf873164df188cf5..0bac21bf8f3d6d5120217f93ff27ac229f6df848 100644
--- a/src/lagrangian/intermediate/parcels/Templates/MPPICParcel/MPPICParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/MPPICParcel/MPPICParcel.C
@@ -76,6 +76,8 @@ bool Foam::MPPICParcel<ParcelType>::move
             p.UCorrect() =
                 td.cloud().dampingModel().velocityCorrection(p, trackTime);
 
+            td.keepParticle = true;
+
             break;
         }
         case TrackData::tpPackingNoTrack:
@@ -83,6 +85,8 @@ bool Foam::MPPICParcel<ParcelType>::move
             p.UCorrect() =
                 td.cloud().packingModel().velocityCorrection(p, trackTime);
 
+            td.keepParticle = true;
+
             break;
         }
         case TrackData::tpCorrectTrack: