diff --git a/applications/solvers/multiphase/interFoam/interMixingFoam/threePhaseInterfaceProperties/threePhaseInterfaceProperties.C b/applications/solvers/multiphase/interFoam/interMixingFoam/threePhaseInterfaceProperties/threePhaseInterfaceProperties.C
index 9a1c7bc215..dce55234e7 100644
--- a/applications/solvers/multiphase/interFoam/interMixingFoam/threePhaseInterfaceProperties/threePhaseInterfaceProperties.C
+++ b/applications/solvers/multiphase/interFoam/interMixingFoam/threePhaseInterfaceProperties/threePhaseInterfaceProperties.C
@@ -145,7 +145,7 @@ void Foam::threePhaseInterfaceProperties::calculateK()
     // Complex expression for curvature.
     // Correction is formally zero but numerically non-zero.
     // volVectorField nHat = gradAlpha/(mag(gradAlpha) + deltaN_);
-    // nHat.boundaryField() = nHatfv.boundaryField();
+    // nHat.boundaryFieldRef() = nHatfv.boundaryField();
     // K_ = -fvc::div(nHatf_) + (nHat & fvc::grad(nHatfv) & nHat);
 }
 
diff --git a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyLayerDriver.C b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyLayerDriver.C
index 54052a3078..8b626896ea 100644
--- a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyLayerDriver.C
+++ b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyLayerDriver.C
@@ -2824,7 +2824,7 @@ bool Foam::snappyLayerDriver::writeLayerData
                     faceRealThickness
                 );
 
-                // Convert patchReal to relavtive thickness
+                // Convert patchReal to relative thickness
                 scalarField pfld(patchReal.size(), 0.0);
                 forAll(patchReal, i)
                 {
@@ -2834,7 +2834,7 @@ bool Foam::snappyLayerDriver::writeLayerData
                     }
                 }
 
-                fld.boundaryField()[patchi] == pfld;
+                fld.boundaryFieldRef()[patchi] == pfld;
             }
             Info<< indent << fld.name()
                 << " : overall layer thickness (fraction"
diff --git a/src/transportModels/interfaceProperties/interfaceProperties.C b/src/transportModels/interfaceProperties/interfaceProperties.C
index 95ddb6aac4..40761a3a59 100644
--- a/src/transportModels/interfaceProperties/interfaceProperties.C
+++ b/src/transportModels/interfaceProperties/interfaceProperties.C
@@ -142,7 +142,7 @@ void Foam::interfaceProperties::calculateK()
     volVectorField nHat(gradAlpha/(mag(gradAlpha) + deltaN_));
     forAll(nHat.boundaryField(), patchi)
     {
-        nHat.boundaryField()[patchi] = nHatfv.boundaryField()[patchi];
+        nHat.boundaryFieldRef()[patchi] = nHatfv.boundaryField()[patchi];
     }
 
     K_ = -fvc::div(nHatf_) + (nHat & fvc::grad(nHatfv) & nHat);