diff --git a/src/lagrangian/dsmc/clouds/Templates/DsmcCloud/DsmcCloud.C b/src/lagrangian/dsmc/clouds/Templates/DsmcCloud/DsmcCloud.C
index e748f72077..80e0dea37d 100644
--- a/src/lagrangian/dsmc/clouds/Templates/DsmcCloud/DsmcCloud.C
+++ b/src/lagrangian/dsmc/clouds/Templates/DsmcCloud/DsmcCloud.C
@@ -28,6 +28,7 @@ License
 #include "WallInteractionModel.H"
 #include "InflowBoundaryModel.H"
 #include "constants.H"
+#include "zeroGradientFvPatchFields.H"
 
 using namespace Foam::constant;
 
@@ -459,6 +460,8 @@ void Foam::DsmcCloud<ParcelType>::collisions()
 
     reduce(collisionCandidates, sumOp<label>());
 
+    sigmaTcRMax_.correctBoundaryConditions();
+
     if (collisionCandidates)
     {
         Info<< "    Collisions                      = "
@@ -550,6 +553,8 @@ void Foam::DsmcCloud<ParcelType>::calculateFields()
     rhoM *= nParticle_/mesh().cellVolumes();
     rhoM_.correctBoundaryConditions();
 
+    dsmcRhoN_.correctBoundaryConditions();
+
     linearKE *= nParticle_/mesh().cellVolumes();
     linearKE_.correctBoundaryConditions();
 
@@ -851,7 +856,8 @@ Foam::DsmcCloud<ParcelType>::DsmcCloud
             IOobject::AUTO_WRITE
         ),
         mesh_,
-        dimensionedScalar("zero",  dimensionSet(0, 3, -1, 0, 0), 0.0)
+        dimensionedScalar("zero",  dimensionSet(0, 3, -1, 0, 0), 0.0),
+        zeroGradientFvPatchScalarField::typeName
     ),
     collisionSelectionRemainder_(),
     q_
diff --git a/src/lagrangian/dsmc/clouds/Templates/DsmcCloud/DsmcCloud.H b/src/lagrangian/dsmc/clouds/Templates/DsmcCloud/DsmcCloud.H
index 1df4935f68..aba26e7bf8 100644
--- a/src/lagrangian/dsmc/clouds/Templates/DsmcCloud/DsmcCloud.H
+++ b/src/lagrangian/dsmc/clouds/Templates/DsmcCloud/DsmcCloud.H
@@ -83,10 +83,10 @@ class DsmcCloud
         //- Dictionary of particle properties
         IOdictionary particleProperties_;
 
-        //- A list of unique instances of molecule types in the simulation.
-        // The position of an entry in the list maps to the label identifying
-        // the typeId, i.e. where typeIdList_ = (N2 O2 CO2)
-        // N2 has typeId label = 0, O2 = 1, CO2 = 2.
+        //- A list of unique instances of molecule types in the
+        //  simulation.  The position of an entry in the list maps to
+        //  the label identifying the typeId, i.e. where typeIdList_ =
+        //  (N2 O2 CO2) N2 has typeId label = 0, O2 = 1, CO2 = 2.
         List<word> typeIdList_;
 
         //- Number of real atoms/molecules represented by a parcel
@@ -95,9 +95,9 @@ class DsmcCloud
         //- A data structure holding which particles are in which cell
         List<DynamicList<ParcelType*> > cellOccupancy_;
 
-        //- An IOField holding the value of (sigmaT * cR)max for each cell (see
-        // Bird p220). Initialised with the parcels, updated as required, and
-        // read in on start/restart.
+        //- A field holding the value of (sigmaT * cR)max for each
+        //  cell (see Bird p220). Initialised with the parcels,
+        //  updated as required, and read in on start/restart.
         volScalarField sigmaTcRMax_;
 
         //- A field holding the remainder from the previous collision selections