From 1528e8f20fad1bfca2f4a362ec8a8c01eb023058 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Thu, 3 Nov 2011 16:34:44 +0000
Subject: [PATCH] adjointShapeOptimizationFoam: Corrected global adjoint
 continuity error

---
 .../adjointContinuityErrs.H                   |  4 +-
 .../adjointShapeOptimizationFoam.C            |  1 +
 .../initAdjointContinuityErrs.H               | 43 +++++++++++++++++++
 3 files changed, 46 insertions(+), 2 deletions(-)
 create mode 100644 applications/solvers/incompressible/adjointShapeOptimizationFoam/initAdjointContinuityErrs.H

diff --git a/applications/solvers/incompressible/adjointShapeOptimizationFoam/adjointContinuityErrs.H b/applications/solvers/incompressible/adjointShapeOptimizationFoam/adjointContinuityErrs.H
index b8d6fd0687..8315aa3770 100644
--- a/applications/solvers/incompressible/adjointShapeOptimizationFoam/adjointContinuityErrs.H
+++ b/applications/solvers/incompressible/adjointShapeOptimizationFoam/adjointContinuityErrs.H
@@ -35,11 +35,11 @@ Description
 
     scalar globalContErr = runTime.deltaTValue()*
         fvc::div(phia)().weightedAverage(mesh.V()).value();
-    cumulativeContErr += globalContErr;
+    cumulativeAdjointContErr += globalContErr;
 
     Info<< "Adjoint continuity errors : sum local = " << sumLocalContErr
         << ", global = " << globalContErr
-        << ", cumulative = " << cumulativeContErr
+        << ", cumulative = " << cumulativeAdjointContErr
         << endl;
 }
 
diff --git a/applications/solvers/incompressible/adjointShapeOptimizationFoam/adjointShapeOptimizationFoam.C b/applications/solvers/incompressible/adjointShapeOptimizationFoam/adjointShapeOptimizationFoam.C
index bedc12b3ba..b5a699cc87 100644
--- a/applications/solvers/incompressible/adjointShapeOptimizationFoam/adjointShapeOptimizationFoam.C
+++ b/applications/solvers/incompressible/adjointShapeOptimizationFoam/adjointShapeOptimizationFoam.C
@@ -74,6 +74,7 @@ int main(int argc, char *argv[])
     #include "createMesh.H"
     #include "createFields.H"
     #include "initContinuityErrs.H"
+    #include "initAdjointContinuityErrs.H"
 
     simpleControl simple(mesh);
 
diff --git a/applications/solvers/incompressible/adjointShapeOptimizationFoam/initAdjointContinuityErrs.H b/applications/solvers/incompressible/adjointShapeOptimizationFoam/initAdjointContinuityErrs.H
new file mode 100644
index 0000000000..9d6c3b55be
--- /dev/null
+++ b/applications/solvers/incompressible/adjointShapeOptimizationFoam/initAdjointContinuityErrs.H
@@ -0,0 +1,43 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+     \\/     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/>.
+
+Global
+    cumulativeAdjointContErr
+
+Description
+    Declare and initialise the cumulative ddjoint continuity error.
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef initAdjointContinuityErrs_H
+#define initAdjointContinuityErrs_H
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+scalar cumulativeAdjointContErr = 0;
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //