diff --git a/applications/solvers/lagrangian/icoUncoupledKinematicParcelDyMFoam/Make/files b/applications/solvers/lagrangian/icoUncoupledKinematicParcelDyMFoam/Make/files
new file mode 100644
index 0000000000..2cce02e105
--- /dev/null
+++ b/applications/solvers/lagrangian/icoUncoupledKinematicParcelDyMFoam/Make/files
@@ -0,0 +1,3 @@
+icoUncoupledKinematicParcelDyMFoam.C
+
+EXE = $(FOAM_APPBIN)/icoUncoupledKinematicParcelDyMFoam
diff --git a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelDyMFoam/Make/options b/applications/solvers/lagrangian/icoUncoupledKinematicParcelDyMFoam/Make/options
similarity index 96%
rename from applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelDyMFoam/Make/options
rename to applications/solvers/lagrangian/icoUncoupledKinematicParcelDyMFoam/Make/options
index 34ba95f0dc..7b471b791a 100644
--- a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelDyMFoam/Make/options
+++ b/applications/solvers/lagrangian/icoUncoupledKinematicParcelDyMFoam/Make/options
@@ -1,4 +1,5 @@
 EXE_INC = \
+    -I../icoUncoupledKinematicParcelFoam \
     -I$(LIB_SRC)/lagrangian/basic/lnInclude \
     -I$(LIB_SRC)/lagrangian/intermediate/lnInclude \
     -I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
diff --git a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelDyMFoam/incompressibleUncoupledKinematicParcelDyMFoam.C b/applications/solvers/lagrangian/icoUncoupledKinematicParcelDyMFoam/icoUncoupledKinematicParcelDyMFoam.C
similarity index 96%
rename from applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelDyMFoam/incompressibleUncoupledKinematicParcelDyMFoam.C
rename to applications/solvers/lagrangian/icoUncoupledKinematicParcelDyMFoam/icoUncoupledKinematicParcelDyMFoam.C
index 40b642a284..e6a33f6c8d 100644
--- a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelDyMFoam/incompressibleUncoupledKinematicParcelDyMFoam.C
+++ b/applications/solvers/lagrangian/icoUncoupledKinematicParcelDyMFoam/icoUncoupledKinematicParcelDyMFoam.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2008-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2008-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -36,7 +36,7 @@ Description
 #include "dynamicFvMesh.H"
 #include "singlePhaseTransportModel.H"
 #include "turbulenceModel.H"
-#include "basicKinematicCloud.H"
+#include "basicKinematicCollidingCloud.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
diff --git a/applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/Make/files b/applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/Make/files
new file mode 100644
index 0000000000..95c98aff1d
--- /dev/null
+++ b/applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/Make/files
@@ -0,0 +1,3 @@
+icoUncoupledKinematicParcelFoam.C
+
+EXE = $(FOAM_APPBIN)/icoUncoupledKinematicParcelFoam
diff --git a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/Make/options b/applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/Make/options
similarity index 100%
rename from applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/Make/options
rename to applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/Make/options
diff --git a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/createFields.H b/applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/createFields.H
similarity index 70%
rename from applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/createFields.H
rename to applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/createFields.H
index 31be0cf8ee..9ef9d6a1a0 100644
--- a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/createFields.H
+++ b/applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/createFields.H
@@ -75,7 +75,7 @@
     args.optionReadIfPresent("cloudName", kinematicCloudName);
 
     Info<< "Constructing kinematicCloud " << kinematicCloudName << endl;
-    basicKinematicCloud kinematicCloud
+    basicKinematicCollidingCloud kinematicCloud
     (
         kinematicCloudName,
         rhoInf,
@@ -89,30 +89,17 @@
         "H",
         runTime.timeName(),
         mesh,
-        IOobject::NO_READ
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
     );
 
-    autoPtr<volVectorField> HPtr_;
+    autoPtr<volVectorField> HPtr;
 
     if (Hheader.headerOk())
     {
         Info<< "\nReading field H\n" << endl;
 
-        HPtr_.reset
-        (
-            new volVectorField
-            (
-                IOobject
-                (
-                    "H",
-                    runTime.timeName(),
-                    mesh,
-                    IOobject::MUST_READ,
-                    IOobject::AUTO_WRITE
-                ),
-                mesh
-            )
-        );
+        HPtr.reset(new volVectorField (Hheader, mesh));
     }
 
     IOobject HdotGradHheader
