diff --git a/src/OpenFOAM/db/Time/Time.C b/src/OpenFOAM/db/Time/Time.C
index e5b8e5f2fd..9d21d32b35 100644
--- a/src/OpenFOAM/db/Time/Time.C
+++ b/src/OpenFOAM/db/Time/Time.C
@@ -950,6 +950,12 @@ bool Foam::Time::stopAt(const stopAtControls stopCtrl) const
 }
 
 
+bool Foam::Time::isAdjustTimeStep() const
+{
+    return controlDict_.lookupOrDefault("adjustTimeStep", false) ? true : false;
+}
+
+
 void Foam::Time::setTime(const Time& t)
 {
     value() = t.value();
diff --git a/src/OpenFOAM/db/Time/Time.H b/src/OpenFOAM/db/Time/Time.H
index 7ab154f604..9b8f37e252 100644
--- a/src/OpenFOAM/db/Time/Time.H
+++ b/src/OpenFOAM/db/Time/Time.H
@@ -422,6 +422,9 @@ public:
             //- Return the stop control information
             virtual stopAtControls stopAt() const;
 
+            //- Return true if adjustTimeStep is true
+            virtual bool isAdjustTimeStep() const;
+
             //- Return the list of function objects
             const functionObjectList& functionObjects() const
             {
diff --git a/src/finiteVolume/Make/files b/src/finiteVolume/Make/files
index 791ef91c32..81500b3e13 100644
--- a/src/finiteVolume/Make/files
+++ b/src/finiteVolume/Make/files
@@ -200,6 +200,7 @@ $(derivedFvPatchFields)/translatingWallVelocity/translatingWallVelocityFvPatchVe
 $(derivedFvPatchFields)/turbulentDFSEMInlet/turbulentDFSEMInletFvPatchVectorField.C
 $(derivedFvPatchFields)/turbulentDFSEMInlet/eddy/eddy.C
 $(derivedFvPatchFields)/turbulentDFSEMInlet/eddy/eddyIO.C
+$(derivedFvPatchFields)/turbulentDigitalFilterInlet/turbulentDigitalFilterInletFvPatchVectorField.C
 $(derivedFvPatchFields)/turbulentInlet/turbulentInletFvPatchFields.C
 $(derivedFvPatchFields)/turbulentIntensityKineticEnergyInlet/turbulentIntensityKineticEnergyInletFvPatchScalarField.C
 $(derivedFvPatchFields)/uniformDensityHydrostaticPressure/uniformDensityHydrostaticPressureFvPatchScalarField.C
diff --git a/src/finiteVolume/fields/fvPatchFields/derived/turbulentDigitalFilterInlet/turbulentDigitalFilterInletFvPatchVectorField.C b/src/finiteVolume/fields/fvPatchFields/derived/turbulentDigitalFilterInlet/turbulentDigitalFilterInletFvPatchVectorField.C
new file mode 100644
index 0000000000..c1a2b841ab
--- /dev/null
+++ b/src/finiteVolume/fields/fvPatchFields/derived/turbulentDigitalFilterInlet/turbulentDigitalFilterInletFvPatchVectorField.C
@@ -0,0 +1,1261 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2019 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 "turbulentDigitalFilterInletFvPatchVectorField.H"
+#include "volFields.H"
+#include "mathematicalConstants.H"
+#include "addToRunTimeSelectionTable.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+const Foam::Enum
+<
+    Foam::turbulentDigitalFilterInletFvPatchVectorField::variantType
+>
+Foam::turbulentDigitalFilterInletFvPatchVectorField::variantNames
+({
+    { variantType::DIGITAL_FILTER, "digitalFilter" },
+    { variantType::DIGITAL_FILTER, "DFM" },
+    { variantType::FORWARD_STEPWISE, "forwardStepwise" },
+    { variantType::FORWARD_STEPWISE, "reducedDigitalFilter" },
+    { variantType::FORWARD_STEPWISE, "FSM" }
+});
+
+// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
+
+const Foam::pointToPointPlanarInterpolation&
+Foam::turbulentDigitalFilterInletFvPatchVectorField::patchMapper() const
+{
+    // Initialise interpolation (2D planar interpolation by triangulation)
+    if (mapperPtr_.empty())
+    {
+        // Reread values and interpolate
+        fileName samplePointsFile
+        (
+            this->db().time().path()
+           /this->db().time().caseConstant()
+           /"boundaryData"
+           /this->patch().name()
+           /"points"
+        );
+
+        pointField samplePoints((IFstream(samplePointsFile)()));
+
+        // tbd: run-time selection
+        bool nearestOnly =
+        (
+           !mapMethod_.empty()
+         && mapMethod_ != "planarInterpolation"
+        );
+
+        // Allocate the interpolator
+        mapperPtr_.reset
+        (
+            new pointToPointPlanarInterpolation
+            (
+                samplePoints,
+                this->patch().patch().faceCentres(),
+                perturb_,
+                nearestOnly
+            )
+        );
+    }
+
+    return *mapperPtr_;
+}
+
+
+Foam::List<Foam::Pair<Foam::label>>
+Foam::turbulentDigitalFilterInletFvPatchVectorField::patchIndexPairs()
+{
+    // Get patch normal direction into the domain
+    const vector nf(computePatchNormal());
+
+    // Find the second local coordinate direction
+    direction minCmpt = -1;
+    scalar minMag = VGREAT;
+    for (direction cmpt = 0; cmpt < pTraits<vector>::nComponents; ++cmpt)
+    {
+        scalar s = mag(nf[cmpt]);
+        if (s < minMag)
+        {
+            minMag = s;
+            minCmpt = cmpt;
+        }
+    }
+
+    // Create the second local coordinate direction
+    vector e2(Zero);
+    e2[minCmpt] = 1.0;
+
+    // Remove normal component
+    e2 -= (nf&e2)*nf;
+
+    // Create the local coordinate system
+    coordSystem::cartesian cs
+    (
+        Zero,   // origin
+        nf,     // normal
+        e2      // 0-axis
+    );
+
+    // Convert patch points into local coordinate system
+    const pointField localPos
+    (
+        cs.localPosition
+        (
+            pointField
+            (
+                patch().patch().points(),
+                patch().patch().meshPoints()
+            )
+        )
+    );
+
+    const boundBox bb(localPos);
+
+    // Normalise to (i, j) coordinates
+    const Vector2D<label> n(planeDivisions_.first(), planeDivisions_.second());
+    invDelta_[0] = n[0]/bb.span()[0];
+    invDelta_[1] = n[1]/bb.span()[1];
+    const point localMinPt(bb.min());
+
+    // Compute virtual-actual patch index pairs
+    List<Pair<label>> indexPairs(this->size(), Pair<label>(Zero, Zero));
+
+    // Virtual turbulence plane indices
+    label j = 0;
+    label k = 0;
+
+    forAll(*this, facei)
+    {
+        const scalar& centre0 = localPos[facei][0];
+        const scalar& centre1 = localPos[facei][1];
+
+        j = label((centre0 - localMinPt[0])*invDelta_[0]);
+        k = label((centre1 - localMinPt[1])*invDelta_[1]);
+
+        indexPairs[facei] = Pair<label>(facei, k*n[0] + j);
+    }
+
+    return indexPairs;
+}
+
+
+void Foam::turbulentDigitalFilterInletFvPatchVectorField::
+checkRTensorRealisable() const
+{
+    const vectorField& faceCentres = this->patch().patch().faceCentres();
+
+    forAll(R_, facei)
+    {
+        if (R_[facei].xx() <= 0)
+        {
+            FatalErrorInFunction
+                << "Reynolds stress tensor component Rxx cannot be negative"
+                   "or zero, where Rxx = " << R_[facei].xx() << " at the face "
+                   "centre = " << faceCentres[facei] << exit(FatalError);
+        }
+
+        if (R_[facei].yy() < 0 || R_[facei].zz() < 0)
+        {
+            FatalErrorInFunction
+                << "Reynolds stress tensor components Ryy or Rzz cannot be"
+                << "negative where Ryy = " << R_[facei].yy() << ", and Rzz = "
+                << R_[facei].zz() << " at the face centre = "
+                << faceCentres[facei] << exit(FatalError);
+        }
+
+        scalar term0 = R_[facei].xx()*R_[facei].yy() - sqr(R_[facei].xy());
+
+        if (term0 <= 0)
+        {
+            FatalErrorInFunction
+                << "Reynolds stress tensor component group, Rxx*Ryy - Rxy^2"
+                << "cannot be negative or zero at the face centre = "
+                << faceCentres[facei] << exit(FatalError);
+        }
+
+        scalar term1 = R_[facei].zz() - sqr(R_[facei].xz())/R_[facei].xx();
+        scalar term2 =
+            sqr(R_[facei].yz() - R_[facei].xy()*R_[facei].xz()
+           /(R_[facei].xx()*term0));
+        scalar term3 = term1 - term2;
+
+        if (term3 < 0)
+        {
+            FatalErrorInFunction
+                << "Reynolds stress tensor component group,"
+                << "Rzz - Rxz^2/Rxx - (Ryz - Rxy*Rxz/(Rxx*(Rxx*Ryy - Rxy^2)))^2"
+                << "cannot be negative at the face centre = "
+                << faceCentres[facei] << exit(FatalError);
+        }
+    }
+
+    #ifdef FULLDEBUG
+    Info<< "Ends: checkRTensorRealisable()."
+        << " Reynolds tensor (on patch) is consistent." << nl;
+    #endif
+}
+
+
+Foam::symmTensorField Foam::turbulentDigitalFilterInletFvPatchVectorField::
+computeLundWuSquires() const
+{
+    checkRTensorRealisable();
+
+    symmTensorField LundWuSquires(symmTensorField(R_.size()));
+
+    forAll(R_, facei)
+    {
+        const symmTensor& R = R_[facei];
+        symmTensor& lws = LundWuSquires[facei];
+
+        // (Klein et al., 2003, Eq. 5)
+        lws.xx() = Foam::sqrt(R.xx());
+        lws.xy() = R.xy()/lws.xx();
+        lws.xz() = R.xz()/lws.xx();
+        lws.yy() = Foam::sqrt(R.yy() - sqr(lws.xy()));
+        lws.yz() = (R.yz() - lws.xy()*lws.xz())/lws.yy();
+        lws.zz() = Foam::sqrt(R.zz() - sqr(lws.xz()) - sqr(lws.yz()));
+    }
+
+    #ifdef FULLDEBUG
+    Info<< "Ends: computeLundWuSquires()." << nl;
+    #endif
+
+    return LundWuSquires;
+}
+
+
+Foam::vector Foam::turbulentDigitalFilterInletFvPatchVectorField::
+computePatchNormal() const
+{
+    vector patchNormal(-gAverage(patch().nf()));
+    return patchNormal.normalise();
+}
+
+
+Foam::scalar Foam::turbulentDigitalFilterInletFvPatchVectorField::
+computeInitialFlowRate() const
+{
+    const vector patchNormal(computePatchNormal());
+    return gSum((UMean_ & patchNormal)*patch().magSf());
+}
+
+
+void Foam::turbulentDigitalFilterInletFvPatchVectorField::convertToTimeScale
+(
+    tensor& L
+) const
+{
+    if (isTaylorHypot_)
+    {
+        forAll(L.x(), i)
+        {
+            L[i] /= patchNormalSpeed_;
+        }
+
+    #ifdef FULLDEBUG
+    Info<< "Ends: convertToTimeScale()."
+        << "Streamwise integral length scales are converted to time scales via"
+        << "Taylor's frozen turbulence hypothesis" << nl;
+    #endif
+    }
+}
+
+
+Foam::tensor Foam::turbulentDigitalFilterInletFvPatchVectorField::
+convertScalesToGridUnits
+(
+    const tensor& L
+) const
+{
+    tensor Ls(L);
+    convertToTimeScale(Ls);
+    const scalar invDeltaT = 1.0/db().time().deltaTValue();
+
+    Ls.row(0, Ls.x()*invDeltaT);
+    Ls.row(1, Ls.y()*invDelta_[0]);
+    Ls.row(2, Ls.z()*invDelta_[1]);
+
+    #ifdef FULLDEBUG
+    Info<< "Ends: convertScalesToGridUnits()."
+        << "Units of input length scales are converted from metres to"
+        << "virtual-patch cell size/time-step" << nl;
+    #endif
+
+    return Ls;
+}
+
+
+Foam::List<Foam::label> Foam::turbulentDigitalFilterInletFvPatchVectorField::
+initLenRandomBox() const
+{
+    label initValue = 0;
+    label rangeModifier = 0;
+
+    if (variant_ == variantType::FORWARD_STEPWISE)
+    {
+        // Initialise with 1 since x-dir possess 1 node with this variant
+        initValue = pTraits<label>::nComponents;
+        rangeModifier = pTraits<vector>::nComponents;
+    }
+
+    List<label> lenRandomBox(pTraits<tensor>::nComponents, initValue);
+    Vector<label> lenGrid
+    (
+        pTraits<label>::nComponents,
+        planeDivisions_.first(),
+        planeDivisions_.second()
+    );
+
+    // Initialise: Main convenience factor, lenRandomBox_
+    for
+    (
+        const label& i
+      : labelRange(rangeModifier, pTraits<tensor>::nComponents - rangeModifier)
+    )
+    {
+        // Slicing index
+        const label sliceI = label(i/pTraits<vector>::nComponents);
+
+        // Refer to 'computeFilterCoeffs()'
+        const label n = ceil(L_[i]);
+        const label twiceN = 4*n;       
+
+        // Initialise: Random-number set sizes
+        lenRandomBox[i] = lenGrid[sliceI] + twiceN;
+    }
+
+    return lenRandomBox;
+}
+
+
+Foam::List<Foam::label> Foam::turbulentDigitalFilterInletFvPatchVectorField::
+initBoxFactors2D() const
+{
+    List<label> boxFactors2D(pTraits<vector>::nComponents);
+
+    for (direction dir = 0; dir < pTraits<vector>::nComponents; ++dir)
+    {
+        boxFactors2D[dir] =
+            lenRandomBox_[pTraits<vector>::nComponents + dir]
+           *lenRandomBox_[pTraits<symmTensor>::nComponents + dir];
+    }
+
+    return boxFactors2D;
+}
+
+
+Foam::List<Foam::label> Foam::turbulentDigitalFilterInletFvPatchVectorField::
+initBoxFactors3D() const
+{
+    List<label> boxFactors3D(pTraits<vector>::nComponents);
+
+    for (direction dir = 0; dir < pTraits<vector>::nComponents; ++dir)
+    {
+        boxFactors3D[dir] = randomBoxFactors2D_[dir]*lenRandomBox_[dir];
+    }
+
+    return boxFactors3D;
+}
+
+
+Foam::List<Foam::label> Foam::turbulentDigitalFilterInletFvPatchVectorField::
+initBoxPlaneFactors() const
+{
+    List<label> planeFactors(pTraits<vector>::nComponents);
+
+    for (direction dir = 0; dir < pTraits<vector>::nComponents; ++dir)
+    {
+        planeFactors[dir] =
+            randomBoxFactors2D_[dir]*(lenRandomBox_[dir] - pTraits<label>::one);
+    }
+
+    return planeFactors;
+}
+
+
+Foam::List<Foam::List<Foam::scalar>>
+Foam::turbulentDigitalFilterInletFvPatchVectorField::fillRandomBox()
+{
+    List<List<scalar>> randomBox(pTraits<vector>::nComponents, List<scalar>());
+
+    // Initialise: Remaining convenience factors for (e1 e2 e3)
+    for (direction dir = 0; dir < pTraits<vector>::nComponents; ++dir)
+    {
+        // Initialise: Random-box content with random-number sets
+        randomBox[dir] =
+            generateRandomSet<List<scalar>, scalar>(randomBoxFactors3D_[dir]);
+    }
+
+    #ifdef FULLDEBUG
+    Info<< "Ends: fillRandomBox()."
+        << "Random-number box is created." << nl;
+    #endif
+
+    return randomBox;
+}
+
+
+Foam::List<Foam::List<Foam::scalar>>
+Foam::turbulentDigitalFilterInletFvPatchVectorField::computeFilterCoeffs() const
+{
+    List<List<scalar>> filterCoeffs
+    (
+        pTraits<tensor>::nComponents,
+        List<scalar>(1, 1.0)
+    );
+
+    label initValue = 0;
+
+    if(variant_ == variantType::FORWARD_STEPWISE)
+    {
+        initValue = pTraits<vector>::nComponents;
+    }
+
+    for (direction dir = initValue; dir < pTraits<tensor>::nComponents; ++dir)
+    {
+        // The smallest filter width larger than length scale
+        // (Klein et al., 2003, 'n' in Eq. 15)
+        const label n  = ceil(L_[dir]);
+
+        // Full filter-width (except mid-zero) according to the condition
+        // (Klein et al., 2003, Eq. 15 whereat N is minimum =2n)
+        const label twiceN = 4*n;
+
+        // Initialise filter-coeffs containers with full filter-width size
+        // Extra elem is mid-zero within [-N, N]
+        filterCoeffs[dir] = List<scalar>(twiceN + 1, Zero);
+
+        // First element: -N within [-N, N]
+        const scalar initElem = -2*scalar(n);
+
+        // Container initialised with [-N, N]
+        // (Klein et al., 2003, p. 658, Item-b)
+        std::iota
+        (
+            filterCoeffs[dir].begin(),
+            filterCoeffs[dir].end(),
+            initElem
+        );
+
+        // Compute filter-coeffs
+        // (Klein et al., 2003, Eq. 14 (Gaussian))
+        // (Bercin et al., 2018, Fig. 9 (Exponential))
+        List<scalar> fTemp(filterCoeffs[dir]);
+        scalar fSum = 0;
+        const scalar nSqr = n*n;
+
+        if (isGaussian_)
+        {
+            fTemp = sqr(exp(modelConst_*sqr(fTemp)/nSqr));
+            fSum = Foam::sqrt(sum(fTemp));
+
+            filterCoeffs[dir] =
+                exp(modelConst_*sqr(filterCoeffs[dir])/nSqr)/fSum;
+        }
+        else
+        {
+            fTemp = sqr(exp(modelConst_*mag(fTemp)/n));
+            fSum = Foam::sqrt(sum(fTemp));
+
+            filterCoeffs[dir] = exp(modelConst_*mag(filterCoeffs[dir])/n)/fSum;
+        }
+    }
+
+    #ifdef FULLDEBUG
+    Info<< "Ends: computeFilterCoeffs()."
+        << " Filter coefficients are computed." << nl;
+    #endif
+
+    return filterCoeffs;
+}
+
+
+void Foam::turbulentDigitalFilterInletFvPatchVectorField::rndShiftRefill()
+{
+    forAll(randomBox_, dir)
+    {
+        List<scalar>& r = randomBox_[dir];
+
+        // Shift forward from the back to the front / First Out
+        inplaceRotateList(r, randomBoxFactors2D_[dir]);
+
+        // Refill the back with a new random-number set / First In
+        for (label i = 0; i < randomBoxFactors2D_[dir]; ++i)
+        {
+            r[i] = rndGen_.GaussNormal<scalar>();
+        }
+    }
+}
+
+
+void Foam::turbulentDigitalFilterInletFvPatchVectorField::mapFilteredRandomBox
+(
+    vectorField& U
+)
+{
+    for (const auto& x : indexPairs_)
+    {
+        const label xf = x.first();
+        const label xs = x.second();
+        U[xf][0] = filteredRandomBox_[0][xs];
+        U[xf][1] = filteredRandomBox_[1][xs];
+        U[xf][2] = filteredRandomBox_[2][xs];
+    }
+}
+
+
+void Foam::turbulentDigitalFilterInletFvPatchVectorField::embedOnePointCorrs
+(
+    vectorField& U
+) const
+{
+    forAll(LundWuSquires_, facei)
+    {
+        vector& Us = U[facei];
+        const symmTensor& lws = LundWuSquires_[facei];
+
+        // (Klein et al. p. 658, Item-e)
+        Us.z() = Us.x()*lws.xz() + Us.y()*lws.yz() + Us.z()*lws.zz();
+        Us.y() = Us.x()*lws.xy() + Us.y()*lws.yy();
+        Us.x() = Us.x()*lws.xx();
+    }
+}
+
+
+void Foam::turbulentDigitalFilterInletFvPatchVectorField::embedMeanVelocity
+(
+    vectorField& U
+) const
+{
+    U += UMean_;
+}
+
+
+void Foam::turbulentDigitalFilterInletFvPatchVectorField::correctFlowRate
+(
+    vectorField& U
+) const
+{
+    U *= (initialFlowRate_/gSum(U & -patch().Sf()));
+}
+
+
+void Foam::turbulentDigitalFilterInletFvPatchVectorField::embedTwoPointCorrs()
+{
+    for (direction dir = 0; dir < pTraits<vector>::nComponents; ++dir)
+    {
+        List<scalar>& in = randomBox_[dir];
+        List<scalar>& out = filteredRandomBox_[dir];
+        const List<scalar>& filter1 = filterCoeffs_[dir];
+        const List<scalar>& filter2 = filterCoeffs_[3 + dir];
+        const List<scalar>& filter3 = filterCoeffs_[6 + dir];
+        
+        const label sz1 = lenRandomBox_[dir];
+        const label sz2 = lenRandomBox_[3 + dir];
+        const label sz3 = lenRandomBox_[6 + dir];
+        const label szfilter1 = filterCoeffs_[dir].size();
+        const label szfilter2 = filterCoeffs_[3 + dir].size();
+        const label szfilter3 = filterCoeffs_[6 + dir].size();
+        const label sz23 = randomBoxFactors2D_[dir];
+        const label sz123 = randomBoxFactors3D_[dir];
+        const label validSlice2 = sz2 - (szfilter2 - label(1));
+        const label validSlice3 = sz3 - (szfilter3 - label(1));
+
+        // Convolution summation - Along 1st direction
+        scalarField tmp(sz123);
+        label filterCentre = label(szfilter2/label(2));
+        label endIndex = sz2 - filterCentre;
+        label i0 = 0;
+        label i1 = 0;
+        label i2 = 0;
+
+        for (label i = 0; i < sz1; ++i)
+        {
+            for (label j = 0; j < sz3; ++j)
+            {
+                i1 += filterCentre;
+
+                for (label k = filterCentre; k < endIndex; ++k)
+                {
+                    tmp[i1] = 0.0;
+                    label q = 0;
+
+                    for (label p = szfilter2 - 1; p >= 0; --p, ++q)
+                    {
+                        tmp[i1] += in[i0 + q]*filter2[p];
+                    }
+                    ++i0;
+                    ++i1;
+
+                }
+                i0 += (filterCentre + filterCentre);
+                i1 += filterCentre;
+            }
+        }
+
+        // Convolution summation - Along 2nd direction
+        scalarField tmp2(tmp);
+        filterCentre = label(szfilter3/label(2));
+        endIndex = sz3 - filterCentre;
+        i1 = 0;
+        i2 = 0;
+
+        for (label i = 0; i < sz1; ++i)
+        {
+           label sl = i*sz23;
+
+            for (label j = 0; j < sz2; ++j)
+            {
+                i1 = j + sl;
+                i2 = i1;
+
+                for (label k = 0; k < endIndex - filterCentre; ++k)
+                {
+                    tmp[i1] = 0.0;
+                    label q = 0;
+
+                    for (label p = szfilter3 - 1; p >= 0; --p, ++q)
+                    {
+                        tmp[i1] += tmp2[i2 + q*sz2]*filter3[p];
+                    }
+                    i1 += sz2;
+                    i2 += sz2;
+
+                }
+                i1 += (sz2 + filterCentre);
+                i2 += (sz2 + filterCentre);
+            }
+        }
+
+        // Convolution summation - Along 3rd direction
+        filterCentre = label(szfilter1/label(2));
+        endIndex = sz1 - filterCentre;
+        i1 = (szfilter2 - label(1))/label(2);
+        i2 = (szfilter2 - label(1))/label(2);
+        label i3 = 0;
+
+        for (label i = 0; i < validSlice3; ++i)
+        {
+            for (label j = 0; j < validSlice2; ++j)
+            {
+                scalar sum = 0.0;
+                i1 = i2 + j;
+
+                for (label k = szfilter1 - 1; k >= 0; --k)
+                {
+                    sum += tmp[i1]*filter1[k];
+                    i1 += sz23;
+                }
+                out[i3] = sum;
+                ++i3;
+
+            }
+            i2 += sz2;
+        }
+    }
+}
+
+
+void Foam::turbulentDigitalFilterInletFvPatchVectorField::computeDFM
+(
+    vectorField& U
+)
+{
+    #ifdef FULLDEBUG
+    Info<< "Starts: computeDFM()" << nl;
+    #endif
+
+    if (Pstream::master())
+    {
+        embedTwoPointCorrs();
+        rndShiftRefill();
+    }
+
+    Pstream::scatter(filteredRandomBox_);
+
+    mapFilteredRandomBox(U);
+
+    embedOnePointCorrs(U);
+
+    embedMeanVelocity(U);
+
+    if (isCorrectedFlowRate_)
+    {
+        correctFlowRate(U);
+    }
+
+    #ifdef FULLDEBUG
+    Info<< "Ends: computeDFM()" << nl;
+    #endif
+}
+
+
+void Foam::turbulentDigitalFilterInletFvPatchVectorField::computeReducedDFM
+(
+    vectorField& U
+)
+{
+    #ifdef FULLDEBUG
+    Info<< "Starts: computeReducedDFM()" << nl;
+    #endif
+
+    if (Pstream::master())
+    {
+        embedTwoPointCorrs();
+        rndShiftRefill();
+    }
+
+    Pstream::scatter(filteredRandomBox_);
+
+    mapFilteredRandomBox(U);
+
+    computeFSM(U);
+
+    embedOnePointCorrs(U);
+
+    embedMeanVelocity(U);
+
+    if (isCorrectedFlowRate_)
+    {
+        correctFlowRate(U);
+    }
+
+    #ifdef FULLDEBUG
+    Info<< "Ends: computeReducedDFM()" << nl;
+    #endif
+}
+
+
+Foam::List<Foam::scalar> Foam::turbulentDigitalFilterInletFvPatchVectorField::
+computeConstList1FSM() const
+{
+    List<scalar> constList1FSM(pTraits<vector>::nComponents);
+
+    forAll(L_.x(), i)
+    {
+        constList1FSM[i] = exp(const1FSM_/L_[i]);
+    }
+
+    return constList1FSM;
+}
+
+
+Foam::List<Foam::scalar> Foam::turbulentDigitalFilterInletFvPatchVectorField::
+computeConstList2FSM() const
+{
+    List<scalar> constList2FSM(pTraits<vector>::nComponents);
+
+    forAll(L_.x(), i)
+    {
+        constList2FSM[i] = Foam::sqrt(1.0 - exp(const2FSM_/L_[i]));
+    }
+
+    return constList2FSM;
+}
+
+
+void Foam::turbulentDigitalFilterInletFvPatchVectorField::computeFSM
+(
+    vectorField& U
+)
+{
+    for (direction dir = 0; dir < pTraits<vector>::nComponents; ++dir)
+    {
+        U.component(dir) = U0_.component(dir)*constList1FSM_[dir]
+                          +U.component(dir)*constList2FSM_[dir];
+    }
+
+    U0_ = U;
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::turbulentDigitalFilterInletFvPatchVectorField::
+turbulentDigitalFilterInletFvPatchVectorField
+(
+    const fvPatch& p,
+    const DimensionedField<vector, volMesh>& iF
+)
+:
+    fixedValueFvPatchField<vector>(p, iF),
+    mapperPtr_(nullptr),
+    variant_(variantType::DIGITAL_FILTER),
+    isGaussian_(true),
+    isFixedSeed_(true),
+    isContinuous_(false),
+    isCorrectedFlowRate_(true),
+    interpolateR_(false),
+    interpolateUMean_(false),
+    isInsideMesh_(false),
+    isTaylorHypot_(true),
+    mapMethod_("nearestCell"),
+    curTimeIndex_(-1),
+    tiny_(1e-8),
+    patchNormalSpeed_(Zero),
+    modelConst_(-0.5*constant::mathematical::pi),
+    perturb_(1e-5),
+    initialFlowRate_(pTraits<scalar>::one),
+    rndGen_(1234),
+    planeDivisions_(Zero, Zero),
+    invDelta_(Zero),
+    indexPairs_(Zero),
+    R_(Zero),
+    LundWuSquires_(Zero),
+    UMean_(Zero),
+    Lbak_(Zero),
+    L_(Zero),
+    const1FSM_(Zero),
+    const2FSM_(Zero),
+    constList1FSM_(Zero),
+    constList2FSM_(Zero),
+    lenRandomBox_(Zero),
+    randomBoxFactors2D_(Zero),
+    randomBoxFactors3D_(Zero),
+    iNextToLastPlane_(Zero),
+    randomBox_(Zero),
+    filterCoeffs_(Zero),
+    filteredRandomBox_(Zero),
+    U0_(Zero),
+    computeVariant(nullptr)
+{}
+
+
+Foam::turbulentDigitalFilterInletFvPatchVectorField::
+turbulentDigitalFilterInletFvPatchVectorField
+(
+    const turbulentDigitalFilterInletFvPatchVectorField& ptf,
+    const fvPatch& p,
+    const DimensionedField<vector, volMesh>& iF,
+    const fvPatchFieldMapper& mapper
+)
+:
+    fixedValueFvPatchField<vector>(ptf, p, iF, mapper),
+    mapperPtr_(nullptr),
+    variant_(ptf.variant_),
+    isGaussian_(ptf.isGaussian_),
+    isFixedSeed_(ptf.isFixedSeed_),
+    isContinuous_(ptf.isContinuous_),
+    isCorrectedFlowRate_(ptf.isCorrectedFlowRate_),
+    interpolateR_(ptf.interpolateR_),
+    interpolateUMean_(ptf.interpolateUMean_),
+    isInsideMesh_(ptf.isInsideMesh_),
+    isTaylorHypot_(ptf.isTaylorHypot_),
+    mapMethod_(ptf.mapMethod_),
+    curTimeIndex_(-1),
+    tiny_(ptf.tiny_),
+    patchNormalSpeed_(ptf.patchNormalSpeed_),
+    modelConst_(ptf.modelConst_),
+    perturb_(ptf.perturb_),
+    initialFlowRate_(ptf.initialFlowRate_),
+    rndGen_(ptf.rndGen_),
+    planeDivisions_(ptf.planeDivisions_),
+    invDelta_(ptf.invDelta_),
+    indexPairs_(ptf.indexPairs_),
+    R_(ptf.R_),
+    LundWuSquires_(ptf.LundWuSquires_),
+    UMean_(ptf.UMean_),
+    Lbak_(ptf.Lbak_),
+    L_(ptf.L_),
+    const1FSM_(ptf.const1FSM_),
+    const2FSM_(ptf.const2FSM_),
+    constList1FSM_(ptf.constList1FSM_),
+    constList2FSM_(ptf.constList2FSM_),
+    lenRandomBox_(ptf.lenRandomBox_),
+    randomBoxFactors2D_(ptf.randomBoxFactors2D_),
+    randomBoxFactors3D_(ptf.randomBoxFactors3D_),
+    iNextToLastPlane_(ptf.iNextToLastPlane_),
+    randomBox_(ptf.randomBox_),
+    filterCoeffs_(ptf.filterCoeffs_),
+    filteredRandomBox_(ptf.filteredRandomBox_),
+    U0_(ptf.U0_),
+    computeVariant(ptf.computeVariant)
+{}
+
+
+Foam::turbulentDigitalFilterInletFvPatchVectorField::
+turbulentDigitalFilterInletFvPatchVectorField
+(
+    const fvPatch& p,
+    const DimensionedField<vector, volMesh>& iF,
+    const dictionary& dict
+)
+:
+    fixedValueFvPatchField<vector>(p, iF, dict),
+    mapperPtr_(nullptr),
+    variant_
+    (
+        variantNames.getOrDefault
+        (
+            "variant",
+            dict,
+            variantType::DIGITAL_FILTER
+        )
+    ),
+    isGaussian_
+    (
+        (variant_ == variantType::FORWARD_STEPWISE)
+      ? false
+      : dict.getOrDefault("isGaussian", true)
+    ),
+    isFixedSeed_(dict.getOrDefault("isFixedSeed", true)),
+    isContinuous_(dict.getOrDefault("isContinuous", false)),
+    isCorrectedFlowRate_(dict.getOrDefault("isCorrectedFlowRate", true)),
+    interpolateR_(dict.getOrDefault("interpolateR", false)),
+    interpolateUMean_(dict.getOrDefault("interpolateUMean", false)),
+    isInsideMesh_(dict.getOrDefault("isInsideMesh", false)),
+    isTaylorHypot_(dict.getOrDefault("isTaylorHypot", true)),
+    mapMethod_(dict.getOrDefault<word>("mapMethod", "nearestCell")),
+    curTimeIndex_(-1),
+    tiny_(dict.getOrDefault<scalar>("threshold", 1e-8)),
+    patchNormalSpeed_
+    (
+        dict.getCheck<scalar>
+        (
+            "patchNormalSpeed",
+            scalarMinMax::ge(tiny_)
+        )
+    ),
+    modelConst_
+    (
+        dict.getOrDefault<scalar>
+        (
+            "modelConst",
+           -0.5*constant::mathematical::pi
+        )
+    ),
+    perturb_(dict.getOrDefault<scalar>("perturb", 1e-5)),
+    initialFlowRate_(pTraits<scalar>::one),
+    rndGen_
+    (
+        isFixedSeed_
+      ? 1234
+      : time(0)
+    ),
+    planeDivisions_
+    (
+        dict.getCheck<Tuple2<label, label>>
+        (
+            "planeDivisions",
+            [&](const Tuple2<label, label>& len)
+            {
+                return tiny_ < min(len.first(), len.second()) ? true : false;
+            }
+        )
+    ),
+    invDelta_(),
+    indexPairs_(patchIndexPairs()),
+    R_(interpolateOrRead<symmTensor>("R", dict, interpolateR_)),
+    LundWuSquires_(computeLundWuSquires()),
+    UMean_(interpolateOrRead<vector>("UMean", dict, interpolateUMean_)),
+    Lbak_
+    (
+        dict.getCheck<tensor>
+        (
+            "L",
+            [&](const tensor& l){return tiny_ < cmptMin(l) ? true : false;}
+        )
+    ),
+    L_(convertScalesToGridUnits(Lbak_)),
+    const1FSM_
+    (
+        (variant_ == variantType::FORWARD_STEPWISE)
+      ? dict.getOrDefault<scalar>
+        (
+            "const1FSM",
+           -0.25*constant::mathematical::pi
+        )
+      : Zero
+    ),
+    const2FSM_
+    (
+        (variant_ == variantType::FORWARD_STEPWISE)
+      ? dict.getOrDefault<scalar>
+        (
+            "const2FSM",
+           -0.5*constant::mathematical::pi
+        )
+      : Zero
+    ),
+    constList1FSM_
+    (
+        (variant_ == variantType::FORWARD_STEPWISE)
+      ? computeConstList1FSM()
+      : List<scalar>()
+    ),
+    constList2FSM_
+    (
+        (variant_ == variantType::FORWARD_STEPWISE)
+      ? computeConstList2FSM()
+      : List<scalar>()
+    ),
+    lenRandomBox_(initLenRandomBox()),
+    randomBoxFactors2D_(initBoxFactors2D()),
+    randomBoxFactors3D_(initBoxFactors3D()),
+    iNextToLastPlane_(initBoxPlaneFactors()),
+    randomBox_
+    (
+        (isContinuous_ && Pstream::master())
+      ? dict.getOrDefault<List<List<scalar>>>
+        (
+            "randomBox",
+            fillRandomBox() // First time-step
+        )
+      :
+        fillRandomBox()
+    ),
+    filterCoeffs_(computeFilterCoeffs()),
+    filteredRandomBox_
+    (
+        pTraits<vector>::nComponents,
+        List<scalar>(planeDivisions_.first()*planeDivisions_.second(), Zero)
+    ),
+    U0_
+    (
+        (variant_ == variantType::FORWARD_STEPWISE)
+      ? generateRandomSet<vectorField, vector>(patch().size())
+      : vectorField()
+    ),
+    computeVariant
+    (
+        (variant_ == variantType::FORWARD_STEPWISE)
+      ? &turbulentDigitalFilterInletFvPatchVectorField::computeReducedDFM
+      : &turbulentDigitalFilterInletFvPatchVectorField::computeDFM
+    )
+{
+    if (isCorrectedFlowRate_)
+    {
+        initialFlowRate_ = computeInitialFlowRate();
+    }
+
+    // Check if varying or fixed time-step computation
+    if (db().time().isAdjustTimeStep())
+    {
+        WarningInFunction
+            << "Varying time-step computations are not fully supported"
+            << " for the moment."<< nl << nl;
+    }
+
+    #ifdef FULLDEBUG
+    Info<< "Ends: Resource acquisition/initialisation for the"
+        << " synthetic turbulence boundary condition." << nl;
+    #endif
+}
+
+
+Foam::turbulentDigitalFilterInletFvPatchVectorField::
+turbulentDigitalFilterInletFvPatchVectorField
+(
+    const turbulentDigitalFilterInletFvPatchVectorField& ptf
+)
+:
+    fixedValueFvPatchField<vector>(ptf),
+    mapperPtr_(nullptr),
+    variant_(ptf.variant_),
+    isGaussian_(ptf.isGaussian_),
+    isFixedSeed_(ptf.isFixedSeed_),
+    isContinuous_(ptf.isContinuous_),
+    isCorrectedFlowRate_(ptf.isCorrectedFlowRate_),
+    interpolateR_(ptf.interpolateR_),
+    interpolateUMean_(ptf.interpolateUMean_),
+    isInsideMesh_(ptf.isInsideMesh_),
+    isTaylorHypot_(ptf.isTaylorHypot_),
+    mapMethod_(ptf.mapMethod_),
+    curTimeIndex_(ptf.curTimeIndex_),
+    tiny_(ptf.tiny_),
+    patchNormalSpeed_(ptf.patchNormalSpeed_),
+    modelConst_(ptf.modelConst_),
+    perturb_(ptf.perturb_),
+    initialFlowRate_(ptf.initialFlowRate_),
+    rndGen_(ptf.rndGen_),
+    planeDivisions_(ptf.planeDivisions_),
+    invDelta_(ptf.invDelta_),
+    indexPairs_(ptf.indexPairs_),
+    R_(ptf.R_),
+    LundWuSquires_(ptf.LundWuSquires_),
+    UMean_(ptf.UMean_),
+    Lbak_(ptf.Lbak_),
+    L_(ptf.L_),
+    const1FSM_(ptf.const1FSM_),
+    const2FSM_(ptf.const2FSM_),
+    constList1FSM_(ptf.constList1FSM_),
+    constList2FSM_(ptf.constList2FSM_),
+    lenRandomBox_(ptf.lenRandomBox_),
+    randomBoxFactors2D_(ptf.randomBoxFactors2D_),
+    randomBoxFactors3D_(ptf.randomBoxFactors3D_),
+    iNextToLastPlane_(ptf.iNextToLastPlane_),
+    randomBox_(ptf.randomBox_),
+    filterCoeffs_(ptf.filterCoeffs_),
+    filteredRandomBox_(ptf.filteredRandomBox_),
+    U0_(ptf.U0_),
+    computeVariant(ptf.computeVariant)
+{}
+
+
+Foam::turbulentDigitalFilterInletFvPatchVectorField::
+turbulentDigitalFilterInletFvPatchVectorField
+(
+    const turbulentDigitalFilterInletFvPatchVectorField& ptf,
+    const DimensionedField<vector, volMesh>& iF
+)
+:
+    fixedValueFvPatchField<vector>(ptf, iF),
+    mapperPtr_(nullptr),
+    variant_(ptf.variant_),
+    isGaussian_(ptf.isGaussian_),
+    isFixedSeed_(ptf.isFixedSeed_),
+    isContinuous_(ptf.isContinuous_),
+    isCorrectedFlowRate_(ptf.isCorrectedFlowRate_),
+    interpolateR_(ptf.interpolateR_),
+    interpolateUMean_(ptf.interpolateUMean_),
+    isInsideMesh_(ptf.isInsideMesh_),
+    isTaylorHypot_(ptf.isTaylorHypot_),
+    mapMethod_(ptf.mapMethod_),
+    curTimeIndex_(-1),
+    tiny_(ptf.tiny_),
+    patchNormalSpeed_(ptf.patchNormalSpeed_),
+    modelConst_(ptf.modelConst_),
+    perturb_(ptf.perturb_),
+    initialFlowRate_(ptf.initialFlowRate_),
+    rndGen_(ptf.rndGen_),
+    planeDivisions_(ptf.planeDivisions_),
+    invDelta_(ptf.invDelta_),
+    indexPairs_(ptf.indexPairs_),
+    R_(ptf.R_),
+    LundWuSquires_(ptf.LundWuSquires_),
+    UMean_(ptf.UMean_),
+    Lbak_(ptf.Lbak_),
+    L_(ptf.L_),
+    const1FSM_(ptf.const1FSM_),
+    const2FSM_(ptf.const2FSM_),
+    constList1FSM_(ptf.constList1FSM_),
+    constList2FSM_(ptf.constList2FSM_),
+    lenRandomBox_(ptf.lenRandomBox_),
+    randomBoxFactors2D_(ptf.randomBoxFactors2D_),
+    randomBoxFactors3D_(ptf.randomBoxFactors3D_),
+    iNextToLastPlane_(ptf.iNextToLastPlane_),
+    randomBox_(ptf.randomBox_),
+    filterCoeffs_(ptf.filterCoeffs_),
+    filteredRandomBox_(ptf.filteredRandomBox_),
+    U0_(ptf.U0_),
+    computeVariant(ptf.computeVariant)
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+void Foam::turbulentDigitalFilterInletFvPatchVectorField::updateCoeffs()
+{
+    if (updated())
+    {
+        return;
+    }
+
+    if (curTimeIndex_ != db().time().timeIndex())
+    {
+        vectorField& U = *this;
+
+        computeVariant(this, U);
+
+        curTimeIndex_ = db().time().timeIndex();
+    }
+
+    fixedValueFvPatchVectorField::updateCoeffs();
+}
+
+
+void Foam::turbulentDigitalFilterInletFvPatchVectorField::write
+(
+    Ostream& os
+) const
+{
+    fixedValueFvPatchField<vector>::write(os),
+    os.writeEntry("variant", variantNames[variant_]);
+    os.writeEntry("planeDivisions", planeDivisions_);
+    os.writeEntry("L", Lbak_);
+
+    if (!interpolateR_)
+    {
+        R_.writeEntry("R", os);
+    }
+
+    if (!interpolateUMean_)
+    {
+        UMean_.writeEntry("UMean", os);
+    }
+
+    os.writeEntryIfDifferent<bool>("isGaussian", true, isGaussian_);
+    os.writeEntryIfDifferent<bool>("isFixedSeed", true, isFixedSeed_);
+    os.writeEntryIfDifferent<bool>("isContinuous", false, isContinuous_);
+    os.writeEntryIfDifferent<bool>
+        ("isCorrectedFlowRate", true, isCorrectedFlowRate_);
+    os.writeEntryIfDifferent<bool>("isInsideMesh", false, isInsideMesh_);
+    os.writeEntryIfDifferent<bool>("isTaylorHypot", true, isTaylorHypot_);
+
+    if (!mapMethod_.empty())
+    {
+        os.writeEntryIfDifferent<word>
+        (
+            "mapMethod",
+            "nearestCell",
+            mapMethod_
+        );
+    }
+
+    os.writeEntry("threshold", tiny_);
+    os.writeEntry("patchNormalSpeed", patchNormalSpeed_);
+    os.writeEntry("modelConst", modelConst_);
+    os.writeEntry("perturb", perturb_);
+
+    if (variant_ == variantType::FORWARD_STEPWISE)
+    {
+        os.writeEntry("const1FSM", const1FSM_);
+        os.writeEntry("const2FSM", const2FSM_);
+    }
+
+    if (isContinuous_ && Pstream::master())
+    {
+        os.writeEntry("randomBox", randomBox_);
+    }
+}
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+   makePatchTypeField
+   (
+       fvPatchVectorField,
+       turbulentDigitalFilterInletFvPatchVectorField
+   );
+}
+
+
+// ************************************************************************* //
diff --git a/src/finiteVolume/fields/fvPatchFields/derived/turbulentDigitalFilterInlet/turbulentDigitalFilterInletFvPatchVectorField.H b/src/finiteVolume/fields/fvPatchFields/derived/turbulentDigitalFilterInlet/turbulentDigitalFilterInletFvPatchVectorField.H
new file mode 100644
index 0000000000..2ec5ef8c7e
--- /dev/null
+++ b/src/finiteVolume/fields/fvPatchFields/derived/turbulentDigitalFilterInlet/turbulentDigitalFilterInletFvPatchVectorField.H
@@ -0,0 +1,569 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2019 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::turbulentDigitalFilterFvPatchVectorField
+
+Group
+    grpInletBoundaryConditions
+
+Description
+    Velocity boundary condition generating synthetic turbulence-alike
+    time-series for LES and DES turbulent flow computations.
+
+    To this end, two synthetic turbulence generators can be chosen:
+
+    - Digital-filter method-based generator (DFM)
+
+    \verbatim
+    Klein, M., Sadiki, A., and Janicka, J.
+        A digital filter based generation of inflow data for spatially 
+        developing direct numerical or large eddy simulations,
+        Journal of Computational Physics (2003) 186(2):652-665.
+        doi:10.1016/S0021-9991(03)00090-1
+    \endverbatim
+
+    - Forward-stepwise method-based generator (FSM)
+
+    \verbatim
+    Xie, Z.-T., and Castro, I.
+        Efficient generation of inflow conditions for large eddy simulation of 
+        street-scale flows, Flow, Turbulence and Combustion (2008) 81(3):449-470
+        doi:10.1007/s10494-008-9151-5
+    \endverbatim
+
+    In DFM or FSM, a random number set (mostly white noise), and a group
+    of target statistics (mostly mean flow, Reynolds stress tensor profiles and
+    length-scale sets) are fused into a new number set (stochastic time-series,
+    yet consisting of the statistics) by a chain of mathematical operations
+    whose characteristics are designated by the target statistics, so that the
+    realised statistics of the new sets could match the target.
+
+    Random number sets ---->-|
+                             |
+                         DFM or FSM ---> New stochastic time-series consisting
+                             |           turbulence statistics
+    Turbulence statistics ->-|
+
+    The main difference between DFM and FSM is that the latter replaces the
+    streamwise convolution summation in DFM by a simpler and a quantitatively
+    justified equivalent procedure in order to reduce computational costs.
+    Accordingly, the latter potentially brings resource advantages for
+    computations involving relatively large length-scale sets and small
+    time-steps.
+
+    Synthetic turbulence is produced on a virtual rectangular structured-mesh
+    turbulence plane, which is parallel to the actual patch, and is mapped onto
+    the chosen patch by the selected mapping method.
+
+Usage
+    \table
+     Property    | Description                     | Required    | Default value
+     planeDivisions   | Number of nodes on turbulence plane (e2, e3) [-] | yes |
+     L        | Integral length-scale set (9-comp):{e1,e2,e3}{u,v,w} [m] | yes |
+     R         | Reynolds stress tensor set (xx xy xz yy yz zz) [m2/s2]  | yes |
+     patchNormalSpeed    | Characteristic mean flow speed  [m/s] | yes |
+     isGaussian  | Autocorrelation function form                 | no | true
+     isFixedSeed | Flag to identify random-number seed is fixed  | no | true
+     isContinuous   | Flag for random-number restart behaviour   | no | false
+     isCorrectedFlowRate | Flag for mass flow rate correction    | no | true
+     interpolateR     | Placeholder flag: interpolate R field    | no | false
+     interpolateUMean | Placeholder flag: interpolate UMean field   | no | false
+     isInsideMesh | Placeholder flag: TP is inside mesh or on patch | no | false
+     isTaylorHypot | Placeholder flag: Taylor's hypothesis is on | no | true
+     mapMethod   | Method to map reference values             | no | nearestCell
+     threshold   | Threshold to avoid unintentional 'tiny' input | no | 1e-8
+     modelConst  | Model constant (Klein et al., Eq. 14)         | no | -0.5*PI
+     perturb     | Point perturbation for interpolation          | no | 1e-5
+     const1FSM | A model coefficient in FSM (Xie-Castro, Eq. 14) | no | -0.25*PI
+     const2FSM | A model coefficient in FSM (Xie-Castro, Eq. 14) | no | -0.5*PI
+    \endtable
+
+    Minimal example of the boundary condition specification with commented
+    options:
+    \verbatim
+    <patchName>
+    {
+        type                turbulentDigitalFilter;
+        variant             digitalFilter;          // reducedDigitalFilter;
+        planeDivisions      (<planeDivisionsHeight> <planeDivisionsWidth>);
+        L               (<Lxu> <Lxv> <Lxw> <Lyu> <Lyv> <Lyw> <Lzu> <Lzv> <Lzw>);
+        R                   (<Rxx> <Rxy> <Rxz> <Ryy> <Ryz> <Rzz>);
+        patchNormalSpeed    <characteristic flow speed>;
+        value               uniform (0 0 0);        // mandatory placeholder
+
+        // Optional entries with default input
+        isGaussian          true;           // false    // always false for FSM
+        isFixedSeed         true;           // false
+        isContinuous        false;          // true
+        isCorrectedFlowRate true;           // false
+        interpolateR        false;          // placeholder
+        interpolateUMean    false;          // placeholder
+        isInsideMesh        false;          // placeholder
+        isTaylorHypot       true;           // placeholder
+        mapMethod           nearestCell;    // planarInterpolation
+        threshold           1e-8;
+        modelConst          -1.5707;        //-0.5*PI;
+        perturb             1e-5;
+
+        // Optional entries for only FSM with default input
+        const1FSM           -0.7854         //-0.25*PI;
+        const2FSM           -1.5707;        //-0.5*PI;
+    }
+    \endverbatim 
+
+    Among the dictionary entries, two entries can be input as patch profiles:
+
+    \verbatim
+     - Reynolds stress tensor, R
+     - Mean velocity, UMean
+    \endverbatim
+
+    Profile data and corresponding coordinates are then input in the following
+    directories:
+
+    \verbatim
+     - $FOAM_CASE/constant/boundaryData/\<patchName\>/points
+     - $FOAM_CASE/constant/boundaryData/\<patchName\>/0/\{R|UMean\}
+    \endverbatim
+
+    The profile data and corresponding coordinates take the same form used by
+    the \c timeVaryingMappedFixedValue and \c turbulentDFSEMInlet boundary
+    conditions, consisting of a \c points file containing a list of 3-D
+    coordinates, and profile data files providing a value per coordinate.
+
+    Reynolds stress tensor and mean velocity input are in the global coordinate
+    system whereas integral length scale set input is in the local patch
+    coordinate system.
+
+    It is assumed that the patch normal direction is \c e1, and the remaining
+    patch plane directions are \c e2 and \c e3 following the right-handed
+    coordinate system, which should be taken into consideration when the
+    integral length scale set is input. The first three integral scale entries,
+    i.e. L_e1u, Le1v, Le2w, should always correspond to the length scales
+    that are in association with the convective mean flow direction.
+
+Note
+    - \c mapMethod \c planarInterpolation option requires point coordinates
+    which can form a plane, thus input point coordinates varying only in a
+    single direction will trigger error.
+    - \c adjustTimeStep = true option is not fully supported at the moment.
+
+SeeAlso
+    turbulentDFSEMInletFvPatchVectorField.C
+
+SourceFiles
+    turbulentDigitalFilterInletFvPatchVectorField.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef turbulentDigitalFilterInletFvPatchVectorField_H
+#define turbulentDigitalFilterInletFvPatchVectorField_H
+
+#include "fixedValueFvPatchFields.H"
+#include "Random.H"
+#include <functional>
+#include "fieldTypes.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+class pointToPointPlanarInterpolation;
+
+/*---------------------------------------------------------------------------*\
+        Class turbulentDigitalFilterInletFvPatchVectorField Declaration
+\*---------------------------------------------------------------------------*/
+
+class turbulentDigitalFilterInletFvPatchVectorField
+:
+    public fixedValueFvPatchVectorField
+{
+
+    // Private Enumerations
+
+        //- Options for the synthetic turbulence generator variant
+        enum variantType : uint8_t
+        {
+            DIGITAL_FILTER = 1,    //!< Digital-filter method (Klein et al.)
+            FORWARD_STEPWISE = 2,  //!< Forward-stepwise method (Xie-Castro)
+        };
+
+        //- Names for variant types
+        static const Enum<variantType> variantNames;
+
+
+    // Private Data
+
+        //- 2D interpolation (for 'planarInterpolation' mapMethod)
+        mutable autoPtr<pointToPointPlanarInterpolation> mapperPtr_;
+
+        //- Selected option for the synthetic turbulence generator variant
+        const enum variantType variant_;
+
+        //- Flag: correlation function form is Gaussian or Exponential
+        //  for variantType::DIGITAL_FILTER, default=Gaussian
+        //  for variantType::FORWARD_STEPWISE, default=Exponential (only option)
+        const bool isGaussian_;
+
+        //- Flag: random-number generator seed is fixed (default=true) or
+        //- generated pseudo-randomly based on clock-time per simulation
+        const bool isFixedSeed_;
+
+        //- Flag: write non-manipulated random-number sets at output time, and
+        //- to read them on restart. Otherwise, generate new random-number sets
+        //- on restart. (default=false)
+        //  true: deterministic & statistically consistent, more expensive
+        //  false: deterministic discontinuity & statistically consistent, cheaper
+        const bool isContinuous_;
+
+        //- Flag: mass flow rate is corrected on turbulence plane (default=true)
+        const bool isCorrectedFlowRate_;
+
+        //- Flag: interpolate R field (default=false)
+        bool interpolateR_;
+
+        //- Flag: interpolate UMean field (default=false)
+        bool interpolateUMean_;
+
+        //- Flag: turbulence plane is inside mesh or on a patch (default=false)
+        //  Currently, true option is not available.
+        const bool isInsideMesh_;
+
+        //- Flag: convert streamwise (x) length scales to time scales by
+        //- Taylor's 'frozen turbulence' hypothesis (default=true)
+        //  Currently, false option is not available.
+        const bool isTaylorHypot_;
+
+        //- Method for interpolation between a patch and turbulence plane
+        //- (default=nearestCell)
+        //  Options:
+        //   - nearestCell: one-to-one direct map, no interpolation
+        //   - planarInterpolation: bilinear interpolation
+        const word mapMethod_;
+
+        //- Current time index
+        label curTimeIndex_;
+
+        //- Threshold to avoid unintentional 'tiny' input scalars (default=1e-8)
+        const scalar tiny_;
+
+        //- Characteristic (e.g. bulk) mean speed of flow in the patch normal
+        //- direction [m/s]
+        const scalar patchNormalSpeed_;
+
+        //- Model constant shaping autocorr function [-] (default='-0.5*pi')
+        //  (Klein et al., 2003, Eq. 14)
+        const scalar modelConst_;
+
+        //- Fraction of perturbation (fraction of bounding box) to add (for
+        //- 'planarInterpolation' mapMethod)
+        const scalar perturb_;
+
+        //- Initial (first time-step) mass/vol flow rate [m^3/s]
+        scalar initialFlowRate_;
+
+        //- Random number generator
+        Random rndGen_;
+
+        //- Number of nodes on turbulence plane (e2 e3) [-]
+        const Tuple2<label, label> planeDivisions_;
+
+        //- Turbulence plane mesh size (reversed) (e2 e3) [1/m]
+        Vector2D<scalar> invDelta_;
+
+        //- Nearest cell mapping: Index pairs between patch and turbulence plane
+        const List<Pair<label>> indexPairs_;
+
+        //- Reynolds stress tensor profile (xx xy xz yy yz zz) in global
+        //- coordinates [m^2/s^2]
+        symmTensorField R_;
+
+        //- Lund-Wu-Squires transformation (Cholesky decomp.) [m/s]
+        //- Mapped onto actual mesh patch rather than turbulence plane
+        //  (Klein et al., 2003, Eq. 5)
+        symmTensorField LundWuSquires_;
+
+        //- Mean inlet velocity profile in global coordinates [m/s]
+        vectorField UMean_;
+
+        //- Integral length-scale set per turbulence plane section in local 
+        //- coordinates (e1u, e1v, e1w, e2u, e2v, e2w, e3u, e3v, e3w) [m]
+        //  First three entries should always correspond to the length scales
+        //  in association with the convective mean flow direction
+        //  Backup of L_ for restart purposes
+        const tensor Lbak_;
+
+        //- Integral length-scale set in mesh units [node]
+        const tensor L_;
+
+        //- One of the two model coefficients in FSM
+        //  (Xie-Castro, 2008, the argument of the first exp func in Eq. 14)
+        const scalar const1FSM_;
+
+        //- One of the two model coefficients in FSM
+        //  (Xie-Castro, 2008, the argument of the second exp func in Eq. 14)
+        const scalar const2FSM_;
+
+        //- One of the two exponential functions in FSM
+        //  (Xie-Castro, 2008, the first exponential function in Eq. 14)
+        const List<scalar> constList1FSM_;
+
+        //- One of the two exponential functions in FSM
+        //  (Xie-Castro, 2008, the first exponential function in Eq. 14)
+        const List<scalar> constList2FSM_;
+
+        //- Number of nodes in random-number box [node]
+        //- Random-number sets within box are filtered with filterCoeffs_
+        //- (e1u, e1v, e1w, e2u, e2v, e2w, e3u, e3v, e3w)
+        const List<label> lenRandomBox_;
+
+        //- Convenience factors for 2-D random-number box [-]
+        const List<label> randomBoxFactors2D_;
+
+        //- Convenience factors for 3-D randomNum box [-]
+        const List<label> randomBoxFactors3D_;
+
+        //- Index to the first elem of last plane of random-number box [-]
+        const List<label> iNextToLastPlane_;
+
+        //- Random-number sets distributed over a 3-D box (u, v, w)
+        List<List<scalar>> randomBox_;
+
+        //- Filter coefficients corresponding to L_ [-]
+        const List<List<scalar>> filterCoeffs_;
+
+        //- Filter-applied random-number sets [m/s] (effectively turb plane)
+        List<List<scalar>> filteredRandomBox_;
+
+        //- Filter-applied previous-time-step velocity field [m/s] used in FSM
+        vectorField U0_;
+
+        //- Run-time function selector between the original and reduced methods
+        const std::function
+        <
+            void(turbulentDigitalFilterInletFvPatchVectorField*, vectorField&)
+        > computeVariant;
+
+
+    // Private Member Functions
+
+        //- Return a reference to the patch mapper object
+        const pointToPointPlanarInterpolation& patchMapper() const;
+
+        //- Helper function to interpolate values from the boundary data or
+        //- read from dictionary
+        template<class Type>
+        tmp<Field<Type>> interpolateOrRead
+        (
+            const word& fieldName,
+            const dictionary& dict,
+            bool& interpolateField
+        ) const;
+
+        //- Helper function to interpolate values from the boundary data
+        template<class Type>
+        tmp<Field<Type>> interpolateBoundaryData
+        (
+            const word& fieldName
+        ) const;
+
+        //- Generate random-number sets obeying the standard normal distribution
+        template<class Form, class Type>
+        Form generateRandomSet(const label len);
+
+        //- Compute nearest cell index-pairs between turbulence plane and patch
+        List<Pair<label>> patchIndexPairs();
+
+        //- Check R_ (mapped on actual mesh) for mathematical domain errors
+        void checkRTensorRealisable() const;
+
+        //- Compute Lund-Wu-Squires transformation
+        //  (Klein et al., 2003, Eq. 5)
+        symmTensorField computeLundWuSquires() const;
+
+        //- Compute patch-normal into the domain
+        vector computePatchNormal() const;
+
+        //- Compute initial (first time-step) mass/vol flow rate based on UMean_
+        scalar computeInitialFlowRate() const;
+
+        //- Convert streamwise integral length scales to integral time scales
+        //- via Taylor's frozen turbulence hypothesis
+        void convertToTimeScale(tensor& L) const;
+
+        //- Convert length scale phys. unit to turbulence plane mesh-size unit
+        //- (Klein et al., 2003, Eq. 13)
+        tensor convertScalesToGridUnits(const tensor& L) const;
+
+        //- Resource allocation functions for the convenience factors
+        List<label> initLenRandomBox() const;
+        List<label> initBoxFactors2D() const;
+        List<label> initBoxFactors3D() const;
+        List<label> initBoxPlaneFactors() const;
+
+        //- Compute various convenience factors for random-number box
+        List<List<scalar>> fillRandomBox();
+
+        //- Compute filter coeffs once per simulation
+        //  (Klein et al., 2003, Eq. 14)
+        List<List<scalar>> computeFilterCoeffs() const;
+
+        //- Discard current time-step random-box plane (closest to patch) by
+        //- shifting from the back to the front, and dd new plane to the back
+        void rndShiftRefill();
+
+        //- Map two-point correlated random-number sets on patch based on chosen
+        //- mapping method
+        void mapFilteredRandomBox(vectorField& U);
+
+        //- Map R_ on patch
+        void embedOnePointCorrs(vectorField& U) const;
+
+        //- Map UMean_ on patch
+        void embedMeanVelocity(vectorField& U) const;
+
+        //- Correct mass/vol flow rate in (only) streamwise direction
+        void correctFlowRate(vectorField& U) const;
+
+        //- 3-D 'valid'-type 'separable' convolution summation algorithm
+        //  'Inspired' from Song Ho Ahn's 2-D 'full'-type convolution algorithm
+        //  with his permission
+        void embedTwoPointCorrs();
+
+        //- Compute the DFM (Klein et al., 2003)
+        void computeDFM(vectorField& U);
+
+        //- Compute the reduced DFM (i.e. hybrid FSM-DFM) (Xie-Castro, 2008)
+        void computeReducedDFM(vectorField& U);
+
+        //- Compute constList1FSM_ once per simulation
+        List<scalar> computeConstList1FSM() const;
+
+        //- Compute constList2FSM_ once per simulation
+        List<scalar> computeConstList2FSM() const;
+
+        //- Compute the forward-stepwise method
+        //  (Xie-Castro, 2008, Eq. 14)
+        void computeFSM(vectorField& U);
+
+
+public:
+
+   //- Runtime type information
+   TypeName("turbulentDigitalFilterInlet");
+
+    // Constructors
+
+        //- Construct from patch and internal field
+        turbulentDigitalFilterInletFvPatchVectorField
+        (
+            const fvPatch&,
+            const DimensionedField<vector, volMesh>&
+        );
+
+        //- Construct from patch, internal field and dictionary
+        turbulentDigitalFilterInletFvPatchVectorField
+        (
+            const fvPatch&,
+            const DimensionedField<vector, volMesh>&,
+            const dictionary&
+        );
+
+        //- Construct by mapping given
+        //- turbulentDigitalFilterInletFvPatchVectorField onto a new patch
+        turbulentDigitalFilterInletFvPatchVectorField
+        (
+            const turbulentDigitalFilterInletFvPatchVectorField&,
+            const fvPatch&,
+            const DimensionedField<vector, volMesh>&,
+            const fvPatchFieldMapper&
+        );
+
+        //- Construct as copy
+        turbulentDigitalFilterInletFvPatchVectorField
+        (
+            const turbulentDigitalFilterInletFvPatchVectorField&
+        );
+
+        //- Construct and return a clone
+        virtual tmp<fvPatchVectorField> clone() const
+        {
+            return tmp<fvPatchVectorField>
+            (
+                new turbulentDigitalFilterInletFvPatchVectorField(*this)
+            );
+        }
+
+        //- Construct as copy setting internal field reference
+        turbulentDigitalFilterInletFvPatchVectorField
+        (
+            const turbulentDigitalFilterInletFvPatchVectorField&,
+            const DimensionedField<vector, volMesh>&
+        );
+
+        //- Construct and return a clone setting internal field reference
+        virtual tmp<fvPatchVectorField> clone
+        (
+            const DimensionedField<vector, volMesh>& iF
+        ) const
+        {
+            return tmp<fvPatchVectorField>
+            (
+                new turbulentDigitalFilterInletFvPatchVectorField(*this, iF)
+            );
+        }
+
+
+    //- Destructor
+    virtual ~turbulentDigitalFilterInletFvPatchVectorField() = default;
+
+
+    // Member Functions
+
+        // Evaluation functions
+
+            //- Update the coefficients associated with the patch field
+            virtual void updateCoeffs();
+
+
+       //- Write
+       virtual void write(Ostream&) const;
+};
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#ifdef NoRepository
+    #include "turbulentDigitalFilterInletFvPatchVectorFieldTemplates.C"
+#endif
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/finiteVolume/fields/fvPatchFields/derived/turbulentDigitalFilterInlet/turbulentDigitalFilterInletFvPatchVectorFieldTemplates.C b/src/finiteVolume/fields/fvPatchFields/derived/turbulentDigitalFilterInlet/turbulentDigitalFilterInletFvPatchVectorFieldTemplates.C
new file mode 100644
index 0000000000..f43b2a2d9f
--- /dev/null
+++ b/src/finiteVolume/fields/fvPatchFields/derived/turbulentDigitalFilterInlet/turbulentDigitalFilterInletFvPatchVectorFieldTemplates.C
@@ -0,0 +1,127 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2016 - 2019 OpenCFD Ltd.
+     \\/     M anipulation  | 
+-------------------------------------------------------------------------------
+                              Copyright (C) 2015 OpenFOAM Foundation
+-------------------------------------------------------------------------------
+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 "pointToPointPlanarInterpolation.H"
+#include "Time.H"
+#include "IFstream.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+template<class Type>
+Foam::tmp<Foam::Field<Type>>
+Foam::turbulentDigitalFilterInletFvPatchVectorField::interpolateOrRead
+(
+    const word& fieldName,
+    const dictionary& dict,
+    bool& interpolateField
+) const
+{
+    if (dict.found(fieldName))
+    {
+        tmp<Field<Type>> tFld
+        (
+            new Field<Type>
+            (
+                fieldName,
+                dict,
+                this->patch().size()
+            )
+        );
+
+        interpolateField = false;
+        return tFld;
+    }
+    else
+    {
+        interpolateField = true;
+        return interpolateBoundaryData<Type>(fieldName);
+    }
+}
+
+
+template<class Type>
+Foam::tmp<Foam::Field<Type>>
+Foam::turbulentDigitalFilterInletFvPatchVectorField::interpolateBoundaryData
+(
+    const word& fieldName
+) const
+{
+    const word& patchName = this->patch().name();
+
+    fileName valsFile
+    (
+        fileHandler().filePath
+        (
+            fileName
+            (
+                this->db().time().path()
+               /this->db().time().caseConstant()
+               /"boundaryData"
+               /patchName
+               /"0"
+               /fieldName
+            )
+        )
+    );
+
+    autoPtr<ISstream> isPtr
+    (
+        fileHandler().NewIFstream
+        (
+            valsFile
+        )
+    );
+
+    Field<Type> vals(isPtr());
+
+    Info<< "Turbulent DFM/FSM patch " << patchName
+        << ": Interpolating field " << fieldName
+        << " from " << valsFile << endl;
+
+    return patchMapper().interpolate(vals);
+}
+
+
+template<class Form, class Type>
+Form Foam::turbulentDigitalFilterInletFvPatchVectorField::generateRandomSet
+(
+    const label len
+)
+{
+    Form randomSet(len);
+
+    std::generate
+    (
+        randomSet.begin(),
+        randomSet.end(),
+        [&]{return rndGen_.GaussNormal<Type>();}
+    );
+
+    return randomSet;
+}
+
+// ************************************************************************* //
diff --git a/tutorials/verificationAndValidation/turbulentInflow/0.orig/U.DFSEM b/tutorials/verificationAndValidation/turbulentInflow/0.orig/U.DFSEM
new file mode 100644
index 0000000000..e3b5a0ad38
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/0.orig/U.DFSEM
@@ -0,0 +1,59 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    bottomWall
+    {
+        type            fixedValue;
+        value           uniform (0 0 0);
+    }
+    topWall
+    {
+        type            fixedValue;
+        value           uniform (0 0 0);
+    }
+    sides_half0
+    {
+        type            cyclic;
+    }
+    sides_half1
+    {
+        type            cyclic;
+    }
+    inlet
+    {
+        type            turbulentDFSEMInlet;
+        delta           2;
+        nCellPerEddy    1;
+        mapMethod       nearestCell;
+        value           uniform (0 0 0);
+    }
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      uniform (0 0 0);
+        value           uniform (0 0 0);
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/verificationAndValidation/turbulentInflow/0.orig/U.digitalFilter b/tutorials/verificationAndValidation/turbulentInflow/0.orig/U.digitalFilter
new file mode 100644
index 0000000000..876710fae8
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/0.orig/U.digitalFilter
@@ -0,0 +1,61 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    bottomWall
+    {
+        type            fixedValue;
+        value           uniform (0 0 0);
+    }
+    topWall
+    {
+        type            fixedValue;
+        value           uniform (0 0 0);
+    }
+    sides_half0
+    {
+        type            cyclic;
+    }
+    sides_half1
+    {
+        type            cyclic;
+    }
+    inlet
+    {
+        type                turbulentDigitalFilterInlet;
+        variant             digitalFilter;
+        planeDivisions      ( 64 70 );
+        L                   ( 0.78035508 0.31085352 0.342261 0.1728125 0.171875
+                              0.22459375 0.172787596 0.171889998 0.224578995 );
+        patchNormalSpeed    20.133;
+        value               uniform (0 0 0);
+    }
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      uniform (0 0 0);
+        value           uniform (0 0 0);
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/verificationAndValidation/turbulentInflow/0.orig/U.reducedDigitalFilter b/tutorials/verificationAndValidation/turbulentInflow/0.orig/U.reducedDigitalFilter
new file mode 100644
index 0000000000..e2a00e5816
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/0.orig/U.reducedDigitalFilter
@@ -0,0 +1,61 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    bottomWall
+    {
+        type            fixedValue;
+        value           uniform (0 0 0);
+    }
+    topWall
+    {
+        type            fixedValue;
+        value           uniform (0 0 0);
+    }
+    sides_half0
+    {
+        type            cyclic;
+    }
+    sides_half1
+    {
+        type            cyclic;
+    }
+    inlet
+    {
+        type                turbulentDigitalFilterInlet;
+        variant             reducedDigitalFilter;
+        planeDivisions      ( 64 70 );
+        L                   ( 0.78035508 0.31085352 0.342261 0.1728125 0.171875
+                              0.22459375 0.172787596 0.171889998 0.224578995 );
+        patchNormalSpeed    20.133;
+        value               uniform (0 0 0);
+    }
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      uniform (0 0 0);
+        value           uniform (0 0 0);
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/verificationAndValidation/turbulentInflow/0.orig/nut b/tutorials/verificationAndValidation/turbulentInflow/0.orig/nut
new file mode 100644
index 0000000000..2496e7ae86
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/0.orig/nut
@@ -0,0 +1,53 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      nut;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    bottomWall
+    {
+        type            zeroGradient;
+    }
+    topWall
+    {
+        type            zeroGradient;
+    }
+    sides_half0
+    {
+        type            cyclic;
+    }
+    sides_half1
+    {
+        type            cyclic;
+    }
+    inlet
+    {
+        type            calculated;
+        value           uniform 1e-08;
+    }
+    outlet
+    {
+        type            calculated;
+        value           uniform 1e-08;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/verificationAndValidation/turbulentInflow/0.orig/p b/tutorials/verificationAndValidation/turbulentInflow/0.orig/p
new file mode 100644
index 0000000000..73a85d2823
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/0.orig/p
@@ -0,0 +1,52 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      p;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -2 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    bottomWall
+    {
+        type            zeroGradient;
+    }
+    topWall
+    {
+        type            zeroGradient;
+    }
+    sides_half0
+    {
+        type            cyclic;
+    }
+    sides_half1
+    {
+        type            cyclic;
+    }
+    inlet
+    {
+        type            zeroGradient;
+    }
+    outlet
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/verificationAndValidation/turbulentInflow/Allclean b/tutorials/verificationAndValidation/turbulentInflow/Allclean
new file mode 100755
index 0000000000..213d289a42
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/Allclean
@@ -0,0 +1,8 @@
+#!/bin/sh
+
+cd ${0%/*} || exit 1    # run from this directory
+
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions  # Tutorial clean functions
+
+cleanCase0
+\rm -r results
diff --git a/tutorials/verificationAndValidation/turbulentInflow/Allrun b/tutorials/verificationAndValidation/turbulentInflow/Allrun
new file mode 100755
index 0000000000..e2679d94ec
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/Allrun
@@ -0,0 +1,34 @@
+#!/bin/sh
+
+cd ${0%/*} || exit 1    # run from this directory
+
+. $WM_PROJECT_DIR/bin/tools/RunFunctions  # Tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions  # Tutorial clean functions
+
+# Compute test case with different synthetic inflow models
+declare -a Models=(
+    "reducedDigitalFilter"
+    "digitalFilter"
+    "DFSEM"
+)
+
+# Collect data in 'results' directory
+[ -d "results" ] || mkdir results
+restore0Dir
+runApplication blockMesh
+
+# Compute and collect data
+for val in ${Models[@]}; do
+  echo "Running the case with the model: $val"
+  \cp 0/U.$val 0/U
+  \rm -r constant/boundaryData/inlet
+  \cp -r constant/boundaryData/inlet.$val constant/boundaryData/inlet
+  runApplication pimpleFoam
+  gnuplot plot.patch
+  gnuplot plot.cell
+  mv postProcessing results/postProcessing.$val
+  mv stress* results/postProcessing.$val/.
+  mv log* results/postProcessing.$val/.
+  cleanTimeDirectories
+done
+
diff --git a/tutorials/verificationAndValidation/turbulentInflow/Allrunparallel b/tutorials/verificationAndValidation/turbulentInflow/Allrunparallel
new file mode 100755
index 0000000000..78f16766c6
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/Allrunparallel
@@ -0,0 +1,37 @@
+#!/bin/sh
+
+cd ${0%/*} || exit 1    # run from this directory
+
+. $WM_PROJECT_DIR/bin/tools/RunFunctions  # Tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions  # Tutorial clean functions
+
+# Compute test case with different synthetic inflow models
+declare -a Models=(
+    "reducedDigitalFilter"
+    "digitalFilter"
+    "DFSEM"
+)
+
+# Collect data in 'results' directory
+[ -d "results" ] || mkdir results
+restore0Dir
+runApplication blockMesh
+
+# Compute parallel and collect data
+for val in ${Models[@]}; do
+  echo "Running the case with the model: $val"
+  \rm 0/U*
+  \cp 0.orig/U.$val 0/U
+  \rm -r constant/boundaryData/inlet
+  \cp -r constant/boundaryData/inlet.$val constant/boundaryData/inlet
+  runApplication decomposePar
+  runParallel pimpleFoam
+  gnuplot plot.patch
+  gnuplot plot.cell
+  mv postProcessing results/postProcessing.$val
+  mv stress* results/postProcessing.$val/.
+  mv log* results/postProcessing.$val/.
+  cleanTimeDirectories
+  \rm -rf processor* > /dev/null 2>&1
+done
+
diff --git a/tutorials/verificationAndValidation/turbulentInflow/README b/tutorials/verificationAndValidation/turbulentInflow/README
new file mode 100644
index 0000000000..e7d02e3af7
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/README
@@ -0,0 +1,34 @@
+Synthetic turbulence inflow tests
+======================
+
+The following three synthetic turbulence inflow boundary conditions are
+examined through a single-cell-domain smooth-wall plane channel flow setup:
+
+- turbulentDFSEMInlet
+- turbulentDigitalFilterInlet variant=digitalFilter
+- turbulentDigitalFilterInlet variant=reducedDigitalFilter
+
+The input statistics are obtained from:
+
+    Moser, Kim & Mansour
+    "DNS of Turbulent Channel Flow up to Re_tau=590",
+    Physics of Fluids (1999) vol 11, 943-945.
+    doi.org/10.1063/1.869966
+
+from which the input first-/second-order turbulence statistics data for the 
+smooth-wall plane channel flow at Re=395 were used.
+
+The data is available online from (Retrieved: 21-06-2019):
+
+    http://turbulence.ices.utexas.edu/MKM_1999.html
+
+Serial executing:
+
+./Allrun
+
+Parallel (decompositionMethod=scotch) executing:
+
+./Allrunparallel
+
+Both scripts run the test case, and collect the OpenFOAM samples, and plots into
+the 'results' directory. 
diff --git a/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.DFSEM/0/L b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.DFSEM/0/L
new file mode 100644
index 0000000000..9c6e3ab67d
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.DFSEM/0/L
@@ -0,0 +1,259 @@
+(
+0
+1.76693e-06
+1.06447e-07
+2.85582e-08
+1.73393e-08
+3.56243e-08
+3.79499e-08
+1.50479e-07
+1.55614e-07
+1.93408e-07
+4.50118e-08
+2.6295e-10
+6.18128e-08
+4.4709e-08
+5.98657e-07
+2.25711e-06
+4.20946e-06
+6.292e-06
+7.5414e-06
+9.02814e-06
+1.04737e-05
+1.25187e-05
+1.37559e-05
+1.40338e-05
+1.24681e-05
+1.08272e-05
+9.10144e-06
+7.81661e-06
+6.08627e-06
+4.68729e-06
+3.01241e-06
+1.93855e-06
+1.03945e-06
+2.49601e-07
+2.33919e-07
+3.47823e-08
+3.05604e-06
+8.25413e-06
+1.04748e-05
+1.36651e-05
+2.35252e-05
+3.87371e-05
+5.21582e-05
+6.71737e-05
+6.89399e-05
+4.57512e-05
+1.94891e-05
+9.7778e-06
+1.47858e-05
+3.40409e-05
+6.06823e-05
+8.59513e-05
+8.29571e-05
+5.42278e-05
+3.46015e-05
+2.20432e-05
+1.73878e-05
+1.05749e-05
+4.67814e-06
+8.66793e-06
+1.79361e-05
+1.22464e-05
+3.93103e-06
+1.00778e-06
+1.81983e-06
+2.70815e-05
+9.56468e-05
+0.000139304
+0.00012493
+9.50459e-05
+4.78307e-05
+8.71962e-06
+1.24275e-05
+4.91383e-05
+8.77264e-05
+0.000114449
+0.00014567
+0.000201758
+0.000251863
+0.000272594
+0.000213599
+0.000145126
+0.000115603
+0.000122779
+0.000116865
+8.36843e-05
+4.06009e-05
+3.49149e-05
+5.22978e-05
+5.07525e-05
+2.5995e-05
+2.6683e-06
+9.29144e-06
+1.14821e-05
+9.41939e-06
+1.01946e-05
+6.64024e-06
+1.1913e-06
+3.25066e-06
+8.47834e-06
+1.42023e-06
+4.30742e-05
+0.000106228
+0.000115468
+9.07632e-05
+7.04511e-05
+9.21776e-05
+0.000104486
+0.000111678
+0.000109852
+0.000111867
+9.10747e-05
+6.82206e-05
+4.09085e-05
+2.3961e-05
+2.0467e-06
+7.74345e-06
+1.66716e-05
+1.54967e-05
+2.95089e-05
+4.82299e-05
+6.99781e-05
+7.16947e-05
+7.33475e-05
+7.40551e-05
+9.45846e-05
+0.000107202
+0.000120068
+0.000122517
+0.000120068
+0.000107202
+9.45846e-05
+7.40551e-05
+7.33475e-05
+7.16947e-05
+6.99781e-05
+4.82299e-05
+2.95089e-05
+1.54967e-05
+1.66716e-05
+7.74345e-06
+2.0467e-06
+2.3961e-05
+4.09085e-05
+6.82206e-05
+9.10747e-05
+0.000111867
+0.000109852
+0.000111678
+0.000104486
+9.21776e-05
+7.04511e-05
+9.07632e-05
+0.000115468
+0.000106228
+4.30742e-05
+1.42023e-06
+8.47834e-06
+3.25066e-06
+1.1913e-06
+6.64024e-06
+1.01946e-05
+9.41939e-06
+1.14821e-05
+9.29144e-06
+2.6683e-06
+2.5995e-05
+5.07525e-05
+5.22978e-05
+3.49149e-05
+4.06009e-05
+8.36843e-05
+0.000116865
+0.000122779
+0.000115603
+0.000145126
+0.000213599
+0.000272594
+0.000251863
+0.000201758
+0.00014567
+0.000114449
+8.77264e-05
+4.91383e-05
+1.24275e-05
+8.71962e-06
+4.78307e-05
+9.50459e-05
+0.00012493
+0.000139304
+9.56468e-05
+2.70815e-05
+1.81983e-06
+1.00778e-06
+3.93103e-06
+1.22464e-05
+1.79361e-05
+8.66793e-06
+4.67814e-06
+1.05749e-05
+1.73878e-05
+2.20432e-05
+3.46015e-05
+5.42278e-05
+8.29571e-05
+8.59513e-05
+6.06823e-05
+3.40409e-05
+1.47858e-05
+9.7778e-06
+1.94891e-05
+4.57512e-05
+6.89399e-05
+6.71737e-05
+5.21582e-05
+3.87371e-05
+2.35252e-05
+1.36651e-05
+1.04748e-05
+8.25413e-06
+3.05604e-06
+3.47823e-08
+2.33919e-07
+2.49601e-07
+1.03945e-06
+1.93855e-06
+3.01241e-06
+4.68729e-06
+6.08627e-06
+7.81661e-06
+9.10144e-06
+1.08272e-05
+1.24681e-05
+1.40338e-05
+1.37559e-05
+1.25187e-05
+1.04737e-05
+9.02814e-06
+7.5414e-06
+6.292e-06
+4.20946e-06
+2.25711e-06
+5.98657e-07
+4.4709e-08
+6.18128e-08
+2.6295e-10
+4.50118e-08
+1.93408e-07
+1.55614e-07
+1.50479e-07
+3.79499e-08
+3.56243e-08
+1.73393e-08
+2.85582e-08
+1.06447e-07
+1.76693e-06
+0
+)
diff --git a/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.DFSEM/0/R b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.DFSEM/0/R
new file mode 100644
index 0000000000..084bdd9f8f
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.DFSEM/0/R
@@ -0,0 +1,259 @@
+(
+(5.3376e-28 1.1206e-30 -6.1282e-30 4.9018e-29 6.2793e-32 5.9575e-28)
+(1.3684e-04 -2.4414e-08 5.9464e-08 9.4734e-11 -2.8525e-11 5.2943e-05)
+(2.1828e-03 -1.5728e-06 8.3228e-07 2.3603e-08 -1.7141e-09 8.2715e-04)
+(1.1020e-02 -1.8113e-05 3.3616e-06 5.7878e-07 -1.9051e-08 4.0392e-03)
+(3.4718e-02 -1.0324e-04 1.0732e-05 5.4400e-06 -1.0411e-07 1.2147e-02)
+(8.4454e-02 -4.0027e-04 2.8297e-05 3.0014e-05 -3.8272e-07 2.7852e-02)
+(1.7438e-01 -1.2143e-03 6.6865e-05 1.1756e-04 -1.0969e-06 5.3564e-02)
+(3.2123e-01 -3.1009e-03 1.4506e-04 3.6204e-04 -2.6288e-06 9.0947e-02)
+(5.4336e-01 -6.9486e-03 2.9097e-04 9.3203e-04 -5.4596e-06 1.4062e-01)
+(8.5851e-01 -1.4015e-02 5.4723e-04 2.0879e-03 -1.0091e-05 2.0209e-01)
+(1.2802e+00 -2.5858e-02 9.7373e-04 4.1852e-03 -1.7119e-05 2.7384e-01)
+(1.8132e+00 -4.4135e-02 1.6529e-03 7.6610e-03 -2.7764e-05 3.5369e-01)
+(2.4490e+00 -7.0276e-02 2.6847e-03 1.3005e-02 -4.4587e-05 4.3919e-01)
+(3.1644e+00 -1.0513e-01 4.1721e-03 2.0724e-02 -7.1044e-05 5.2804e-01)
+(3.9229e+00 -1.4869e-01 6.1889e-03 3.1300e-02 -1.0893e-04 6.1834e-01)
+(4.6811e+00 -1.9999e-01 8.7562e-03 4.5157e-02 -1.5460e-04 7.0871e-01)
+(5.3953e+00 -2.5724e-01 1.1834e-02 6.2636e-02 -1.9552e-04 7.9819e-01)
+(6.0286e+00 -3.1808e-01 1.5387e-02 8.3971e-02 -2.1029e-04 8.8610e-01)
+(6.5558e+00 -3.8001e-01 1.9449e-02 1.0928e-01 -1.7262e-04 9.7178e-01)
+(6.9642e+00 -4.4073e-01 2.4133e-02 1.3854e-01 -6.0490e-05 1.0545e+00)
+(7.2524e+00 -4.9837e-01 2.9501e-02 1.7162e-01 1.3268e-04 1.1333e+00)
+(7.4274e+00 -5.5162e-01 3.5418e-02 2.0824e-01 3.9086e-04 1.2075e+00)
+(7.5016e+00 -5.9969e-01 4.1511e-02 2.4803e-01 6.7818e-04 1.2762e+00)
+(7.4897e+00 -6.4224e-01 4.7330e-02 2.9050e-01 9.5511e-04 1.3387e+00)
+(7.4069e+00 -6.7924e-01 5.2512e-02 3.3510e-01 1.1989e-03 1.3948e+00)
+(7.2679e+00 -7.1090e-01 5.6926e-02 3.8121e-01 1.4016e-03 1.4444e+00)
+(7.0865e+00 -7.3759e-01 6.0659e-02 4.2818e-01 1.5532e-03 1.4878e+00)
+(6.8750e+00 -7.5980e-01 6.3835e-02 4.7535e-01 1.6360e-03 1.5251e+00)
+(6.6441e+00 -7.7810e-01 6.6418e-02 5.2211e-01 1.6465e-03 1.5568e+00)
+(6.4025e+00 -7.9302e-01 6.8207e-02 5.6790e-01 1.6026e-03 1.5833e+00)
+(6.1573e+00 -8.0503e-01 6.9001e-02 6.1229e-01 1.5073e-03 1.6047e+00)
+(5.9141e+00 -8.1452e-01 6.8896e-02 6.5493e-01 1.3268e-03 1.6216e+00)
+(5.6767e+00 -8.2184e-01 6.8285e-02 6.9554e-01 1.0303e-03 1.6344e+00)
+(5.4477e+00 -8.2733e-01 6.7481e-02 7.3388e-01 6.4119e-04 1.6439e+00)
+(5.2291e+00 -8.3133e-01 6.6450e-02 7.6971e-01 2.3642e-04 1.6512e+00)
+(5.0223e+00 -8.3412e-01 6.4937e-02 8.0285e-01 -1.3612e-04 1.6567e+00)
+(4.8279e+00 -8.3590e-01 6.2875e-02 8.3316e-01 -4.7569e-04 1.6606e+00)
+(4.6462e+00 -8.3683e-01 6.0710e-02 8.6061e-01 -7.8412e-04 1.6623e+00)
+(4.4768e+00 -8.3700e-01 5.9182e-02 8.8520e-01 -1.0539e-03 1.6618e+00)
+(4.3194e+00 -8.3649e-01 5.8527e-02 9.0700e-01 -1.2655e-03 1.6594e+00)
+(4.1730e+00 -8.3523e-01 5.8437e-02 9.2607e-01 -1.4126e-03 1.6553e+00)
+(4.0370e+00 -8.3321e-01 5.8512e-02 9.4247e-01 -1.4930e-03 1.6500e+00)
+(3.9106e+00 -8.3044e-01 5.8376e-02 9.5626e-01 -1.5118e-03 1.6439e+00)
+(3.7930e+00 -8.2687e-01 5.8084e-02 9.6754e-01 -1.5196e-03 1.6373e+00)
+(3.6834e+00 -8.2252e-01 5.8125e-02 9.7651e-01 -1.6228e-03 1.6299e+00)
+(3.5821e+00 -8.1762e-01 5.8537e-02 9.8341e-01 -1.8295e-03 1.6209e+00)
+(3.4896e+00 -8.1249e-01 5.8809e-02 9.8850e-01 -2.0145e-03 1.6105e+00)
+(3.4053e+00 -8.0726e-01 5.8494e-02 9.9197e-01 -2.0826e-03 1.5991e+00)
+(3.3278e+00 -8.0180e-01 5.7876e-02 9.9396e-01 -2.1101e-03 1.5872e+00)
+(3.2557e+00 -7.9592e-01 5.7034e-02 9.9455e-01 -2.1628e-03 1.5754e+00)
+(3.1886e+00 -7.8961e-01 5.5746e-02 9.9379e-01 -2.2441e-03 1.5637e+00)
+(3.1269e+00 -7.8296e-01 5.4259e-02 9.9183e-01 -2.3220e-03 1.5510e+00)
+(3.0702e+00 -7.7613e-01 5.2957e-02 9.8887e-01 -2.3519e-03 1.5367e+00)
+(3.0178e+00 -7.6929e-01 5.1819e-02 9.8512e-01 -2.3018e-03 1.5212e+00)
+(2.9678e+00 -7.6243e-01 5.0872e-02 9.8075e-01 -2.0851e-03 1.5052e+00)
+(2.9191e+00 -7.5552e-01 5.0514e-02 9.7583e-01 -1.6209e-03 1.4889e+00)
+(2.8712e+00 -7.4851e-01 5.0564e-02 9.7044e-01 -8.5802e-04 1.4721e+00)
+(2.8239e+00 -7.4146e-01 5.0229e-02 9.6465e-01 1.1783e-04 1.4551e+00)
+(2.7776e+00 -7.3445e-01 4.9171e-02 9.5846e-01 1.0337e-03 1.4379e+00)
+(2.7324e+00 -7.2737e-01 4.7565e-02 9.5179e-01 1.6440e-03 1.4206e+00)
+(2.6883e+00 -7.2012e-01 4.5716e-02 9.4459e-01 1.9851e-03 1.4035e+00)
+(2.6453e+00 -7.1271e-01 4.3916e-02 9.3679e-01 2.2291e-03 1.3867e+00)
+(2.6031e+00 -7.0521e-01 4.2134e-02 9.2845e-01 2.4143e-03 1.3696e+00)
+(2.5612e+00 -6.9761e-01 4.0077e-02 9.1976e-01 2.5128e-03 1.3515e+00)
+(2.5196e+00 -6.8985e-01 3.7520e-02 9.1087e-01 2.5006e-03 1.3322e+00)
+(2.4793e+00 -6.8192e-01 3.4358e-02 9.0185e-01 2.3726e-03 1.3118e+00)
+(2.4411e+00 -6.7396e-01 3.0858e-02 8.9268e-01 2.1952e-03 1.2910e+00)
+(2.4055e+00 -6.6598e-01 2.7094e-02 8.8335e-01 2.0517e-03 1.2700e+00)
+(2.3714e+00 -6.5794e-01 2.3371e-02 8.7385e-01 1.9999e-03 1.2488e+00)
+(2.3380e+00 -6.4977e-01 2.0500e-02 8.6418e-01 1.9700e-03 1.2276e+00)
+(2.3046e+00 -6.4148e-01 1.8795e-02 8.5436e-01 1.8189e-03 1.2066e+00)
+(2.2708e+00 -6.3298e-01 1.7891e-02 8.4439e-01 1.4336e-03 1.1855e+00)
+(2.2365e+00 -6.2415e-01 1.7572e-02 8.3421e-01 7.5375e-04 1.1641e+00)
+(2.2023e+00 -6.1500e-01 1.7517e-02 8.2369e-01 -7.9290e-05 1.1423e+00)
+(2.1682e+00 -6.0551e-01 1.7095e-02 8.1285e-01 -8.3720e-04 1.1202e+00)
+(2.1335e+00 -5.9561e-01 1.5913e-02 8.0181e-01 -1.3413e-03 1.0983e+00)
+(2.0983e+00 -5.8530e-01 1.4306e-02 7.9074e-01 -1.5771e-03 1.0767e+00)
+(2.0631e+00 -5.7460e-01 1.2812e-02 7.7968e-01 -1.5179e-03 1.0551e+00)
+(2.0277e+00 -5.6350e-01 1.1907e-02 7.6856e-01 -1.1762e-03 1.0339e+00)
+(1.9920e+00 -5.5204e-01 1.1640e-02 7.5731e-01 -7.0498e-04 1.0131e+00)
+(1.9568e+00 -5.4040e-01 1.1503e-02 7.4588e-01 -2.0689e-04 9.9279e-01)
+(1.9229e+00 -5.2883e-01 1.1097e-02 7.3425e-01 2.0529e-04 9.7254e-01)
+(1.8896e+00 -5.1747e-01 1.0218e-02 7.2249e-01 5.0971e-04 9.5187e-01)
+(1.8562e+00 -5.0626e-01 9.0004e-03 7.1064e-01 7.6211e-04 9.3057e-01)
+(1.8232e+00 -4.9507e-01 7.8753e-03 6.9862e-01 9.2292e-04 9.0936e-01)
+(1.7908e+00 -4.8386e-01 7.1617e-03 6.8642e-01 9.3017e-04 8.8926e-01)
+(1.7591e+00 -4.7269e-01 7.0121e-03 6.7417e-01 7.2814e-04 8.7042e-01)
+(1.7274e+00 -4.6164e-01 7.3802e-03 6.6201e-01 2.7352e-04 8.5243e-01)
+(1.6948e+00 -4.5072e-01 8.2299e-03 6.5007e-01 -2.9303e-04 8.3469e-01)
+(1.6610e+00 -4.3999e-01 9.3954e-03 6.3843e-01 -7.4565e-04 8.1690e-01)
+(1.6264e+00 -4.2949e-01 1.0585e-02 6.2712e-01 -9.5384e-04 7.9939e-01)
+(1.5910e+00 -4.1910e-01 1.1834e-02 6.1613e-01 -9.5389e-04 7.8241e-01)
+(1.5550e+00 -4.0876e-01 1.3728e-02 6.0545e-01 -1.0075e-03 7.6587e-01)
+(1.5184e+00 -3.9853e-01 1.6563e-02 5.9512e-01 -1.2970e-03 7.4931e-01)
+(1.4808e+00 -3.8843e-01 1.9835e-02 5.8523e-01 -1.7620e-03 7.3258e-01)
+(1.4425e+00 -3.7835e-01 2.2555e-02 5.7589e-01 -2.2305e-03 7.1572e-01)
+(1.4037e+00 -3.6810e-01 2.3965e-02 5.6715e-01 -2.5637e-03 6.9868e-01)
+(1.3647e+00 -3.5761e-01 2.4099e-02 5.5896e-01 -2.7183e-03 6.8159e-01)
+(1.3260e+00 -3.4695e-01 2.3378e-02 5.5120e-01 -2.7095e-03 6.6492e-01)
+(1.2881e+00 -3.3621e-01 2.2236e-02 5.4372e-01 -2.5722e-03 6.4916e-01)
+(1.2513e+00 -3.2544e-01 2.0871e-02 5.3642e-01 -2.4646e-03 6.3424e-01)
+(1.2162e+00 -3.1468e-01 1.9224e-02 5.2919e-01 -2.5472e-03 6.2001e-01)
+(1.1827e+00 -3.0393e-01 1.7273e-02 5.2197e-01 -2.7871e-03 6.0660e-01)
+(1.1500e+00 -2.9315e-01 1.5116e-02 5.1473e-01 -3.0094e-03 5.9377e-01)
+(1.1175e+00 -2.8231e-01 1.3055e-02 5.0750e-01 -3.1202e-03 5.8101e-01)
+(1.0843e+00 -2.7133e-01 1.1360e-02 5.0029e-01 -3.1772e-03 5.6815e-01)
+(1.0507e+00 -2.6019e-01 9.9787e-03 4.9312e-01 -3.2410e-03 5.5541e-01)
+(1.0173e+00 -2.4893e-01 8.7677e-03 4.8605e-01 -3.3171e-03 5.4304e-01)
+(9.8478e-01 -2.3753e-01 7.8868e-03 4.7922e-01 -3.3738e-03 5.3113e-01)
+(9.5377e-01 -2.2605e-01 7.6156e-03 4.7268e-01 -3.4135e-03 5.1971e-01)
+(9.2454e-01 -2.1462e-01 7.9424e-03 4.6644e-01 -3.4897e-03 5.0882e-01)
+(8.9706e-01 -2.0326e-01 8.6379e-03 4.6047e-01 -3.6720e-03 4.9853e-01)
+(8.7113e-01 -1.9188e-01 9.5840e-03 4.5478e-01 -3.9882e-03 4.8883e-01)
+(8.4645e-01 -1.8043e-01 1.0540e-02 4.4932e-01 -4.3896e-03 4.7982e-01)
+(8.2296e-01 -1.6898e-01 1.1159e-02 4.4410e-01 -4.7685e-03 4.7164e-01)
+(8.0062e-01 -1.5752e-01 1.1209e-02 4.3915e-01 -5.0417e-03 4.6392e-01)
+(7.7929e-01 -1.4601e-01 1.0617e-02 4.3458e-01 -5.1687e-03 4.5618e-01)
+(7.5889e-01 -1.3442e-01 9.5260e-03 4.3044e-01 -5.1040e-03 4.4842e-01)
+(7.3932e-01 -1.2272e-01 8.2998e-03 4.2671e-01 -4.8287e-03 4.4071e-01)
+(7.2051e-01 -1.1094e-01 7.2287e-03 4.2339e-01 -4.4056e-03 4.3345e-01)
+(7.0261e-01 -9.9120e-02 6.4502e-03 4.2046e-01 -3.8984e-03 4.2708e-01)
+(6.8576e-01 -8.7301e-02 6.1696e-03 4.1796e-01 -3.3462e-03 4.2169e-01)
+(6.7023e-01 -7.5441e-02 6.4277e-03 4.1591e-01 -2.7927e-03 4.1709e-01)
+(6.5620e-01 -6.3452e-02 7.0238e-03 4.1431e-01 -2.2376e-03 4.1319e-01)
+(6.4389e-01 -5.1252e-02 7.7818e-03 4.1313e-01 -1.6855e-03 4.0984e-01)
+(6.3375e-01 -3.8772e-02 8.7489e-03 4.1234e-01 -1.1858e-03 4.0689e-01)
+(6.2626e-01 -2.6001e-02 9.8627e-03 4.1188e-01 -7.6554e-04 4.0431e-01)
+(6.2174e-01 -1.3031e-02 1.0731e-02 4.1166e-01 -3.8835e-04 4.0247e-01)
+(6.2024e-01 0.0000e+00 1.1032e-02 4.1160e-01 0.0000e+00 4.0179e-01)
+(6.2174e-01 1.3031e-02 1.0731e-02 4.1166e-01 3.8835e-04 4.0247e-01)
+(6.2626e-01 2.6001e-02 9.8627e-03 4.1188e-01 7.6554e-04 4.0431e-01)
+(6.3375e-01 3.8772e-02 8.7489e-03 4.1234e-01 1.1858e-03 4.0689e-01)
+(6.4389e-01 5.1252e-02 7.7818e-03 4.1313e-01 1.6855e-03 4.0984e-01)
+(6.5620e-01 6.3452e-02 7.0238e-03 4.1431e-01 2.2376e-03 4.1319e-01)
+(6.7023e-01 7.5441e-02 6.4277e-03 4.1591e-01 2.7927e-03 4.1709e-01)
+(6.8576e-01 8.7301e-02 6.1696e-03 4.1796e-01 3.3462e-03 4.2169e-01)
+(7.0261e-01 9.9120e-02 6.4502e-03 4.2046e-01 3.8984e-03 4.2708e-01)
+(7.2051e-01 1.1094e-01 7.2287e-03 4.2339e-01 4.4056e-03 4.3345e-01)
+(7.3932e-01 1.2272e-01 8.2998e-03 4.2671e-01 4.8287e-03 4.4071e-01)
+(7.5889e-01 1.3442e-01 9.5260e-03 4.3044e-01 5.1040e-03 4.4842e-01)
+(7.7929e-01 1.4601e-01 1.0617e-02 4.3458e-01 5.1687e-03 4.5618e-01)
+(8.0062e-01 1.5752e-01 1.1209e-02 4.3915e-01 5.0417e-03 4.6392e-01)
+(8.2296e-01 1.6898e-01 1.1159e-02 4.4410e-01 4.7685e-03 4.7164e-01)
+(8.4645e-01 1.8043e-01 1.0540e-02 4.4932e-01 4.3896e-03 4.7982e-01)
+(8.7113e-01 1.9188e-01 9.5840e-03 4.5478e-01 3.9882e-03 4.8883e-01)
+(8.9706e-01 2.0326e-01 8.6379e-03 4.6047e-01 3.6720e-03 4.9853e-01)
+(9.2454e-01 2.1462e-01 7.9424e-03 4.6644e-01 3.4897e-03 5.0882e-01)
+(9.5377e-01 2.2605e-01 7.6156e-03 4.7268e-01 3.4135e-03 5.1971e-01)
+(9.8478e-01 2.3753e-01 7.8868e-03 4.7922e-01 3.3738e-03 5.3113e-01)
+(1.0173e+00 2.4893e-01 8.7677e-03 4.8605e-01 3.3171e-03 5.4304e-01)
+(1.0507e+00 2.6019e-01 9.9787e-03 4.9312e-01 3.2410e-03 5.5541e-01)
+(1.0843e+00 2.7133e-01 1.1360e-02 5.0029e-01 3.1772e-03 5.6815e-01)
+(1.1175e+00 2.8231e-01 1.3055e-02 5.0750e-01 3.1202e-03 5.8101e-01)
+(1.1500e+00 2.9315e-01 1.5116e-02 5.1473e-01 3.0094e-03 5.9377e-01)
+(1.1827e+00 3.0393e-01 1.7273e-02 5.2197e-01 2.7871e-03 6.0660e-01)
+(1.2162e+00 3.1468e-01 1.9224e-02 5.2919e-01 2.5472e-03 6.2001e-01)
+(1.2513e+00 3.2544e-01 2.0871e-02 5.3642e-01 2.4646e-03 6.3424e-01)
+(1.2881e+00 3.3621e-01 2.2236e-02 5.4372e-01 2.5722e-03 6.4916e-01)
+(1.3260e+00 3.4695e-01 2.3378e-02 5.5120e-01 2.7095e-03 6.6492e-01)
+(1.3647e+00 3.5761e-01 2.4099e-02 5.5896e-01 2.7183e-03 6.8159e-01)
+(1.4037e+00 3.6810e-01 2.3965e-02 5.6715e-01 2.5637e-03 6.9868e-01)
+(1.4425e+00 3.7835e-01 2.2555e-02 5.7589e-01 2.2305e-03 7.1572e-01)
+(1.4808e+00 3.8843e-01 1.9835e-02 5.8523e-01 1.7620e-03 7.3258e-01)
+(1.5184e+00 3.9853e-01 1.6563e-02 5.9512e-01 1.2970e-03 7.4931e-01)
+(1.5550e+00 4.0876e-01 1.3728e-02 6.0545e-01 1.0075e-03 7.6587e-01)
+(1.5910e+00 4.1910e-01 1.1834e-02 6.1613e-01 9.5389e-04 7.8241e-01)
+(1.6264e+00 4.2949e-01 1.0585e-02 6.2712e-01 9.5384e-04 7.9939e-01)
+(1.6610e+00 4.3999e-01 9.3954e-03 6.3843e-01 7.4565e-04 8.1690e-01)
+(1.6948e+00 4.5072e-01 8.2299e-03 6.5007e-01 2.9303e-04 8.3469e-01)
+(1.7274e+00 4.6164e-01 7.3802e-03 6.6201e-01 -2.7352e-04 8.5243e-01)
+(1.7591e+00 4.7269e-01 7.0121e-03 6.7417e-01 -7.2814e-04 8.7042e-01)
+(1.7908e+00 4.8386e-01 7.1617e-03 6.8642e-01 -9.3017e-04 8.8926e-01)
+(1.8232e+00 4.9507e-01 7.8753e-03 6.9862e-01 -9.2292e-04 9.0936e-01)
+(1.8562e+00 5.0626e-01 9.0004e-03 7.1064e-01 -7.6211e-04 9.3057e-01)
+(1.8896e+00 5.1747e-01 1.0218e-02 7.2249e-01 -5.0971e-04 9.5187e-01)
+(1.9229e+00 5.2883e-01 1.1097e-02 7.3425e-01 -2.0529e-04 9.7254e-01)
+(1.9568e+00 5.4040e-01 1.1503e-02 7.4588e-01 2.0689e-04 9.9279e-01)
+(1.9920e+00 5.5204e-01 1.1640e-02 7.5731e-01 7.0498e-04 1.0131e+00)
+(2.0277e+00 5.6350e-01 1.1907e-02 7.6856e-01 1.1762e-03 1.0339e+00)
+(2.0631e+00 5.7460e-01 1.2812e-02 7.7968e-01 1.5179e-03 1.0551e+00)
+(2.0983e+00 5.8530e-01 1.4306e-02 7.9074e-01 1.5771e-03 1.0767e+00)
+(2.1335e+00 5.9561e-01 1.5913e-02 8.0181e-01 1.3413e-03 1.0983e+00)
+(2.1682e+00 6.0551e-01 1.7095e-02 8.1285e-01 8.3720e-04 1.1202e+00)
+(2.2023e+00 6.1500e-01 1.7517e-02 8.2369e-01 7.9290e-05 1.1423e+00)
+(2.2365e+00 6.2415e-01 1.7572e-02 8.3421e-01 -7.5375e-04 1.1641e+00)
+(2.2708e+00 6.3298e-01 1.7891e-02 8.4439e-01 -1.4336e-03 1.1855e+00)
+(2.3046e+00 6.4148e-01 1.8795e-02 8.5436e-01 -1.8189e-03 1.2066e+00)
+(2.3380e+00 6.4977e-01 2.0500e-02 8.6418e-01 -1.9700e-03 1.2276e+00)
+(2.3714e+00 6.5794e-01 2.3371e-02 8.7385e-01 -1.9999e-03 1.2488e+00)
+(2.4055e+00 6.6598e-01 2.7094e-02 8.8335e-01 -2.0517e-03 1.2700e+00)
+(2.4411e+00 6.7396e-01 3.0858e-02 8.9268e-01 -2.1952e-03 1.2910e+00)
+(2.4793e+00 6.8192e-01 3.4358e-02 9.0185e-01 -2.3726e-03 1.3118e+00)
+(2.5196e+00 6.8985e-01 3.7520e-02 9.1087e-01 -2.5006e-03 1.3322e+00)
+(2.5612e+00 6.9761e-01 4.0077e-02 9.1976e-01 -2.5128e-03 1.3515e+00)
+(2.6031e+00 7.0521e-01 4.2134e-02 9.2845e-01 -2.4143e-03 1.3696e+00)
+(2.6453e+00 7.1271e-01 4.3916e-02 9.3679e-01 -2.2291e-03 1.3867e+00)
+(2.6883e+00 7.2012e-01 4.5716e-02 9.4459e-01 -1.9851e-03 1.4035e+00)
+(2.7324e+00 7.2737e-01 4.7565e-02 9.5179e-01 -1.6440e-03 1.4206e+00)
+(2.7776e+00 7.3445e-01 4.9171e-02 9.5846e-01 -1.0337e-03 1.4379e+00)
+(2.8239e+00 7.4146e-01 5.0229e-02 9.6465e-01 -1.1783e-04 1.4551e+00)
+(2.8712e+00 7.4851e-01 5.0564e-02 9.7044e-01 8.5802e-04 1.4721e+00)
+(2.9191e+00 7.5552e-01 5.0514e-02 9.7583e-01 1.6209e-03 1.4889e+00)
+(2.9678e+00 7.6243e-01 5.0872e-02 9.8075e-01 2.0851e-03 1.5052e+00)
+(3.0178e+00 7.6929e-01 5.1819e-02 9.8512e-01 2.3018e-03 1.5212e+00)
+(3.0702e+00 7.7613e-01 5.2957e-02 9.8887e-01 2.3519e-03 1.5367e+00)
+(3.1269e+00 7.8296e-01 5.4259e-02 9.9183e-01 2.3220e-03 1.5510e+00)
+(3.1886e+00 7.8961e-01 5.5746e-02 9.9379e-01 2.2441e-03 1.5637e+00)
+(3.2557e+00 7.9592e-01 5.7034e-02 9.9455e-01 2.1628e-03 1.5754e+00)
+(3.3278e+00 8.0180e-01 5.7876e-02 9.9396e-01 2.1101e-03 1.5872e+00)
+(3.4053e+00 8.0726e-01 5.8494e-02 9.9197e-01 2.0826e-03 1.5991e+00)
+(3.4896e+00 8.1249e-01 5.8809e-02 9.8850e-01 2.0145e-03 1.6105e+00)
+(3.5821e+00 8.1762e-01 5.8537e-02 9.8341e-01 1.8295e-03 1.6209e+00)
+(3.6834e+00 8.2252e-01 5.8125e-02 9.7651e-01 1.6228e-03 1.6299e+00)
+(3.7930e+00 8.2687e-01 5.8084e-02 9.6754e-01 1.5196e-03 1.6373e+00)
+(3.9106e+00 8.3044e-01 5.8376e-02 9.5626e-01 1.5118e-03 1.6439e+00)
+(4.0370e+00 8.3321e-01 5.8512e-02 9.4247e-01 1.4930e-03 1.6500e+00)
+(4.1730e+00 8.3523e-01 5.8437e-02 9.2607e-01 1.4126e-03 1.6553e+00)
+(4.3194e+00 8.3649e-01 5.8527e-02 9.0700e-01 1.2655e-03 1.6594e+00)
+(4.4768e+00 8.3700e-01 5.9182e-02 8.8520e-01 1.0539e-03 1.6618e+00)
+(4.6462e+00 8.3683e-01 6.0710e-02 8.6061e-01 7.8412e-04 1.6623e+00)
+(4.8279e+00 8.3590e-01 6.2875e-02 8.3316e-01 4.7569e-04 1.6606e+00)
+(5.0223e+00 8.3412e-01 6.4937e-02 8.0285e-01 1.3612e-04 1.6567e+00)
+(5.2291e+00 8.3133e-01 6.6450e-02 7.6971e-01 -2.3642e-04 1.6512e+00)
+(5.4477e+00 8.2733e-01 6.7481e-02 7.3388e-01 -6.4119e-04 1.6439e+00)
+(5.6767e+00 8.2184e-01 6.8285e-02 6.9554e-01 -1.0303e-03 1.6344e+00)
+(5.9141e+00 8.1452e-01 6.8896e-02 6.5493e-01 -1.3268e-03 1.6216e+00)
+(6.1573e+00 8.0503e-01 6.9001e-02 6.1229e-01 -1.5073e-03 1.6047e+00)
+(6.4025e+00 7.9302e-01 6.8207e-02 5.6790e-01 -1.6026e-03 1.5833e+00)
+(6.6441e+00 7.7810e-01 6.6418e-02 5.2211e-01 -1.6465e-03 1.5568e+00)
+(6.8750e+00 7.5980e-01 6.3835e-02 4.7535e-01 -1.6360e-03 1.5251e+00)
+(7.0865e+00 7.3759e-01 6.0659e-02 4.2818e-01 -1.5532e-03 1.4878e+00)
+(7.2679e+00 7.1090e-01 5.6926e-02 3.8121e-01 -1.4016e-03 1.4444e+00)
+(7.4069e+00 6.7924e-01 5.2512e-02 3.3510e-01 -1.1989e-03 1.3948e+00)
+(7.4897e+00 6.4224e-01 4.7330e-02 2.9050e-01 -9.5511e-04 1.3387e+00)
+(7.5016e+00 5.9969e-01 4.1511e-02 2.4803e-01 -6.7818e-04 1.2762e+00)
+(7.4274e+00 5.5162e-01 3.5418e-02 2.0824e-01 -3.9086e-04 1.2075e+00)
+(7.2524e+00 4.9837e-01 2.9501e-02 1.7162e-01 -1.3268e-04 1.1333e+00)
+(6.9642e+00 4.4073e-01 2.4133e-02 1.3854e-01 6.0490e-05 1.0545e+00)
+(6.5558e+00 3.8001e-01 1.9449e-02 1.0928e-01 1.7262e-04 9.7178e-01)
+(6.0286e+00 3.1808e-01 1.5387e-02 8.3971e-02 2.1029e-04 8.8610e-01)
+(5.3953e+00 2.5724e-01 1.1834e-02 6.2636e-02 1.9552e-04 7.9819e-01)
+(4.6811e+00 1.9999e-01 8.7562e-03 4.5157e-02 1.5460e-04 7.0871e-01)
+(3.9229e+00 1.4869e-01 6.1889e-03 3.1300e-02 1.0893e-04 6.1834e-01)
+(3.1644e+00 1.0513e-01 4.1721e-03 2.0724e-02 7.1044e-05 5.2804e-01)
+(2.4490e+00 7.0276e-02 2.6847e-03 1.3005e-02 4.4587e-05 4.3919e-01)
+(1.8132e+00 4.4135e-02 1.6529e-03 7.6610e-03 2.7764e-05 3.5369e-01)
+(1.2802e+00 2.5858e-02 9.7373e-04 4.1852e-03 1.7119e-05 2.7384e-01)
+(8.5851e-01 1.4015e-02 5.4723e-04 2.0879e-03 1.0091e-05 2.0209e-01)
+(5.4336e-01 6.9486e-03 2.9097e-04 9.3203e-04 5.4596e-06 1.4062e-01)
+(3.2123e-01 3.1009e-03 1.4506e-04 3.6204e-04 2.6288e-06 9.0947e-02)
+(1.7438e-01 1.2143e-03 6.6865e-05 1.1756e-04 1.0969e-06 5.3564e-02)
+(8.4454e-02 4.0027e-04 2.8297e-05 3.0014e-05 3.8272e-07 2.7852e-02)
+(3.4718e-02 1.0324e-04 1.0732e-05 5.4400e-06 1.0411e-07 1.2147e-02)
+(1.1020e-02 1.8113e-05 3.3616e-06 5.7878e-07 1.9051e-08 4.0392e-03)
+(2.1828e-03 1.5728e-06 8.3228e-07 2.3603e-08 1.7141e-09 8.2715e-04)
+(1.3684e-04 2.4414e-08 5.9464e-08 9.4734e-11 2.8525e-11 5.2943e-05)
+(5.3376e-28 -1.1206e-30 -6.1282e-30 4.9018e-29 -6.2793e-32 5.9575e-28)
+)
diff --git a/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.DFSEM/0/U b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.DFSEM/0/U
new file mode 100644
index 0000000000..2ecd2a6ad4
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.DFSEM/0/U
@@ -0,0 +1,259 @@
+(
+( 0.0000e+00 0 0 )
+( 2.9538e-02 0 0 )
+( 1.1811e-01 0 0 )
+( 2.6562e-01 0 0 )
+( 4.7198e-01 0 0 )
+( 7.3701e-01 0 0 )
+( 1.0605e+00 0 0 )
+( 1.4417e+00 0 0 )
+( 1.8799e+00 0 0 )
+( 2.3729e+00 0 0 )
+( 2.9177e+00 0 0 )
+( 3.5093e+00 0 0 )
+( 4.1409e+00 0 0 )
+( 4.8032e+00 0 0 )
+( 5.4854e+00 0 0 )
+( 6.1754e+00 0 0 )
+( 6.8611e+00 0 0 )
+( 7.5309e+00 0 0 )
+( 8.1754e+00 0 0 )
+( 8.7870e+00 0 0 )
+( 9.3607e+00 0 0 )
+( 9.8937e+00 0 0 )
+( 1.0385e+01 0 0 )
+( 1.0836e+01 0 0 )
+( 1.1248e+01 0 0 )
+( 1.1624e+01 0 0 )
+( 1.1966e+01 0 0 )
+( 1.2278e+01 0 0 )
+( 1.2563e+01 0 0 )
+( 1.2822e+01 0 0 )
+( 1.3060e+01 0 0 )
+( 1.3278e+01 0 0 )
+( 1.3479e+01 0 0 )
+( 1.3664e+01 0 0 )
+( 1.3837e+01 0 0 )
+( 1.3998e+01 0 0 )
+( 1.4148e+01 0 0 )
+( 1.4290e+01 0 0 )
+( 1.4425e+01 0 0 )
+( 1.4552e+01 0 0 )
+( 1.4673e+01 0 0 )
+( 1.4790e+01 0 0 )
+( 1.4902e+01 0 0 )
+( 1.5011e+01 0 0 )
+( 1.5117e+01 0 0 )
+( 1.5221e+01 0 0 )
+( 1.5322e+01 0 0 )
+( 1.5421e+01 0 0 )
+( 1.5518e+01 0 0 )
+( 1.5614e+01 0 0 )
+( 1.5707e+01 0 0 )
+( 1.5799e+01 0 0 )
+( 1.5890e+01 0 0 )
+( 1.5979e+01 0 0 )
+( 1.6067e+01 0 0 )
+( 1.6153e+01 0 0 )
+( 1.6239e+01 0 0 )
+( 1.6324e+01 0 0 )
+( 1.6409e+01 0 0 )
+( 1.6493e+01 0 0 )
+( 1.6576e+01 0 0 )
+( 1.6659e+01 0 0 )
+( 1.6741e+01 0 0 )
+( 1.6823e+01 0 0 )
+( 1.6903e+01 0 0 )
+( 1.6984e+01 0 0 )
+( 1.7063e+01 0 0 )
+( 1.7141e+01 0 0 )
+( 1.7218e+01 0 0 )
+( 1.7294e+01 0 0 )
+( 1.7369e+01 0 0 )
+( 1.7443e+01 0 0 )
+( 1.7517e+01 0 0 )
+( 1.7590e+01 0 0 )
+( 1.7664e+01 0 0 )
+( 1.7738e+01 0 0 )
+( 1.7812e+01 0 0 )
+( 1.7886e+01 0 0 )
+( 1.7960e+01 0 0 )
+( 1.8034e+01 0 0 )
+( 1.8108e+01 0 0 )
+( 1.8182e+01 0 0 )
+( 1.8254e+01 0 0 )
+( 1.8326e+01 0 0 )
+( 1.8396e+01 0 0 )
+( 1.8466e+01 0 0 )
+( 1.8535e+01 0 0 )
+( 1.8603e+01 0 0 )
+( 1.8669e+01 0 0 )
+( 1.8734e+01 0 0 )
+( 1.8797e+01 0 0 )
+( 1.8859e+01 0 0 )
+( 1.8919e+01 0 0 )
+( 1.8978e+01 0 0 )
+( 1.9035e+01 0 0 )
+( 1.9092e+01 0 0 )
+( 1.9148e+01 0 0 )
+( 1.9202e+01 0 0 )
+( 1.9256e+01 0 0 )
+( 1.9308e+01 0 0 )
+( 1.9359e+01 0 0 )
+( 1.9408e+01 0 0 )
+( 1.9456e+01 0 0 )
+( 1.9503e+01 0 0 )
+( 1.9548e+01 0 0 )
+( 1.9593e+01 0 0 )
+( 1.9636e+01 0 0 )
+( 1.9678e+01 0 0 )
+( 1.9719e+01 0 0 )
+( 1.9758e+01 0 0 )
+( 1.9796e+01 0 0 )
+( 1.9832e+01 0 0 )
+( 1.9865e+01 0 0 )
+( 1.9897e+01 0 0 )
+( 1.9927e+01 0 0 )
+( 1.9955e+01 0 0 )
+( 1.9981e+01 0 0 )
+( 2.0004e+01 0 0 )
+( 2.0026e+01 0 0 )
+( 2.0046e+01 0 0 )
+( 2.0064e+01 0 0 )
+( 2.0080e+01 0 0 )
+( 2.0094e+01 0 0 )
+( 2.0106e+01 0 0 )
+( 2.0116e+01 0 0 )
+( 2.0123e+01 0 0 )
+( 2.0129e+01 0 0 )
+( 2.0132e+01 0 0 )
+( 2.0133e+01 0 0 )
+( 2.0132e+01 0 0 )
+( 2.0129e+01 0 0 )
+( 2.0123e+01 0 0 )
+( 2.0116e+01 0 0 )
+( 2.0106e+01 0 0 )
+( 2.0094e+01 0 0 )
+( 2.0080e+01 0 0 )
+( 2.0064e+01 0 0 )
+( 2.0046e+01 0 0 )
+( 2.0026e+01 0 0 )
+( 2.0004e+01 0 0 )
+( 1.9981e+01 0 0 )
+( 1.9955e+01 0 0 )
+( 1.9927e+01 0 0 )
+( 1.9897e+01 0 0 )
+( 1.9865e+01 0 0 )
+( 1.9832e+01 0 0 )
+( 1.9796e+01 0 0 )
+( 1.9758e+01 0 0 )
+( 1.9719e+01 0 0 )
+( 1.9678e+01 0 0 )
+( 1.9636e+01 0 0 )
+( 1.9593e+01 0 0 )
+( 1.9548e+01 0 0 )
+( 1.9503e+01 0 0 )
+( 1.9456e+01 0 0 )
+( 1.9408e+01 0 0 )
+( 1.9359e+01 0 0 )
+( 1.9308e+01 0 0 )
+( 1.9256e+01 0 0 )
+( 1.9202e+01 0 0 )
+( 1.9148e+01 0 0 )
+( 1.9092e+01 0 0 )
+( 1.9035e+01 0 0 )
+( 1.8978e+01 0 0 )
+( 1.8919e+01 0 0 )
+( 1.8859e+01 0 0 )
+( 1.8797e+01 0 0 )
+( 1.8734e+01 0 0 )
+( 1.8669e+01 0 0 )
+( 1.8603e+01 0 0 )
+( 1.8535e+01 0 0 )
+( 1.8466e+01 0 0 )
+( 1.8396e+01 0 0 )
+( 1.8326e+01 0 0 )
+( 1.8254e+01 0 0 )
+( 1.8182e+01 0 0 )
+( 1.8108e+01 0 0 )
+( 1.8034e+01 0 0 )
+( 1.7960e+01 0 0 )
+( 1.7886e+01 0 0 )
+( 1.7812e+01 0 0 )
+( 1.7738e+01 0 0 )
+( 1.7664e+01 0 0 )
+( 1.7590e+01 0 0 )
+( 1.7517e+01 0 0 )
+( 1.7443e+01 0 0 )
+( 1.7369e+01 0 0 )
+( 1.7294e+01 0 0 )
+( 1.7218e+01 0 0 )
+( 1.7141e+01 0 0 )
+( 1.7063e+01 0 0 )
+( 1.6984e+01 0 0 )
+( 1.6903e+01 0 0 )
+( 1.6823e+01 0 0 )
+( 1.6741e+01 0 0 )
+( 1.6659e+01 0 0 )
+( 1.6576e+01 0 0 )
+( 1.6493e+01 0 0 )
+( 1.6409e+01 0 0 )
+( 1.6324e+01 0 0 )
+( 1.6239e+01 0 0 )
+( 1.6153e+01 0 0 )
+( 1.6067e+01 0 0 )
+( 1.5979e+01 0 0 )
+( 1.5890e+01 0 0 )
+( 1.5799e+01 0 0 )
+( 1.5707e+01 0 0 )
+( 1.5614e+01 0 0 )
+( 1.5518e+01 0 0 )
+( 1.5421e+01 0 0 )
+( 1.5322e+01 0 0 )
+( 1.5221e+01 0 0 )
+( 1.5117e+01 0 0 )
+( 1.5011e+01 0 0 )
+( 1.4902e+01 0 0 )
+( 1.4790e+01 0 0 )
+( 1.4673e+01 0 0 )
+( 1.4552e+01 0 0 )
+( 1.4425e+01 0 0 )
+( 1.4290e+01 0 0 )
+( 1.4148e+01 0 0 )
+( 1.3998e+01 0 0 )
+( 1.3837e+01 0 0 )
+( 1.3664e+01 0 0 )
+( 1.3479e+01 0 0 )
+( 1.3278e+01 0 0 )
+( 1.3060e+01 0 0 )
+( 1.2822e+01 0 0 )
+( 1.2563e+01 0 0 )
+( 1.2278e+01 0 0 )
+( 1.1966e+01 0 0 )
+( 1.1624e+01 0 0 )
+( 1.1248e+01 0 0 )
+( 1.0836e+01 0 0 )
+( 1.0385e+01 0 0 )
+( 9.8937e+00 0 0 )
+( 9.3607e+00 0 0 )
+( 8.7870e+00 0 0 )
+( 8.1754e+00 0 0 )
+( 7.5309e+00 0 0 )
+( 6.8611e+00 0 0 )
+( 6.1754e+00 0 0 )
+( 5.4854e+00 0 0 )
+( 4.8032e+00 0 0 )
+( 4.1409e+00 0 0 )
+( 3.5093e+00 0 0 )
+( 2.9177e+00 0 0 )
+( 2.3729e+00 0 0 )
+( 1.8799e+00 0 0 )
+( 1.4417e+00 0 0 )
+( 1.0605e+00 0 0 )
+( 7.3701e-01 0 0 )
+( 4.7198e-01 0 0 )
+( 2.6562e-01 0 0 )
+( 1.1811e-01 0 0 )
+( 2.9538e-02 0 0 )
+( 0.0000e+00 0 0 )
+)
diff --git a/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.DFSEM/points b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.DFSEM/points
new file mode 100644
index 0000000000..5474eb4eae
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.DFSEM/points
@@ -0,0 +1,259 @@
+(
+(0 0.0000e+00 0)
+(0 7.5298e-05 0)
+(0 3.0118e-04 0)
+(0 6.7762e-04 0)
+(0 1.2045e-03 0)
+(0 1.8819e-03 0)
+(0 2.7095e-03 0)
+(0 3.6874e-03 0)
+(0 4.8153e-03 0)
+(0 6.0930e-03 0)
+(0 7.5205e-03 0)
+(0 9.0974e-03 0)
+(0 1.0823e-02 0)
+(0 1.2699e-02 0)
+(0 1.4722e-02 0)
+(0 1.6895e-02 0)
+(0 1.9215e-02 0)
+(0 2.1683e-02 0)
+(0 2.4298e-02 0)
+(0 2.7060e-02 0)
+(0 2.9969e-02 0)
+(0 3.3024e-02 0)
+(0 3.6224e-02 0)
+(0 3.9569e-02 0)
+(0 4.3060e-02 0)
+(0 4.6694e-02 0)
+(0 5.0472e-02 0)
+(0 5.4393e-02 0)
+(0 5.8456e-02 0)
+(0 6.2661e-02 0)
+(0 6.7007e-02 0)
+(0 7.1494e-02 0)
+(0 7.6120e-02 0)
+(0 8.0886e-02 0)
+(0 8.5790e-02 0)
+(0 9.0832e-02 0)
+(0 9.6011e-02 0)
+(0 1.0133e-01 0)
+(0 1.0678e-01 0)
+(0 1.1236e-01 0)
+(0 1.1808e-01 0)
+(0 1.2393e-01 0)
+(0 1.2991e-01 0)
+(0 1.3603e-01 0)
+(0 1.4227e-01 0)
+(0 1.4864e-01 0)
+(0 1.5515e-01 0)
+(0 1.6178e-01 0)
+(0 1.6853e-01 0)
+(0 1.7541e-01 0)
+(0 1.8242e-01 0)
+(0 1.8954e-01 0)
+(0 1.9679e-01 0)
+(0 2.0416e-01 0)
+(0 2.1165e-01 0)
+(0 2.1926e-01 0)
+(0 2.2699e-01 0)
+(0 2.3483e-01 0)
+(0 2.4279e-01 0)
+(0 2.5086e-01 0)
+(0 2.5905e-01 0)
+(0 2.6735e-01 0)
+(0 2.7575e-01 0)
+(0 2.8427e-01 0)
+(0 2.9289e-01 0)
+(0 3.0162e-01 0)
+(0 3.1046e-01 0)
+(0 3.1940e-01 0)
+(0 3.2844e-01 0)
+(0 3.3758e-01 0)
+(0 3.4683e-01 0)
+(0 3.5617e-01 0)
+(0 3.6561e-01 0)
+(0 3.7514e-01 0)
+(0 3.8477e-01 0)
+(0 3.9449e-01 0)
+(0 4.0430e-01 0)
+(0 4.1420e-01 0)
+(0 4.2419e-01 0)
+(0 4.3427e-01 0)
+(0 4.4443e-01 0)
+(0 4.5467e-01 0)
+(0 4.6500e-01 0)
+(0 4.7541e-01 0)
+(0 4.8590e-01 0)
+(0 4.9646e-01 0)
+(0 5.0710e-01 0)
+(0 5.1782e-01 0)
+(0 5.2860e-01 0)
+(0 5.3946e-01 0)
+(0 5.5039e-01 0)
+(0 5.6138e-01 0)
+(0 5.7244e-01 0)
+(0 5.8357e-01 0)
+(0 5.9476e-01 0)
+(0 6.0601e-01 0)
+(0 6.1732e-01 0)
+(0 6.2868e-01 0)
+(0 6.4010e-01 0)
+(0 6.5158e-01 0)
+(0 6.6311e-01 0)
+(0 6.7469e-01 0)
+(0 6.8632e-01 0)
+(0 6.9799e-01 0)
+(0 7.0972e-01 0)
+(0 7.2148e-01 0)
+(0 7.3329e-01 0)
+(0 7.4513e-01 0)
+(0 7.5702e-01 0)
+(0 7.6894e-01 0)
+(0 7.8090e-01 0)
+(0 7.9289e-01 0)
+(0 8.0491e-01 0)
+(0 8.1696e-01 0)
+(0 8.2904e-01 0)
+(0 8.4114e-01 0)
+(0 8.5327e-01 0)
+(0 8.6542e-01 0)
+(0 8.7759e-01 0)
+(0 8.8978e-01 0)
+(0 9.0198e-01 0)
+(0 9.1420e-01 0)
+(0 9.2644e-01 0)
+(0 9.3868e-01 0)
+(0 9.5093e-01 0)
+(0 9.6319e-01 0)
+(0 9.7546e-01 0)
+(0 9.8773e-01 0)
+(0 1.0000e-00 0)
+(0 1.01227 0)
+(0 1.02454 0)
+(0 1.03681 0)
+(0 1.04907 0)
+(0 1.06132 0)
+(0 1.07356 0)
+(0 1.0858 0)
+(0 1.09802 0)
+(0 1.11022 0)
+(0 1.12241 0)
+(0 1.13458 0)
+(0 1.14673 0)
+(0 1.15886 0)
+(0 1.17096 0)
+(0 1.18304 0)
+(0 1.19509 0)
+(0 1.20711 0)
+(0 1.2191 0)
+(0 1.23106 0)
+(0 1.24298 0)
+(0 1.25487 0)
+(0 1.26671 0)
+(0 1.27852 0)
+(0 1.29028 0)
+(0 1.30201 0)
+(0 1.31368 0)
+(0 1.32531 0)
+(0 1.33689 0)
+(0 1.34842 0)
+(0 1.3599 0)
+(0 1.37132 0)
+(0 1.38268 0)
+(0 1.39399 0)
+(0 1.40524 0)
+(0 1.41643 0)
+(0 1.42756 0)
+(0 1.43862 0)
+(0 1.44961 0)
+(0 1.46054 0)
+(0 1.4714 0)
+(0 1.48218 0)
+(0 1.4929 0)
+(0 1.50354 0)
+(0 1.5141 0)
+(0 1.52459 0)
+(0 1.535 0)
+(0 1.54533 0)
+(0 1.55557 0)
+(0 1.56573 0)
+(0 1.57581 0)
+(0 1.5858 0)
+(0 1.5957 0)
+(0 1.60551 0)
+(0 1.61523 0)
+(0 1.62486 0)
+(0 1.63439 0)
+(0 1.64383 0)
+(0 1.65317 0)
+(0 1.66242 0)
+(0 1.67156 0)
+(0 1.6806 0)
+(0 1.68954 0)
+(0 1.69838 0)
+(0 1.70711 0)
+(0 1.71573 0)
+(0 1.72425 0)
+(0 1.73265 0)
+(0 1.74095 0)
+(0 1.74914 0)
+(0 1.75721 0)
+(0 1.76517 0)
+(0 1.77301 0)
+(0 1.78074 0)
+(0 1.78835 0)
+(0 1.79584 0)
+(0 1.80321 0)
+(0 1.81046 0)
+(0 1.81758 0)
+(0 1.82459 0)
+(0 1.83147 0)
+(0 1.83822 0)
+(0 1.84485 0)
+(0 1.85136 0)
+(0 1.85773 0)
+(0 1.86397 0)
+(0 1.87009 0)
+(0 1.87607 0)
+(0 1.88192 0)
+(0 1.88764 0)
+(0 1.89322 0)
+(0 1.89867 0)
+(0 1.90399 0)
+(0 1.90917 0)
+(0 1.91421 0)
+(0 1.91911 0)
+(0 1.92388 0)
+(0 1.92851 0)
+(0 1.93299 0)
+(0 1.93734 0)
+(0 1.94154 0)
+(0 1.94561 0)
+(0 1.94953 0)
+(0 1.95331 0)
+(0 1.95694 0)
+(0 1.96043 0)
+(0 1.96378 0)
+(0 1.96698 0)
+(0 1.97003 0)
+(0 1.97294 0)
+(0 1.9757 0)
+(0 1.97832 0)
+(0 1.98079 0)
+(0 1.9831 0)
+(0 1.98528 0)
+(0 1.9873 0)
+(0 1.98918 0)
+(0 1.9909 0)
+(0 1.99248 0)
+(0 1.99391 0)
+(0 1.99518 0)
+(0 1.99631 0)
+(0 1.99729 0)
+(0 1.99812 0)
+(0 1.9988 0)
+(0 1.99932 0)
+(0 1.9997 0)
+(0 1.99992 0)
+(0 2 0)
+)
diff --git a/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.digitalFilter/0/R b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.digitalFilter/0/R
new file mode 100644
index 0000000000..084bdd9f8f
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.digitalFilter/0/R
@@ -0,0 +1,259 @@
+(
+(5.3376e-28 1.1206e-30 -6.1282e-30 4.9018e-29 6.2793e-32 5.9575e-28)
+(1.3684e-04 -2.4414e-08 5.9464e-08 9.4734e-11 -2.8525e-11 5.2943e-05)
+(2.1828e-03 -1.5728e-06 8.3228e-07 2.3603e-08 -1.7141e-09 8.2715e-04)
+(1.1020e-02 -1.8113e-05 3.3616e-06 5.7878e-07 -1.9051e-08 4.0392e-03)
+(3.4718e-02 -1.0324e-04 1.0732e-05 5.4400e-06 -1.0411e-07 1.2147e-02)
+(8.4454e-02 -4.0027e-04 2.8297e-05 3.0014e-05 -3.8272e-07 2.7852e-02)
+(1.7438e-01 -1.2143e-03 6.6865e-05 1.1756e-04 -1.0969e-06 5.3564e-02)
+(3.2123e-01 -3.1009e-03 1.4506e-04 3.6204e-04 -2.6288e-06 9.0947e-02)
+(5.4336e-01 -6.9486e-03 2.9097e-04 9.3203e-04 -5.4596e-06 1.4062e-01)
+(8.5851e-01 -1.4015e-02 5.4723e-04 2.0879e-03 -1.0091e-05 2.0209e-01)
+(1.2802e+00 -2.5858e-02 9.7373e-04 4.1852e-03 -1.7119e-05 2.7384e-01)
+(1.8132e+00 -4.4135e-02 1.6529e-03 7.6610e-03 -2.7764e-05 3.5369e-01)
+(2.4490e+00 -7.0276e-02 2.6847e-03 1.3005e-02 -4.4587e-05 4.3919e-01)
+(3.1644e+00 -1.0513e-01 4.1721e-03 2.0724e-02 -7.1044e-05 5.2804e-01)
+(3.9229e+00 -1.4869e-01 6.1889e-03 3.1300e-02 -1.0893e-04 6.1834e-01)
+(4.6811e+00 -1.9999e-01 8.7562e-03 4.5157e-02 -1.5460e-04 7.0871e-01)
+(5.3953e+00 -2.5724e-01 1.1834e-02 6.2636e-02 -1.9552e-04 7.9819e-01)
+(6.0286e+00 -3.1808e-01 1.5387e-02 8.3971e-02 -2.1029e-04 8.8610e-01)
+(6.5558e+00 -3.8001e-01 1.9449e-02 1.0928e-01 -1.7262e-04 9.7178e-01)
+(6.9642e+00 -4.4073e-01 2.4133e-02 1.3854e-01 -6.0490e-05 1.0545e+00)
+(7.2524e+00 -4.9837e-01 2.9501e-02 1.7162e-01 1.3268e-04 1.1333e+00)
+(7.4274e+00 -5.5162e-01 3.5418e-02 2.0824e-01 3.9086e-04 1.2075e+00)
+(7.5016e+00 -5.9969e-01 4.1511e-02 2.4803e-01 6.7818e-04 1.2762e+00)
+(7.4897e+00 -6.4224e-01 4.7330e-02 2.9050e-01 9.5511e-04 1.3387e+00)
+(7.4069e+00 -6.7924e-01 5.2512e-02 3.3510e-01 1.1989e-03 1.3948e+00)
+(7.2679e+00 -7.1090e-01 5.6926e-02 3.8121e-01 1.4016e-03 1.4444e+00)
+(7.0865e+00 -7.3759e-01 6.0659e-02 4.2818e-01 1.5532e-03 1.4878e+00)
+(6.8750e+00 -7.5980e-01 6.3835e-02 4.7535e-01 1.6360e-03 1.5251e+00)
+(6.6441e+00 -7.7810e-01 6.6418e-02 5.2211e-01 1.6465e-03 1.5568e+00)
+(6.4025e+00 -7.9302e-01 6.8207e-02 5.6790e-01 1.6026e-03 1.5833e+00)
+(6.1573e+00 -8.0503e-01 6.9001e-02 6.1229e-01 1.5073e-03 1.6047e+00)
+(5.9141e+00 -8.1452e-01 6.8896e-02 6.5493e-01 1.3268e-03 1.6216e+00)
+(5.6767e+00 -8.2184e-01 6.8285e-02 6.9554e-01 1.0303e-03 1.6344e+00)
+(5.4477e+00 -8.2733e-01 6.7481e-02 7.3388e-01 6.4119e-04 1.6439e+00)
+(5.2291e+00 -8.3133e-01 6.6450e-02 7.6971e-01 2.3642e-04 1.6512e+00)
+(5.0223e+00 -8.3412e-01 6.4937e-02 8.0285e-01 -1.3612e-04 1.6567e+00)
+(4.8279e+00 -8.3590e-01 6.2875e-02 8.3316e-01 -4.7569e-04 1.6606e+00)
+(4.6462e+00 -8.3683e-01 6.0710e-02 8.6061e-01 -7.8412e-04 1.6623e+00)
+(4.4768e+00 -8.3700e-01 5.9182e-02 8.8520e-01 -1.0539e-03 1.6618e+00)
+(4.3194e+00 -8.3649e-01 5.8527e-02 9.0700e-01 -1.2655e-03 1.6594e+00)
+(4.1730e+00 -8.3523e-01 5.8437e-02 9.2607e-01 -1.4126e-03 1.6553e+00)
+(4.0370e+00 -8.3321e-01 5.8512e-02 9.4247e-01 -1.4930e-03 1.6500e+00)
+(3.9106e+00 -8.3044e-01 5.8376e-02 9.5626e-01 -1.5118e-03 1.6439e+00)
+(3.7930e+00 -8.2687e-01 5.8084e-02 9.6754e-01 -1.5196e-03 1.6373e+00)
+(3.6834e+00 -8.2252e-01 5.8125e-02 9.7651e-01 -1.6228e-03 1.6299e+00)
+(3.5821e+00 -8.1762e-01 5.8537e-02 9.8341e-01 -1.8295e-03 1.6209e+00)
+(3.4896e+00 -8.1249e-01 5.8809e-02 9.8850e-01 -2.0145e-03 1.6105e+00)
+(3.4053e+00 -8.0726e-01 5.8494e-02 9.9197e-01 -2.0826e-03 1.5991e+00)
+(3.3278e+00 -8.0180e-01 5.7876e-02 9.9396e-01 -2.1101e-03 1.5872e+00)
+(3.2557e+00 -7.9592e-01 5.7034e-02 9.9455e-01 -2.1628e-03 1.5754e+00)
+(3.1886e+00 -7.8961e-01 5.5746e-02 9.9379e-01 -2.2441e-03 1.5637e+00)
+(3.1269e+00 -7.8296e-01 5.4259e-02 9.9183e-01 -2.3220e-03 1.5510e+00)
+(3.0702e+00 -7.7613e-01 5.2957e-02 9.8887e-01 -2.3519e-03 1.5367e+00)
+(3.0178e+00 -7.6929e-01 5.1819e-02 9.8512e-01 -2.3018e-03 1.5212e+00)
+(2.9678e+00 -7.6243e-01 5.0872e-02 9.8075e-01 -2.0851e-03 1.5052e+00)
+(2.9191e+00 -7.5552e-01 5.0514e-02 9.7583e-01 -1.6209e-03 1.4889e+00)
+(2.8712e+00 -7.4851e-01 5.0564e-02 9.7044e-01 -8.5802e-04 1.4721e+00)
+(2.8239e+00 -7.4146e-01 5.0229e-02 9.6465e-01 1.1783e-04 1.4551e+00)
+(2.7776e+00 -7.3445e-01 4.9171e-02 9.5846e-01 1.0337e-03 1.4379e+00)
+(2.7324e+00 -7.2737e-01 4.7565e-02 9.5179e-01 1.6440e-03 1.4206e+00)
+(2.6883e+00 -7.2012e-01 4.5716e-02 9.4459e-01 1.9851e-03 1.4035e+00)
+(2.6453e+00 -7.1271e-01 4.3916e-02 9.3679e-01 2.2291e-03 1.3867e+00)
+(2.6031e+00 -7.0521e-01 4.2134e-02 9.2845e-01 2.4143e-03 1.3696e+00)
+(2.5612e+00 -6.9761e-01 4.0077e-02 9.1976e-01 2.5128e-03 1.3515e+00)
+(2.5196e+00 -6.8985e-01 3.7520e-02 9.1087e-01 2.5006e-03 1.3322e+00)
+(2.4793e+00 -6.8192e-01 3.4358e-02 9.0185e-01 2.3726e-03 1.3118e+00)
+(2.4411e+00 -6.7396e-01 3.0858e-02 8.9268e-01 2.1952e-03 1.2910e+00)
+(2.4055e+00 -6.6598e-01 2.7094e-02 8.8335e-01 2.0517e-03 1.2700e+00)
+(2.3714e+00 -6.5794e-01 2.3371e-02 8.7385e-01 1.9999e-03 1.2488e+00)
+(2.3380e+00 -6.4977e-01 2.0500e-02 8.6418e-01 1.9700e-03 1.2276e+00)
+(2.3046e+00 -6.4148e-01 1.8795e-02 8.5436e-01 1.8189e-03 1.2066e+00)
+(2.2708e+00 -6.3298e-01 1.7891e-02 8.4439e-01 1.4336e-03 1.1855e+00)
+(2.2365e+00 -6.2415e-01 1.7572e-02 8.3421e-01 7.5375e-04 1.1641e+00)
+(2.2023e+00 -6.1500e-01 1.7517e-02 8.2369e-01 -7.9290e-05 1.1423e+00)
+(2.1682e+00 -6.0551e-01 1.7095e-02 8.1285e-01 -8.3720e-04 1.1202e+00)
+(2.1335e+00 -5.9561e-01 1.5913e-02 8.0181e-01 -1.3413e-03 1.0983e+00)
+(2.0983e+00 -5.8530e-01 1.4306e-02 7.9074e-01 -1.5771e-03 1.0767e+00)
+(2.0631e+00 -5.7460e-01 1.2812e-02 7.7968e-01 -1.5179e-03 1.0551e+00)
+(2.0277e+00 -5.6350e-01 1.1907e-02 7.6856e-01 -1.1762e-03 1.0339e+00)
+(1.9920e+00 -5.5204e-01 1.1640e-02 7.5731e-01 -7.0498e-04 1.0131e+00)
+(1.9568e+00 -5.4040e-01 1.1503e-02 7.4588e-01 -2.0689e-04 9.9279e-01)
+(1.9229e+00 -5.2883e-01 1.1097e-02 7.3425e-01 2.0529e-04 9.7254e-01)
+(1.8896e+00 -5.1747e-01 1.0218e-02 7.2249e-01 5.0971e-04 9.5187e-01)
+(1.8562e+00 -5.0626e-01 9.0004e-03 7.1064e-01 7.6211e-04 9.3057e-01)
+(1.8232e+00 -4.9507e-01 7.8753e-03 6.9862e-01 9.2292e-04 9.0936e-01)
+(1.7908e+00 -4.8386e-01 7.1617e-03 6.8642e-01 9.3017e-04 8.8926e-01)
+(1.7591e+00 -4.7269e-01 7.0121e-03 6.7417e-01 7.2814e-04 8.7042e-01)
+(1.7274e+00 -4.6164e-01 7.3802e-03 6.6201e-01 2.7352e-04 8.5243e-01)
+(1.6948e+00 -4.5072e-01 8.2299e-03 6.5007e-01 -2.9303e-04 8.3469e-01)
+(1.6610e+00 -4.3999e-01 9.3954e-03 6.3843e-01 -7.4565e-04 8.1690e-01)
+(1.6264e+00 -4.2949e-01 1.0585e-02 6.2712e-01 -9.5384e-04 7.9939e-01)
+(1.5910e+00 -4.1910e-01 1.1834e-02 6.1613e-01 -9.5389e-04 7.8241e-01)
+(1.5550e+00 -4.0876e-01 1.3728e-02 6.0545e-01 -1.0075e-03 7.6587e-01)
+(1.5184e+00 -3.9853e-01 1.6563e-02 5.9512e-01 -1.2970e-03 7.4931e-01)
+(1.4808e+00 -3.8843e-01 1.9835e-02 5.8523e-01 -1.7620e-03 7.3258e-01)
+(1.4425e+00 -3.7835e-01 2.2555e-02 5.7589e-01 -2.2305e-03 7.1572e-01)
+(1.4037e+00 -3.6810e-01 2.3965e-02 5.6715e-01 -2.5637e-03 6.9868e-01)
+(1.3647e+00 -3.5761e-01 2.4099e-02 5.5896e-01 -2.7183e-03 6.8159e-01)
+(1.3260e+00 -3.4695e-01 2.3378e-02 5.5120e-01 -2.7095e-03 6.6492e-01)
+(1.2881e+00 -3.3621e-01 2.2236e-02 5.4372e-01 -2.5722e-03 6.4916e-01)
+(1.2513e+00 -3.2544e-01 2.0871e-02 5.3642e-01 -2.4646e-03 6.3424e-01)
+(1.2162e+00 -3.1468e-01 1.9224e-02 5.2919e-01 -2.5472e-03 6.2001e-01)
+(1.1827e+00 -3.0393e-01 1.7273e-02 5.2197e-01 -2.7871e-03 6.0660e-01)
+(1.1500e+00 -2.9315e-01 1.5116e-02 5.1473e-01 -3.0094e-03 5.9377e-01)
+(1.1175e+00 -2.8231e-01 1.3055e-02 5.0750e-01 -3.1202e-03 5.8101e-01)
+(1.0843e+00 -2.7133e-01 1.1360e-02 5.0029e-01 -3.1772e-03 5.6815e-01)
+(1.0507e+00 -2.6019e-01 9.9787e-03 4.9312e-01 -3.2410e-03 5.5541e-01)
+(1.0173e+00 -2.4893e-01 8.7677e-03 4.8605e-01 -3.3171e-03 5.4304e-01)
+(9.8478e-01 -2.3753e-01 7.8868e-03 4.7922e-01 -3.3738e-03 5.3113e-01)
+(9.5377e-01 -2.2605e-01 7.6156e-03 4.7268e-01 -3.4135e-03 5.1971e-01)
+(9.2454e-01 -2.1462e-01 7.9424e-03 4.6644e-01 -3.4897e-03 5.0882e-01)
+(8.9706e-01 -2.0326e-01 8.6379e-03 4.6047e-01 -3.6720e-03 4.9853e-01)
+(8.7113e-01 -1.9188e-01 9.5840e-03 4.5478e-01 -3.9882e-03 4.8883e-01)
+(8.4645e-01 -1.8043e-01 1.0540e-02 4.4932e-01 -4.3896e-03 4.7982e-01)
+(8.2296e-01 -1.6898e-01 1.1159e-02 4.4410e-01 -4.7685e-03 4.7164e-01)
+(8.0062e-01 -1.5752e-01 1.1209e-02 4.3915e-01 -5.0417e-03 4.6392e-01)
+(7.7929e-01 -1.4601e-01 1.0617e-02 4.3458e-01 -5.1687e-03 4.5618e-01)
+(7.5889e-01 -1.3442e-01 9.5260e-03 4.3044e-01 -5.1040e-03 4.4842e-01)
+(7.3932e-01 -1.2272e-01 8.2998e-03 4.2671e-01 -4.8287e-03 4.4071e-01)
+(7.2051e-01 -1.1094e-01 7.2287e-03 4.2339e-01 -4.4056e-03 4.3345e-01)
+(7.0261e-01 -9.9120e-02 6.4502e-03 4.2046e-01 -3.8984e-03 4.2708e-01)
+(6.8576e-01 -8.7301e-02 6.1696e-03 4.1796e-01 -3.3462e-03 4.2169e-01)
+(6.7023e-01 -7.5441e-02 6.4277e-03 4.1591e-01 -2.7927e-03 4.1709e-01)
+(6.5620e-01 -6.3452e-02 7.0238e-03 4.1431e-01 -2.2376e-03 4.1319e-01)
+(6.4389e-01 -5.1252e-02 7.7818e-03 4.1313e-01 -1.6855e-03 4.0984e-01)
+(6.3375e-01 -3.8772e-02 8.7489e-03 4.1234e-01 -1.1858e-03 4.0689e-01)
+(6.2626e-01 -2.6001e-02 9.8627e-03 4.1188e-01 -7.6554e-04 4.0431e-01)
+(6.2174e-01 -1.3031e-02 1.0731e-02 4.1166e-01 -3.8835e-04 4.0247e-01)
+(6.2024e-01 0.0000e+00 1.1032e-02 4.1160e-01 0.0000e+00 4.0179e-01)
+(6.2174e-01 1.3031e-02 1.0731e-02 4.1166e-01 3.8835e-04 4.0247e-01)
+(6.2626e-01 2.6001e-02 9.8627e-03 4.1188e-01 7.6554e-04 4.0431e-01)
+(6.3375e-01 3.8772e-02 8.7489e-03 4.1234e-01 1.1858e-03 4.0689e-01)
+(6.4389e-01 5.1252e-02 7.7818e-03 4.1313e-01 1.6855e-03 4.0984e-01)
+(6.5620e-01 6.3452e-02 7.0238e-03 4.1431e-01 2.2376e-03 4.1319e-01)
+(6.7023e-01 7.5441e-02 6.4277e-03 4.1591e-01 2.7927e-03 4.1709e-01)
+(6.8576e-01 8.7301e-02 6.1696e-03 4.1796e-01 3.3462e-03 4.2169e-01)
+(7.0261e-01 9.9120e-02 6.4502e-03 4.2046e-01 3.8984e-03 4.2708e-01)
+(7.2051e-01 1.1094e-01 7.2287e-03 4.2339e-01 4.4056e-03 4.3345e-01)
+(7.3932e-01 1.2272e-01 8.2998e-03 4.2671e-01 4.8287e-03 4.4071e-01)
+(7.5889e-01 1.3442e-01 9.5260e-03 4.3044e-01 5.1040e-03 4.4842e-01)
+(7.7929e-01 1.4601e-01 1.0617e-02 4.3458e-01 5.1687e-03 4.5618e-01)
+(8.0062e-01 1.5752e-01 1.1209e-02 4.3915e-01 5.0417e-03 4.6392e-01)
+(8.2296e-01 1.6898e-01 1.1159e-02 4.4410e-01 4.7685e-03 4.7164e-01)
+(8.4645e-01 1.8043e-01 1.0540e-02 4.4932e-01 4.3896e-03 4.7982e-01)
+(8.7113e-01 1.9188e-01 9.5840e-03 4.5478e-01 3.9882e-03 4.8883e-01)
+(8.9706e-01 2.0326e-01 8.6379e-03 4.6047e-01 3.6720e-03 4.9853e-01)
+(9.2454e-01 2.1462e-01 7.9424e-03 4.6644e-01 3.4897e-03 5.0882e-01)
+(9.5377e-01 2.2605e-01 7.6156e-03 4.7268e-01 3.4135e-03 5.1971e-01)
+(9.8478e-01 2.3753e-01 7.8868e-03 4.7922e-01 3.3738e-03 5.3113e-01)
+(1.0173e+00 2.4893e-01 8.7677e-03 4.8605e-01 3.3171e-03 5.4304e-01)
+(1.0507e+00 2.6019e-01 9.9787e-03 4.9312e-01 3.2410e-03 5.5541e-01)
+(1.0843e+00 2.7133e-01 1.1360e-02 5.0029e-01 3.1772e-03 5.6815e-01)
+(1.1175e+00 2.8231e-01 1.3055e-02 5.0750e-01 3.1202e-03 5.8101e-01)
+(1.1500e+00 2.9315e-01 1.5116e-02 5.1473e-01 3.0094e-03 5.9377e-01)
+(1.1827e+00 3.0393e-01 1.7273e-02 5.2197e-01 2.7871e-03 6.0660e-01)
+(1.2162e+00 3.1468e-01 1.9224e-02 5.2919e-01 2.5472e-03 6.2001e-01)
+(1.2513e+00 3.2544e-01 2.0871e-02 5.3642e-01 2.4646e-03 6.3424e-01)
+(1.2881e+00 3.3621e-01 2.2236e-02 5.4372e-01 2.5722e-03 6.4916e-01)
+(1.3260e+00 3.4695e-01 2.3378e-02 5.5120e-01 2.7095e-03 6.6492e-01)
+(1.3647e+00 3.5761e-01 2.4099e-02 5.5896e-01 2.7183e-03 6.8159e-01)
+(1.4037e+00 3.6810e-01 2.3965e-02 5.6715e-01 2.5637e-03 6.9868e-01)
+(1.4425e+00 3.7835e-01 2.2555e-02 5.7589e-01 2.2305e-03 7.1572e-01)
+(1.4808e+00 3.8843e-01 1.9835e-02 5.8523e-01 1.7620e-03 7.3258e-01)
+(1.5184e+00 3.9853e-01 1.6563e-02 5.9512e-01 1.2970e-03 7.4931e-01)
+(1.5550e+00 4.0876e-01 1.3728e-02 6.0545e-01 1.0075e-03 7.6587e-01)
+(1.5910e+00 4.1910e-01 1.1834e-02 6.1613e-01 9.5389e-04 7.8241e-01)
+(1.6264e+00 4.2949e-01 1.0585e-02 6.2712e-01 9.5384e-04 7.9939e-01)
+(1.6610e+00 4.3999e-01 9.3954e-03 6.3843e-01 7.4565e-04 8.1690e-01)
+(1.6948e+00 4.5072e-01 8.2299e-03 6.5007e-01 2.9303e-04 8.3469e-01)
+(1.7274e+00 4.6164e-01 7.3802e-03 6.6201e-01 -2.7352e-04 8.5243e-01)
+(1.7591e+00 4.7269e-01 7.0121e-03 6.7417e-01 -7.2814e-04 8.7042e-01)
+(1.7908e+00 4.8386e-01 7.1617e-03 6.8642e-01 -9.3017e-04 8.8926e-01)
+(1.8232e+00 4.9507e-01 7.8753e-03 6.9862e-01 -9.2292e-04 9.0936e-01)
+(1.8562e+00 5.0626e-01 9.0004e-03 7.1064e-01 -7.6211e-04 9.3057e-01)
+(1.8896e+00 5.1747e-01 1.0218e-02 7.2249e-01 -5.0971e-04 9.5187e-01)
+(1.9229e+00 5.2883e-01 1.1097e-02 7.3425e-01 -2.0529e-04 9.7254e-01)
+(1.9568e+00 5.4040e-01 1.1503e-02 7.4588e-01 2.0689e-04 9.9279e-01)
+(1.9920e+00 5.5204e-01 1.1640e-02 7.5731e-01 7.0498e-04 1.0131e+00)
+(2.0277e+00 5.6350e-01 1.1907e-02 7.6856e-01 1.1762e-03 1.0339e+00)
+(2.0631e+00 5.7460e-01 1.2812e-02 7.7968e-01 1.5179e-03 1.0551e+00)
+(2.0983e+00 5.8530e-01 1.4306e-02 7.9074e-01 1.5771e-03 1.0767e+00)
+(2.1335e+00 5.9561e-01 1.5913e-02 8.0181e-01 1.3413e-03 1.0983e+00)
+(2.1682e+00 6.0551e-01 1.7095e-02 8.1285e-01 8.3720e-04 1.1202e+00)
+(2.2023e+00 6.1500e-01 1.7517e-02 8.2369e-01 7.9290e-05 1.1423e+00)
+(2.2365e+00 6.2415e-01 1.7572e-02 8.3421e-01 -7.5375e-04 1.1641e+00)
+(2.2708e+00 6.3298e-01 1.7891e-02 8.4439e-01 -1.4336e-03 1.1855e+00)
+(2.3046e+00 6.4148e-01 1.8795e-02 8.5436e-01 -1.8189e-03 1.2066e+00)
+(2.3380e+00 6.4977e-01 2.0500e-02 8.6418e-01 -1.9700e-03 1.2276e+00)
+(2.3714e+00 6.5794e-01 2.3371e-02 8.7385e-01 -1.9999e-03 1.2488e+00)
+(2.4055e+00 6.6598e-01 2.7094e-02 8.8335e-01 -2.0517e-03 1.2700e+00)
+(2.4411e+00 6.7396e-01 3.0858e-02 8.9268e-01 -2.1952e-03 1.2910e+00)
+(2.4793e+00 6.8192e-01 3.4358e-02 9.0185e-01 -2.3726e-03 1.3118e+00)
+(2.5196e+00 6.8985e-01 3.7520e-02 9.1087e-01 -2.5006e-03 1.3322e+00)
+(2.5612e+00 6.9761e-01 4.0077e-02 9.1976e-01 -2.5128e-03 1.3515e+00)
+(2.6031e+00 7.0521e-01 4.2134e-02 9.2845e-01 -2.4143e-03 1.3696e+00)
+(2.6453e+00 7.1271e-01 4.3916e-02 9.3679e-01 -2.2291e-03 1.3867e+00)
+(2.6883e+00 7.2012e-01 4.5716e-02 9.4459e-01 -1.9851e-03 1.4035e+00)
+(2.7324e+00 7.2737e-01 4.7565e-02 9.5179e-01 -1.6440e-03 1.4206e+00)
+(2.7776e+00 7.3445e-01 4.9171e-02 9.5846e-01 -1.0337e-03 1.4379e+00)
+(2.8239e+00 7.4146e-01 5.0229e-02 9.6465e-01 -1.1783e-04 1.4551e+00)
+(2.8712e+00 7.4851e-01 5.0564e-02 9.7044e-01 8.5802e-04 1.4721e+00)
+(2.9191e+00 7.5552e-01 5.0514e-02 9.7583e-01 1.6209e-03 1.4889e+00)
+(2.9678e+00 7.6243e-01 5.0872e-02 9.8075e-01 2.0851e-03 1.5052e+00)
+(3.0178e+00 7.6929e-01 5.1819e-02 9.8512e-01 2.3018e-03 1.5212e+00)
+(3.0702e+00 7.7613e-01 5.2957e-02 9.8887e-01 2.3519e-03 1.5367e+00)
+(3.1269e+00 7.8296e-01 5.4259e-02 9.9183e-01 2.3220e-03 1.5510e+00)
+(3.1886e+00 7.8961e-01 5.5746e-02 9.9379e-01 2.2441e-03 1.5637e+00)
+(3.2557e+00 7.9592e-01 5.7034e-02 9.9455e-01 2.1628e-03 1.5754e+00)
+(3.3278e+00 8.0180e-01 5.7876e-02 9.9396e-01 2.1101e-03 1.5872e+00)
+(3.4053e+00 8.0726e-01 5.8494e-02 9.9197e-01 2.0826e-03 1.5991e+00)
+(3.4896e+00 8.1249e-01 5.8809e-02 9.8850e-01 2.0145e-03 1.6105e+00)
+(3.5821e+00 8.1762e-01 5.8537e-02 9.8341e-01 1.8295e-03 1.6209e+00)
+(3.6834e+00 8.2252e-01 5.8125e-02 9.7651e-01 1.6228e-03 1.6299e+00)
+(3.7930e+00 8.2687e-01 5.8084e-02 9.6754e-01 1.5196e-03 1.6373e+00)
+(3.9106e+00 8.3044e-01 5.8376e-02 9.5626e-01 1.5118e-03 1.6439e+00)
+(4.0370e+00 8.3321e-01 5.8512e-02 9.4247e-01 1.4930e-03 1.6500e+00)
+(4.1730e+00 8.3523e-01 5.8437e-02 9.2607e-01 1.4126e-03 1.6553e+00)
+(4.3194e+00 8.3649e-01 5.8527e-02 9.0700e-01 1.2655e-03 1.6594e+00)
+(4.4768e+00 8.3700e-01 5.9182e-02 8.8520e-01 1.0539e-03 1.6618e+00)
+(4.6462e+00 8.3683e-01 6.0710e-02 8.6061e-01 7.8412e-04 1.6623e+00)
+(4.8279e+00 8.3590e-01 6.2875e-02 8.3316e-01 4.7569e-04 1.6606e+00)
+(5.0223e+00 8.3412e-01 6.4937e-02 8.0285e-01 1.3612e-04 1.6567e+00)
+(5.2291e+00 8.3133e-01 6.6450e-02 7.6971e-01 -2.3642e-04 1.6512e+00)
+(5.4477e+00 8.2733e-01 6.7481e-02 7.3388e-01 -6.4119e-04 1.6439e+00)
+(5.6767e+00 8.2184e-01 6.8285e-02 6.9554e-01 -1.0303e-03 1.6344e+00)
+(5.9141e+00 8.1452e-01 6.8896e-02 6.5493e-01 -1.3268e-03 1.6216e+00)
+(6.1573e+00 8.0503e-01 6.9001e-02 6.1229e-01 -1.5073e-03 1.6047e+00)
+(6.4025e+00 7.9302e-01 6.8207e-02 5.6790e-01 -1.6026e-03 1.5833e+00)
+(6.6441e+00 7.7810e-01 6.6418e-02 5.2211e-01 -1.6465e-03 1.5568e+00)
+(6.8750e+00 7.5980e-01 6.3835e-02 4.7535e-01 -1.6360e-03 1.5251e+00)
+(7.0865e+00 7.3759e-01 6.0659e-02 4.2818e-01 -1.5532e-03 1.4878e+00)
+(7.2679e+00 7.1090e-01 5.6926e-02 3.8121e-01 -1.4016e-03 1.4444e+00)
+(7.4069e+00 6.7924e-01 5.2512e-02 3.3510e-01 -1.1989e-03 1.3948e+00)
+(7.4897e+00 6.4224e-01 4.7330e-02 2.9050e-01 -9.5511e-04 1.3387e+00)
+(7.5016e+00 5.9969e-01 4.1511e-02 2.4803e-01 -6.7818e-04 1.2762e+00)
+(7.4274e+00 5.5162e-01 3.5418e-02 2.0824e-01 -3.9086e-04 1.2075e+00)
+(7.2524e+00 4.9837e-01 2.9501e-02 1.7162e-01 -1.3268e-04 1.1333e+00)
+(6.9642e+00 4.4073e-01 2.4133e-02 1.3854e-01 6.0490e-05 1.0545e+00)
+(6.5558e+00 3.8001e-01 1.9449e-02 1.0928e-01 1.7262e-04 9.7178e-01)
+(6.0286e+00 3.1808e-01 1.5387e-02 8.3971e-02 2.1029e-04 8.8610e-01)
+(5.3953e+00 2.5724e-01 1.1834e-02 6.2636e-02 1.9552e-04 7.9819e-01)
+(4.6811e+00 1.9999e-01 8.7562e-03 4.5157e-02 1.5460e-04 7.0871e-01)
+(3.9229e+00 1.4869e-01 6.1889e-03 3.1300e-02 1.0893e-04 6.1834e-01)
+(3.1644e+00 1.0513e-01 4.1721e-03 2.0724e-02 7.1044e-05 5.2804e-01)
+(2.4490e+00 7.0276e-02 2.6847e-03 1.3005e-02 4.4587e-05 4.3919e-01)
+(1.8132e+00 4.4135e-02 1.6529e-03 7.6610e-03 2.7764e-05 3.5369e-01)
+(1.2802e+00 2.5858e-02 9.7373e-04 4.1852e-03 1.7119e-05 2.7384e-01)
+(8.5851e-01 1.4015e-02 5.4723e-04 2.0879e-03 1.0091e-05 2.0209e-01)
+(5.4336e-01 6.9486e-03 2.9097e-04 9.3203e-04 5.4596e-06 1.4062e-01)
+(3.2123e-01 3.1009e-03 1.4506e-04 3.6204e-04 2.6288e-06 9.0947e-02)
+(1.7438e-01 1.2143e-03 6.6865e-05 1.1756e-04 1.0969e-06 5.3564e-02)
+(8.4454e-02 4.0027e-04 2.8297e-05 3.0014e-05 3.8272e-07 2.7852e-02)
+(3.4718e-02 1.0324e-04 1.0732e-05 5.4400e-06 1.0411e-07 1.2147e-02)
+(1.1020e-02 1.8113e-05 3.3616e-06 5.7878e-07 1.9051e-08 4.0392e-03)
+(2.1828e-03 1.5728e-06 8.3228e-07 2.3603e-08 1.7141e-09 8.2715e-04)
+(1.3684e-04 2.4414e-08 5.9464e-08 9.4734e-11 2.8525e-11 5.2943e-05)
+(5.3376e-28 -1.1206e-30 -6.1282e-30 4.9018e-29 -6.2793e-32 5.9575e-28)
+)
diff --git a/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.digitalFilter/0/UMean b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.digitalFilter/0/UMean
new file mode 100644
index 0000000000..2ecd2a6ad4
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.digitalFilter/0/UMean
@@ -0,0 +1,259 @@
+(
+( 0.0000e+00 0 0 )
+( 2.9538e-02 0 0 )
+( 1.1811e-01 0 0 )
+( 2.6562e-01 0 0 )
+( 4.7198e-01 0 0 )
+( 7.3701e-01 0 0 )
+( 1.0605e+00 0 0 )
+( 1.4417e+00 0 0 )
+( 1.8799e+00 0 0 )
+( 2.3729e+00 0 0 )
+( 2.9177e+00 0 0 )
+( 3.5093e+00 0 0 )
+( 4.1409e+00 0 0 )
+( 4.8032e+00 0 0 )
+( 5.4854e+00 0 0 )
+( 6.1754e+00 0 0 )
+( 6.8611e+00 0 0 )
+( 7.5309e+00 0 0 )
+( 8.1754e+00 0 0 )
+( 8.7870e+00 0 0 )
+( 9.3607e+00 0 0 )
+( 9.8937e+00 0 0 )
+( 1.0385e+01 0 0 )
+( 1.0836e+01 0 0 )
+( 1.1248e+01 0 0 )
+( 1.1624e+01 0 0 )
+( 1.1966e+01 0 0 )
+( 1.2278e+01 0 0 )
+( 1.2563e+01 0 0 )
+( 1.2822e+01 0 0 )
+( 1.3060e+01 0 0 )
+( 1.3278e+01 0 0 )
+( 1.3479e+01 0 0 )
+( 1.3664e+01 0 0 )
+( 1.3837e+01 0 0 )
+( 1.3998e+01 0 0 )
+( 1.4148e+01 0 0 )
+( 1.4290e+01 0 0 )
+( 1.4425e+01 0 0 )
+( 1.4552e+01 0 0 )
+( 1.4673e+01 0 0 )
+( 1.4790e+01 0 0 )
+( 1.4902e+01 0 0 )
+( 1.5011e+01 0 0 )
+( 1.5117e+01 0 0 )
+( 1.5221e+01 0 0 )
+( 1.5322e+01 0 0 )
+( 1.5421e+01 0 0 )
+( 1.5518e+01 0 0 )
+( 1.5614e+01 0 0 )
+( 1.5707e+01 0 0 )
+( 1.5799e+01 0 0 )
+( 1.5890e+01 0 0 )
+( 1.5979e+01 0 0 )
+( 1.6067e+01 0 0 )
+( 1.6153e+01 0 0 )
+( 1.6239e+01 0 0 )
+( 1.6324e+01 0 0 )
+( 1.6409e+01 0 0 )
+( 1.6493e+01 0 0 )
+( 1.6576e+01 0 0 )
+( 1.6659e+01 0 0 )
+( 1.6741e+01 0 0 )
+( 1.6823e+01 0 0 )
+( 1.6903e+01 0 0 )
+( 1.6984e+01 0 0 )
+( 1.7063e+01 0 0 )
+( 1.7141e+01 0 0 )
+( 1.7218e+01 0 0 )
+( 1.7294e+01 0 0 )
+( 1.7369e+01 0 0 )
+( 1.7443e+01 0 0 )
+( 1.7517e+01 0 0 )
+( 1.7590e+01 0 0 )
+( 1.7664e+01 0 0 )
+( 1.7738e+01 0 0 )
+( 1.7812e+01 0 0 )
+( 1.7886e+01 0 0 )
+( 1.7960e+01 0 0 )
+( 1.8034e+01 0 0 )
+( 1.8108e+01 0 0 )
+( 1.8182e+01 0 0 )
+( 1.8254e+01 0 0 )
+( 1.8326e+01 0 0 )
+( 1.8396e+01 0 0 )
+( 1.8466e+01 0 0 )
+( 1.8535e+01 0 0 )
+( 1.8603e+01 0 0 )
+( 1.8669e+01 0 0 )
+( 1.8734e+01 0 0 )
+( 1.8797e+01 0 0 )
+( 1.8859e+01 0 0 )
+( 1.8919e+01 0 0 )
+( 1.8978e+01 0 0 )
+( 1.9035e+01 0 0 )
+( 1.9092e+01 0 0 )
+( 1.9148e+01 0 0 )
+( 1.9202e+01 0 0 )
+( 1.9256e+01 0 0 )
+( 1.9308e+01 0 0 )
+( 1.9359e+01 0 0 )
+( 1.9408e+01 0 0 )
+( 1.9456e+01 0 0 )
+( 1.9503e+01 0 0 )
+( 1.9548e+01 0 0 )
+( 1.9593e+01 0 0 )
+( 1.9636e+01 0 0 )
+( 1.9678e+01 0 0 )
+( 1.9719e+01 0 0 )
+( 1.9758e+01 0 0 )
+( 1.9796e+01 0 0 )
+( 1.9832e+01 0 0 )
+( 1.9865e+01 0 0 )
+( 1.9897e+01 0 0 )
+( 1.9927e+01 0 0 )
+( 1.9955e+01 0 0 )
+( 1.9981e+01 0 0 )
+( 2.0004e+01 0 0 )
+( 2.0026e+01 0 0 )
+( 2.0046e+01 0 0 )
+( 2.0064e+01 0 0 )
+( 2.0080e+01 0 0 )
+( 2.0094e+01 0 0 )
+( 2.0106e+01 0 0 )
+( 2.0116e+01 0 0 )
+( 2.0123e+01 0 0 )
+( 2.0129e+01 0 0 )
+( 2.0132e+01 0 0 )
+( 2.0133e+01 0 0 )
+( 2.0132e+01 0 0 )
+( 2.0129e+01 0 0 )
+( 2.0123e+01 0 0 )
+( 2.0116e+01 0 0 )
+( 2.0106e+01 0 0 )
+( 2.0094e+01 0 0 )
+( 2.0080e+01 0 0 )
+( 2.0064e+01 0 0 )
+( 2.0046e+01 0 0 )
+( 2.0026e+01 0 0 )
+( 2.0004e+01 0 0 )
+( 1.9981e+01 0 0 )
+( 1.9955e+01 0 0 )
+( 1.9927e+01 0 0 )
+( 1.9897e+01 0 0 )
+( 1.9865e+01 0 0 )
+( 1.9832e+01 0 0 )
+( 1.9796e+01 0 0 )
+( 1.9758e+01 0 0 )
+( 1.9719e+01 0 0 )
+( 1.9678e+01 0 0 )
+( 1.9636e+01 0 0 )
+( 1.9593e+01 0 0 )
+( 1.9548e+01 0 0 )
+( 1.9503e+01 0 0 )
+( 1.9456e+01 0 0 )
+( 1.9408e+01 0 0 )
+( 1.9359e+01 0 0 )
+( 1.9308e+01 0 0 )
+( 1.9256e+01 0 0 )
+( 1.9202e+01 0 0 )
+( 1.9148e+01 0 0 )
+( 1.9092e+01 0 0 )
+( 1.9035e+01 0 0 )
+( 1.8978e+01 0 0 )
+( 1.8919e+01 0 0 )
+( 1.8859e+01 0 0 )
+( 1.8797e+01 0 0 )
+( 1.8734e+01 0 0 )
+( 1.8669e+01 0 0 )
+( 1.8603e+01 0 0 )
+( 1.8535e+01 0 0 )
+( 1.8466e+01 0 0 )
+( 1.8396e+01 0 0 )
+( 1.8326e+01 0 0 )
+( 1.8254e+01 0 0 )
+( 1.8182e+01 0 0 )
+( 1.8108e+01 0 0 )
+( 1.8034e+01 0 0 )
+( 1.7960e+01 0 0 )
+( 1.7886e+01 0 0 )
+( 1.7812e+01 0 0 )
+( 1.7738e+01 0 0 )
+( 1.7664e+01 0 0 )
+( 1.7590e+01 0 0 )
+( 1.7517e+01 0 0 )
+( 1.7443e+01 0 0 )
+( 1.7369e+01 0 0 )
+( 1.7294e+01 0 0 )
+( 1.7218e+01 0 0 )
+( 1.7141e+01 0 0 )
+( 1.7063e+01 0 0 )
+( 1.6984e+01 0 0 )
+( 1.6903e+01 0 0 )
+( 1.6823e+01 0 0 )
+( 1.6741e+01 0 0 )
+( 1.6659e+01 0 0 )
+( 1.6576e+01 0 0 )
+( 1.6493e+01 0 0 )
+( 1.6409e+01 0 0 )
+( 1.6324e+01 0 0 )
+( 1.6239e+01 0 0 )
+( 1.6153e+01 0 0 )
+( 1.6067e+01 0 0 )
+( 1.5979e+01 0 0 )
+( 1.5890e+01 0 0 )
+( 1.5799e+01 0 0 )
+( 1.5707e+01 0 0 )
+( 1.5614e+01 0 0 )
+( 1.5518e+01 0 0 )
+( 1.5421e+01 0 0 )
+( 1.5322e+01 0 0 )
+( 1.5221e+01 0 0 )
+( 1.5117e+01 0 0 )
+( 1.5011e+01 0 0 )
+( 1.4902e+01 0 0 )
+( 1.4790e+01 0 0 )
+( 1.4673e+01 0 0 )
+( 1.4552e+01 0 0 )
+( 1.4425e+01 0 0 )
+( 1.4290e+01 0 0 )
+( 1.4148e+01 0 0 )
+( 1.3998e+01 0 0 )
+( 1.3837e+01 0 0 )
+( 1.3664e+01 0 0 )
+( 1.3479e+01 0 0 )
+( 1.3278e+01 0 0 )
+( 1.3060e+01 0 0 )
+( 1.2822e+01 0 0 )
+( 1.2563e+01 0 0 )
+( 1.2278e+01 0 0 )
+( 1.1966e+01 0 0 )
+( 1.1624e+01 0 0 )
+( 1.1248e+01 0 0 )
+( 1.0836e+01 0 0 )
+( 1.0385e+01 0 0 )
+( 9.8937e+00 0 0 )
+( 9.3607e+00 0 0 )
+( 8.7870e+00 0 0 )
+( 8.1754e+00 0 0 )
+( 7.5309e+00 0 0 )
+( 6.8611e+00 0 0 )
+( 6.1754e+00 0 0 )
+( 5.4854e+00 0 0 )
+( 4.8032e+00 0 0 )
+( 4.1409e+00 0 0 )
+( 3.5093e+00 0 0 )
+( 2.9177e+00 0 0 )
+( 2.3729e+00 0 0 )
+( 1.8799e+00 0 0 )
+( 1.4417e+00 0 0 )
+( 1.0605e+00 0 0 )
+( 7.3701e-01 0 0 )
+( 4.7198e-01 0 0 )
+( 2.6562e-01 0 0 )
+( 1.1811e-01 0 0 )
+( 2.9538e-02 0 0 )
+( 0.0000e+00 0 0 )
+)
diff --git a/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.digitalFilter/points b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.digitalFilter/points
new file mode 100644
index 0000000000..5474eb4eae
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.digitalFilter/points
@@ -0,0 +1,259 @@
+(
+(0 0.0000e+00 0)
+(0 7.5298e-05 0)
+(0 3.0118e-04 0)
+(0 6.7762e-04 0)
+(0 1.2045e-03 0)
+(0 1.8819e-03 0)
+(0 2.7095e-03 0)
+(0 3.6874e-03 0)
+(0 4.8153e-03 0)
+(0 6.0930e-03 0)
+(0 7.5205e-03 0)
+(0 9.0974e-03 0)
+(0 1.0823e-02 0)
+(0 1.2699e-02 0)
+(0 1.4722e-02 0)
+(0 1.6895e-02 0)
+(0 1.9215e-02 0)
+(0 2.1683e-02 0)
+(0 2.4298e-02 0)
+(0 2.7060e-02 0)
+(0 2.9969e-02 0)
+(0 3.3024e-02 0)
+(0 3.6224e-02 0)
+(0 3.9569e-02 0)
+(0 4.3060e-02 0)
+(0 4.6694e-02 0)
+(0 5.0472e-02 0)
+(0 5.4393e-02 0)
+(0 5.8456e-02 0)
+(0 6.2661e-02 0)
+(0 6.7007e-02 0)
+(0 7.1494e-02 0)
+(0 7.6120e-02 0)
+(0 8.0886e-02 0)
+(0 8.5790e-02 0)
+(0 9.0832e-02 0)
+(0 9.6011e-02 0)
+(0 1.0133e-01 0)
+(0 1.0678e-01 0)
+(0 1.1236e-01 0)
+(0 1.1808e-01 0)
+(0 1.2393e-01 0)
+(0 1.2991e-01 0)
+(0 1.3603e-01 0)
+(0 1.4227e-01 0)
+(0 1.4864e-01 0)
+(0 1.5515e-01 0)
+(0 1.6178e-01 0)
+(0 1.6853e-01 0)
+(0 1.7541e-01 0)
+(0 1.8242e-01 0)
+(0 1.8954e-01 0)
+(0 1.9679e-01 0)
+(0 2.0416e-01 0)
+(0 2.1165e-01 0)
+(0 2.1926e-01 0)
+(0 2.2699e-01 0)
+(0 2.3483e-01 0)
+(0 2.4279e-01 0)
+(0 2.5086e-01 0)
+(0 2.5905e-01 0)
+(0 2.6735e-01 0)
+(0 2.7575e-01 0)
+(0 2.8427e-01 0)
+(0 2.9289e-01 0)
+(0 3.0162e-01 0)
+(0 3.1046e-01 0)
+(0 3.1940e-01 0)
+(0 3.2844e-01 0)
+(0 3.3758e-01 0)
+(0 3.4683e-01 0)
+(0 3.5617e-01 0)
+(0 3.6561e-01 0)
+(0 3.7514e-01 0)
+(0 3.8477e-01 0)
+(0 3.9449e-01 0)
+(0 4.0430e-01 0)
+(0 4.1420e-01 0)
+(0 4.2419e-01 0)
+(0 4.3427e-01 0)
+(0 4.4443e-01 0)
+(0 4.5467e-01 0)
+(0 4.6500e-01 0)
+(0 4.7541e-01 0)
+(0 4.8590e-01 0)
+(0 4.9646e-01 0)
+(0 5.0710e-01 0)
+(0 5.1782e-01 0)
+(0 5.2860e-01 0)
+(0 5.3946e-01 0)
+(0 5.5039e-01 0)
+(0 5.6138e-01 0)
+(0 5.7244e-01 0)
+(0 5.8357e-01 0)
+(0 5.9476e-01 0)
+(0 6.0601e-01 0)
+(0 6.1732e-01 0)
+(0 6.2868e-01 0)
+(0 6.4010e-01 0)
+(0 6.5158e-01 0)
+(0 6.6311e-01 0)
+(0 6.7469e-01 0)
+(0 6.8632e-01 0)
+(0 6.9799e-01 0)
+(0 7.0972e-01 0)
+(0 7.2148e-01 0)
+(0 7.3329e-01 0)
+(0 7.4513e-01 0)
+(0 7.5702e-01 0)
+(0 7.6894e-01 0)
+(0 7.8090e-01 0)
+(0 7.9289e-01 0)
+(0 8.0491e-01 0)
+(0 8.1696e-01 0)
+(0 8.2904e-01 0)
+(0 8.4114e-01 0)
+(0 8.5327e-01 0)
+(0 8.6542e-01 0)
+(0 8.7759e-01 0)
+(0 8.8978e-01 0)
+(0 9.0198e-01 0)
+(0 9.1420e-01 0)
+(0 9.2644e-01 0)
+(0 9.3868e-01 0)
+(0 9.5093e-01 0)
+(0 9.6319e-01 0)
+(0 9.7546e-01 0)
+(0 9.8773e-01 0)
+(0 1.0000e-00 0)
+(0 1.01227 0)
+(0 1.02454 0)
+(0 1.03681 0)
+(0 1.04907 0)
+(0 1.06132 0)
+(0 1.07356 0)
+(0 1.0858 0)
+(0 1.09802 0)
+(0 1.11022 0)
+(0 1.12241 0)
+(0 1.13458 0)
+(0 1.14673 0)
+(0 1.15886 0)
+(0 1.17096 0)
+(0 1.18304 0)
+(0 1.19509 0)
+(0 1.20711 0)
+(0 1.2191 0)
+(0 1.23106 0)
+(0 1.24298 0)
+(0 1.25487 0)
+(0 1.26671 0)
+(0 1.27852 0)
+(0 1.29028 0)
+(0 1.30201 0)
+(0 1.31368 0)
+(0 1.32531 0)
+(0 1.33689 0)
+(0 1.34842 0)
+(0 1.3599 0)
+(0 1.37132 0)
+(0 1.38268 0)
+(0 1.39399 0)
+(0 1.40524 0)
+(0 1.41643 0)
+(0 1.42756 0)
+(0 1.43862 0)
+(0 1.44961 0)
+(0 1.46054 0)
+(0 1.4714 0)
+(0 1.48218 0)
+(0 1.4929 0)
+(0 1.50354 0)
+(0 1.5141 0)
+(0 1.52459 0)
+(0 1.535 0)
+(0 1.54533 0)
+(0 1.55557 0)
+(0 1.56573 0)
+(0 1.57581 0)
+(0 1.5858 0)
+(0 1.5957 0)
+(0 1.60551 0)
+(0 1.61523 0)
+(0 1.62486 0)
+(0 1.63439 0)
+(0 1.64383 0)
+(0 1.65317 0)
+(0 1.66242 0)
+(0 1.67156 0)
+(0 1.6806 0)
+(0 1.68954 0)
+(0 1.69838 0)
+(0 1.70711 0)
+(0 1.71573 0)
+(0 1.72425 0)
+(0 1.73265 0)
+(0 1.74095 0)
+(0 1.74914 0)
+(0 1.75721 0)
+(0 1.76517 0)
+(0 1.77301 0)
+(0 1.78074 0)
+(0 1.78835 0)
+(0 1.79584 0)
+(0 1.80321 0)
+(0 1.81046 0)
+(0 1.81758 0)
+(0 1.82459 0)
+(0 1.83147 0)
+(0 1.83822 0)
+(0 1.84485 0)
+(0 1.85136 0)
+(0 1.85773 0)
+(0 1.86397 0)
+(0 1.87009 0)
+(0 1.87607 0)
+(0 1.88192 0)
+(0 1.88764 0)
+(0 1.89322 0)
+(0 1.89867 0)
+(0 1.90399 0)
+(0 1.90917 0)
+(0 1.91421 0)
+(0 1.91911 0)
+(0 1.92388 0)
+(0 1.92851 0)
+(0 1.93299 0)
+(0 1.93734 0)
+(0 1.94154 0)
+(0 1.94561 0)
+(0 1.94953 0)
+(0 1.95331 0)
+(0 1.95694 0)
+(0 1.96043 0)
+(0 1.96378 0)
+(0 1.96698 0)
+(0 1.97003 0)
+(0 1.97294 0)
+(0 1.9757 0)
+(0 1.97832 0)
+(0 1.98079 0)
+(0 1.9831 0)
+(0 1.98528 0)
+(0 1.9873 0)
+(0 1.98918 0)
+(0 1.9909 0)
+(0 1.99248 0)
+(0 1.99391 0)
+(0 1.99518 0)
+(0 1.99631 0)
+(0 1.99729 0)
+(0 1.99812 0)
+(0 1.9988 0)
+(0 1.99932 0)
+(0 1.9997 0)
+(0 1.99992 0)
+(0 2 0)
+)
diff --git a/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.reducedDigitalFilter/0/R b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.reducedDigitalFilter/0/R
new file mode 100644
index 0000000000..084bdd9f8f
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.reducedDigitalFilter/0/R
@@ -0,0 +1,259 @@
+(
+(5.3376e-28 1.1206e-30 -6.1282e-30 4.9018e-29 6.2793e-32 5.9575e-28)
+(1.3684e-04 -2.4414e-08 5.9464e-08 9.4734e-11 -2.8525e-11 5.2943e-05)
+(2.1828e-03 -1.5728e-06 8.3228e-07 2.3603e-08 -1.7141e-09 8.2715e-04)
+(1.1020e-02 -1.8113e-05 3.3616e-06 5.7878e-07 -1.9051e-08 4.0392e-03)
+(3.4718e-02 -1.0324e-04 1.0732e-05 5.4400e-06 -1.0411e-07 1.2147e-02)
+(8.4454e-02 -4.0027e-04 2.8297e-05 3.0014e-05 -3.8272e-07 2.7852e-02)
+(1.7438e-01 -1.2143e-03 6.6865e-05 1.1756e-04 -1.0969e-06 5.3564e-02)
+(3.2123e-01 -3.1009e-03 1.4506e-04 3.6204e-04 -2.6288e-06 9.0947e-02)
+(5.4336e-01 -6.9486e-03 2.9097e-04 9.3203e-04 -5.4596e-06 1.4062e-01)
+(8.5851e-01 -1.4015e-02 5.4723e-04 2.0879e-03 -1.0091e-05 2.0209e-01)
+(1.2802e+00 -2.5858e-02 9.7373e-04 4.1852e-03 -1.7119e-05 2.7384e-01)
+(1.8132e+00 -4.4135e-02 1.6529e-03 7.6610e-03 -2.7764e-05 3.5369e-01)
+(2.4490e+00 -7.0276e-02 2.6847e-03 1.3005e-02 -4.4587e-05 4.3919e-01)
+(3.1644e+00 -1.0513e-01 4.1721e-03 2.0724e-02 -7.1044e-05 5.2804e-01)
+(3.9229e+00 -1.4869e-01 6.1889e-03 3.1300e-02 -1.0893e-04 6.1834e-01)
+(4.6811e+00 -1.9999e-01 8.7562e-03 4.5157e-02 -1.5460e-04 7.0871e-01)
+(5.3953e+00 -2.5724e-01 1.1834e-02 6.2636e-02 -1.9552e-04 7.9819e-01)
+(6.0286e+00 -3.1808e-01 1.5387e-02 8.3971e-02 -2.1029e-04 8.8610e-01)
+(6.5558e+00 -3.8001e-01 1.9449e-02 1.0928e-01 -1.7262e-04 9.7178e-01)
+(6.9642e+00 -4.4073e-01 2.4133e-02 1.3854e-01 -6.0490e-05 1.0545e+00)
+(7.2524e+00 -4.9837e-01 2.9501e-02 1.7162e-01 1.3268e-04 1.1333e+00)
+(7.4274e+00 -5.5162e-01 3.5418e-02 2.0824e-01 3.9086e-04 1.2075e+00)
+(7.5016e+00 -5.9969e-01 4.1511e-02 2.4803e-01 6.7818e-04 1.2762e+00)
+(7.4897e+00 -6.4224e-01 4.7330e-02 2.9050e-01 9.5511e-04 1.3387e+00)
+(7.4069e+00 -6.7924e-01 5.2512e-02 3.3510e-01 1.1989e-03 1.3948e+00)
+(7.2679e+00 -7.1090e-01 5.6926e-02 3.8121e-01 1.4016e-03 1.4444e+00)
+(7.0865e+00 -7.3759e-01 6.0659e-02 4.2818e-01 1.5532e-03 1.4878e+00)
+(6.8750e+00 -7.5980e-01 6.3835e-02 4.7535e-01 1.6360e-03 1.5251e+00)
+(6.6441e+00 -7.7810e-01 6.6418e-02 5.2211e-01 1.6465e-03 1.5568e+00)
+(6.4025e+00 -7.9302e-01 6.8207e-02 5.6790e-01 1.6026e-03 1.5833e+00)
+(6.1573e+00 -8.0503e-01 6.9001e-02 6.1229e-01 1.5073e-03 1.6047e+00)
+(5.9141e+00 -8.1452e-01 6.8896e-02 6.5493e-01 1.3268e-03 1.6216e+00)
+(5.6767e+00 -8.2184e-01 6.8285e-02 6.9554e-01 1.0303e-03 1.6344e+00)
+(5.4477e+00 -8.2733e-01 6.7481e-02 7.3388e-01 6.4119e-04 1.6439e+00)
+(5.2291e+00 -8.3133e-01 6.6450e-02 7.6971e-01 2.3642e-04 1.6512e+00)
+(5.0223e+00 -8.3412e-01 6.4937e-02 8.0285e-01 -1.3612e-04 1.6567e+00)
+(4.8279e+00 -8.3590e-01 6.2875e-02 8.3316e-01 -4.7569e-04 1.6606e+00)
+(4.6462e+00 -8.3683e-01 6.0710e-02 8.6061e-01 -7.8412e-04 1.6623e+00)
+(4.4768e+00 -8.3700e-01 5.9182e-02 8.8520e-01 -1.0539e-03 1.6618e+00)
+(4.3194e+00 -8.3649e-01 5.8527e-02 9.0700e-01 -1.2655e-03 1.6594e+00)
+(4.1730e+00 -8.3523e-01 5.8437e-02 9.2607e-01 -1.4126e-03 1.6553e+00)
+(4.0370e+00 -8.3321e-01 5.8512e-02 9.4247e-01 -1.4930e-03 1.6500e+00)
+(3.9106e+00 -8.3044e-01 5.8376e-02 9.5626e-01 -1.5118e-03 1.6439e+00)
+(3.7930e+00 -8.2687e-01 5.8084e-02 9.6754e-01 -1.5196e-03 1.6373e+00)
+(3.6834e+00 -8.2252e-01 5.8125e-02 9.7651e-01 -1.6228e-03 1.6299e+00)
+(3.5821e+00 -8.1762e-01 5.8537e-02 9.8341e-01 -1.8295e-03 1.6209e+00)
+(3.4896e+00 -8.1249e-01 5.8809e-02 9.8850e-01 -2.0145e-03 1.6105e+00)
+(3.4053e+00 -8.0726e-01 5.8494e-02 9.9197e-01 -2.0826e-03 1.5991e+00)
+(3.3278e+00 -8.0180e-01 5.7876e-02 9.9396e-01 -2.1101e-03 1.5872e+00)
+(3.2557e+00 -7.9592e-01 5.7034e-02 9.9455e-01 -2.1628e-03 1.5754e+00)
+(3.1886e+00 -7.8961e-01 5.5746e-02 9.9379e-01 -2.2441e-03 1.5637e+00)
+(3.1269e+00 -7.8296e-01 5.4259e-02 9.9183e-01 -2.3220e-03 1.5510e+00)
+(3.0702e+00 -7.7613e-01 5.2957e-02 9.8887e-01 -2.3519e-03 1.5367e+00)
+(3.0178e+00 -7.6929e-01 5.1819e-02 9.8512e-01 -2.3018e-03 1.5212e+00)
+(2.9678e+00 -7.6243e-01 5.0872e-02 9.8075e-01 -2.0851e-03 1.5052e+00)
+(2.9191e+00 -7.5552e-01 5.0514e-02 9.7583e-01 -1.6209e-03 1.4889e+00)
+(2.8712e+00 -7.4851e-01 5.0564e-02 9.7044e-01 -8.5802e-04 1.4721e+00)
+(2.8239e+00 -7.4146e-01 5.0229e-02 9.6465e-01 1.1783e-04 1.4551e+00)
+(2.7776e+00 -7.3445e-01 4.9171e-02 9.5846e-01 1.0337e-03 1.4379e+00)
+(2.7324e+00 -7.2737e-01 4.7565e-02 9.5179e-01 1.6440e-03 1.4206e+00)
+(2.6883e+00 -7.2012e-01 4.5716e-02 9.4459e-01 1.9851e-03 1.4035e+00)
+(2.6453e+00 -7.1271e-01 4.3916e-02 9.3679e-01 2.2291e-03 1.3867e+00)
+(2.6031e+00 -7.0521e-01 4.2134e-02 9.2845e-01 2.4143e-03 1.3696e+00)
+(2.5612e+00 -6.9761e-01 4.0077e-02 9.1976e-01 2.5128e-03 1.3515e+00)
+(2.5196e+00 -6.8985e-01 3.7520e-02 9.1087e-01 2.5006e-03 1.3322e+00)
+(2.4793e+00 -6.8192e-01 3.4358e-02 9.0185e-01 2.3726e-03 1.3118e+00)
+(2.4411e+00 -6.7396e-01 3.0858e-02 8.9268e-01 2.1952e-03 1.2910e+00)
+(2.4055e+00 -6.6598e-01 2.7094e-02 8.8335e-01 2.0517e-03 1.2700e+00)
+(2.3714e+00 -6.5794e-01 2.3371e-02 8.7385e-01 1.9999e-03 1.2488e+00)
+(2.3380e+00 -6.4977e-01 2.0500e-02 8.6418e-01 1.9700e-03 1.2276e+00)
+(2.3046e+00 -6.4148e-01 1.8795e-02 8.5436e-01 1.8189e-03 1.2066e+00)
+(2.2708e+00 -6.3298e-01 1.7891e-02 8.4439e-01 1.4336e-03 1.1855e+00)
+(2.2365e+00 -6.2415e-01 1.7572e-02 8.3421e-01 7.5375e-04 1.1641e+00)
+(2.2023e+00 -6.1500e-01 1.7517e-02 8.2369e-01 -7.9290e-05 1.1423e+00)
+(2.1682e+00 -6.0551e-01 1.7095e-02 8.1285e-01 -8.3720e-04 1.1202e+00)
+(2.1335e+00 -5.9561e-01 1.5913e-02 8.0181e-01 -1.3413e-03 1.0983e+00)
+(2.0983e+00 -5.8530e-01 1.4306e-02 7.9074e-01 -1.5771e-03 1.0767e+00)
+(2.0631e+00 -5.7460e-01 1.2812e-02 7.7968e-01 -1.5179e-03 1.0551e+00)
+(2.0277e+00 -5.6350e-01 1.1907e-02 7.6856e-01 -1.1762e-03 1.0339e+00)
+(1.9920e+00 -5.5204e-01 1.1640e-02 7.5731e-01 -7.0498e-04 1.0131e+00)
+(1.9568e+00 -5.4040e-01 1.1503e-02 7.4588e-01 -2.0689e-04 9.9279e-01)
+(1.9229e+00 -5.2883e-01 1.1097e-02 7.3425e-01 2.0529e-04 9.7254e-01)
+(1.8896e+00 -5.1747e-01 1.0218e-02 7.2249e-01 5.0971e-04 9.5187e-01)
+(1.8562e+00 -5.0626e-01 9.0004e-03 7.1064e-01 7.6211e-04 9.3057e-01)
+(1.8232e+00 -4.9507e-01 7.8753e-03 6.9862e-01 9.2292e-04 9.0936e-01)
+(1.7908e+00 -4.8386e-01 7.1617e-03 6.8642e-01 9.3017e-04 8.8926e-01)
+(1.7591e+00 -4.7269e-01 7.0121e-03 6.7417e-01 7.2814e-04 8.7042e-01)
+(1.7274e+00 -4.6164e-01 7.3802e-03 6.6201e-01 2.7352e-04 8.5243e-01)
+(1.6948e+00 -4.5072e-01 8.2299e-03 6.5007e-01 -2.9303e-04 8.3469e-01)
+(1.6610e+00 -4.3999e-01 9.3954e-03 6.3843e-01 -7.4565e-04 8.1690e-01)
+(1.6264e+00 -4.2949e-01 1.0585e-02 6.2712e-01 -9.5384e-04 7.9939e-01)
+(1.5910e+00 -4.1910e-01 1.1834e-02 6.1613e-01 -9.5389e-04 7.8241e-01)
+(1.5550e+00 -4.0876e-01 1.3728e-02 6.0545e-01 -1.0075e-03 7.6587e-01)
+(1.5184e+00 -3.9853e-01 1.6563e-02 5.9512e-01 -1.2970e-03 7.4931e-01)
+(1.4808e+00 -3.8843e-01 1.9835e-02 5.8523e-01 -1.7620e-03 7.3258e-01)
+(1.4425e+00 -3.7835e-01 2.2555e-02 5.7589e-01 -2.2305e-03 7.1572e-01)
+(1.4037e+00 -3.6810e-01 2.3965e-02 5.6715e-01 -2.5637e-03 6.9868e-01)
+(1.3647e+00 -3.5761e-01 2.4099e-02 5.5896e-01 -2.7183e-03 6.8159e-01)
+(1.3260e+00 -3.4695e-01 2.3378e-02 5.5120e-01 -2.7095e-03 6.6492e-01)
+(1.2881e+00 -3.3621e-01 2.2236e-02 5.4372e-01 -2.5722e-03 6.4916e-01)
+(1.2513e+00 -3.2544e-01 2.0871e-02 5.3642e-01 -2.4646e-03 6.3424e-01)
+(1.2162e+00 -3.1468e-01 1.9224e-02 5.2919e-01 -2.5472e-03 6.2001e-01)
+(1.1827e+00 -3.0393e-01 1.7273e-02 5.2197e-01 -2.7871e-03 6.0660e-01)
+(1.1500e+00 -2.9315e-01 1.5116e-02 5.1473e-01 -3.0094e-03 5.9377e-01)
+(1.1175e+00 -2.8231e-01 1.3055e-02 5.0750e-01 -3.1202e-03 5.8101e-01)
+(1.0843e+00 -2.7133e-01 1.1360e-02 5.0029e-01 -3.1772e-03 5.6815e-01)
+(1.0507e+00 -2.6019e-01 9.9787e-03 4.9312e-01 -3.2410e-03 5.5541e-01)
+(1.0173e+00 -2.4893e-01 8.7677e-03 4.8605e-01 -3.3171e-03 5.4304e-01)
+(9.8478e-01 -2.3753e-01 7.8868e-03 4.7922e-01 -3.3738e-03 5.3113e-01)
+(9.5377e-01 -2.2605e-01 7.6156e-03 4.7268e-01 -3.4135e-03 5.1971e-01)
+(9.2454e-01 -2.1462e-01 7.9424e-03 4.6644e-01 -3.4897e-03 5.0882e-01)
+(8.9706e-01 -2.0326e-01 8.6379e-03 4.6047e-01 -3.6720e-03 4.9853e-01)
+(8.7113e-01 -1.9188e-01 9.5840e-03 4.5478e-01 -3.9882e-03 4.8883e-01)
+(8.4645e-01 -1.8043e-01 1.0540e-02 4.4932e-01 -4.3896e-03 4.7982e-01)
+(8.2296e-01 -1.6898e-01 1.1159e-02 4.4410e-01 -4.7685e-03 4.7164e-01)
+(8.0062e-01 -1.5752e-01 1.1209e-02 4.3915e-01 -5.0417e-03 4.6392e-01)
+(7.7929e-01 -1.4601e-01 1.0617e-02 4.3458e-01 -5.1687e-03 4.5618e-01)
+(7.5889e-01 -1.3442e-01 9.5260e-03 4.3044e-01 -5.1040e-03 4.4842e-01)
+(7.3932e-01 -1.2272e-01 8.2998e-03 4.2671e-01 -4.8287e-03 4.4071e-01)
+(7.2051e-01 -1.1094e-01 7.2287e-03 4.2339e-01 -4.4056e-03 4.3345e-01)
+(7.0261e-01 -9.9120e-02 6.4502e-03 4.2046e-01 -3.8984e-03 4.2708e-01)
+(6.8576e-01 -8.7301e-02 6.1696e-03 4.1796e-01 -3.3462e-03 4.2169e-01)
+(6.7023e-01 -7.5441e-02 6.4277e-03 4.1591e-01 -2.7927e-03 4.1709e-01)
+(6.5620e-01 -6.3452e-02 7.0238e-03 4.1431e-01 -2.2376e-03 4.1319e-01)
+(6.4389e-01 -5.1252e-02 7.7818e-03 4.1313e-01 -1.6855e-03 4.0984e-01)
+(6.3375e-01 -3.8772e-02 8.7489e-03 4.1234e-01 -1.1858e-03 4.0689e-01)
+(6.2626e-01 -2.6001e-02 9.8627e-03 4.1188e-01 -7.6554e-04 4.0431e-01)
+(6.2174e-01 -1.3031e-02 1.0731e-02 4.1166e-01 -3.8835e-04 4.0247e-01)
+(6.2024e-01 0.0000e+00 1.1032e-02 4.1160e-01 0.0000e+00 4.0179e-01)
+(6.2174e-01 1.3031e-02 1.0731e-02 4.1166e-01 3.8835e-04 4.0247e-01)
+(6.2626e-01 2.6001e-02 9.8627e-03 4.1188e-01 7.6554e-04 4.0431e-01)
+(6.3375e-01 3.8772e-02 8.7489e-03 4.1234e-01 1.1858e-03 4.0689e-01)
+(6.4389e-01 5.1252e-02 7.7818e-03 4.1313e-01 1.6855e-03 4.0984e-01)
+(6.5620e-01 6.3452e-02 7.0238e-03 4.1431e-01 2.2376e-03 4.1319e-01)
+(6.7023e-01 7.5441e-02 6.4277e-03 4.1591e-01 2.7927e-03 4.1709e-01)
+(6.8576e-01 8.7301e-02 6.1696e-03 4.1796e-01 3.3462e-03 4.2169e-01)
+(7.0261e-01 9.9120e-02 6.4502e-03 4.2046e-01 3.8984e-03 4.2708e-01)
+(7.2051e-01 1.1094e-01 7.2287e-03 4.2339e-01 4.4056e-03 4.3345e-01)
+(7.3932e-01 1.2272e-01 8.2998e-03 4.2671e-01 4.8287e-03 4.4071e-01)
+(7.5889e-01 1.3442e-01 9.5260e-03 4.3044e-01 5.1040e-03 4.4842e-01)
+(7.7929e-01 1.4601e-01 1.0617e-02 4.3458e-01 5.1687e-03 4.5618e-01)
+(8.0062e-01 1.5752e-01 1.1209e-02 4.3915e-01 5.0417e-03 4.6392e-01)
+(8.2296e-01 1.6898e-01 1.1159e-02 4.4410e-01 4.7685e-03 4.7164e-01)
+(8.4645e-01 1.8043e-01 1.0540e-02 4.4932e-01 4.3896e-03 4.7982e-01)
+(8.7113e-01 1.9188e-01 9.5840e-03 4.5478e-01 3.9882e-03 4.8883e-01)
+(8.9706e-01 2.0326e-01 8.6379e-03 4.6047e-01 3.6720e-03 4.9853e-01)
+(9.2454e-01 2.1462e-01 7.9424e-03 4.6644e-01 3.4897e-03 5.0882e-01)
+(9.5377e-01 2.2605e-01 7.6156e-03 4.7268e-01 3.4135e-03 5.1971e-01)
+(9.8478e-01 2.3753e-01 7.8868e-03 4.7922e-01 3.3738e-03 5.3113e-01)
+(1.0173e+00 2.4893e-01 8.7677e-03 4.8605e-01 3.3171e-03 5.4304e-01)
+(1.0507e+00 2.6019e-01 9.9787e-03 4.9312e-01 3.2410e-03 5.5541e-01)
+(1.0843e+00 2.7133e-01 1.1360e-02 5.0029e-01 3.1772e-03 5.6815e-01)
+(1.1175e+00 2.8231e-01 1.3055e-02 5.0750e-01 3.1202e-03 5.8101e-01)
+(1.1500e+00 2.9315e-01 1.5116e-02 5.1473e-01 3.0094e-03 5.9377e-01)
+(1.1827e+00 3.0393e-01 1.7273e-02 5.2197e-01 2.7871e-03 6.0660e-01)
+(1.2162e+00 3.1468e-01 1.9224e-02 5.2919e-01 2.5472e-03 6.2001e-01)
+(1.2513e+00 3.2544e-01 2.0871e-02 5.3642e-01 2.4646e-03 6.3424e-01)
+(1.2881e+00 3.3621e-01 2.2236e-02 5.4372e-01 2.5722e-03 6.4916e-01)
+(1.3260e+00 3.4695e-01 2.3378e-02 5.5120e-01 2.7095e-03 6.6492e-01)
+(1.3647e+00 3.5761e-01 2.4099e-02 5.5896e-01 2.7183e-03 6.8159e-01)
+(1.4037e+00 3.6810e-01 2.3965e-02 5.6715e-01 2.5637e-03 6.9868e-01)
+(1.4425e+00 3.7835e-01 2.2555e-02 5.7589e-01 2.2305e-03 7.1572e-01)
+(1.4808e+00 3.8843e-01 1.9835e-02 5.8523e-01 1.7620e-03 7.3258e-01)
+(1.5184e+00 3.9853e-01 1.6563e-02 5.9512e-01 1.2970e-03 7.4931e-01)
+(1.5550e+00 4.0876e-01 1.3728e-02 6.0545e-01 1.0075e-03 7.6587e-01)
+(1.5910e+00 4.1910e-01 1.1834e-02 6.1613e-01 9.5389e-04 7.8241e-01)
+(1.6264e+00 4.2949e-01 1.0585e-02 6.2712e-01 9.5384e-04 7.9939e-01)
+(1.6610e+00 4.3999e-01 9.3954e-03 6.3843e-01 7.4565e-04 8.1690e-01)
+(1.6948e+00 4.5072e-01 8.2299e-03 6.5007e-01 2.9303e-04 8.3469e-01)
+(1.7274e+00 4.6164e-01 7.3802e-03 6.6201e-01 -2.7352e-04 8.5243e-01)
+(1.7591e+00 4.7269e-01 7.0121e-03 6.7417e-01 -7.2814e-04 8.7042e-01)
+(1.7908e+00 4.8386e-01 7.1617e-03 6.8642e-01 -9.3017e-04 8.8926e-01)
+(1.8232e+00 4.9507e-01 7.8753e-03 6.9862e-01 -9.2292e-04 9.0936e-01)
+(1.8562e+00 5.0626e-01 9.0004e-03 7.1064e-01 -7.6211e-04 9.3057e-01)
+(1.8896e+00 5.1747e-01 1.0218e-02 7.2249e-01 -5.0971e-04 9.5187e-01)
+(1.9229e+00 5.2883e-01 1.1097e-02 7.3425e-01 -2.0529e-04 9.7254e-01)
+(1.9568e+00 5.4040e-01 1.1503e-02 7.4588e-01 2.0689e-04 9.9279e-01)
+(1.9920e+00 5.5204e-01 1.1640e-02 7.5731e-01 7.0498e-04 1.0131e+00)
+(2.0277e+00 5.6350e-01 1.1907e-02 7.6856e-01 1.1762e-03 1.0339e+00)
+(2.0631e+00 5.7460e-01 1.2812e-02 7.7968e-01 1.5179e-03 1.0551e+00)
+(2.0983e+00 5.8530e-01 1.4306e-02 7.9074e-01 1.5771e-03 1.0767e+00)
+(2.1335e+00 5.9561e-01 1.5913e-02 8.0181e-01 1.3413e-03 1.0983e+00)
+(2.1682e+00 6.0551e-01 1.7095e-02 8.1285e-01 8.3720e-04 1.1202e+00)
+(2.2023e+00 6.1500e-01 1.7517e-02 8.2369e-01 7.9290e-05 1.1423e+00)
+(2.2365e+00 6.2415e-01 1.7572e-02 8.3421e-01 -7.5375e-04 1.1641e+00)
+(2.2708e+00 6.3298e-01 1.7891e-02 8.4439e-01 -1.4336e-03 1.1855e+00)
+(2.3046e+00 6.4148e-01 1.8795e-02 8.5436e-01 -1.8189e-03 1.2066e+00)
+(2.3380e+00 6.4977e-01 2.0500e-02 8.6418e-01 -1.9700e-03 1.2276e+00)
+(2.3714e+00 6.5794e-01 2.3371e-02 8.7385e-01 -1.9999e-03 1.2488e+00)
+(2.4055e+00 6.6598e-01 2.7094e-02 8.8335e-01 -2.0517e-03 1.2700e+00)
+(2.4411e+00 6.7396e-01 3.0858e-02 8.9268e-01 -2.1952e-03 1.2910e+00)
+(2.4793e+00 6.8192e-01 3.4358e-02 9.0185e-01 -2.3726e-03 1.3118e+00)
+(2.5196e+00 6.8985e-01 3.7520e-02 9.1087e-01 -2.5006e-03 1.3322e+00)
+(2.5612e+00 6.9761e-01 4.0077e-02 9.1976e-01 -2.5128e-03 1.3515e+00)
+(2.6031e+00 7.0521e-01 4.2134e-02 9.2845e-01 -2.4143e-03 1.3696e+00)
+(2.6453e+00 7.1271e-01 4.3916e-02 9.3679e-01 -2.2291e-03 1.3867e+00)
+(2.6883e+00 7.2012e-01 4.5716e-02 9.4459e-01 -1.9851e-03 1.4035e+00)
+(2.7324e+00 7.2737e-01 4.7565e-02 9.5179e-01 -1.6440e-03 1.4206e+00)
+(2.7776e+00 7.3445e-01 4.9171e-02 9.5846e-01 -1.0337e-03 1.4379e+00)
+(2.8239e+00 7.4146e-01 5.0229e-02 9.6465e-01 -1.1783e-04 1.4551e+00)
+(2.8712e+00 7.4851e-01 5.0564e-02 9.7044e-01 8.5802e-04 1.4721e+00)
+(2.9191e+00 7.5552e-01 5.0514e-02 9.7583e-01 1.6209e-03 1.4889e+00)
+(2.9678e+00 7.6243e-01 5.0872e-02 9.8075e-01 2.0851e-03 1.5052e+00)
+(3.0178e+00 7.6929e-01 5.1819e-02 9.8512e-01 2.3018e-03 1.5212e+00)
+(3.0702e+00 7.7613e-01 5.2957e-02 9.8887e-01 2.3519e-03 1.5367e+00)
+(3.1269e+00 7.8296e-01 5.4259e-02 9.9183e-01 2.3220e-03 1.5510e+00)
+(3.1886e+00 7.8961e-01 5.5746e-02 9.9379e-01 2.2441e-03 1.5637e+00)
+(3.2557e+00 7.9592e-01 5.7034e-02 9.9455e-01 2.1628e-03 1.5754e+00)
+(3.3278e+00 8.0180e-01 5.7876e-02 9.9396e-01 2.1101e-03 1.5872e+00)
+(3.4053e+00 8.0726e-01 5.8494e-02 9.9197e-01 2.0826e-03 1.5991e+00)
+(3.4896e+00 8.1249e-01 5.8809e-02 9.8850e-01 2.0145e-03 1.6105e+00)
+(3.5821e+00 8.1762e-01 5.8537e-02 9.8341e-01 1.8295e-03 1.6209e+00)
+(3.6834e+00 8.2252e-01 5.8125e-02 9.7651e-01 1.6228e-03 1.6299e+00)
+(3.7930e+00 8.2687e-01 5.8084e-02 9.6754e-01 1.5196e-03 1.6373e+00)
+(3.9106e+00 8.3044e-01 5.8376e-02 9.5626e-01 1.5118e-03 1.6439e+00)
+(4.0370e+00 8.3321e-01 5.8512e-02 9.4247e-01 1.4930e-03 1.6500e+00)
+(4.1730e+00 8.3523e-01 5.8437e-02 9.2607e-01 1.4126e-03 1.6553e+00)
+(4.3194e+00 8.3649e-01 5.8527e-02 9.0700e-01 1.2655e-03 1.6594e+00)
+(4.4768e+00 8.3700e-01 5.9182e-02 8.8520e-01 1.0539e-03 1.6618e+00)
+(4.6462e+00 8.3683e-01 6.0710e-02 8.6061e-01 7.8412e-04 1.6623e+00)
+(4.8279e+00 8.3590e-01 6.2875e-02 8.3316e-01 4.7569e-04 1.6606e+00)
+(5.0223e+00 8.3412e-01 6.4937e-02 8.0285e-01 1.3612e-04 1.6567e+00)
+(5.2291e+00 8.3133e-01 6.6450e-02 7.6971e-01 -2.3642e-04 1.6512e+00)
+(5.4477e+00 8.2733e-01 6.7481e-02 7.3388e-01 -6.4119e-04 1.6439e+00)
+(5.6767e+00 8.2184e-01 6.8285e-02 6.9554e-01 -1.0303e-03 1.6344e+00)
+(5.9141e+00 8.1452e-01 6.8896e-02 6.5493e-01 -1.3268e-03 1.6216e+00)
+(6.1573e+00 8.0503e-01 6.9001e-02 6.1229e-01 -1.5073e-03 1.6047e+00)
+(6.4025e+00 7.9302e-01 6.8207e-02 5.6790e-01 -1.6026e-03 1.5833e+00)
+(6.6441e+00 7.7810e-01 6.6418e-02 5.2211e-01 -1.6465e-03 1.5568e+00)
+(6.8750e+00 7.5980e-01 6.3835e-02 4.7535e-01 -1.6360e-03 1.5251e+00)
+(7.0865e+00 7.3759e-01 6.0659e-02 4.2818e-01 -1.5532e-03 1.4878e+00)
+(7.2679e+00 7.1090e-01 5.6926e-02 3.8121e-01 -1.4016e-03 1.4444e+00)
+(7.4069e+00 6.7924e-01 5.2512e-02 3.3510e-01 -1.1989e-03 1.3948e+00)
+(7.4897e+00 6.4224e-01 4.7330e-02 2.9050e-01 -9.5511e-04 1.3387e+00)
+(7.5016e+00 5.9969e-01 4.1511e-02 2.4803e-01 -6.7818e-04 1.2762e+00)
+(7.4274e+00 5.5162e-01 3.5418e-02 2.0824e-01 -3.9086e-04 1.2075e+00)
+(7.2524e+00 4.9837e-01 2.9501e-02 1.7162e-01 -1.3268e-04 1.1333e+00)
+(6.9642e+00 4.4073e-01 2.4133e-02 1.3854e-01 6.0490e-05 1.0545e+00)
+(6.5558e+00 3.8001e-01 1.9449e-02 1.0928e-01 1.7262e-04 9.7178e-01)
+(6.0286e+00 3.1808e-01 1.5387e-02 8.3971e-02 2.1029e-04 8.8610e-01)
+(5.3953e+00 2.5724e-01 1.1834e-02 6.2636e-02 1.9552e-04 7.9819e-01)
+(4.6811e+00 1.9999e-01 8.7562e-03 4.5157e-02 1.5460e-04 7.0871e-01)
+(3.9229e+00 1.4869e-01 6.1889e-03 3.1300e-02 1.0893e-04 6.1834e-01)
+(3.1644e+00 1.0513e-01 4.1721e-03 2.0724e-02 7.1044e-05 5.2804e-01)
+(2.4490e+00 7.0276e-02 2.6847e-03 1.3005e-02 4.4587e-05 4.3919e-01)
+(1.8132e+00 4.4135e-02 1.6529e-03 7.6610e-03 2.7764e-05 3.5369e-01)
+(1.2802e+00 2.5858e-02 9.7373e-04 4.1852e-03 1.7119e-05 2.7384e-01)
+(8.5851e-01 1.4015e-02 5.4723e-04 2.0879e-03 1.0091e-05 2.0209e-01)
+(5.4336e-01 6.9486e-03 2.9097e-04 9.3203e-04 5.4596e-06 1.4062e-01)
+(3.2123e-01 3.1009e-03 1.4506e-04 3.6204e-04 2.6288e-06 9.0947e-02)
+(1.7438e-01 1.2143e-03 6.6865e-05 1.1756e-04 1.0969e-06 5.3564e-02)
+(8.4454e-02 4.0027e-04 2.8297e-05 3.0014e-05 3.8272e-07 2.7852e-02)
+(3.4718e-02 1.0324e-04 1.0732e-05 5.4400e-06 1.0411e-07 1.2147e-02)
+(1.1020e-02 1.8113e-05 3.3616e-06 5.7878e-07 1.9051e-08 4.0392e-03)
+(2.1828e-03 1.5728e-06 8.3228e-07 2.3603e-08 1.7141e-09 8.2715e-04)
+(1.3684e-04 2.4414e-08 5.9464e-08 9.4734e-11 2.8525e-11 5.2943e-05)
+(5.3376e-28 -1.1206e-30 -6.1282e-30 4.9018e-29 -6.2793e-32 5.9575e-28)
+)
diff --git a/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.reducedDigitalFilter/0/UMean b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.reducedDigitalFilter/0/UMean
new file mode 100644
index 0000000000..2ecd2a6ad4
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.reducedDigitalFilter/0/UMean
@@ -0,0 +1,259 @@
+(
+( 0.0000e+00 0 0 )
+( 2.9538e-02 0 0 )
+( 1.1811e-01 0 0 )
+( 2.6562e-01 0 0 )
+( 4.7198e-01 0 0 )
+( 7.3701e-01 0 0 )
+( 1.0605e+00 0 0 )
+( 1.4417e+00 0 0 )
+( 1.8799e+00 0 0 )
+( 2.3729e+00 0 0 )
+( 2.9177e+00 0 0 )
+( 3.5093e+00 0 0 )
+( 4.1409e+00 0 0 )
+( 4.8032e+00 0 0 )
+( 5.4854e+00 0 0 )
+( 6.1754e+00 0 0 )
+( 6.8611e+00 0 0 )
+( 7.5309e+00 0 0 )
+( 8.1754e+00 0 0 )
+( 8.7870e+00 0 0 )
+( 9.3607e+00 0 0 )
+( 9.8937e+00 0 0 )
+( 1.0385e+01 0 0 )
+( 1.0836e+01 0 0 )
+( 1.1248e+01 0 0 )
+( 1.1624e+01 0 0 )
+( 1.1966e+01 0 0 )
+( 1.2278e+01 0 0 )
+( 1.2563e+01 0 0 )
+( 1.2822e+01 0 0 )
+( 1.3060e+01 0 0 )
+( 1.3278e+01 0 0 )
+( 1.3479e+01 0 0 )
+( 1.3664e+01 0 0 )
+( 1.3837e+01 0 0 )
+( 1.3998e+01 0 0 )
+( 1.4148e+01 0 0 )
+( 1.4290e+01 0 0 )
+( 1.4425e+01 0 0 )
+( 1.4552e+01 0 0 )
+( 1.4673e+01 0 0 )
+( 1.4790e+01 0 0 )
+( 1.4902e+01 0 0 )
+( 1.5011e+01 0 0 )
+( 1.5117e+01 0 0 )
+( 1.5221e+01 0 0 )
+( 1.5322e+01 0 0 )
+( 1.5421e+01 0 0 )
+( 1.5518e+01 0 0 )
+( 1.5614e+01 0 0 )
+( 1.5707e+01 0 0 )
+( 1.5799e+01 0 0 )
+( 1.5890e+01 0 0 )
+( 1.5979e+01 0 0 )
+( 1.6067e+01 0 0 )
+( 1.6153e+01 0 0 )
+( 1.6239e+01 0 0 )
+( 1.6324e+01 0 0 )
+( 1.6409e+01 0 0 )
+( 1.6493e+01 0 0 )
+( 1.6576e+01 0 0 )
+( 1.6659e+01 0 0 )
+( 1.6741e+01 0 0 )
+( 1.6823e+01 0 0 )
+( 1.6903e+01 0 0 )
+( 1.6984e+01 0 0 )
+( 1.7063e+01 0 0 )
+( 1.7141e+01 0 0 )
+( 1.7218e+01 0 0 )
+( 1.7294e+01 0 0 )
+( 1.7369e+01 0 0 )
+( 1.7443e+01 0 0 )
+( 1.7517e+01 0 0 )
+( 1.7590e+01 0 0 )
+( 1.7664e+01 0 0 )
+( 1.7738e+01 0 0 )
+( 1.7812e+01 0 0 )
+( 1.7886e+01 0 0 )
+( 1.7960e+01 0 0 )
+( 1.8034e+01 0 0 )
+( 1.8108e+01 0 0 )
+( 1.8182e+01 0 0 )
+( 1.8254e+01 0 0 )
+( 1.8326e+01 0 0 )
+( 1.8396e+01 0 0 )
+( 1.8466e+01 0 0 )
+( 1.8535e+01 0 0 )
+( 1.8603e+01 0 0 )
+( 1.8669e+01 0 0 )
+( 1.8734e+01 0 0 )
+( 1.8797e+01 0 0 )
+( 1.8859e+01 0 0 )
+( 1.8919e+01 0 0 )
+( 1.8978e+01 0 0 )
+( 1.9035e+01 0 0 )
+( 1.9092e+01 0 0 )
+( 1.9148e+01 0 0 )
+( 1.9202e+01 0 0 )
+( 1.9256e+01 0 0 )
+( 1.9308e+01 0 0 )
+( 1.9359e+01 0 0 )
+( 1.9408e+01 0 0 )
+( 1.9456e+01 0 0 )
+( 1.9503e+01 0 0 )
+( 1.9548e+01 0 0 )
+( 1.9593e+01 0 0 )
+( 1.9636e+01 0 0 )
+( 1.9678e+01 0 0 )
+( 1.9719e+01 0 0 )
+( 1.9758e+01 0 0 )
+( 1.9796e+01 0 0 )
+( 1.9832e+01 0 0 )
+( 1.9865e+01 0 0 )
+( 1.9897e+01 0 0 )
+( 1.9927e+01 0 0 )
+( 1.9955e+01 0 0 )
+( 1.9981e+01 0 0 )
+( 2.0004e+01 0 0 )
+( 2.0026e+01 0 0 )
+( 2.0046e+01 0 0 )
+( 2.0064e+01 0 0 )
+( 2.0080e+01 0 0 )
+( 2.0094e+01 0 0 )
+( 2.0106e+01 0 0 )
+( 2.0116e+01 0 0 )
+( 2.0123e+01 0 0 )
+( 2.0129e+01 0 0 )
+( 2.0132e+01 0 0 )
+( 2.0133e+01 0 0 )
+( 2.0132e+01 0 0 )
+( 2.0129e+01 0 0 )
+( 2.0123e+01 0 0 )
+( 2.0116e+01 0 0 )
+( 2.0106e+01 0 0 )
+( 2.0094e+01 0 0 )
+( 2.0080e+01 0 0 )
+( 2.0064e+01 0 0 )
+( 2.0046e+01 0 0 )
+( 2.0026e+01 0 0 )
+( 2.0004e+01 0 0 )
+( 1.9981e+01 0 0 )
+( 1.9955e+01 0 0 )
+( 1.9927e+01 0 0 )
+( 1.9897e+01 0 0 )
+( 1.9865e+01 0 0 )
+( 1.9832e+01 0 0 )
+( 1.9796e+01 0 0 )
+( 1.9758e+01 0 0 )
+( 1.9719e+01 0 0 )
+( 1.9678e+01 0 0 )
+( 1.9636e+01 0 0 )
+( 1.9593e+01 0 0 )
+( 1.9548e+01 0 0 )
+( 1.9503e+01 0 0 )
+( 1.9456e+01 0 0 )
+( 1.9408e+01 0 0 )
+( 1.9359e+01 0 0 )
+( 1.9308e+01 0 0 )
+( 1.9256e+01 0 0 )
+( 1.9202e+01 0 0 )
+( 1.9148e+01 0 0 )
+( 1.9092e+01 0 0 )
+( 1.9035e+01 0 0 )
+( 1.8978e+01 0 0 )
+( 1.8919e+01 0 0 )
+( 1.8859e+01 0 0 )
+( 1.8797e+01 0 0 )
+( 1.8734e+01 0 0 )
+( 1.8669e+01 0 0 )
+( 1.8603e+01 0 0 )
+( 1.8535e+01 0 0 )
+( 1.8466e+01 0 0 )
+( 1.8396e+01 0 0 )
+( 1.8326e+01 0 0 )
+( 1.8254e+01 0 0 )
+( 1.8182e+01 0 0 )
+( 1.8108e+01 0 0 )
+( 1.8034e+01 0 0 )
+( 1.7960e+01 0 0 )
+( 1.7886e+01 0 0 )
+( 1.7812e+01 0 0 )
+( 1.7738e+01 0 0 )
+( 1.7664e+01 0 0 )
+( 1.7590e+01 0 0 )
+( 1.7517e+01 0 0 )
+( 1.7443e+01 0 0 )
+( 1.7369e+01 0 0 )
+( 1.7294e+01 0 0 )
+( 1.7218e+01 0 0 )
+( 1.7141e+01 0 0 )
+( 1.7063e+01 0 0 )
+( 1.6984e+01 0 0 )
+( 1.6903e+01 0 0 )
+( 1.6823e+01 0 0 )
+( 1.6741e+01 0 0 )
+( 1.6659e+01 0 0 )
+( 1.6576e+01 0 0 )
+( 1.6493e+01 0 0 )
+( 1.6409e+01 0 0 )
+( 1.6324e+01 0 0 )
+( 1.6239e+01 0 0 )
+( 1.6153e+01 0 0 )
+( 1.6067e+01 0 0 )
+( 1.5979e+01 0 0 )
+( 1.5890e+01 0 0 )
+( 1.5799e+01 0 0 )
+( 1.5707e+01 0 0 )
+( 1.5614e+01 0 0 )
+( 1.5518e+01 0 0 )
+( 1.5421e+01 0 0 )
+( 1.5322e+01 0 0 )
+( 1.5221e+01 0 0 )
+( 1.5117e+01 0 0 )
+( 1.5011e+01 0 0 )
+( 1.4902e+01 0 0 )
+( 1.4790e+01 0 0 )
+( 1.4673e+01 0 0 )
+( 1.4552e+01 0 0 )
+( 1.4425e+01 0 0 )
+( 1.4290e+01 0 0 )
+( 1.4148e+01 0 0 )
+( 1.3998e+01 0 0 )
+( 1.3837e+01 0 0 )
+( 1.3664e+01 0 0 )
+( 1.3479e+01 0 0 )
+( 1.3278e+01 0 0 )
+( 1.3060e+01 0 0 )
+( 1.2822e+01 0 0 )
+( 1.2563e+01 0 0 )
+( 1.2278e+01 0 0 )
+( 1.1966e+01 0 0 )
+( 1.1624e+01 0 0 )
+( 1.1248e+01 0 0 )
+( 1.0836e+01 0 0 )
+( 1.0385e+01 0 0 )
+( 9.8937e+00 0 0 )
+( 9.3607e+00 0 0 )
+( 8.7870e+00 0 0 )
+( 8.1754e+00 0 0 )
+( 7.5309e+00 0 0 )
+( 6.8611e+00 0 0 )
+( 6.1754e+00 0 0 )
+( 5.4854e+00 0 0 )
+( 4.8032e+00 0 0 )
+( 4.1409e+00 0 0 )
+( 3.5093e+00 0 0 )
+( 2.9177e+00 0 0 )
+( 2.3729e+00 0 0 )
+( 1.8799e+00 0 0 )
+( 1.4417e+00 0 0 )
+( 1.0605e+00 0 0 )
+( 7.3701e-01 0 0 )
+( 4.7198e-01 0 0 )
+( 2.6562e-01 0 0 )
+( 1.1811e-01 0 0 )
+( 2.9538e-02 0 0 )
+( 0.0000e+00 0 0 )
+)
diff --git a/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.reducedDigitalFilter/points b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.reducedDigitalFilter/points
new file mode 100644
index 0000000000..5474eb4eae
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/constant/boundaryData/inlet.reducedDigitalFilter/points
@@ -0,0 +1,259 @@
+(
+(0 0.0000e+00 0)
+(0 7.5298e-05 0)
+(0 3.0118e-04 0)
+(0 6.7762e-04 0)
+(0 1.2045e-03 0)
+(0 1.8819e-03 0)
+(0 2.7095e-03 0)
+(0 3.6874e-03 0)
+(0 4.8153e-03 0)
+(0 6.0930e-03 0)
+(0 7.5205e-03 0)
+(0 9.0974e-03 0)
+(0 1.0823e-02 0)
+(0 1.2699e-02 0)
+(0 1.4722e-02 0)
+(0 1.6895e-02 0)
+(0 1.9215e-02 0)
+(0 2.1683e-02 0)
+(0 2.4298e-02 0)
+(0 2.7060e-02 0)
+(0 2.9969e-02 0)
+(0 3.3024e-02 0)
+(0 3.6224e-02 0)
+(0 3.9569e-02 0)
+(0 4.3060e-02 0)
+(0 4.6694e-02 0)
+(0 5.0472e-02 0)
+(0 5.4393e-02 0)
+(0 5.8456e-02 0)
+(0 6.2661e-02 0)
+(0 6.7007e-02 0)
+(0 7.1494e-02 0)
+(0 7.6120e-02 0)
+(0 8.0886e-02 0)
+(0 8.5790e-02 0)
+(0 9.0832e-02 0)
+(0 9.6011e-02 0)
+(0 1.0133e-01 0)
+(0 1.0678e-01 0)
+(0 1.1236e-01 0)
+(0 1.1808e-01 0)
+(0 1.2393e-01 0)
+(0 1.2991e-01 0)
+(0 1.3603e-01 0)
+(0 1.4227e-01 0)
+(0 1.4864e-01 0)
+(0 1.5515e-01 0)
+(0 1.6178e-01 0)
+(0 1.6853e-01 0)
+(0 1.7541e-01 0)
+(0 1.8242e-01 0)
+(0 1.8954e-01 0)
+(0 1.9679e-01 0)
+(0 2.0416e-01 0)
+(0 2.1165e-01 0)
+(0 2.1926e-01 0)
+(0 2.2699e-01 0)
+(0 2.3483e-01 0)
+(0 2.4279e-01 0)
+(0 2.5086e-01 0)
+(0 2.5905e-01 0)
+(0 2.6735e-01 0)
+(0 2.7575e-01 0)
+(0 2.8427e-01 0)
+(0 2.9289e-01 0)
+(0 3.0162e-01 0)
+(0 3.1046e-01 0)
+(0 3.1940e-01 0)
+(0 3.2844e-01 0)
+(0 3.3758e-01 0)
+(0 3.4683e-01 0)
+(0 3.5617e-01 0)
+(0 3.6561e-01 0)
+(0 3.7514e-01 0)
+(0 3.8477e-01 0)
+(0 3.9449e-01 0)
+(0 4.0430e-01 0)
+(0 4.1420e-01 0)
+(0 4.2419e-01 0)
+(0 4.3427e-01 0)
+(0 4.4443e-01 0)
+(0 4.5467e-01 0)
+(0 4.6500e-01 0)
+(0 4.7541e-01 0)
+(0 4.8590e-01 0)
+(0 4.9646e-01 0)
+(0 5.0710e-01 0)
+(0 5.1782e-01 0)
+(0 5.2860e-01 0)
+(0 5.3946e-01 0)
+(0 5.5039e-01 0)
+(0 5.6138e-01 0)
+(0 5.7244e-01 0)
+(0 5.8357e-01 0)
+(0 5.9476e-01 0)
+(0 6.0601e-01 0)
+(0 6.1732e-01 0)
+(0 6.2868e-01 0)
+(0 6.4010e-01 0)
+(0 6.5158e-01 0)
+(0 6.6311e-01 0)
+(0 6.7469e-01 0)
+(0 6.8632e-01 0)
+(0 6.9799e-01 0)
+(0 7.0972e-01 0)
+(0 7.2148e-01 0)
+(0 7.3329e-01 0)
+(0 7.4513e-01 0)
+(0 7.5702e-01 0)
+(0 7.6894e-01 0)
+(0 7.8090e-01 0)
+(0 7.9289e-01 0)
+(0 8.0491e-01 0)
+(0 8.1696e-01 0)
+(0 8.2904e-01 0)
+(0 8.4114e-01 0)
+(0 8.5327e-01 0)
+(0 8.6542e-01 0)
+(0 8.7759e-01 0)
+(0 8.8978e-01 0)
+(0 9.0198e-01 0)
+(0 9.1420e-01 0)
+(0 9.2644e-01 0)
+(0 9.3868e-01 0)
+(0 9.5093e-01 0)
+(0 9.6319e-01 0)
+(0 9.7546e-01 0)
+(0 9.8773e-01 0)
+(0 1.0000e-00 0)
+(0 1.01227 0)
+(0 1.02454 0)
+(0 1.03681 0)
+(0 1.04907 0)
+(0 1.06132 0)
+(0 1.07356 0)
+(0 1.0858 0)
+(0 1.09802 0)
+(0 1.11022 0)
+(0 1.12241 0)
+(0 1.13458 0)
+(0 1.14673 0)
+(0 1.15886 0)
+(0 1.17096 0)
+(0 1.18304 0)
+(0 1.19509 0)
+(0 1.20711 0)
+(0 1.2191 0)
+(0 1.23106 0)
+(0 1.24298 0)
+(0 1.25487 0)
+(0 1.26671 0)
+(0 1.27852 0)
+(0 1.29028 0)
+(0 1.30201 0)
+(0 1.31368 0)
+(0 1.32531 0)
+(0 1.33689 0)
+(0 1.34842 0)
+(0 1.3599 0)
+(0 1.37132 0)
+(0 1.38268 0)
+(0 1.39399 0)
+(0 1.40524 0)
+(0 1.41643 0)
+(0 1.42756 0)
+(0 1.43862 0)
+(0 1.44961 0)
+(0 1.46054 0)
+(0 1.4714 0)
+(0 1.48218 0)
+(0 1.4929 0)
+(0 1.50354 0)
+(0 1.5141 0)
+(0 1.52459 0)
+(0 1.535 0)
+(0 1.54533 0)
+(0 1.55557 0)
+(0 1.56573 0)
+(0 1.57581 0)
+(0 1.5858 0)
+(0 1.5957 0)
+(0 1.60551 0)
+(0 1.61523 0)
+(0 1.62486 0)
+(0 1.63439 0)
+(0 1.64383 0)
+(0 1.65317 0)
+(0 1.66242 0)
+(0 1.67156 0)
+(0 1.6806 0)
+(0 1.68954 0)
+(0 1.69838 0)
+(0 1.70711 0)
+(0 1.71573 0)
+(0 1.72425 0)
+(0 1.73265 0)
+(0 1.74095 0)
+(0 1.74914 0)
+(0 1.75721 0)
+(0 1.76517 0)
+(0 1.77301 0)
+(0 1.78074 0)
+(0 1.78835 0)
+(0 1.79584 0)
+(0 1.80321 0)
+(0 1.81046 0)
+(0 1.81758 0)
+(0 1.82459 0)
+(0 1.83147 0)
+(0 1.83822 0)
+(0 1.84485 0)
+(0 1.85136 0)
+(0 1.85773 0)
+(0 1.86397 0)
+(0 1.87009 0)
+(0 1.87607 0)
+(0 1.88192 0)
+(0 1.88764 0)
+(0 1.89322 0)
+(0 1.89867 0)
+(0 1.90399 0)
+(0 1.90917 0)
+(0 1.91421 0)
+(0 1.91911 0)
+(0 1.92388 0)
+(0 1.92851 0)
+(0 1.93299 0)
+(0 1.93734 0)
+(0 1.94154 0)
+(0 1.94561 0)
+(0 1.94953 0)
+(0 1.95331 0)
+(0 1.95694 0)
+(0 1.96043 0)
+(0 1.96378 0)
+(0 1.96698 0)
+(0 1.97003 0)
+(0 1.97294 0)
+(0 1.9757 0)
+(0 1.97832 0)
+(0 1.98079 0)
+(0 1.9831 0)
+(0 1.98528 0)
+(0 1.9873 0)
+(0 1.98918 0)
+(0 1.9909 0)
+(0 1.99248 0)
+(0 1.99391 0)
+(0 1.99518 0)
+(0 1.99631 0)
+(0 1.99729 0)
+(0 1.99812 0)
+(0 1.9988 0)
+(0 1.99932 0)
+(0 1.9997 0)
+(0 1.99992 0)
+(0 2 0)
+)
diff --git a/tutorials/verificationAndValidation/turbulentInflow/constant/transportProperties b/tutorials/verificationAndValidation/turbulentInflow/constant/transportProperties
new file mode 100644
index 0000000000..f97c34ebfc
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/constant/transportProperties
@@ -0,0 +1,31 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      transportProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+transportModel  Newtonian;
+
+// Re_tau = u_tau L / nu
+// Re_tau = 395
+// L = half channel height = 1
+// Ubulk/u_tau = 17.55
+// U_bulk = 17.55 -> u_tau = 1
+// -> nu = 1*1/395 = 2.532e-3
+
+nu              nu [ 0 2 -1 0 0 0 0 ] 2.532e-3;
+
+
+
+// ************************************************************************* //
diff --git a/tutorials/verificationAndValidation/turbulentInflow/constant/turbulenceProperties b/tutorials/verificationAndValidation/turbulentInflow/constant/turbulenceProperties
new file mode 100644
index 0000000000..ded27993b6
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/constant/turbulenceProperties
@@ -0,0 +1,48 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      turbulenceProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+simulationType LES;
+
+LES
+{
+    turbulence      on;
+
+    LESModel        Smagorinsky;
+    SmagorinskyCoeffs
+    {
+        Ce              1.05;
+        Ck              0.0266; // Updated to give Cs = 0.065
+    }
+
+    printCoeffs     on;
+
+    delta           vanDriest;
+
+    vanDriestCoeffs
+    {
+        delta           cubeRootVol;
+        cubeRootVolCoeffs
+        {
+            deltaCoeff      2;
+        }
+        Aplus           26;
+        Cdelta          0.158;
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/tutorials/verificationAndValidation/turbulentInflow/plot.cell b/tutorials/verificationAndValidation/turbulentInflow/plot.cell
new file mode 100755
index 0000000000..d5d6a716a5
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/plot.cell
@@ -0,0 +1,26 @@
+set terminal png size 1000,800 enhanced font "Helvetica,24"
+set output 'stress-cell.png'
+
+set xrange [0:1]
+set yrange [-1:8]
+set xlabel "Channel height"
+set ylabel "<u_i u_i>"
+
+set offset .05, .05
+set style data linespoints
+
+set grid
+set linetype 1 lc rgb 'black' lw 2
+set linetype 2 lc rgb 'red' lw 2
+set linetype 3 lc rgb 'blue' lw 2
+set linetype 4 lc rgb 'green' lw 2
+set linetype 5 lc rgb 'black' pi -8 pt 4 ps 1.5
+set linetype 6 lc rgb 'red' pi -8 pt 4 ps 1.5
+set linetype 7 lc rgb 'blue' pi -8 pt 4 ps 1.5
+set linetype 8 lc rgb 'green' pi -8 pt 4 ps 1.5
+
+plot \
+    "postProcessing/inletSampling/85/inletCell_UPrime2Mean.xy" u 1:2 w lines t "<uu>" lt 1, \
+    "postProcessing/inletSampling/85/inletCell_UPrime2Mean.xy" u 1:5 w lines t "<vv>" lt 2, \
+    "postProcessing/inletSampling/85/inletCell_UPrime2Mean.xy" u 1:7 w lines t "<ww>" lt 3, \
+    "postProcessing/inletSampling/85/inletCell_UPrime2Mean.xy" u 1:3 w lines t "<uv>" lt 4
diff --git a/tutorials/verificationAndValidation/turbulentInflow/plot.patch b/tutorials/verificationAndValidation/turbulentInflow/plot.patch
new file mode 100755
index 0000000000..e209ae43e3
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/plot.patch
@@ -0,0 +1,26 @@
+set terminal png size 1000,800 enhanced font "Helvetica,24"
+set output 'stress-patch.png'
+
+set xrange [0:1]
+set yrange [-1:8]
+set xlabel "Channel height"
+set ylabel "<u_i u_i>"
+
+set offset .05, .05
+set style data linespoints
+
+set grid
+set linetype 1 lc rgb 'black' lw 2
+set linetype 2 lc rgb 'red' lw 2
+set linetype 3 lc rgb 'blue' lw 2
+set linetype 4 lc rgb 'green' lw 2
+set linetype 5 lc rgb 'black' pi -8 pt 4 ps 1.5
+set linetype 6 lc rgb 'red' pi -8 pt 4 ps 1.5
+set linetype 7 lc rgb 'blue' pi -8 pt 4 ps 1.5
+set linetype 8 lc rgb 'green' pi -8 pt 4 ps 1.5
+
+plot \
+    "postProcessing/inletSampling/85/inletPatch_UPrime2Mean.xy" u 1:2 w lines t "<uu>" lt 1, \
+    "postProcessing/inletSampling/85/inletPatch_UPrime2Mean.xy" u 1:5 w lines t "<vv>" lt 2, \
+    "postProcessing/inletSampling/85/inletPatch_UPrime2Mean.xy" u 1:7 w lines t "<ww>" lt 3, \
+    "postProcessing/inletSampling/85/inletPatch_UPrime2Mean.xy" u 1:3 w lines t "<uv>" lt 4
diff --git a/tutorials/verificationAndValidation/turbulentInflow/system/blockMeshDict b/tutorials/verificationAndValidation/turbulentInflow/system/blockMeshDict
new file mode 100644
index 0000000000..ee1f25b003
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/system/blockMeshDict
@@ -0,0 +1,91 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      blockMeshDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+scale 1;
+
+L 0.125664;// "20*M_PI";
+H #calc "1";
+H2 #calc "2*$H";
+W #calc "M_PI";
+
+vertices
+(
+    ( 0   0   0)
+    ($L   0   0)
+    ($L  $H   0)
+    ($L $H2   0)
+    ( 0 $H2   0)
+    ( 0  $H   0)
+
+    ( 0   0  $W)
+    ($L   0  $W)
+    ($L  $H  $W)
+    ($L $H2  $W)
+    ( 0 $H2  $W)
+    ( 0  $H  $W)
+);
+
+blocks
+(
+    hex ( 0  1  2  5  6  7  8 11) (1 23 82) simpleGrading (1 15 1)
+    hex ( 5  2  3  4 11  8  9 10) (1 23 82) simpleGrading (1 0.0666 1)
+);
+
+edges
+(
+);
+
+boundary
+(
+    bottomWall
+    {
+        type            wall;
+        faces           ((0 6 7 1));
+    }
+    topWall
+    {
+        type            wall;
+        faces           ((4 3 9 10));
+    }
+    sides_half0
+    {
+        type            cyclic;
+        neighbourPatch  sides_half1;
+        faces           ((1 2 5 0)(2 3 4 5));
+    }
+    sides_half1
+    {
+        type            cyclic;
+        neighbourPatch  sides_half0;
+        faces           ((6 11 8 7)(11 10 9 8));
+    }
+    inlet
+    {
+        type            patch;
+        faces           ((0 5 11 6)(5 4 10 11));
+    }
+    outlet
+    {
+        type            patch;
+        faces           ((1 7 8 2)(2 8 9 3));
+    }
+);
+
+mergePatchPairs
+(
+);
+
+// ************************************************************************* //
diff --git a/tutorials/verificationAndValidation/turbulentInflow/system/controlDict b/tutorials/verificationAndValidation/turbulentInflow/system/controlDict
new file mode 100644
index 0000000000..498fc8b9f9
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/system/controlDict
@@ -0,0 +1,103 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      controlDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+application     pimpleFoam;
+
+startFrom       latestTime;
+
+startTime       0;
+
+stopAt          endTime;
+
+endTime         85;
+
+deltaT          4e-3;
+
+writeControl    timeStep;
+
+writeInterval   1250;
+
+purgeWrite      0;
+
+writeFormat     ascii;
+
+writePrecision  8;
+
+writeCompression off;
+
+timeFormat      general;
+
+timePrecision   8;
+
+runTimeModifiable false;
+
+adjustTimeStep false;
+
+functions
+{
+    fieldAverage1
+    {
+        type                fieldAverage;
+        libs                ("libfieldFunctionObjects.so");
+        enabled             true;
+        writeControl        writeTime;
+        timeStart           8.5;
+        restartOnRestart    false;
+        fields
+        (
+            U
+            {
+                mean        on;
+                prime2Mean  on;
+                base        time;
+            }
+        );
+    }
+
+    inletSampling
+    {
+        type                sets;
+        libs                ("libsampling.so");
+        enabled             true;
+        writeControl        writeTime;
+        timeStart           8.5;
+        restartOnRestart    false;
+        interpolationScheme cellPoint;
+        setFormat           raw;
+        fields              (UPrime2Mean);
+
+        sets
+        (
+            inletPatch
+            {
+                type        face;
+                axis        y;
+                start       (0.0 0 1.57);
+                end         (0.0 2 1.57);
+            }
+            inletCell
+            {
+                type        face;
+                axis        y;
+                start       (0.062832 0 1.57);
+                end         (0.062832 2 1.57);
+            }
+        );
+    }
+}
+
+// ************************************************************************* //
diff --git a/tutorials/verificationAndValidation/turbulentInflow/system/decomposeParDict b/tutorials/verificationAndValidation/turbulentInflow/system/decomposeParDict
new file mode 100644
index 0000000000..58f3c532be
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/system/decomposeParDict
@@ -0,0 +1,22 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    note        "mesh decomposition control dictionary";
+    object      decomposeParDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+numberOfSubdomains  6;
+
+method         scotch;
+
+// ************************************************************************* //
diff --git a/tutorials/verificationAndValidation/turbulentInflow/system/fvSchemes b/tutorials/verificationAndValidation/turbulentInflow/system/fvSchemes
new file mode 100644
index 0000000000..71463c82da
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/system/fvSchemes
@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSchemes;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+ddtSchemes
+{
+    default         backward;
+}
+
+gradSchemes
+{
+    default         leastSquares;
+}
+
+divSchemes
+{
+    default         none;
+    div(phi,U)      Gauss linear;
+    div(phi,k)      Gauss limitedLinear 0.1;
+    div(phi,B)      Gauss limitedLinear 0.1;
+    div(B)          Gauss linear;
+    div(phi,nuTilda) Gauss limitedLinear 0.1;
+    div((nuEff*dev2(T(grad(U))))) Gauss linear;
+}
+
+laplacianSchemes
+{
+    default         Gauss linear corrected;
+}
+
+interpolationSchemes
+{
+    default         linear;
+}
+
+snGradSchemes
+{
+    default         corrected;
+}
+
+// ************************************************************************* //
diff --git a/tutorials/verificationAndValidation/turbulentInflow/system/fvSolution b/tutorials/verificationAndValidation/turbulentInflow/system/fvSolution
new file mode 100644
index 0000000000..0a8dcda666
--- /dev/null
+++ b/tutorials/verificationAndValidation/turbulentInflow/system/fvSolution
@@ -0,0 +1,66 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  v1906                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSolution;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+solvers
+{
+    p
+    {
+        solver          GAMG;
+        tolerance       0;
+        relTol          0.01;
+        smoother        DICGaussSeidel;
+        nPreSweeps      0;
+        nPostSweeps     2;
+        cacheAgglomeration true;
+        nCellsInCoarsestLevel 10;
+        agglomerator    faceAreaPair;
+        mergeLevels     1;
+    }
+
+    pFinal
+    {
+        $p;
+        smoother        DICGaussSeidel;
+        tolerance       1e-06;
+        relTol          0;
+    }
+
+    "(U|k)"
+    {
+        solver          PBiCG;
+        preconditioner  DILU;
+        tolerance       1e-05;
+        relTol          0.1;
+    }
+
+    "(U|k)Final"
+    {
+        $U;
+        tolerance       1e-06;
+        relTol          0;
+    }
+}
+
+PIMPLE
+{
+    nOuterCorrectors 2;
+    nCorrectors     1;
+    nNonOrthogonalCorrectors 0;
+}
+
+
+// ************************************************************************* //