From 5f21e8ca887d7f6b2bcba6edda35952ff0e37188 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Tue, 26 Nov 2013 12:31:58 +0000
Subject: [PATCH] ENH: force/Coeffs - added writing of bin co-ordinates to
 output

---
 .../forces/forceCoeffs/forceCoeffs.C          | 24 +++++++++++++++++++
 .../functionObjects/forces/forces/forces.C    | 24 +++++++++++++++++++
 2 files changed, 48 insertions(+)

diff --git a/src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C b/src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C
index 1daa0c279c..42ae77e264 100644
--- a/src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C
+++ b/src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C
@@ -72,6 +72,30 @@ void Foam::forceCoeffs::writeFileHeader(const label i)
         writeHeaderValue(file(i), "start", binMin_);
         writeHeaderValue(file(i), "delta", binDx_);
         writeHeaderValue(file(i), "direction", binDir_);
+
+        vectorField binPoints(nBin_);
+        writeCommented(file(i), "x co-ords  :");
+        forAll(binPoints, pointI)
+        {
+            binPoints[pointI] = (binMin_ + (pointI + 1)*binDx_)*binDir_;
+            file(i) << tab << binPoints[pointI].x();
+        }
+        file(i) << nl;
+
+        writeCommented(file(i), "y co-ords  :");
+        forAll(binPoints, pointI)
+        {
+            file(i) << tab << binPoints[pointI].y();
+        }
+        file(i) << nl;
+
+        writeCommented(file(i), "z co-ords  :");
+        forAll(binPoints, pointI)
+        {
+            file(i) << tab << binPoints[pointI].z();
+        }
+        file(i) << nl;
+
         writeCommented(file(i), "Time");
 
         for (label j = 0; j < nBin_; j++)
diff --git a/src/postProcessing/functionObjects/forces/forces/forces.C b/src/postProcessing/functionObjects/forces/forces/forces.C
index 70c1e00e78..c096ad5d4b 100644
--- a/src/postProcessing/functionObjects/forces/forces/forces.C
+++ b/src/postProcessing/functionObjects/forces/forces/forces.C
@@ -98,6 +98,30 @@ void Foam::forces::writeFileHeader(const label i)
         writeHeaderValue(file(i), "start", binMin_);
         writeHeaderValue(file(i), "delta", binDx_);
         writeHeaderValue(file(i), "direction", binDir_);
+
+        vectorField binPoints(nBin_);
+        writeCommented(file(i), "x co-ords  :");
+        forAll(binPoints, pointI)
+        {
+            binPoints[pointI] = (binMin_ + (pointI + 1)*binDx_)*binDir_;
+            file(i) << tab << binPoints[pointI].x();
+        }
+        file(i) << nl;
+
+        writeCommented(file(i), "y co-ords  :");
+        forAll(binPoints, pointI)
+        {
+            file(i) << tab << binPoints[pointI].y();
+        }
+        file(i) << nl;
+
+        writeCommented(file(i), "z co-ords  :");
+        forAll(binPoints, pointI)
+        {
+            file(i) << tab << binPoints[pointI].z();
+        }
+        file(i) << nl;
+
         writeCommented(file(i), "Time");
 
         for (label j = 0; j < nBin_; j++)