@@ -120,28 +107,17 @@
         "HdotGradH",
         runTime.timeName(),
         mesh,
-        IOobject::NO_READ
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
     );
 
-    autoPtr<volVectorField> HdotGradHPtr_;
+    autoPtr<volVectorField> HdotGradHPtr;
 
     if (HdotGradHheader.headerOk())
     {
         Info<< "Reading field HdotGradH" << endl;
 
-        HdotGradHPtr_.reset
-        (
-            new volVectorField
-            (
-                IOobject
-                (
-                    "HdotGradH",
-                    runTime.timeName(),
-                    mesh,
-                    IOobject::MUST_READ,
-                    IOobject::AUTO_WRITE
-                ),
-                mesh
-            )
-        );
+        HdotGradHPtr.reset(new volVectorField(HdotGradHheader, mesh));
     }
+
+    #include "createNonInertialFrameFields.H"
diff --git a/applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/createNonInertialFrameFields.H b/applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/createNonInertialFrameFields.H
new file mode 100644
index 0000000000..abaede9385
--- /dev/null
+++ b/applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/createNonInertialFrameFields.H
@@ -0,0 +1,88 @@
+    Info<< "Reading non-inertial frame fields" << endl;
+
+    IOobject linearAccelerationHeader
+    (
+        "linearAcceleration",
+        runTime.timeName(),
+        mesh,
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
+    );
+
+    autoPtr<uniformDimensionedVectorField> linearAccelerationPtr;
+
+    if (linearAccelerationHeader.headerOk())
+    {
+        Info<< "    Reading " << linearAccelerationHeader.name() << endl;
+
+        linearAccelerationPtr.reset
+        (
+            new uniformDimensionedVectorField(linearAccelerationHeader)
+        );
+    }
+
+
+    IOobject angularVelocityHeader
+    (
+        "angularVelocity",
+        runTime.timeName(),
+        mesh,
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
+    );
+
+    autoPtr<uniformDimensionedVectorField> angularVelocityPtr;
+
+    if (angularVelocityHeader.headerOk())
+    {
+        Info<< "    Reading " << angularVelocityHeader.name() << endl;
+
+        angularVelocityPtr.reset
+        (
+            new uniformDimensionedVectorField(angularVelocityHeader)
+        );
+    }
+
+
+    IOobject angularAccelerationHeader
+    (
+        "angularAcceleration",
+        runTime.timeName(),
+        mesh,
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
+    );
+
+    autoPtr<uniformDimensionedVectorField> angularAccelerationPtr;
+
+    if (angularAccelerationHeader.headerOk())
+    {
+        Info<< "    Reading " << angularAccelerationHeader.name() << endl;
+
+        angularAccelerationPtr.reset
+        (
+            new uniformDimensionedVectorField(angularAccelerationHeader)
+        );
+    }
+
+
+    IOobject centreOfRotationHeader
+    (
+        "centreOfRotation",
+        runTime.timeName(),
+        mesh,
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
+    );
+
+    autoPtr<uniformDimensionedVectorField> centreOfRotationPtr;
+
+    if (centreOfRotationHeader.headerOk())
+    {
+        Info<< "    Reading " << centreOfRotationHeader.name() << endl;
+
+        centreOfRotationPtr.reset
+        (
+            new uniformDimensionedVectorField(centreOfRotationHeader)
+        );
+    }
diff --git a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/incompressibleUncoupledKinematicParcelFoam.C b/applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/icoUncoupledKinematicParcelFoam.C
similarity index 95%
rename from applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/incompressibleUncoupledKinematicParcelFoam.C
rename to applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/icoUncoupledKinematicParcelFoam.C
index 34d1d3adaa..c04ee90a7b 100644
--- a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/incompressibleUncoupledKinematicParcelFoam.C
+++ b/applications/solvers/lagrangian/icoUncoupledKinematicParcelFoam/icoUncoupledKinematicParcelFoam.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2008-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2008-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -35,7 +35,7 @@ Description
 #include "fvCFD.H"
 #include "singlePhaseTransportModel.H"
 #include "turbulenceModel.H"
