diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/solid/solveSolid.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/solid/solveSolid.H
index ca29fcd0d9..c480614510 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/solid/solveSolid.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/solid/solveSolid.H
@@ -3,9 +3,11 @@
     {
         fvScalarMatrix hEqn
         (
-            thermo.isotropic()
-            ? -fvm::laplacian(betav*thermo.alpha(), h, "laplacian(alpha,h)")
-            : -fvm::laplacian(betav*tAnialpha(), h, "laplacian(alpha,h)")
+            (
+                thermo.isotropic()
+                ? -fvm::laplacian(betav*thermo.alpha(), h, "laplacian(alpha,h)")
+                : -fvm::laplacian(betav*tAnialpha(), h, "laplacian(alpha,h)")
+            )
           ==
             fvOptions(rho, h)
         );
diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/createSolidFields.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/createSolidFields.H
index 703ed6e4c2..e444914c02 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/createSolidFields.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/createSolidFields.H
@@ -11,10 +11,10 @@
         Info<< "*** Reading solid mesh thermophysical properties for region "
             << solidRegions[i].name() << nl << endl;
 
-        Info<< "    Adding to thermos\n" << endl;
+        Info<< "    Adding to thermo models\n" << endl;
         thermos.set(i, solidThermo::New(solidRegions[i]));
 
-        Info<< "    Adding to radiations\n" << endl;
+        Info<< "    Adding to radiation models\n" << endl;
         radiations.set(i, radiation::radiationModel::New(thermos[i].T()));
 
         Info<< "    Adding fvOptions\n" << endl;
diff --git a/applications/utilities/postProcessing/miscellaneous/execFlowFunctionObjects/Make/options b/applications/utilities/postProcessing/miscellaneous/execFlowFunctionObjects/Make/options
index 2d8a878038..012600720f 100644
--- a/applications/utilities/postProcessing/miscellaneous/execFlowFunctionObjects/Make/options
+++ b/applications/utilities/postProcessing/miscellaneous/execFlowFunctionObjects/Make/options
@@ -3,7 +3,10 @@ EXE_INC = \
     -I$(LIB_SRC)/turbulenceModels \
     -I$(LIB_SRC)/turbulenceModels/LES/LESdeltas/lnInclude \
     -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
-    -I$(LIB_SRC)/finiteVolume/lnInclude
+    -I$(LIB_SRC)/finiteVolume/lnInclude \
+    -I$(LIB_SRC)/fvOptions/lnInclude \
+    -I$(LIB_SRC)/meshTools/lnInclude \
+    -I$(LIB_SRC)/sampling/lnInclude
 
 EXE_LIBS = \
     -lincompressibleTransportModels \
@@ -14,4 +17,7 @@ EXE_LIBS = \
     -lcompressibleRASModels \
     -lcompressibleLESModels \
     -lfiniteVolume \
+    -lfvOptions \
+    -lmeshTools \
+    -lsampling \
     -lgenericPatchFields
diff --git a/applications/utilities/postProcessing/miscellaneous/execFlowFunctionObjects/execFlowFunctionObjects.C b/applications/utilities/postProcessing/miscellaneous/execFlowFunctionObjects/execFlowFunctionObjects.C
index 3c13617a78..4849232d69 100644
--- a/applications/utilities/postProcessing/miscellaneous/execFlowFunctionObjects/execFlowFunctionObjects.C
+++ b/applications/utilities/postProcessing/miscellaneous/execFlowFunctionObjects/execFlowFunctionObjects.C
@@ -42,6 +42,7 @@ Description
 #include "surfaceFields.H"
 #include "pointFields.H"
 #include "ReadFields.H"
+#include "fvIOoptionList.H"
 
 #include "incompressible/singlePhaseTransportModel/singlePhaseTransportModel.H"
 
@@ -167,6 +168,8 @@ void calc
             mesh
         );
 
