From 0db9fa570258e9d8c6ea3bdae1ef16cbec5aba90 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Tue, 9 Nov 2010 18:04:32 +0000
Subject: [PATCH] ENH: Ensight output for sampledSurfaces

---
 .../postProcessing/sampling/sample/sampleDict |   1 +
 src/sampling/Make/files                       |   1 +
 .../sampledSurfacesTemplates.C                |   6 +-
 .../writers/dx/dxSurfaceWriter.C              |   3 +-
 .../writers/dx/dxSurfaceWriter.H              |   1 +
 .../writers/ensight/ensightSurfaceWriter.C    | 372 ++++++++++++++++++
 .../writers/ensight/ensightSurfaceWriter.H    | 138 +++++++
 .../ensight/ensightSurfaceWriterRunTime.C     |  44 +++
 .../writers/foamFile/foamFileSurfaceWriter.C  |   1 +
 .../writers/foamFile/foamFileSurfaceWriter.H  |   1 +
 .../writers/null/nullSurfaceWriter.C          |   1 +
 .../writers/null/nullSurfaceWriter.H          |   1 +
 .../writers/proxy/proxySurfaceWriter.H        |   1 +
 .../writers/raw/rawSurfaceWriter.C            |  21 +-
 .../writers/raw/rawSurfaceWriter.H            |   6 +
 .../sampledSurface/writers/surfaceWriter.H    |   3 +-
 .../writers/vtk/vtkSurfaceWriter.C            |   3 +-
 .../writers/vtk/vtkSurfaceWriter.H            |   1 +
 18 files changed, 593 insertions(+), 12 deletions(-)
 create mode 100644 src/sampling/sampledSurface/writers/ensight/ensightSurfaceWriter.C
 create mode 100644 src/sampling/sampledSurface/writers/ensight/ensightSurfaceWriter.H
 create mode 100644 src/sampling/sampledSurface/writers/ensight/ensightSurfaceWriterRunTime.C

diff --git a/applications/utilities/postProcessing/sampling/sample/sampleDict b/applications/utilities/postProcessing/sampling/sample/sampleDict
index f48237e9ed..156b7cf54f 100644
--- a/applications/utilities/postProcessing/sampling/sample/sampleDict
+++ b/applications/utilities/postProcessing/sampling/sample/sampleDict
@@ -27,6 +27,7 @@ setFormat raw;
 
 // Surface output format. Choice of
 //      null        : suppress output
+//      ensight     : Ensight Gold format, one field per case file
 //      foamFile    : separate points, faces and values file
 //      dx          : DX scalar or vector format
 //      vtk         : VTK ascii format
diff --git a/src/sampling/Make/files b/src/sampling/Make/files
index c385e9aa93..6fee75b8f7 100644
--- a/src/sampling/Make/files
+++ b/src/sampling/Make/files
@@ -46,6 +46,7 @@ surfWriters = sampledSurface/writers
 
 $(surfWriters)/surfaceWriters.C
 $(surfWriters)/dx/dxSurfaceWriterRunTime.C
+$(surfWriters)/ensight/ensightSurfaceWriterRunTime.C
 $(surfWriters)/foamFile/foamFileSurfaceWriterRunTime.C
 $(surfWriters)/null/nullSurfaceWriterRunTime.C
 $(surfWriters)/proxy/proxySurfaceWriterRunTime.C
diff --git a/src/sampling/sampledSurface/sampledSurfaces/sampledSurfacesTemplates.C b/src/sampling/sampledSurface/sampledSurfaces/sampledSurfacesTemplates.C
index 4dd5deffa5..128f7d2606 100644
--- a/src/sampling/sampledSurface/sampledSurfaces/sampledSurfacesTemplates.C
+++ b/src/sampling/sampledSurface/sampledSurfaces/sampledSurfacesTemplates.C
@@ -103,7 +103,8 @@ void Foam::sampledSurfaces::sampleAndWrite
                         mergeList_[surfI].points,
                         mergeList_[surfI].faces,
                         fieldName,
-                        allValues
+                        allValues,
+                        s.interpolate()
                     );
                 }
             }
@@ -121,7 +122,8 @@ void Foam::sampledSurfaces::sampleAndWrite
                     s.points(),
                     s.faces(),
                     fieldName,