-#include "basicKinematicCloud.H"
+#include "basicKinematicCollidingCloud.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
diff --git a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelDyMFoam/Make/files b/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelDyMFoam/Make/files
deleted file mode 100644
index e233a549b8..0000000000
--- a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelDyMFoam/Make/files
+++ /dev/null
@@ -1,3 +0,0 @@
-incompressibleUncoupledKinematicParcelDyMFoam.C
-
-EXE = $(FOAM_APPBIN)/incompressibleUncoupledKinematicParcelDyMFoam
diff --git a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelDyMFoam/createFields.H b/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelDyMFoam/createFields.H
deleted file mode 100644
index 31be0cf8ee..0000000000
--- a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelDyMFoam/createFields.H
+++ /dev/null
@@ -1,147 +0,0 @@
-    Info<< "\nReading transportProperties\n" << endl;
-
-    IOdictionary transportProperties
-    (
-        IOobject
-        (
-            "transportProperties",
-            runTime.constant(),
-            mesh,
-            IOobject::MUST_READ_IF_MODIFIED,
-            IOobject::NO_WRITE
-        )
-    );
-
-    dimensionedScalar rhoInfValue
-    (
-        transportProperties.lookup("rhoInf")
-    );
-
-    volScalarField rhoInf
-    (
-        IOobject
-        (
-            "rho",
-            runTime.timeName(),
-            mesh,
-            IOobject::NO_READ,
-            IOobject::AUTO_WRITE
-        ),
-        mesh,
-        rhoInfValue
-    );
-
-    Info<< "Reading field U\n" << endl;
-    volVectorField U
-    (
-        IOobject
-        (
-            "U",
-            runTime.timeName(),
-            mesh,
-            IOobject::MUST_READ,
-            IOobject::AUTO_WRITE
-        ),
-        mesh
-    );
-
-    #include "createPhi.H"
-
-    Info<< "Creating turbulence model\n" << endl;
-
-    singlePhaseTransportModel laminarTransport(U, phi);
-
-    const volScalarField nu(laminarTransport.nu());
-
-    autoPtr<incompressible::turbulenceModel> turbulence
-    (
-        incompressible::turbulenceModel::New(U, phi, laminarTransport)
-    );
-
-    volScalarField mu
-    (
-        IOobject
-        (
-            "mu",
-            runTime.timeName(),
-            mesh,
-            IOobject::NO_READ,
-            IOobject::AUTO_WRITE
-        ),
-        nu*rhoInfValue
-    );
-
-    word kinematicCloudName("kinematicCloud");
-    args.optionReadIfPresent("cloudName", kinematicCloudName);
-
-    Info<< "Constructing kinematicCloud " << kinematicCloudName << endl;
-    basicKinematicCloud kinematicCloud
-    (
-        kinematicCloudName,
-        rhoInf,
-        U,
-        mu,
-        g
-    );
-
-    IOobject Hheader
-    (
-        "H",
-        runTime.timeName(),
-        mesh,
-        IOobject::NO_READ
-    );
-
-    autoPtr<volVectorField> HPtr_;
-
-    if (Hheader.headerOk())
-    {
-        Info<< "\nReading field H\n" << endl;
-
-        HPtr_.reset
-        (
-            new volVectorField
-            (
-                IOobject
-                (
-                    "H",
-                    runTime.timeName(),
-                    mesh,
-                    IOobject::MUST_READ,
-                    IOobject::AUTO_WRITE
-                ),
-                mesh
-            )
-        );
-    }
-
-    IOobject HdotGradHheader
-    (
-        "HdotGradH",
-        runTime.timeName(),
-        mesh,
-        IOobject::NO_READ
-    );
-
-    autoPtr<volVectorField> HdotGradHPtr_;
-
-    if (HdotGradHheader.headerOk())
-    {
-        Info<< "Reading field HdotGradH" << endl;
-
-        HdotGradHPtr_.reset
-        (
-            new volVectorField
-            (
-                IOobject
-                (
-                    "HdotGradH",
-                    runTime.timeName(),
-                    mesh,
-                    IOobject::MUST_READ,
-                    IOobject::AUTO_WRITE
-                ),
-                mesh
-            )
-        );
-    }
diff --git a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/Make/files b/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/Make/files
deleted file mode 100644
index 7a27f008c8..0000000000
--- a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/Make/files
+++ /dev/null
@@ -1,3 +0,0 @@
-incompressibleUncoupledKinematicParcelFoam.C
-
-EXE = $(FOAM_APPBIN)/incompressibleUncoupledKinematicParcelFoam
diff --git a/src/OpenFOAM/fields/UniformDimensionedFields/UniformDimensionedField.C b/src/OpenFOAM/fields/UniformDimensionedFields/UniformDimensionedField.C
index 14ecb4bc3f..4f49599796 100644
--- a/src/OpenFOAM/fields/UniformDimensionedFields/UniformDimensionedField.C
+++ b/src/OpenFOAM/fields/UniformDimensionedFields/UniformDimensionedField.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2009-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2009-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -61,7 +61,7 @@ Foam::UniformDimensionedField<Type>::UniformDimensionedField
 {
     dictionary dict(readStream(typeName));
     this->dimensions().reset(dict.lookup("dimensions"));
-    this->value() = dict.lookup("value");
+    dict.lookup("value") >> this->value();
 }
 
 