+        #include "createFvOptions.H"
+
         if (phi.dimensions() == dimVolume/dimTime)
         {
             IOobject RASPropertiesHeader
diff --git a/src/OpenFOAM/primitives/functions/DataEntry/Table/Table.C b/src/OpenFOAM/primitives/functions/DataEntry/Table/Table.C
index 56c7a2b669..0ee107e079 100644
--- a/src/OpenFOAM/primitives/functions/DataEntry/Table/Table.C
+++ b/src/OpenFOAM/primitives/functions/DataEntry/Table/Table.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -31,7 +31,7 @@ template<class Type>
 Foam::Table<Type>::Table(const word& entryName, const dictionary& dict)
 :
     DataEntry<Type>(entryName),
-    TableBase<Type>(entryName, dictionary::null)
+    TableBase<Type>(entryName, dict)
 {
     Istream& is(dict.lookup(entryName));
     word entryType(is);
diff --git a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObject/CloudFunctionObject.C b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObject/CloudFunctionObject.C
index e67b0b938c..b0399c5acf 100644
--- a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObject/CloudFunctionObject.C
+++ b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObject/CloudFunctionObject.C
@@ -113,7 +113,7 @@ void Foam::CloudFunctionObject<CloudType>::postEvolve()
 template<class CloudType>
 void Foam::CloudFunctionObject<CloudType>::postMove
 (
-    const typename CloudType::parcelType&,
+    typename CloudType::parcelType&,
     const label,
     const scalar,
     const point&,
diff --git a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObject/CloudFunctionObject.H b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObject/CloudFunctionObject.H
index 4fe6798676..d5b9369df7 100644
--- a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObject/CloudFunctionObject.H
+++ b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObject/CloudFunctionObject.H
@@ -144,7 +144,7 @@ public:
             //- Post-move hook
             virtual void postMove
             (
-                const typename CloudType::parcelType& p,
+                typename CloudType::parcelType& p,
                 const label cellI,
                 const scalar dt,
                 const point& position0,
diff --git a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObjectList/CloudFunctionObjectList.C b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObjectList/CloudFunctionObjectList.C
index ed3d5534d0..84d2c7d41e 100644
--- a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObjectList/CloudFunctionObjectList.C
+++ b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObjectList/CloudFunctionObjectList.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -136,7 +136,7 @@ void Foam::CloudFunctionObjectList<CloudType>::postEvolve()
 template<class CloudType>
 void Foam::CloudFunctionObjectList<CloudType>::postMove
 (
-    const typename CloudType::parcelType& p,
+    typename CloudType::parcelType& p,
     const label cellI,
     const scalar dt,
     const point& position0,
diff --git a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObjectList/CloudFunctionObjectList.H b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObjectList/CloudFunctionObjectList.H
index 477afe37a5..9448d2559c 100644
--- a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObjectList/CloudFunctionObjectList.H
+++ b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/CloudFunctionObjectList/CloudFunctionObjectList.H
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -112,7 +112,7 @@ public:
             //- Post-move hook
             virtual void postMove
             (
-                const typename CloudType::parcelType& p,
+                typename CloudType::parcelType& p,
                 const label cellI,
                 const scalar dt,
                 const point& position0,
diff --git a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleCollector/ParticleCollector.C b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleCollector/ParticleCollector.C
index bc2c559aa7..bf347c9f30 100644
--- a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleCollector/ParticleCollector.C
+++ b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleCollector/ParticleCollector.C
@@ -106,11 +106,13 @@ void Foam::ParticleCollector<CloudType>::makeLogFile
 
 
 template<class CloudType>
-void Foam::ParticleCollector<CloudType>::initPolygons()
+void Foam::ParticleCollector<CloudType>::initPolygons
+(
+    const List<Field<point> >& polygons
+)
 {
     mode_ = mtPolygon;
 
-    List<Field<point> > polygons(this->coeffDict().lookup("polygons"));
     label nPoints = 0;
     forAll(polygons, polyI)
     {
@@ -168,7 +170,7 @@ void Foam::ParticleCollector<CloudType>::initConcentricCircles()
     if (nSector_ > 1)
     {
         refDir = this->coeffDict().lookup("refDir");
-        refDir -= normal_*(normal_ & refDir);
+        refDir -= normal_[0]*(normal_[0] & refDir);
         refDir /= mag(refDir);
     }
     else
@@ -184,7 +186,7 @@ void Foam::ParticleCollector<CloudType>::initConcentricCircles()
         {
             vector v = rnd.vector01();
 
-            tangent = v - (v & normal_)*normal_;
+            tangent = v - (v & normal_[0])*normal_[0];
             magTangent = mag(tangent);
         }
 
@@ -209,7 +211,7 @@ void Foam::ParticleCollector<CloudType>::initConcentricCircles()
     faces_.setSize(nFace);
     area_.setSize(nFace);
 
-    coordSys_ = cylindricalCS("coordSys", origin, normal_, refDir, false);
+    coordSys_ = cylindricalCS("coordSys", origin, normal_[0], refDir, false);
 
     List<label> ptIDs(identity(nPointPerRadius));
 
@@ -301,8 +303,8 @@ void Foam::ParticleCollector<CloudType>::collectParcelPolygon
 
         const point& pf = points_[facePoint0];
 
-        const scalar d1 = normal_ & (p1 - pf);
-        const scalar d2 = normal_ & (p2 - pf);
+        const scalar d1 = normal_[faceI] & (p1 - pf);
+        const scalar d2 = normal_[faceI] & (p2 - pf);
 
         if (sign(d1) == sign(d2))
         {
@@ -344,8 +346,8 @@ void Foam::ParticleCollector<CloudType>::collectParcelConcentricCircles
 {
     label secI = -1;
 
-    const scalar d1 = normal_ & (p1 - coordSys_.origin());
-    const scalar d2 = normal_ & (p2 - coordSys_.origin());
+    const scalar d1 = normal_[0] & (p1 - coordSys_.origin());
+    const scalar d2 = normal_[0] & (p2 - coordSys_.origin());
 
     if (sign(d1) == sign(d2))
     {
@@ -526,7 +528,7 @@ Foam::ParticleCollector<CloudType>::ParticleCollector
     nSector_(0),
     radius_(),
     coordSys_(false),
-    normal_(this->coeffDict().lookup("normal")),
+    normal_(),
     negateParcelsOppositeNormal_
     (
         readBool(this->coeffDict().lookup("negateParcelsOppositeNormal"))
@@ -547,10 +549,37 @@ Foam::ParticleCollector<CloudType>::ParticleCollector
     word mode(this->coeffDict().lookup("mode"));
     if (mode == "polygon")
     {
-        initPolygons();
+        List<Field<point> > polygons(this->coeffDict().lookup("polygons"));
+
+        initPolygons(polygons);
+
+        vector n0(this->coeffDict().lookup("normal"));
+        normal_ = vectorField(faces_.size(), n0);
+    }
+    else if (mode == "polygonWithNormal")
+    {
+        List<Tuple2<Field<point>, vector> > polygonAndNormal
+        (
+            this->coeffDict().lookup("polygons")
+        );
+
+        List<Field<point> > polygons(polygonAndNormal.size());
+        normal_.setSize(polygonAndNormal.size());
+
+        forAll(polygons, polyI)
+        {
+            polygons[polyI] = polygonAndNormal[polyI].first();
+            normal_[polyI] = polygonAndNormal[polyI].second();
+            normal_[polyI] /= mag(normal_[polyI]) + ROOTVSMALL;
+        }
+
+        initPolygons(polygons);
     }
     else if (mode == "concentricCircle")
     {
+        vector n0(this->coeffDict().lookup("normal"));
+        normal_ = vectorField(1, n0);
+
         initConcentricCircles();
     }
     else
@@ -560,12 +589,14 @@ Foam::ParticleCollector<CloudType>::ParticleCollector
             "Foam::ParticleCollector<CloudType>::ParticleCollector"
             "("
                 "const dictionary&,"
-                "CloudType&"
+                "CloudType&, "
+                "const word&"
             ")",
             this->coeffDict()
         )
             << "Unknown mode " << mode << ".  Available options are "
-            << "polygon and concentricCircle" << exit(FatalIOError);
+            << "polygon, polygonWithNormal and concentricCircle"
+            << exit(FatalIOError);
     }
 
     mass_.setSize(faces_.size(), 0.0);
@@ -619,7 +650,7 @@ Foam::ParticleCollector<CloudType>::~ParticleCollector()
 template<class CloudType>
 void Foam::ParticleCollector<CloudType>::postMove
 (
-    const parcelType& p,
+    parcelType& p,
     const label cellI,
     const scalar dt,
     const point& position0,
@@ -663,7 +694,7 @@ void Foam::ParticleCollector<CloudType>::postMove
         {
             vector Uhat = p.U();
             Uhat /= mag(Uhat) + ROOTVSMALL;
-            if ((Uhat & normal_) < 0)
+            if ((Uhat & normal_[faceI]) < 0)
             {
                 m *= -1.0;
             }
diff --git a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleCollector/ParticleCollector.H b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleCollector/ParticleCollector.H
index 2f6faea6b5..e7ed170e97 100644
--- a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleCollector/ParticleCollector.H
+++ b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleCollector/ParticleCollector.H
@@ -164,8 +164,8 @@ private:
         //- Face areas
         Field<scalar> area_;
 
-        //- Polygon normal vector
-        vector normal_;
+        //- Polygon normal vector per face
+        Field<vector> normal_;
 
         //- Remove mass of parcel travelling in opposite direction to normal_
         bool negateParcelsOppositeNormal_;
@@ -212,7 +212,7 @@ private:
         );
 
         //- Initialise polygon collectors
-        void initPolygons();
+        void initPolygons(const List<Field<point> >& polygons);
 
         //- Initialise concentric circle collectors
         void initConcentricCircles();
@@ -286,7 +286,7 @@ public:
             //- Post-move hook
             virtual void postMove
             (
-                const parcelType& p,
+                parcelType& p,
                 const label cellI,
                 const scalar dt,
                 const point& position0,
diff --git a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleTrap/ParticleTrap.C b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleTrap/ParticleTrap.C
index 317a033844..b26d7766e5 100644
--- a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleTrap/ParticleTrap.C
+++ b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleTrap/ParticleTrap.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2012-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2012-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -106,6 +106,7 @@ void Foam::ParticleTrap<CloudType>::postMove
     parcelType& p,
     const label cellI,
     const scalar,
+    const point&,
     bool&
 )
 {
diff --git a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleTrap/ParticleTrap.H b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleTrap/ParticleTrap.H
index 81ccb68e3a..fc03b0c91a 100644
--- a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleTrap/ParticleTrap.H
+++ b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/ParticleTrap/ParticleTrap.H
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2012-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2012-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -132,6 +132,7 @@ public:
                 typename CloudType::parcelType& p,
                 const label cellI,
                 const scalar dt,
+                const point& position0,
                 bool& keepParticle
             );
 };
diff --git a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/VoidFraction/VoidFraction.C b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/VoidFraction/VoidFraction.C
index 12188730c3..7314d531ba 100644
--- a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/VoidFraction/VoidFraction.C
+++ b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/VoidFraction/VoidFraction.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -124,9 +124,10 @@ void Foam::VoidFraction<CloudType>::postEvolve()
 template<class CloudType>
 void Foam::VoidFraction<CloudType>::postMove
 (
-    const parcelType& p,
+    parcelType& p,
     const label cellI,
     const scalar dt,
+    const point&,
     bool&
 )
 {
diff --git a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/VoidFraction/VoidFraction.H b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/VoidFraction/VoidFraction.H
index 502e4aa3ec..a5fb87265f 100644
--- a/src/lagrangian/intermediate/submodels/CloudFunctionObjects/VoidFraction/VoidFraction.H
+++ b/src/lagrangian/intermediate/submodels/CloudFunctionObjects/VoidFraction/VoidFraction.H
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -118,9 +118,10 @@ public:
             //- Post-move hook
             virtual void postMove
             (
-                const parcelType& p,
+                parcelType& p,
                 const label cellI,
                 const scalar dt,
+                const point& position0,
                 bool& keepParticle
             );
 };
diff --git a/src/lagrangian/spray/Make/options b/src/lagrangian/spray/Make/options
index d1deb55ae6..ea1cc90a54 100644
--- a/src/lagrangian/spray/Make/options
+++ b/src/lagrangian/spray/Make/options
@@ -4,6 +4,7 @@ EXE_INC = \
     -I$(LIB_SRC)/lagrangian/basic/lnInclude \
     -I$(LIB_SRC)/lagrangian/intermediate/lnInclude \
     -I$(LIB_SRC)/lagrangian/distributionModels/lnInclude \
+    -I$(LIB_SRC)/lagrangian/turbulence/lnInclude \
     -I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
     -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
     -I$(LIB_SRC)/thermophysicalModels/properties/liquidProperties/lnInclude \
diff --git a/src/lagrangian/spray/parcels/derived/basicSprayParcel/makeBasicSprayParcelSubmodels.C b/src/lagrangian/spray/parcels/derived/basicSprayParcel/makeBasicSprayParcelSubmodels.C
index 3af2751d2f..3f92b79a57 100644
--- a/src/lagrangian/spray/parcels/derived/basicSprayParcel/makeBasicSprayParcelSubmodels.C
+++ b/src/lagrangian/spray/parcels/derived/basicSprayParcel/makeBasicSprayParcelSubmodels.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -29,7 +29,9 @@ License
 
 // Kinematic
 #include "makeThermoParcelForces.H" // thermo variant
+#include "makeThermoParcelTurbulenceForces.H" // add turbulence variant
 #include "makeParcelDispersionModels.H"
+#include "makeParcelTurbulenceDispersionModels.H" // add turbulence variant
 #include "makeSprayParcelInjectionModels.H" // Spray variant
 #include "makeParcelPatchInteractionModels.H"
 #include "makeSprayParcelStochasticCollisionModels.H" // Spray variant
@@ -54,7 +56,9 @@ namespace Foam
 
     // Kinematic sub-models
     makeThermoParcelForces(basicSprayCloud);
+    makeThermoParcelTurbulenceForces(basicSprayCloud);
     makeParcelDispersionModels(basicSprayCloud);
+    makeParcelTurbulenceDispersionModels(basicSprayCloud);
     makeSprayParcelInjectionModels(basicSprayCloud);
     makeParcelPatchInteractionModels(basicSprayCloud);
     makeSprayParcelStochasticCollisionModels(basicSprayCloud);
diff --git a/src/lagrangian/turbulence/parcels/derived/basicKinematicCollidingParcel/makeBasicKinematicCollidingParcelSubmodels.C b/src/lagrangian/turbulence/parcels/derived/basicKinematicCollidingParcel/makeBasicKinematicCollidingParcelSubmodels.C
index 584dc18630..c0d2f76c47 100644
--- a/src/lagrangian/turbulence/parcels/derived/basicKinematicCollidingParcel/makeBasicKinematicCollidingParcelSubmodels.C
+++ b/src/lagrangian/turbulence/parcels/derived/basicKinematicCollidingParcel/makeBasicKinematicCollidingParcelSubmodels.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -24,33 +24,13 @@ License
 \*---------------------------------------------------------------------------*/
 
 #include "basicKinematicCollidingCloud.H"
-
-#include "makeParcelCloudFunctionObjects.H"
-
-#include "GradientDispersionRAS.H"
-#include "StochasticDispersionRAS.H"
+#include "makeParcelTurbulenceDispersionModels.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 namespace Foam
 {
-    typedef basicKinematicCollidingCloud::kinematicCloudType
-        kinematicCloudType_K;
-    defineNamedTemplateTypeNameAndDebug
-    (
-        DispersionRASModel<kinematicCloudType_K>,
-        0
-    );
-    makeDispersionModelType
-    (
-        GradientDispersionRAS,
-        basicKinematicCollidingCloud
-    );
-    makeDispersionModelType
-    (
-        StochasticDispersionRAS,
-        basicKinematicCollidingCloud
-    );
+    makeParcelTurbulenceDispersionModels(basicKinematicCollidingCloud);
 }
 
 
diff --git a/src/lagrangian/turbulence/parcels/derived/basicKinematicParcel/makeBasicKinematicParcelSubmodels.C b/src/lagrangian/turbulence/parcels/derived/basicKinematicParcel/makeBasicKinematicParcelSubmodels.C
index 119de86cc6..4d7c1e908d 100644
--- a/src/lagrangian/turbulence/parcels/derived/basicKinematicParcel/makeBasicKinematicParcelSubmodels.C
+++ b/src/lagrangian/turbulence/parcels/derived/basicKinematicParcel/makeBasicKinematicParcelSubmodels.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -24,24 +24,13 @@ License
 \*---------------------------------------------------------------------------*/
 
 #include "basicKinematicCloud.H"
-
-#include "makeParcelCloudFunctionObjects.H"
-
-#include "GradientDispersionRAS.H"
-#include "StochasticDispersionRAS.H"
+#include "makeParcelTurbulenceDispersionModels.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 namespace Foam
 {
-    typedef basicKinematicCloud::kinematicCloudType kinematicCloudType_K;
-    defineNamedTemplateTypeNameAndDebug
-    (
-        DispersionRASModel<kinematicCloudType_K>,
-        0
-    );
-    makeDispersionModelType(GradientDispersionRAS, basicKinematicCloud);
-    makeDispersionModelType(StochasticDispersionRAS, basicKinematicCloud);
+    makeParcelTurbulenceDispersionModels(basicKinematicCloud);
 }
 
 
diff --git a/src/lagrangian/turbulence/parcels/derived/basicReactingMultiphaseParcel/makeBasicReactingMultiphaseParcelSubmodels.C b/src/lagrangian/turbulence/parcels/derived/basicReactingMultiphaseParcel/makeBasicReactingMultiphaseParcelSubmodels.C
index 08720bbc76..341bf7f349 100644
--- a/src/lagrangian/turbulence/parcels/derived/basicReactingMultiphaseParcel/makeBasicReactingMultiphaseParcelSubmodels.C
+++ b/src/lagrangian/turbulence/parcels/derived/basicReactingMultiphaseParcel/makeBasicReactingMultiphaseParcelSubmodels.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -25,40 +25,15 @@ License
 
 #include "basicReactingMultiphaseCloud.H"
 
-#include "makeParcelCloudFunctionObjects.H"
-
-#include "GradientDispersionRAS.H"
-#include "StochasticDispersionRAS.H"
-
-#include "BrownianMotionForce.H"
+#include "makeParcelTurbulenceDispersionModels.H"
+#include "makeThermoParcelTurbulenceForces.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 namespace Foam
 {
-    typedef basicReactingMultiphaseCloud::kinematicCloudType
-        kinematicCloudType_K;
-    defineNamedTemplateTypeNameAndDebug
-    (
-        DispersionRASModel<kinematicCloudType_K>,
-        0
-    );
-    makeDispersionModelType
-    (
-        GradientDispersionRAS,
-        basicReactingMultiphaseCloud
-    );
-    makeDispersionModelType
-    (
-        StochasticDispersionRAS,
-        basicReactingMultiphaseCloud
-    );
-
-    makeParticleForceModelType
-    (
-        BrownianMotionForce,
-        basicReactingMultiphaseCloud
-    );
+    makeThermoParcelTurbulenceForces(basicReactingMultiphaseCloud);
+    makeParcelTurbulenceDispersionModels(basicReactingMultiphaseCloud);
 }
 
 
diff --git a/src/lagrangian/turbulence/parcels/derived/basicReactingParcel/makeBasicReactingParcelSubmodels.C b/src/lagrangian/turbulence/parcels/derived/basicReactingParcel/makeBasicReactingParcelSubmodels.C
index da13ed2e1e..c259d405fd 100644
--- a/src/lagrangian/turbulence/parcels/derived/basicReactingParcel/makeBasicReactingParcelSubmodels.C
+++ b/src/lagrangian/turbulence/parcels/derived/basicReactingParcel/makeBasicReactingParcelSubmodels.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -25,27 +25,15 @@ License
 
 #include "basicReactingCloud.H"
 
-#include "makeParcelCloudFunctionObjects.H"
-
-#include "GradientDispersionRAS.H"
-#include "StochasticDispersionRAS.H"
-
-#include "BrownianMotionForce.H"
+#include "makeParcelTurbulenceDispersionModels.H"
+#include "makeThermoParcelTurbulenceForces.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 namespace Foam
 {
-    typedef basicReactingCloud::kinematicCloudType kinematicCloudType_K;
-    defineNamedTemplateTypeNameAndDebug
-    (
-        DispersionRASModel<kinematicCloudType_K>,
-        0
-    );
-    makeDispersionModelType(GradientDispersionRAS, basicReactingCloud);
-    makeDispersionModelType(StochasticDispersionRAS, basicReactingCloud);
-
-    makeParticleForceModelType(BrownianMotionForce, basicReactingCloud);
+    makeThermoParcelTurbulenceForces(basicReactingCloud);
+    makeParcelTurbulenceDispersionModels(basicReactingCloud);
 }
 
 
diff --git a/src/lagrangian/turbulence/parcels/derived/basicThermoParcel/makeBasicThermoParcelSubmodels.C b/src/lagrangian/turbulence/parcels/derived/basicThermoParcel/makeBasicThermoParcelSubmodels.C
index 5c71da29f5..d5c398dd8c 100644
--- a/src/lagrangian/turbulence/parcels/derived/basicThermoParcel/makeBasicThermoParcelSubmodels.C
+++ b/src/lagrangian/turbulence/parcels/derived/basicThermoParcel/makeBasicThermoParcelSubmodels.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -25,27 +25,15 @@ License
 
 #include "basicThermoCloud.H"
 
-#include "makeParcelCloudFunctionObjects.H"
-
-#include "GradientDispersionRAS.H"
-#include "StochasticDispersionRAS.H"
-
-#include "BrownianMotionForce.H"
+#include "makeParcelTurbulenceDispersionModels.H"
+#include "makeThermoParcelTurbulenceForces.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 namespace Foam
 {
-    typedef basicThermoCloud::kinematicCloudType kinematicCloudType_K;
-    defineNamedTemplateTypeNameAndDebug
-    (
-        DispersionRASModel<kinematicCloudType_K>,
-        0
-    );
-    makeDispersionModelType(GradientDispersionRAS, basicThermoCloud);
-    makeDispersionModelType(StochasticDispersionRAS, basicThermoCloud);
-
-    makeParticleForceModelType(BrownianMotionForce, basicThermoCloud);
+    makeThermoParcelTurbulenceForces(basicThermoCloud);
+    makeParcelTurbulenceDispersionModels(basicThermoCloud);
 }
 
 
diff --git a/src/lagrangian/turbulence/parcels/include/makeParcelDispersionModels.H b/src/lagrangian/turbulence/parcels/include/makeParcelTurbulenceDispersionModels.H
similarity index 90%
rename from src/lagrangian/turbulence/parcels/include/makeParcelDispersionModels.H
rename to src/lagrangian/turbulence/parcels/include/makeParcelTurbulenceDispersionModels.H
index be281ecb1e..d2428ba4b2 100644
--- a/src/lagrangian/turbulence/parcels/include/makeParcelDispersionModels.H
+++ b/src/lagrangian/turbulence/parcels/include/makeParcelTurbulenceDispersionModels.H
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -23,8 +23,8 @@ License
 
 \*---------------------------------------------------------------------------*/
 
-#ifndef makeParcelDispersionModels_h
-#define makeParcelDispersionModels_h
+#ifndef makeParcelTurbulenceDispersionModels_h
+#define makeParcelTurbulenceDispersionModels_h
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -33,7 +33,7 @@ License
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-#define makeParcelDispersionModels(CloudType)                                 \
+#define makeParcelTurbulenceDispersionModels(CloudType)                       \
                                                                               \
     typedef CloudType::kinematicCloudType kinematicCloudType;                 \
     defineNamedTemplateTypeNameAndDebug                                       \
diff --git a/src/lagrangian/turbulence/parcels/include/makeThermoParcelTurbulenceForces.H b/src/lagrangian/turbulence/parcels/include/makeThermoParcelTurbulenceForces.H
new file mode 100644
index 0000000000..d7581110b5
--- /dev/null
+++ b/src/lagrangian/turbulence/parcels/include/makeThermoParcelTurbulenceForces.H
@@ -0,0 +1,44 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2014 OpenFOAM Foundation
+     \\/     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/>.
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef makeThermoParcelTurbulenceForces_h
+#define makeThermoParcelTurbulenceForces_h
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#include "BrownianMotionForce.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#define makeThermoParcelTurbulenceForces(CloudType)                           \
+                                                                              \
+    makeParticleForceModelType(BrownianMotionForce, CloudType);
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/lagrangian/turbulence/submodels/Kinematic/DispersionModel/GradientDispersionRAS/GradientDispersionRAS.C b/src/lagrangian/turbulence/submodels/Kinematic/DispersionModel/GradientDispersionRAS/GradientDispersionRAS.C
index c96a5c5b3d..06fa810ac9 100644
--- a/src/lagrangian/turbulence/submodels/Kinematic/DispersionModel/GradientDispersionRAS/GradientDispersionRAS.C
+++ b/src/lagrangian/turbulence/submodels/Kinematic/DispersionModel/GradientDispersionRAS/GradientDispersionRAS.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -25,6 +25,7 @@ License
 
 #include "GradientDispersionRAS.H"
 #include "demandDrivenData.H"
+#include "fvcGrad.H"
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
diff --git a/src/meshTools/AMIInterpolation/AMIInterpolation/AMIMethod/mapNearestAMI/mapNearestAMI.C b/src/meshTools/AMIInterpolation/AMIInterpolation/AMIMethod/mapNearestAMI/mapNearestAMI.C
index 61ed01727e..be6f590fa0 100644
--- a/src/meshTools/AMIInterpolation/AMIInterpolation/AMIMethod/mapNearestAMI/mapNearestAMI.C
+++ b/src/meshTools/AMIInterpolation/AMIInterpolation/AMIMethod/mapNearestAMI/mapNearestAMI.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2013-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -310,16 +310,19 @@ void Foam::mapNearestAMI<SourcePatch, TargetPatch>::calculate
         {
             label srcFaceI = findMappedSrcFace(tgtFaceI, tgtAddr);
 
-            // note - reversed search from src->tgt to tgt->src
-            findNearestFace
-            (
-                this->tgtPatch_,
-                this->srcPatch_,
-                tgtFaceI,
-                srcFaceI
-            );
+            if (srcFaceI >= 0)
+            {
+                // note - reversed search from src->tgt to tgt->src
+                findNearestFace
+                (
+                    this->tgtPatch_,
+                    this->srcPatch_,
+                    tgtFaceI,
+                    srcFaceI
+                );
 
-            tgtAddr[tgtFaceI].append(srcFaceI);
+                tgtAddr[tgtFaceI].append(srcFaceI);
+            }
         }
     }
 
diff --git a/src/postProcessing/functionObjects/forces/forces/forces.C b/src/postProcessing/functionObjects/forces/forces/forces.C
index 84863f7f7b..488a71a881 100644
--- a/src/postProcessing/functionObjects/forces/forces/forces.C
+++ b/src/postProcessing/functionObjects/forces/forces/forces.C
@@ -912,6 +912,14 @@ void Foam::forces::calcForcesMoment()
         const HashTable<const porosityModel*> models =
             obr_.lookupClass<porosityModel>();
 
+        if (models.empty())
+        {
+            WarningIn("void Foam::forces::calcForcesMoment()")
+                << "Porosity effects requested, but no porosity models found "
+                << "in the database"
+                << endl;
+        }
+
         forAllConstIter(HashTable<const porosityModel*>, models, iter)
         {
             const porosityModel& pm = *iter();
diff --git a/src/sampling/meshToMeshInterpolation/meshToMesh/meshToMeshTemplates.C b/src/sampling/meshToMeshInterpolation/meshToMesh/meshToMeshTemplates.C
index 5a9926046f..10a1c507c4 100644
--- a/src/sampling/meshToMeshInterpolation/meshToMesh/meshToMeshTemplates.C
+++ b/src/sampling/meshToMeshInterpolation/meshToMesh/meshToMeshTemplates.C
@@ -26,6 +26,7 @@ License
 #include "fvMesh.H"
 #include "volFields.H"
 #include "directFvPatchFieldMapper.H"
+#include "calculatedFvPatchField.H"
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
@@ -381,30 +382,21 @@ Foam::meshToMesh::mapSrcToTgt
 
     const fvMesh& tgtMesh = static_cast<const fvMesh&>(tgtRegion_);
 
-    const typename fieldType::GeometricBoundaryField& bfld =
+    const fvBoundaryMesh& tgtBm = tgtMesh.boundary();
+    const typename fieldType::GeometricBoundaryField& srcBfld =
         field.boundaryField();
 
-    PtrList<fvPatchField<Type> > patchFields(bfld.size());
+    wordList patchTypes(tgtBm.size(), calculatedFvPatchField<Type>::typeName);
 
-    // constuct tgt boundary patch types as copy of field boundary types
+    // constuct tgt boundary patch types as copy of 'field' boundary types
     // note: this will provide place holders for fields with additional
     // entries, but these values will need to be reset
-    forAll(bfld, patchI)
+    forAll(tgtPatchID_, i)
     {
-        patchFields.set
-        (
-            patchI,
-            fvPatchField<Type>::New
-            (
-                bfld[patchI],
-                tgtMesh.boundary()[patchI],
-                DimensionedField<Type, volMesh>::null(),
-                directFvPatchFieldMapper
-                (
-                    labelList(tgtMesh.boundary()[patchI].size(), -1)
-                )
-            )
-        );
+        label srcPatchI = srcPatchID_[i];
+        label tgtPatchI = tgtPatchID_[i];
+
+        patchTypes[tgtPatchI] = srcBfld[srcPatchI].type();
     }
 
     tmp<fieldType> tresult
@@ -420,9 +412,8 @@ Foam::meshToMesh::mapSrcToTgt
                 IOobject::NO_WRITE
             ),
             tgtMesh,
-            field.dimensions(),
-            Field<Type>(tgtMesh.nCells(), pTraits<Type>::zero),
-            patchFields
+            dimensioned<Type>("0", field.dimensions(), pTraits<Type>::zero),
+            patchTypes
         )
     );
 
@@ -517,30 +508,21 @@ Foam::meshToMesh::mapTgtToSrc
 
     const fvMesh& srcMesh = static_cast<const fvMesh&>(srcRegion_);
 
-    const typename fieldType::GeometricBoundaryField& bfld =
+    const fvBoundaryMesh& srcBm = srcMesh.boundary();
+    const typename fieldType::GeometricBoundaryField& tgtBfld =
         field.boundaryField();
 
-    PtrList<fvPatchField<Type> > patchFields(bfld.size());
+    wordList patchTypes(srcBm.size(), calculatedFvPatchField<Type>::typeName);
 
-    // constuct src boundary patch types as copy of field boundary types
+    // constuct src boundary patch types as copy of 'field' boundary types
     // note: this will provide place holders for fields with additional
     // entries, but these values will need to be reset
-    forAll(bfld, patchI)
+    forAll(srcPatchID_, i)
     {
-        patchFields.set
-        (
-            patchI,
-            fvPatchField<Type>::New
-            (
-                bfld[patchI],
-                srcMesh.boundary()[patchI],
-                DimensionedField<Type, volMesh>::null(),
-                directFvPatchFieldMapper
-                (
-                    labelList(srcMesh.boundary()[patchI].size(), -1)
-                )
-            )
-        );
+        label srcPatchI = srcPatchID_[i];
+        label tgtPatchI = tgtPatchID_[i];
+
+        patchTypes[srcPatchI] = tgtBfld[tgtPatchI].type();
     }
 
     tmp<fieldType> tresult
@@ -556,9 +538,8 @@ Foam::meshToMesh::mapTgtToSrc
                 IOobject::NO_WRITE
             ),
             srcMesh,
-            field.dimensions(),
-            Field<Type>(srcMesh.nCells(), pTraits<Type>::zero),
-            patchFields
+            dimensioned<Type>("0", field.dimensions(), pTraits<Type>::zero),
+            patchTypes
         )
     );
 
diff --git a/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/system/snappyHexMeshDict b/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/system/snappyHexMeshDict
index 89829fcc32..1faa6e4939 100644
--- a/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/system/snappyHexMeshDict
+++ b/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/system/snappyHexMeshDict
@@ -275,7 +275,7 @@ addLayersControls
         {
             nSurfaceLayers 3;
         }
-        "igloo_.*"
+        igloo
         {
             nSurfaceLayers 1;
         }