-                    values
+                    values,
+                    s.interpolate()
                 );
             }
         }
diff --git a/src/sampling/sampledSurface/writers/dx/dxSurfaceWriter.C b/src/sampling/sampledSurface/writers/dx/dxSurfaceWriter.C
index 4950e37664..f4e7ebe216 100644
--- a/src/sampling/sampledSurface/writers/dx/dxSurfaceWriter.C
+++ b/src/sampling/sampledSurface/writers/dx/dxSurfaceWriter.C
@@ -250,6 +250,7 @@ void Foam::dxSurfaceWriter<Type>::write
     const faceList& faces,
     const word& fieldName,
     const Field<Type>& values,
+    const bool isNodeValues,
     const bool verbose
 ) const
 {
@@ -272,7 +273,7 @@ void Foam::dxSurfaceWriter<Type>::write
 
     writeData(os, values);
 
-    if (values.size() == points.size())
+    if (isNodeValues)
     {
         os  << nl << "attribute \"dep\" string \"positions\""
             << nl << nl;
diff --git a/src/sampling/sampledSurface/writers/dx/dxSurfaceWriter.H b/src/sampling/sampledSurface/writers/dx/dxSurfaceWriter.H
index 897133a005..fc996ee3d9 100644
--- a/src/sampling/sampledSurface/writers/dx/dxSurfaceWriter.H
+++ b/src/sampling/sampledSurface/writers/dx/dxSurfaceWriter.H
@@ -88,6 +88,7 @@ public:
             const faceList& faces,
             const word& fieldName,
             const Field<Type>& values,
+            const bool isNodeValues,
             const bool verbose = false
         ) const;
 };