diff --git a/src/lagrangian/basic/Cloud/Cloud.C b/src/lagrangian/basic/Cloud/Cloud.C
index 425c8221cd..6706a47f0e 100644
--- a/src/lagrangian/basic/Cloud/Cloud.C
+++ b/src/lagrangian/basic/Cloud/Cloud.C
@@ -334,11 +334,14 @@ void Foam::Cloud<ParticleType>::move(TrackData& td, const scalar trackTime)
         }
     }
 
-    reduce(nTrackingRescues_, sumOp<label>());
-
-    if (nTrackingRescues_ > 0)
+    if (cloud::debug)
     {
-        Info<< nTrackingRescues_ << " tracking rescue corrections" << endl;
+        reduce(nTrackingRescues_, sumOp<label>());
+
+        if (nTrackingRescues_ > 0)
+        {
+            Info<< nTrackingRescues_ << " tracking rescue corrections" << endl;
+        }
     }
 }
 
diff --git a/src/lagrangian/basic/particle/particleI.H b/src/lagrangian/basic/particle/particleI.H
index ebb72c7cb5..a2d974f52f 100644
--- a/src/lagrangian/basic/particle/particleI.H
+++ b/src/lagrangian/basic/particle/particleI.H
@@ -784,23 +784,26 @@ inline void Foam::particle::initCellFacePt()
                         << position_ << abort(FatalError);
                 }
 
-                WarningIn("void Foam::particle::initCellFacePt()")
-                    << "Particle moved from " << position_
-                    << " to " << newPosition
-                    << " in cell " << cellI_
-                    << " tetFace " << tetFaceI_
-                    << " tetPt " << tetPtI_ << nl
-                    << "    (A fraction of "
-                    << 1.0 - mag(cC - newPosition)/mag(cC - position_)
-                    << " of the distance to the cell centre)"
-                    << " because a decomposition tetFace and tetPt "
-                    << "could not be found."
-                    << endl;
+                if (debug)
+                {
+                    WarningIn("void Foam::particle::initCellFacePt()")
+                        << "Particle moved from " << position_
+                        << " to " << newPosition
+                        << " in cell " << cellI_
+                        << " tetFace " << tetFaceI_
+                        << " tetPt " << tetPtI_ << nl
+                        << "    (A fraction of "
+                        << 1.0 - mag(cC - newPosition)/mag(cC - position_)
+                        << " of the distance to the cell centre)"
+                        << " because a decomposition tetFace and tetPt "
+                        << "could not be found."
+                        << endl;
+                }
 
                 position_ = newPosition;
             }
 
-            if (cellI_ != oldCellI)
+            if (debug && cellI_ != oldCellI)
             {
                 WarningIn("void Foam::particle::initCellFacePt()")
                     << "Particle at position " << position_
diff --git a/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloud.H b/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloud.H
index 5a124609f6..1ab7054e77 100644
--- a/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloud.H
+++ b/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloud.H
@@ -304,6 +304,8 @@ public:
             //- Return a reference to the cloud copy
             inline const KinematicCloud& cloudCopy() const;
 
+            inline bool hasWallImpactDistance() const;
+
 
             // References to the mesh and databases
 
diff --git a/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloudI.H b/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloudI.H
index 5987192bfd..05bd862d8a 100644
--- a/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloudI.H
+++ b/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloudI.H
@@ -35,6 +35,13 @@ Foam::KinematicCloud<CloudType>::cloudCopy() const
 }
 
 
