From b0276efa48a57642b31cdcbbeedb839176cdd042 Mon Sep 17 00:00:00 2001
From: andy <a.heather@opencfd.co.uk>
Date: Fri, 14 Jan 2011 15:43:36 +0000
Subject: [PATCH] ENH: Update to parcel move()
---
.../parcels/Templates/KinematicParcel/KinematicParcel.C | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C
index 81ba4c10231..d21aad68b11 100644
--- a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C
@@ -338,12 +338,12 @@ bool Foam::KinematicParcel<ParcelType>::move
// Remember which cell the Parcel is in since this
// will change if a face is hit
- label cellI = p.cell();
+ const label cellI = p.cell();
const scalar magU = mag(U_);
if (p.active() && magU > ROOTVSMALL)
{
- const scalar d = mag(dt*U_) + ROOTVSMALL;
+ const scalar d = dt*magU;
const vector n = U_/magU;
const scalar dCorr = min(d, mag(n & mesh.bounds().span()));
dt *= dCorr/d*p.trackToFace(p.position() + dCorr*n, td);
--
GitLab