diff --git a/src/sampling/sampledSurface/writers/ensight/ensightSurfaceWriter.C b/src/sampling/sampledSurface/writers/ensight/ensightSurfaceWriter.C
new file mode 100644
index 0000000000..81eb8510fa
--- /dev/null
+++ b/src/sampling/sampledSurface/writers/ensight/ensightSurfaceWriter.C
@@ -0,0 +1,372 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2010-2010 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 "ensightSurfaceWriter.H"
+
+#include "OFstream.H"
+#include "OSspecific.H"
+#include "IOmanip.H"
+
+// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
+
+template<class Type>
+void Foam::ensightSurfaceWriter<Type>::binShapes
+(
+    const pointField& points,
+    const faceList& faces,
+    DynamicList<label>& tris,
+    DynamicList<label>& quads,
+    DynamicList<label>& polys
+)
+{
+    tris.setCapacity(faces.size());
+    quads.setCapacity(faces.size());
+    polys.setCapacity(faces.size());
+
+    forAll(faces, i)
+    {
+        const face& f = faces[i];
+        if (f.size() == 3)
+        {
+            tris.append(i);
+        }
+        else if (f.size() == 4)
+        {
+            quads.append(i);
+        }
+        else
+        {
+            polys.append(i);
+        }
+    }
+}
+
+
+template<class Type>
+void Foam::ensightSurfaceWriter<Type>::writeGeometry
+(
+    const fileName& surfaceName,
+    Ostream& os,
+    const pointField& points,
+    const faceList& faces,
+    const DynamicList<label>& tris,
+    const DynamicList<label>& quads,
+    const DynamicList<label>& polys
+)
+{
+    os  << "EnSight Geometry File" << nl
+        << (string("written by OpenFOAM-") + Foam::FOAMversion).c_str() << nl
+        << "node id assign" << nl
+        << "element id assign" << nl
+        << "part" << nl
+        << setw(10) << 1 << nl
+        << surfaceName.c_str() << nl;
+
+    os  << "coordinates" << nl
+        << setw(10) << points.size() << nl;
+    for (direction cmpt = 0; cmpt < vector::nComponents; cmpt++)
+    {
+        scalarField v = points.component(cmpt);
+        forAll(v, i)
+        {
+            os  << v[i] << nl;
+        }
+    }
+
+
+    if (tris.size())
+    {
+        os  << "tria3" << nl
+            << setw(10) << tris.size() << nl;
+        forAll(tris, i)
+        {
+            const face& f = faces[tris[i]];
+            forAll(f, fp)
+            {
+                os  << setw(10) << f[fp]+1;
+            }
+            os  << nl;
+        }
+    }
+    if (quads.size())
+    {
+        os  << "quad4" << nl
+            << setw(10) << quads.size() << nl;
+        forAll(quads, i)
+        {
+            const face& f = faces[quads[i]];
+            forAll(f, fp)
+            {
+                os  << setw(10) << f[fp]+1;
+            }
+            os  << nl;
+        }
+    }
+    if (polys.size())
+    {
+        os  << "nsided" << nl
+            << setw(10) << polys.size() << nl;
+        forAll(polys, i)
+        {
+            const face& f = faces[polys[i]];
+            forAll(f, fp)
+            {
+                os  << setw(10) << f[fp]+1;
+            }
+            os  << nl;
+        }
+    }
+}
+
+
+namespace Foam
+{
+
+    // Write scalarField in ensight format
+    template<>
+    void Foam::ensightSurfaceWriter<Foam::scalar>::writeData
+    (
+        Ostream& os,
+        const Field<Foam::scalar>& values
+    )
+    {
+        forAll(values, i)
+        {
+            os << values[i] << nl;
+        }
+    }
+
+
+    // Write booField in ensight format
+    template<>
+    void Foam::ensightSurfaceWriter<bool>::writeData
+    (
+        Ostream& os,
+        const Field<bool>& values
+    )
+    {
+        forAll(values, i)
+        {
+            os << values[i] << nl;
+        }
+    }
+}
+
+
+// Write generic field in ensight format
+template<class Type>
+void Foam::ensightSurfaceWriter<Type>::writeData
+(
+    Ostream& os,
+    const Field<Type>& values
+)
+{
+    for (direction cmpt = 0; cmpt < vector::nComponents; cmpt++)
+    {
+        scalarField v = values.component(cmpt);
+        forAll(v, i)
+        {
+            os << v[i] << nl;
+        }
+    }
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+// Construct from components
+template<class Type>
+Foam::ensightSurfaceWriter<Type>::ensightSurfaceWriter()
+:
+    surfaceWriter<Type>()
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+template<class Type>
+Foam::ensightSurfaceWriter<Type>::~ensightSurfaceWriter()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+template<class Type>
+void Foam::ensightSurfaceWriter<Type>::write
+(
+    const fileName& outputDir,
+    const fileName& surfaceName,
+    const pointField& points,
+    const faceList& faces,
+    const bool verbose
+) const
+{
+    if (!isDir(outputDir))
+    {
+        mkDir(outputDir);
+    }
+
+    //const scalar timeValue = Foam::name(this->mesh().time().timeValue());
+    const scalar timeValue = 0.0;
+
+
+    OFstream caseStr(outputDir/surfaceName + ".case");
+    OFstream geomStr(outputDir/surfaceName + ".***.mesh");
+
+    if (verbose)
+    {
+        Info<< "Writing case file to " << caseStr.name() << endl;
+    }
+
+    caseStr
+        << "FORMAT" << nl
+        << "type: ensight gold" << nl
+        << nl
+        << "GEOMETRY" << nl
+        << "model:        1     " << geomStr.name().name() << nl
+        << nl
+        << "TIME" << nl
+        << "time set:                      1" << nl
+        << "number of steps:               1" << nl
+        << "filename start number:         0" << nl
+        << "filename increment:            1" << nl
+        << "time values:" << nl
+        << timeValue << nl
+        << nl;
+
+    DynamicList<label> tris;
+    DynamicList<label> quads;
+    DynamicList<label> polys;
+    binShapes(points, faces, tris, quads, polys);
+
+    writeGeometry(surfaceName, geomStr, points, faces, tris, quads, polys);
+}
+
+
+template<class Type>
+void Foam::ensightSurfaceWriter<Type>::write
+(
+    const fileName& outputDir,
+    const fileName& surfaceName,
+    const pointField& points,
+    const faceList& faces,
+    const word& fieldName,
+    const Field<Type>& values,
+    const bool isNodeValues,
+    const bool verbose
+) const
+{
+    if (!isDir(outputDir/fieldName))
+    {
+        mkDir(outputDir/fieldName);
+    }
+
+    //const scalar timeValue = Foam::name(this->mesh().time().timeValue());
+    const scalar timeValue = 0.0;
+
+
+    OFstream caseStr(outputDir/fieldName/surfaceName + ".case");
+    OFstream geomStr(outputDir/fieldName/surfaceName + ".000.mesh");
+    OFstream fieldStr(outputDir/fieldName/surfaceName + ".000." + fieldName);
+
+    if (verbose)
+    {
+        Info<< "Writing case file to " << caseStr.name() << endl;
+    }
+
+    caseStr
+        << "FORMAT" << nl
+        << "type: ensight gold" << nl
+        << nl
+        << "GEOMETRY" << nl
+        << "model:        1     " << geomStr.name().name() << nl
+        << nl
+        << "VARIABLE" << nl;
+    if (isNodeValues)
+    {
+        caseStr
+            << pTraits<Type>::typeName << " per node:" << setw(10) << 1
+        << "       " << fieldName
+        << "       " << surfaceName.c_str() << ".***." << fieldName << nl;
+    }
+    else
+    {
+        caseStr
+            << pTraits<Type>::typeName << " per element:" << setw(10) << 1
+        << "       " << fieldName
+        << "       " << surfaceName.c_str() << ".***." << fieldName << nl;
+    }
+
+    caseStr
+        << nl
+        << "TIME" << nl
+        << "time set:                      1" << nl
+        << "number of steps:               1" << nl
+        << "filename start number:         0" << nl
+        << "filename increment:            1" << nl
+        << "time values:" << nl
+        << timeValue << nl
+        << nl;
+
+    DynamicList<label> tris;
+    DynamicList<label> quads;
+    DynamicList<label> polys;
+    binShapes(points, faces, tris, quads, polys);
+
+    writeGeometry(surfaceName, geomStr, points, faces, tris, quads, polys);
+
+    // Write field
+    fieldStr
+        << pTraits<Type>::typeName << nl
+        << "part" << nl
+        << setw(10) << 1 << nl;
+
+    if (isNodeValues)
+    {
+        fieldStr << "coordinates" << nl;
+        writeData(fieldStr, values);
+    }
+    else
+    {
+        if (tris.size())
+        {
+            fieldStr << "tria3" << nl;
+            writeData(fieldStr, Field<Type>(values, tris));
+        }
+        if (quads.size())
+        {
+            fieldStr << "quad4" << nl;
+            writeData(fieldStr, Field<Type>(values, quads));
+        }
+        if (polys.size())
+        {
+            fieldStr << "nsided" << nl;
+            writeData(fieldStr, Field<Type>(values, polys));
+        }
+    }
+}
+
+
+// ************************************************************************* //
diff --git a/src/sampling/sampledSurface/writers/ensight/ensightSurfaceWriter.H b/src/sampling/sampledSurface/writers/ensight/ensightSurfaceWriter.H
new file mode 100644
index 0000000000..adb974f133
--- /dev/null
+++ b/src/sampling/sampledSurface/writers/ensight/ensightSurfaceWriter.H
@@ -0,0 +1,138 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2010-2010 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::ensightSurfaceWriter
+
+Description
+
+SourceFiles
+    ensightSurfaceWriter.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef ensightSurfaceWriter_H
+#define ensightSurfaceWriter_H
+
+#include "surfaceWriter.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+/*---------------------------------------------------------------------------*\
+                     Class ensightSurfaceWriter Declaration
+\*---------------------------------------------------------------------------*/
+
+template<class Type>
+class ensightSurfaceWriter
+:
+    public surfaceWriter<Type>
+{
+    // Private Member Functions
+
+        static void binShapes
+        (
+            const pointField& points,
+            const faceList& faces,
+            DynamicList<label>& tris,
+            DynamicList<label>& quads,
+            DynamicList<label>& polys
+        );
+
+        static void writeGeometry
+        (
+            const fileName&,
+            Ostream&,
+            const pointField&,
+            const faceList&,
+            const DynamicList<label>&,
+            const DynamicList<label>&,
+            const DynamicList<label>&
+        );
+
+        static void writeData(Ostream&, const Field<Type>& values);
+
+
+public:
+
+    //- Runtime type information
+    TypeName("ensight");
+
+
+    // Constructors
+
+        //- Construct null
+        ensightSurfaceWriter();
+
+
+    //- Destructor
+    virtual ~ensightSurfaceWriter();
+
+
+    // Member Functions
+
+        // Write
+
+        //- Write geometry to file.
+        virtual void write
+        (
+            const fileName& outputDir,
+            const fileName& surfaceName,
+            const pointField& points,
+            const faceList& faces,
+            const bool verbose = false
+        ) const;
+
+
+        //- Writes single surface to file.
+        virtual void write
+        (
+            const fileName& outputDir,
+            const fileName& surfaceName,
+            const pointField& points,
+            const faceList& faces,
+            const word& fieldName,
+            const Field<Type>& values,
+            const bool isNodeValues,
+            const bool verbose = false
+        ) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#ifdef NoRepository
+#   include "ensightSurfaceWriter.C"
+#endif
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/sampling/sampledSurface/writers/ensight/ensightSurfaceWriterRunTime.C b/src/sampling/sampledSurface/writers/ensight/ensightSurfaceWriterRunTime.C
new file mode 100644
index 0000000000..8355a6d991
--- /dev/null
+++ b/src/sampling/sampledSurface/writers/ensight/ensightSurfaceWriterRunTime.C
@@ -0,0 +1,44 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2010-2010 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 "ensightSurfaceWriter.H"
+#include "surfaceWriters.H"
+#include "addToRunTimeSelectionTable.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+makeSurfaceWriterType(ensightSurfaceWriter, bool);
+makeSurfaceWriters(ensightSurfaceWriter);
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// ************************************************************************* //
diff --git a/src/sampling/sampledSurface/writers/foamFile/foamFileSurfaceWriter.C b/src/sampling/sampledSurface/writers/foamFile/foamFileSurfaceWriter.C
index 887abdd7e4..ea121695bb 100644
--- a/src/sampling/sampledSurface/writers/foamFile/foamFileSurfaceWriter.C
+++ b/src/sampling/sampledSurface/writers/foamFile/foamFileSurfaceWriter.C
@@ -96,6 +96,7 @@ void Foam::foamFileSurfaceWriter<Type>::write
     const faceList& faces,
     const word& fieldName,
     const Field<Type>& values,
+    const bool isNodeValues,
     const bool verbose
 ) const
 {
diff --git a/src/sampling/sampledSurface/writers/foamFile/foamFileSurfaceWriter.H b/src/sampling/sampledSurface/writers/foamFile/foamFileSurfaceWriter.H
index 8e1fb57062..b77c91388b 100644
--- a/src/sampling/sampledSurface/writers/foamFile/foamFileSurfaceWriter.H
+++ b/src/sampling/sampledSurface/writers/foamFile/foamFileSurfaceWriter.H
@@ -95,6 +95,7 @@ public:
             const faceList& faces,
             const word& fieldName,
             const Field<Type>& values,
+            const bool isNodeValues,
             const bool verbose = false
         ) const;
 };