+template<class CloudType>
+inline bool Foam::KinematicCloud<CloudType>::hasWallImpactDistance() const
+{
+    return true;
+}
+
+
 template<class CloudType>
 inline const Foam::fvMesh& Foam::KinematicCloud<CloudType>::mesh() const
 {
diff --git a/src/lagrangian/intermediate/parcels/include/makeParcelForces.H b/src/lagrangian/intermediate/parcels/include/makeParcelForces.H
index 9ae02038aa..cc824fd038 100644
--- a/src/lagrangian/intermediate/parcels/include/makeParcelForces.H
+++ b/src/lagrangian/intermediate/parcels/include/makeParcelForces.H
@@ -32,6 +32,7 @@ License
 #include "NonSphereDragForce.H"
 
 #include "GravityForce.H"
+#include "NonInertialFrameForce.H"
 #include "ParamagneticForce.H"
 #include "PressureGradientForce.H"
 #include "SRFForce.H"
@@ -44,6 +45,7 @@ License
     makeParticleForceModelType(SphereDragForce, CloudType);                   \
     makeParticleForceModelType(NonSphereDragForce, CloudType);                \
     makeParticleForceModelType(GravityForce, CloudType);                      \
+    makeParticleForceModelType(NonInertialFrameForce, CloudType);             \
     makeParticleForceModelType(ParamagneticForce, CloudType);                 \
     makeParticleForceModelType(PressureGradientForce, CloudType);             \
     makeParticleForceModelType(SRFForce, CloudType);
diff --git a/src/lagrangian/intermediate/parcels/include/makeThermoParcelForces.H b/src/lagrangian/intermediate/parcels/include/makeThermoParcelForces.H
index 45b4a86215..fd20d0191d 100644
--- a/src/lagrangian/intermediate/parcels/include/makeThermoParcelForces.H
+++ b/src/lagrangian/intermediate/parcels/include/makeThermoParcelForces.H
@@ -33,6 +33,7 @@ License
 
 #include "BrownianMotionForce.H"
 #include "GravityForce.H"
+#include "NonInertialFrameForce.H"
 #include "ParamagneticForce.H"
 #include "PressureGradientForce.H"
 #include "SRFForce.H"
@@ -46,6 +47,7 @@ License
     makeParticleForceModelType(NonSphereDragForce, CloudType);                \
     makeParticleForceModelType(BrownianMotionForce, CloudType);               \
     makeParticleForceModelType(GravityForce, CloudType);                      \
+    makeParticleForceModelType(NonInertialFrameForce, CloudType);             \
     makeParticleForceModelType(ParamagneticForce, CloudType);                 \
     makeParticleForceModelType(PressureGradientForce, CloudType);             \
     makeParticleForceModelType(SRFForce, CloudType);
diff --git a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/Gravity/GravityForce.H b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/Gravity/GravityForce.H
index a215d0ef66..1bc573df04 100644
--- a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/Gravity/GravityForce.H
+++ b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/Gravity/GravityForce.H
@@ -98,7 +98,7 @@ public:
 
         // Access
 
-            //- Return the the acceleration due to gravity
+            //- Return the acceleration due to gravity
             inline const vector& g() const;
 
 
diff --git a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/NonInertialFrame/NonInertialFrameForce.C b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/NonInertialFrame/NonInertialFrameForce.C
new file mode 100644
index 0000000000..77d8c3e037
--- /dev/null
+++ b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/NonInertialFrame/NonInertialFrameForce.C
@@ -0,0 +1,206 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2011-2011 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "NonInertialFrameForce.H"
+#include "uniformDimensionedFields.H"
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+template<class CloudType>
+Foam::NonInertialFrameForce<CloudType>::NonInertialFrameForce
+(
+    CloudType& owner,
+    const fvMesh& mesh,
+    const dictionary& dict,
+    const word& forceType
+)
+:
+    ParticleForce<CloudType>(owner, mesh, dict),
+    WName_
+    (
+        this->coeffs().template lookupOrDefault<word>
+        (
+            "linearAccelerationName",
+            "linearAcceleration"
+        )
+    ),
+    W_(vector::zero),
+    omegaName_
+    (
+        this->coeffs().template lookupOrDefault<word>
+        (
+            "angularVelocityName",
+            "angularVelocity"
+        )
+    ),
+    omega_(vector::zero),
+    omegaDotName_
+    (
+        this->coeffs().template lookupOrDefault<word>
+        (
+            "angularAccelerationName",
+            "angularAcceleration"
+        )
+    ),
+    omegaDot_(vector::zero),
+    centreOfRotationName_
+    (
+        this->coeffs().template lookupOrDefault<word>
+        (
+            "centreOfRotationName",
+            "centreOfRotation"
+        )
+    ),
+    centreOfRotation_(vector::zero)
+{}
+
+
+template<class CloudType>
+Foam::NonInertialFrameForce<CloudType>::NonInertialFrameForce
+(
+    const NonInertialFrameForce& niff
+)
+:
+    ParticleForce<CloudType>(niff),
+    WName_(niff.WName_),
+    W_(niff.W_),
+    omegaName_(niff.omegaName_),
+    omega_(niff.omega_),
+    omegaDotName_(niff.omegaDotName_),
+    omegaDot_(niff.omegaDot_),
+    centreOfRotationName_(niff.centreOfRotationName_),
+    centreOfRotation_(niff.centreOfRotation_)
+{}
+
+
+// * * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * //
+
+template<class CloudType>
+Foam::NonInertialFrameForce<CloudType>::~NonInertialFrameForce()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+template<class CloudType>
+void Foam::NonInertialFrameForce<CloudType>::cacheFields(const bool store)
+{
+    W_ = vector::zero;
+    omega_ = vector::zero;
+    omegaDot_ = vector::zero;
+    centreOfRotation_ = vector::zero;
+
+    if (store)
+    {
+        if
+        (
+            this->mesh().template foundObject<uniformDimensionedVectorField>
+            (
+                WName_
+            )
+        )
+        {
+            uniformDimensionedVectorField W = this->mesh().template
+                lookupObject<uniformDimensionedVectorField>(WName_);
+
+            W_ = W.value();
+        }
+
+        if
+        (
+            this->mesh().template foundObject<uniformDimensionedVectorField>
+            (
+                omegaName_
+            )
+        )
+        {
+            uniformDimensionedVectorField omega = this->mesh().template
+                lookupObject<uniformDimensionedVectorField>(omegaName_);
+
+            omega_ = omega.value();
+        }
+
+        if
+        (
+            this->mesh().template foundObject<uniformDimensionedVectorField>
+            (
+                omegaDotName_
+            )
+        )
+        {
+            uniformDimensionedVectorField omegaDot = this->mesh().template
+                lookupObject<uniformDimensionedVectorField>(omegaDotName_);
+
+            omegaDot_ = omegaDot.value();
+        }
+
+        if
+        (
+            this->mesh().template foundObject<uniformDimensionedVectorField>
+            (
+                centreOfRotationName_
+            )
+        )
+        {
+            uniformDimensionedVectorField omegaDot = this->mesh().template
+                lookupObject<uniformDimensionedVectorField>
+                (
+                    centreOfRotationName_
+                );
+
+            centreOfRotation_ = omegaDot.value();
+        }
+    }
+}
+
+
+template<class CloudType>
+Foam::forceSuSp Foam::NonInertialFrameForce<CloudType>::calcNonCoupled
+(
+    const typename CloudType::parcelType& p,
+    const scalar dt,
+    const scalar mass,
+    const scalar Re,
+    const scalar muc
+) const
+{
+    forceSuSp value(vector::zero, 0.0);
+
+    const vector& r = p.position() - centreOfRotation_;
+
+    value.Su() =
+        mass
+       *(
+           -W_
+          + (r ^ omegaDot_)
+          + 2.0*(p.U() ^ omega_)
+          + (omega_ ^ (r ^ omega_))
+        );
+
+    return value;
+}
+
+
+// ************************************************************************* //
diff --git a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/NonInertialFrame/NonInertialFrameForce.H b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/NonInertialFrame/NonInertialFrameForce.H
new file mode 100644
index 0000000000..aaf8f1065a
--- /dev/null
+++ b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/NonInertialFrame/NonInertialFrameForce.H
@@ -0,0 +1,168 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2011-2011 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::NonInertialFrameForce
+
+Description
+    Calculates particle non-inertial reference frame force.  Variable names as
+    from Landau and Lifshitz, Mechanics, 3rd Ed, p126-129.
+
+SourceFiles
+    NonInertialFrameForce.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef NonInertialFrameForce_H
+#define NonInertialFrameForce_H
+
+#include "ParticleForce.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+class fvMesh;
+
+/*---------------------------------------------------------------------------*\
+                  Class NonInertialFrameForce Declaration
+\*---------------------------------------------------------------------------*/
+
+template<class CloudType>
+class NonInertialFrameForce
+:
+    public ParticleForce<CloudType>
+{
+    // Private data
+
+        //- Name of the linear acceleration field
+        word WName_;
+
+        //- The linear acceleration of the reference frame
+        vector W_;
+
+        //- Name of the angular velocity field
+        word omegaName_;
+
+        //- The angular velocity of the reference frame
+        vector omega_;
+
+        //- Name of the angular acceleration field
+        word omegaDotName_;
+
+        //- The angular acceleration of the reference frame
+        vector omegaDot_;
+
+        //- Name of the centre of rotation field
+        word centreOfRotationName_;
+
+        //- The centre of rotation of the reference frame
+        vector centreOfRotation_;
+
+
+public:
+
+    //- Runtime type information
+    TypeName("nonInertialFrame");
+
+
+    // Constructors
+
+        //- Construct from mesh
+        NonInertialFrameForce
+        (
+            CloudType& owner,
+            const fvMesh& mesh,
+            const dictionary& dict,
+            const word& forceType
+        );
+
+        //- Construct copy
+        NonInertialFrameForce(const NonInertialFrameForce& niff);
+
+        //- Construct and return a clone
+        virtual autoPtr<ParticleForce<CloudType> > clone() const
+        {
+            return autoPtr<ParticleForce<CloudType> >
+            (
+                new ParticleForce<CloudType>(*this)
+            );
+        }
+
+
+    //- Destructor
+    virtual ~NonInertialFrameForce();
+
+
+    // Member Functions
+
+        // Access
+
+            //- Return the linear acceleration of the reference frame
+            inline const vector& W() const;
+
+            //- Return the angular velocity of the reference frame
+            inline const vector& omega() const;
+
+            //- Return the angular acceleration of the reference frame
+            inline const vector& omegaDot() const;
+
+            //- Return the centre of rotation of the reference frame
+            inline const vector& centreOfRotation() const;
+
+
+        // Evaluation
+
+            //- Cache fields
+            virtual void cacheFields(const bool store);
+
+            //- Calculate the non-coupled force
+            virtual forceSuSp calcNonCoupled
+            (
+                const typename CloudType::parcelType& p,
+                const scalar dt,
+                const scalar mass,
+                const scalar Re,
+                const scalar muc
+            ) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#include "NonInertialFrameForceI.H"
+
+#ifdef NoRepository
+    #include "NonInertialFrameForce.C"
+#endif
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/NonInertialFrame/NonInertialFrameForceI.H b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/NonInertialFrame/NonInertialFrameForceI.H
new file mode 100644
index 0000000000..fb59297222
--- /dev/null
+++ b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/NonInertialFrame/NonInertialFrameForceI.H
@@ -0,0 +1,58 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2011-2011 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+template<class CloudType>
+inline const Foam::vector& Foam::NonInertialFrameForce<CloudType>::W() const
+{
+    return W_;
+}
+
+
+template<class CloudType>
+inline const Foam::vector& Foam::NonInertialFrameForce<CloudType>::omega() const
+{
+    return omega_;
+}
+
+
+template<class CloudType>
+inline const Foam::vector&
+Foam::NonInertialFrameForce<CloudType>::omegaDot() const
+{
+    return omegaDot_;
+}
+
+
+template<class CloudType>
+inline const Foam::vector&
+Foam::NonInertialFrameForce<CloudType>::centreOfRotation() const
+{
+    return centreOfRotation_;
+}
+
+
+// ************************************************************************* //
diff --git a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.C b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.C
index 2f9ca986b4..e001d29f19 100644
--- a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.C
+++ b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.C
@@ -92,12 +92,11 @@ Foam::forceSuSp Foam::SRFForce<CloudType>::calcNonCoupled
     const typename SRF::SRFModel& srf = *srfPtr_;
 
     const vector& omega = srf.omega().value();
-    const vector& axis = srf.axis();
 
-    const vector r = p.position() - axis*(axis & p.position());
+    const vector& r = p.position();
 
     // Coriolis and centrifugal acceleration terms
-    value.Su() = mass*2.0*(p.U() ^ omega) + (omega ^ (r ^ omega));
+    value.Su() = mass*(2.0*(p.U() ^ omega) + (omega ^ (r ^ omega)));
 
     return value;
 }
diff --git a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.H b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.H
index 2e646a37d7..b8f10d0ea9 100644
--- a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.H
+++ b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.H
@@ -28,7 +28,6 @@ Description
     Calculates particle SRF reference frame force
 
 SourceFiles
-    SRFForceI.H
     SRFForce.C
 
 \*---------------------------------------------------------------------------*/