diff --git a/src/sampling/sampledSurface/writers/null/nullSurfaceWriter.C b/src/sampling/sampledSurface/writers/null/nullSurfaceWriter.C
index 6ae7a20933..bfda994af5 100644
--- a/src/sampling/sampledSurface/writers/null/nullSurfaceWriter.C
+++ b/src/sampling/sampledSurface/writers/null/nullSurfaceWriter.C
@@ -54,6 +54,7 @@ void Foam::nullSurfaceWriter<Type>::write
     const faceList& faces,
     const word& fieldName,
     const Field<Type>& values,
+    const bool isNodeValues,
     const bool verbose
 ) const
 {}
diff --git a/src/sampling/sampledSurface/writers/null/nullSurfaceWriter.H b/src/sampling/sampledSurface/writers/null/nullSurfaceWriter.H
index 6b8b3f1fde..8ea0f6a847 100644
--- a/src/sampling/sampledSurface/writers/null/nullSurfaceWriter.H
+++ b/src/sampling/sampledSurface/writers/null/nullSurfaceWriter.H
@@ -80,6 +80,7 @@ public:
             const faceList& faces,
             const word& fieldName,
             const Field<Type>& values,
+            const bool isNodeValues,
             const bool verbose = false
         ) const;
 };
diff --git a/src/sampling/sampledSurface/writers/proxy/proxySurfaceWriter.H b/src/sampling/sampledSurface/writers/proxy/proxySurfaceWriter.H
index 3fdc991d54..426816b32c 100644
--- a/src/sampling/sampledSurface/writers/proxy/proxySurfaceWriter.H
+++ b/src/sampling/sampledSurface/writers/proxy/proxySurfaceWriter.H
@@ -102,6 +102,7 @@ public:
             const faceList& faces,
             const word& fieldName,
             const Field<Type>& values,
+            const bool isNodeValues,
             const bool verbose = false
         ) const
         {}
diff --git a/src/sampling/sampledSurface/writers/raw/rawSurfaceWriter.C b/src/sampling/sampledSurface/writers/raw/rawSurfaceWriter.C
index 1dc8ad80a5..a3a7832e4f 100644
--- a/src/sampling/sampledSurface/writers/raw/rawSurfaceWriter.C
+++ b/src/sampling/sampledSurface/writers/raw/rawSurfaceWriter.C
@@ -68,6 +68,7 @@ void Foam::rawSurfaceWriter<Type>::writeData
     const pointField& points,
     const faceList& faces,
     const scalarField& values,
+    const bool isNodeValues,
     Ostream& os
 )
 {
@@ -75,7 +76,7 @@ void Foam::rawSurfaceWriter<Type>::writeData
     os  << "#  x  y  z  " << fieldName << endl;
 
     // Write data
-    if (values.size() == points.size())
+    if (isNodeValues)
     {
         forAll(values, elemI)
         {
@@ -104,6 +105,7 @@ void Foam::rawSurfaceWriter<Type>::writeData
     const pointField& points,
     const faceList& faces,
     const vectorField& values,
+    const bool isNodeValues,
     Ostream& os
 )
 {
@@ -115,7 +117,7 @@ void Foam::rawSurfaceWriter<Type>::writeData
         << endl;
 
     // Write data
-    if (values.size() == points.size())
+    if (isNodeValues)
     {
         forAll(values, elemI)
         {
@@ -147,6 +149,7 @@ void Foam::rawSurfaceWriter<Type>::writeData
     const pointField& points,
     const faceList& faces,
     const sphericalTensorField& values,
+    const bool isNodeValues,
     Ostream& os
 )
 {
@@ -155,7 +158,7 @@ void Foam::rawSurfaceWriter<Type>::writeData
     os << fieldName << "_ii" << endl;
 
     // Write data
-    if (values.size() == points.size())
+    if (isNodeValues)
     {
         forAll(values, elemI)
         {
@@ -186,6 +189,7 @@ void Foam::rawSurfaceWriter<Type>::writeData
     const pointField& points,
     const faceList& faces,
     const symmTensorField& values,
+    const bool isNodeValues,
     Ostream& os
 )
 {
@@ -198,7 +202,7 @@ void Foam::rawSurfaceWriter<Type>::writeData
     os << endl;
 
     // Write data
-    if (values.size() == points.size())
+    if (isNodeValues)
     {
         forAll(values, elemI)
         {
@@ -233,6 +237,7 @@ void Foam::rawSurfaceWriter<Type>::writeData
     const pointField& points,
     const faceList& faces,
     const tensorField& values,
+    const bool isNodeValues,
     Ostream& os
 )
 {
@@ -245,7 +250,7 @@ void Foam::rawSurfaceWriter<Type>::writeData
     os << endl;
 
     // Write data
-    if (values.size() == points.size())
+    if (isNodeValues)
     {
         forAll(values, elemI)
         {
@@ -343,6 +348,7 @@ namespace Foam
         const faceList& faces,
         const word& fieldName,
         const Field<bool>& values,
+        const bool isNodeValues,
         const bool verbose
     ) const
     {}
@@ -358,6 +364,7 @@ void Foam::rawSurfaceWriter<Type>::write
     const faceList& faces,
     const word& fieldName,
     const Field<Type>& values,
+    const bool isNodeValues,
     const bool verbose
 ) const
 {
@@ -379,7 +386,7 @@ void Foam::rawSurfaceWriter<Type>::write
 
     // header
     os  << "# " << fieldName;
-    if (values.size() == points.size())
+    if (isNodeValues)
     {
         os  << "  POINT_DATA ";
     }
@@ -390,7 +397,7 @@ void Foam::rawSurfaceWriter<Type>::write
 
     os  << values.size() << nl;
 
-    writeData(fieldName, points, faces, values, os);
+    writeData(fieldName, points, faces, values, isNodeValues, os);
 }
 
 
diff --git a/src/sampling/sampledSurface/writers/raw/rawSurfaceWriter.H b/src/sampling/sampledSurface/writers/raw/rawSurfaceWriter.H
index 0fc4478f42..a7dee21262 100644
--- a/src/sampling/sampledSurface/writers/raw/rawSurfaceWriter.H
+++ b/src/sampling/sampledSurface/writers/raw/rawSurfaceWriter.H
@@ -73,6 +73,7 @@ class rawSurfaceWriter
             const pointField& points,
             const faceList& faces,
             const scalarField& values,
+            const bool isNodeValues,
             Ostream& os
         );
 
@@ -82,6 +83,7 @@ class rawSurfaceWriter
             const pointField& points,
             const faceList& faces,
             const vectorField& values,
+            const bool isNodeValues,
             Ostream& os
         );
 
@@ -91,6 +93,7 @@ class rawSurfaceWriter
             const pointField& points,
             const faceList& faces,
             const sphericalTensorField& values,
+            const bool isNodeValues,
             Ostream& os
         );
 
@@ -100,6 +103,7 @@ class rawSurfaceWriter
             const pointField& points,
             const faceList& faces,
             const symmTensorField& values,
+            const bool isNodeValues,
             Ostream& os
         );
 
@@ -109,6 +113,7 @@ class rawSurfaceWriter
             const pointField& points,
             const faceList& faces,
             const tensorField& values,
+            const bool isNodeValues,
             Ostream& os
         );
 
@@ -152,6 +157,7 @@ public:
             const faceList& faces,
             const word& fieldName,
             const Field<Type>& values,
+            const bool isNodeValues,
             const bool verbose = false
         ) const;
 };
diff --git a/src/sampling/sampledSurface/writers/surfaceWriter.H b/src/sampling/sampledSurface/writers/surfaceWriter.H
index 6fcf0f0133..c3f138ad0d 100644
--- a/src/sampling/sampledSurface/writers/surfaceWriter.H
+++ b/src/sampling/sampledSurface/writers/surfaceWriter.H
@@ -119,7 +119,7 @@ public:
 
 
         //- Writes single surface to file. Either one value per vertex or
-        //  one value per face (detected by values.size()==faces.size())
+        //  one value per face (isNodeValues = false)
         virtual void write
         (
             const fileName& outputDir,      // <root>/<case>/surface/TIME
@@ -128,6 +128,7 @@ public:
             const faceList& faces,
             const word& fieldName,          // name of field
             const Field<Type>& values,
+            const bool isNodeValues,
             const bool verbose = false
         ) const = 0;
 };
diff --git a/src/sampling/sampledSurface/writers/vtk/vtkSurfaceWriter.C b/src/sampling/sampledSurface/writers/vtk/vtkSurfaceWriter.C
index 61da515b40..c1aa725c13 100644
--- a/src/sampling/sampledSurface/writers/vtk/vtkSurfaceWriter.C
+++ b/src/sampling/sampledSurface/writers/vtk/vtkSurfaceWriter.C
@@ -267,6 +267,7 @@ void Foam::vtkSurfaceWriter<Type>::write
     const faceList& faces,
     const word& fieldName,
     const Field<Type>& values,
+    const bool isNodeValues,
     const bool verbose
 ) const
 {
@@ -288,7 +289,7 @@ void Foam::vtkSurfaceWriter<Type>::write
     writeGeometry(os, points, faces);
 
     // start writing data
-    if (values.size() == points.size())
+    if (isNodeValues)
     {
         os  << "POINT_DATA ";
     }
diff --git a/src/sampling/sampledSurface/writers/vtk/vtkSurfaceWriter.H b/src/sampling/sampledSurface/writers/vtk/vtkSurfaceWriter.H
index 7524cb963c..1d7b005d5d 100644
--- a/src/sampling/sampledSurface/writers/vtk/vtkSurfaceWriter.H
+++ b/src/sampling/sampledSurface/writers/vtk/vtkSurfaceWriter.H
@@ -97,6 +97,7 @@ public:
             const faceList& faces,
             const word& fieldName,
             const Field<Type>& values,
+            const bool isNodeValues,
             const bool verbose = false
         ) const;
 };