From 7c80ada340a393d183ecf4bdc20bba5a34553c62 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Tue, 12 Oct 2010 18:46:57 +0100
Subject: [PATCH 1/9] Fixed bug in usage of dynamic list

---
 .../directTopoChange/directTopoChange.C       | 36 ++++++++++---------
 .../repatchPolyTopoChanger.H                  |  4 +--
 2 files changed, 21 insertions(+), 19 deletions(-)

diff --git a/src/dynamicMesh/directTopoChange/directTopoChange/directTopoChange.C b/src/dynamicMesh/directTopoChange/directTopoChange/directTopoChange.C
index b0aa049a2..fe3e727d7 100644
--- a/src/dynamicMesh/directTopoChange/directTopoChange/directTopoChange.C
+++ b/src/dynamicMesh/directTopoChange/directTopoChange/directTopoChange.C
@@ -577,7 +577,7 @@ Foam::label Foam::directTopoChange::getCellOrder
     // Now we have new-to-old in newOrder.
     newOrder.setSize(cellInOrder);
 
-    // Invert to get old-to-new. Make sure removed (i.e. unmapped) cells are -1.
+    // Invert to get old-to-new. Make sure removed (i.e. unmapped) cells are -1
     oldToNew = invert(cellCellAddressing.size(), newOrder);
 
     return cellInOrder;
@@ -904,7 +904,7 @@ void Foam::directTopoChange::compact
                                     << " at position " << faceI << endl
                                     << "Probably face has not been adapted for"
                                     << " removed points." << abort(FatalError);
-                            } 
+                            }
                             localPointMap[pointI] = newPointI++;
                         }
                     }
@@ -1487,7 +1487,8 @@ void Foam::directTopoChange::resetZones
             {
                 if (newZoneAddr[i] < pointMap_.size())
                 {
-                    curPzRnb[i] = oldZone.whichPoint(pointMap_[newZoneAddr[i]]);
+                    curPzRnb[i] =
+                        oldZone.whichPoint(pointMap_[newZoneAddr[i]]);
                 }
                 else
                 {
@@ -1719,13 +1720,13 @@ void Foam::directTopoChange::calcFaceZonePointMap
 
         curFzPointRnb.setSize(newZoneMeshPoints.size());
 
-	if( zoneI < oldFaceZoneMeshPointMaps.size() )
+        if( zoneI < oldFaceZoneMeshPointMaps.size() )
         {
-     	    //HR 22.11.09: cannot take reference to zero element of 
+            //HR 22.11.09: cannot take reference to zero element of
             //  oldFaceZoneMeshPointMaps. It may be empty. Hence this
             //  if-statement needs to move out of the loop
-	    const Map<label>& oldZoneMeshPointMap =
-	        oldFaceZoneMeshPointMaps[zoneI];
+            const Map<label>& oldZoneMeshPointMap =
+                oldFaceZoneMeshPointMaps[zoneI];
 
             forAll (newZoneMeshPoints, pointI)
             {
@@ -2034,7 +2035,8 @@ void Foam::directTopoChange::compactAndReorder
     {
         //HR 22.11.09: Fixes error when faceZone is empty
         label maxFaceIndex = -1;
-	if ( mesh.faceZones()[zoneI].size() > 0 )
+
+        if ( mesh.faceZones()[zoneI].size() > 0 )
         {
             maxFaceIndex = max(mesh.faceZones()[zoneI]);
         }
@@ -2211,8 +2213,8 @@ void Foam::directTopoChange::addMesh
         const pointZoneMesh& pointZones = mesh.pointZones();
 
         // Resize
-        points_.setSize(points_.size() + points.size());
-        pointMap_.setSize(pointMap_.size() + points.size());
+        points_.setCapacity(points_.size() + points.size());
+        pointMap_.setCapacity(pointMap_.size() + points.size());
         pointZone_.resize(pointZone_.size() + points.size()/100);
 
         // Precalc offset zones
@@ -2251,11 +2253,11 @@ void Foam::directTopoChange::addMesh
         // always equals nCells
         label nAllCells = mesh.nCells();
 
-        cellMap_.setSize(cellMap_.size() + nAllCells);
+        cellMap_.setCapacity(cellMap_.size() + nAllCells);
         cellFromPoint_.resize(cellFromPoint_.size() + nAllCells/100);
         cellFromEdge_.resize(cellFromEdge_.size() + nAllCells/100);
         cellFromFace_.resize(cellFromFace_.size() + nAllCells/100);
-        cellZone_.setSize(cellZone_.size() + nAllCells);
+        cellZone_.setCapacity(cellZone_.size() + nAllCells);
 
 
         // Precalc offset zones
@@ -2311,11 +2313,11 @@ void Foam::directTopoChange::addMesh
         // Resize
         label nAllFaces = mesh.faces().size();
 
-        faces_.setSize(faces_.size() + nAllFaces);
-        region_.setSize(region_.size() + nAllFaces);
-        faceOwner_.setSize(faceOwner_.size() + nAllFaces);
-        faceNeighbour_.setSize(faceNeighbour_.size() + nAllFaces);
-        faceMap_.setSize(faceMap_.size() + nAllFaces);
+        faces_.setCapacity(faces_.size() + nAllFaces);
+        region_.setCapacity(region_.size() + nAllFaces);
+        faceOwner_.setCapacity(faceOwner_.size() + nAllFaces);
+        faceNeighbour_.setCapacity(faceNeighbour_.size() + nAllFaces);
+        faceMap_.setCapacity(faceMap_.size() + nAllFaces);
         faceFromPoint_.resize(faceFromPoint_.size() + nAllFaces/100);
         faceFromEdge_.resize(faceFromEdge_.size() + nAllFaces/100);
         flipFaceFlux_.resize(flipFaceFlux_.size() + nAllFaces/100);
diff --git a/src/dynamicMesh/directTopoChange/repatchPolyTopoChanger/repatchPolyTopoChanger.H b/src/dynamicMesh/directTopoChange/repatchPolyTopoChanger/repatchPolyTopoChanger.H
index 53e382b3d..778c0977a 100644
--- a/src/dynamicMesh/directTopoChange/repatchPolyTopoChanger/repatchPolyTopoChanger.H
+++ b/src/dynamicMesh/directTopoChange/repatchPolyTopoChanger/repatchPolyTopoChanger.H
@@ -28,7 +28,7 @@ Class
 Description
     A mesh which allows changes in the patch distribution of the
     boundary faces.  The change in patching is set using changePatchID. For a
-    boundary face, a new patch ID is given. 
+    boundary face, a new patch ID is given.
 
 SourceFiles
     repatchPolyTopoChanger.C
@@ -47,7 +47,7 @@ namespace Foam
 {
 
 /*---------------------------------------------------------------------------*\
-                       Class repatchPolyTopoChanger Declaration
+                  Class repatchPolyTopoChanger Declaration
 \*---------------------------------------------------------------------------*/
 
 class repatchPolyTopoChanger

From 801cfd90e9a24fc1e4986a25f2e26f4e6ca7a5ac Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Tue, 12 Oct 2010 23:05:08 +0100
Subject: [PATCH 2/9] Resolving Janaf lookup unboundedness problems

---
 etc/controlDict                               |  4 ++-
 .../specie/thermo/janaf/janafThermo.H         | 18 ++++++++----
 .../specie/thermo/janaf/janafThermoI.H        | 29 ++++++++++++-------
 .../specie/thermo/specieThermo/specieThermo.C |  7 +++++
 .../specie/thermo/specieThermo/specieThermo.H |  7 +++--
 .../thermo/specieThermo/specieThermoI.H       | 20 +++++++++----
 6 files changed, 61 insertions(+), 24 deletions(-)

diff --git a/etc/controlDict b/etc/controlDict
index 4812cbc73..a589f4c9b 100644
--- a/etc/controlDict
+++ b/etc/controlDict
@@ -911,7 +911,9 @@ Tolerances
     patchToPatchProjectionTol   0.05;
 
     // Thermophysical models
-    specieThermoTol     1e-4;
+    specieThermoTol      1e-4;
+    speciesThermoTJump     20;
+    speciesThermoMaxIter  100;
 
     // Intersection tolerance
     intersectionPlanarTol      0.2;
diff --git a/src/thermophysicalModels/specie/thermo/janaf/janafThermo.H b/src/thermophysicalModels/specie/thermo/janaf/janafThermo.H
index 10df8f5b6..03fd65867 100644
--- a/src/thermophysicalModels/specie/thermo/janaf/janafThermo.H
+++ b/src/thermophysicalModels/specie/thermo/janaf/janafThermo.H
@@ -114,10 +114,12 @@ private:
     // Private member functions
 
         //- Check given temperature is within the range of the fitted coeffs
-        inline void checkT(const scalar T) const;
+        //  Note: bounding T within range.  HJ, 12/Oct/2010
+        inline void checkT(scalar& T) const;
 
         //- Return the coefficients corresponding to the given temperature
-        inline const coeffArray& coeffs(const scalar T) const;
+        //  Note: bounding T within range.  HJ, 12/Oct/2010
+        inline const coeffArray& coeffs(scalar& T) const;
 
 
 public:
@@ -145,19 +147,23 @@ public:
     // Member Functions
 
         //- Heat capacity at constant pressure [J/(kmol K)]
-        inline scalar cp(const scalar T) const;
+        //  Note: bounding T within range.  HJ, 12/Oct/2010
+        inline scalar cp(scalar T) const;
 
         //- Enthalpy [J/kmol]
-        inline scalar h(const scalar T) const;
+        //  Note: bounding T within range.  HJ, 12/Oct/2010
+        inline scalar h(scalar T) const;
 
         //- Sensible enthalpy [J/kmol]
-        inline scalar hs(const scalar T) const;
+        //  Note: bounding T within range.  HJ, 12/Oct/2010
+        inline scalar hs(scalar T) const;
 
         //- Chemical enthalpy [J/kmol]
         inline scalar hc() const;
 
         //- Entropy [J/(kmol K)]
-        inline scalar s(const scalar T) const;
+        //  Note: bounding T within range.  HJ, 12/Oct/2010
+        inline scalar s(scalar T) const;
 
 
     // Member operators
diff --git a/src/thermophysicalModels/specie/thermo/janaf/janafThermoI.H b/src/thermophysicalModels/specie/thermo/janaf/janafThermoI.H
index 249960df3..853d90711 100644
--- a/src/thermophysicalModels/specie/thermo/janaf/janafThermoI.H
+++ b/src/thermophysicalModels/specie/thermo/janaf/janafThermoI.H
@@ -54,17 +54,21 @@ inline Foam::janafThermo<equationOfState>::janafThermo
 
 
 template<class equationOfState>
-inline void Foam::janafThermo<equationOfState>::checkT(const scalar T) const
+inline void Foam::janafThermo<equationOfState>::checkT(scalar& T) const
 {
     if (T <  Tlow_ || T > Thigh_)
     {
-        FatalErrorIn
+        // Improvements: graceful exit with recovery.  HJ, 11/Oct/2010
+        InfoIn
         (
-            "janafThermo<equationOfState>::checkT(const scalar T) const"
+            "janafThermo<equationOfState>::checkT(scalar& T) const"
         )   << "attempt to use janafThermo<equationOfState>"
                " out of temperature range "
             << Tlow_ << " -> " << Thigh_ << ";  T = " << T
-            << abort(FatalError);
+            << endl;
+
+        // Bracket T to avoid out-of-range error
+        T = Foam::min(Thigh_, Foam::max(T, Tlow_));
     }
 }
 
@@ -73,9 +77,11 @@ template<class equationOfState>
 inline const typename Foam::janafThermo<equationOfState>::coeffArray&
 Foam::janafThermo<equationOfState>::coeffs
 (
-    const scalar T
+    scalar& T
 ) const
 {
+    // Note: T will be bounded by checkT in coeffs(T).  No longer const
+    // HJ, 12/Oct/2010
     checkT(T);
 
     if (T < Tcommon_)
@@ -103,7 +109,7 @@ inline Foam::janafThermo<equationOfState>::janafThermo
     Thigh_(jt.Thigh_),
     Tcommon_(jt.Tcommon_)
 {
-    for (register label coefLabel=0; coefLabel<nCoeffs_; coefLabel++)
+    for (register label coefLabel = 0; coefLabel < nCoeffs_; coefLabel++)
     {
         highCpCoeffs_[coefLabel] = jt.highCpCoeffs_[coefLabel];
         lowCpCoeffs_[coefLabel] = jt.lowCpCoeffs_[coefLabel];
@@ -116,7 +122,7 @@ inline Foam::janafThermo<equationOfState>::janafThermo
 template<class equationOfState>
 inline Foam::scalar Foam::janafThermo<equationOfState>::cp
 (
-    const scalar T
+    scalar T
 ) const
 {
     const coeffArray& a = coeffs(T);
@@ -127,7 +133,7 @@ inline Foam::scalar Foam::janafThermo<equationOfState>::cp
 template<class equationOfState>
 inline Foam::scalar Foam::janafThermo<equationOfState>::h
 (
-    const scalar T
+    scalar T
 ) const
 {
     const coeffArray& a = coeffs(T);
@@ -142,7 +148,7 @@ inline Foam::scalar Foam::janafThermo<equationOfState>::h
 template<class equationOfState>
 inline Foam::scalar Foam::janafThermo<equationOfState>::hs
 (
-    const scalar T
+    scalar T
 ) const
 {
     return h(T) - hc();
@@ -167,9 +173,12 @@ inline Foam::scalar Foam::janafThermo<equationOfState>::hc() const
 template<class equationOfState>
 inline Foam::scalar Foam::janafThermo<equationOfState>::s
 (
-    const scalar T
+    scalar T
 ) const
 {
+    // Note: T will be bounded by checkT in coeffs(T).  No longer const
+    // HJ, 12/Oct/2010
+
     const coeffArray& a = coeffs(T);
     return
     this->RR*
diff --git a/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.C b/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.C
index 3336c5687..b7fd075d5 100644
--- a/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.C
+++ b/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.C
@@ -41,6 +41,13 @@ const Foam::scalar Foam::specieThermo<thermo>::tol_
 );
 
 
+template<class thermo>
+const Foam::scalar Foam::specieThermo<thermo>::TJump_
+(
+    debug::tolerances("speciesThermoTJump", 20)
+);
+
+
 template<class thermo>
 const int Foam::specieThermo<thermo>::maxIter_
 (
diff --git a/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.H b/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.H
index 074a13bc9..c5b66278a 100644
--- a/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.H
+++ b/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.H
@@ -98,13 +98,16 @@ class specieThermo
         //- Convergence tolerance of energy -> temperature inversion functions
         static const scalar tol_;
 
+        //- Max temperature jump of energy -> temperature inversion functions
+        static const scalar TJump_;
+
         //- Max number of iterations in energy->temperature inversion functions
         static const int maxIter_;
 
 
     // Private member functions
 
-        //- return the temperature corresponding to the value of the
+        //- Return the temperature corresponding to the value of the
         //  thermodynamic property f, given the function f = F(T) and dF(T)/dT
         inline scalar T
         (
@@ -131,7 +134,7 @@ public:
 
     // Member Functions
 
-        // Fundamaental properties
+        // Fundamental properties
         // (These functions must be provided in derived types)
 
             // Heat capacity at constant pressure [J/(kmol K)]
diff --git a/src/thermophysicalModels/specie/thermo/specieThermo/specieThermoI.H b/src/thermophysicalModels/specie/thermo/specieThermo/specieThermoI.H
index 50b9409e2..41dac7833 100644
--- a/src/thermophysicalModels/specie/thermo/specieThermo/specieThermoI.H
+++ b/src/thermophysicalModels/specie/thermo/specieThermo/specieThermoI.H
@@ -50,23 +50,33 @@ inline Foam::scalar Foam::specieThermo<thermo>::T
     scalar Test = T0;
     scalar Tnew = T0;
     scalar Ttol = T0*tol_;
-    int    iter = 0;
+    int iter = 0;
 
     do
     {
+        // Limit the temperature jump in a single corrector to TJump_
+        // HJ, 12/Oct/2010
+
         Test = Tnew;
-        Tnew = Test - ((this->*F)(Test) - f)/(this->*dFdT)(Test);
+        Tnew = Test
+            - Foam::min(((this->*F)(Test) - f)/(this->*dFdT)(Test), TJump_);
 
         if (iter++ > maxIter_)
         {
-            FatalErrorIn
+            // Improvements: graceful exit with recovery.  HJ, 11/Oct/2010
+            InfoIn
             (
                 "specieThermo<thermo>::T(scalar f, scalar T0, "
                 "scalar (specieThermo<thermo>::*F)(const scalar) const, "
                 "scalar (specieThermo<thermo>::*dFdT)(const scalar) const"
                 ") const"
-            )   << "Maximum number of iterations exceeded"
-                << abort(FatalError);
+            )   << "Maximum number of iterations exceeded.  Rescue by HJ"
+                << endl;
+
+            // Use value where dFdT is calculated using T0.  HJ, 11/Oct/2010
+            Tnew = f/(this->*dFdT)(T0);
+
+            return Tnew;
         }
 
     } while (mag(Tnew - Test) > Ttol);

From 69479b46f814bf819abff06d8c406871fec75386 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Wed, 13 Oct 2010 09:43:57 +0100
Subject: [PATCH 3/9] Fixed lduSolvers control format: dictionary

---
 src/lduSolvers/lduSolver/bicgSolver/bicgSolver.C           | 2 +-
 src/lduSolvers/lduSolver/bicgStabSolver/bicgStabSolver.C   | 2 +-
 src/lduSolvers/lduSolver/cgSolver/cgSolver.C               | 2 +-
 src/lduSolvers/lduSolver/deflationSolver/deflationSolver.C | 2 +-
 src/lduSolvers/lduSolver/gmresSolver/gmresSolver.C         | 2 +-
 5 files changed, 5 insertions(+), 5 deletions(-)

diff --git a/src/lduSolvers/lduSolver/bicgSolver/bicgSolver.C b/src/lduSolvers/lduSolver/bicgSolver/bicgSolver.C
index 58c4801ee..d3c618a0e 100644
--- a/src/lduSolvers/lduSolver/bicgSolver/bicgSolver.C
+++ b/src/lduSolvers/lduSolver/bicgSolver/bicgSolver.C
@@ -74,7 +74,7 @@ Foam::bicgSolver::bicgSolver
             coupleBouCoeffs,
             coupleIntCoeffs,
             interfaces,
-            dict.subDict("preconditioner")
+            dict
         )
     )
 {}
diff --git a/src/lduSolvers/lduSolver/bicgStabSolver/bicgStabSolver.C b/src/lduSolvers/lduSolver/bicgStabSolver/bicgStabSolver.C
index 7d8e1f155..eb8066528 100644
--- a/src/lduSolvers/lduSolver/bicgStabSolver/bicgStabSolver.C
+++ b/src/lduSolvers/lduSolver/bicgStabSolver/bicgStabSolver.C
@@ -78,7 +78,7 @@ Foam::bicgStabSolver::bicgStabSolver
             coupleBouCoeffs,
             coupleIntCoeffs,
             interfaces,
-            dict.subDict("preconditioner")
+            dict
         )
     )
 {}
diff --git a/src/lduSolvers/lduSolver/cgSolver/cgSolver.C b/src/lduSolvers/lduSolver/cgSolver/cgSolver.C
index 457dd8295..ce01e3d38 100644
--- a/src/lduSolvers/lduSolver/cgSolver/cgSolver.C
+++ b/src/lduSolvers/lduSolver/cgSolver/cgSolver.C
@@ -74,7 +74,7 @@ Foam::cgSolver::cgSolver
             coupleBouCoeffs,
             coupleIntCoeffs,
             interfaces,
-            dict.subDict("preconditioner")
+            dict
         )
     )
 {}
diff --git a/src/lduSolvers/lduSolver/deflationSolver/deflationSolver.C b/src/lduSolvers/lduSolver/deflationSolver/deflationSolver.C
index 1107c1ffa..207d00d4d 100644
--- a/src/lduSolvers/lduSolver/deflationSolver/deflationSolver.C
+++ b/src/lduSolvers/lduSolver/deflationSolver/deflationSolver.C
@@ -146,7 +146,7 @@ Foam::deflationSolver::deflationSolver
             coupleBouCoeffs,
             coupleIntCoeffs,
             interfaces,
-            dict.subDict("preconditioner")
+            dict
         )
     ),
     rpmOrder_(readLabel(dict.lookup("rpmOrder"))),
diff --git a/src/lduSolvers/lduSolver/gmresSolver/gmresSolver.C b/src/lduSolvers/lduSolver/gmresSolver/gmresSolver.C
index 479f6335e..e5aff27b6 100644
--- a/src/lduSolvers/lduSolver/gmresSolver/gmresSolver.C
+++ b/src/lduSolvers/lduSolver/gmresSolver/gmresSolver.C
@@ -109,7 +109,7 @@ Foam::gmresSolver::gmresSolver
             coupleBouCoeffs,
             coupleIntCoeffs,
             interfaces,
-            dict.subDict("preconditioner")
+            dict
         )
     ),
     nDirs_(readLabel(dict.lookup("nDirections")))

From 974746246f34787bcc1eb7e5c281a96c6ba2abb6 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Wed, 13 Oct 2010 09:45:45 +0100
Subject: [PATCH 4/9] Added demonstration code for dynamic mesh

---
 .../icoDyMFoam/.svn/all-wcprops               |  35 +++
 .../icoDyMFoam/.svn/dir-prop-base             |   7 +
 .../incompressible/icoDyMFoam/.svn/entries    | 201 ++++++++++++++++++
 .../incompressible/icoDyMFoam/.svn/format     |   1 +
 .../icoDyMFoam/.svn/prop-base/UEqn.H.svn-base |   5 +
 .../.svn/prop-base/correctPhi.H.svn-base      |   5 +
 .../.svn/prop-base/createFields.H.svn-base    |   5 +
 .../.svn/prop-base/icoDyMFoam.C.svn-base      |   5 +
 .../icoDyMFoam/.svn/text-base/UEqn.H.svn-base |  11 +
 .../.svn/text-base/correctPhi.H.svn-base      |  47 ++++
 .../.svn/text-base/createFields.H.svn-base    |  90 ++++++++
 .../.svn/text-base/icoDyMFoam.C.svn-base      | 144 +++++++++++++
 .../.svn/text-base/readControls.H.svn-base    |  14 ++
 .../icoDyMFoam/Make/.svn/all-wcprops          |  17 ++
 .../icoDyMFoam/Make/.svn/dir-prop-base        |   9 +
 .../icoDyMFoam/Make/.svn/entries              |  96 +++++++++
 .../icoDyMFoam/Make/.svn/format               |   1 +
 .../Make/.svn/text-base/files.svn-base        |   3 +
 .../Make/.svn/text-base/options.svn-base      |  13 ++
 .../incompressible/icoDyMFoam/Make/files      |   3 +
 .../incompressible/icoDyMFoam/Make/options    |  13 ++
 .../solvers/incompressible/icoDyMFoam/UEqn.H  |  11 +
 .../incompressible/icoDyMFoam/correctPhi.H    |  47 ++++
 .../incompressible/icoDyMFoam/createFields.H  |  90 ++++++++
 .../incompressible/icoDyMFoam/icoDyMFoam.C    | 149 +++++++++++++
 .../incompressible/icoDyMFoam/readControls.H  |  14 ++
 26 files changed, 1036 insertions(+)
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/all-wcprops
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/dir-prop-base
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/entries
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/format
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/UEqn.H.svn-base
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/correctPhi.H.svn-base
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/createFields.H.svn-base
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/icoDyMFoam.C.svn-base
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/text-base/UEqn.H.svn-base
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/text-base/correctPhi.H.svn-base
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/text-base/createFields.H.svn-base
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/text-base/icoDyMFoam.C.svn-base
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/text-base/readControls.H.svn-base
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/Make/.svn/all-wcprops
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/Make/.svn/dir-prop-base
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/Make/.svn/entries
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/Make/.svn/format
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/files.svn-base
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/options.svn-base
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/Make/files
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/Make/options
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/UEqn.H
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/correctPhi.H
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/createFields.H
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C
 create mode 100644 applications/solvers/incompressible/icoDyMFoam/readControls.H

diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/all-wcprops b/applications/solvers/incompressible/icoDyMFoam/.svn/all-wcprops
new file mode 100644
index 000000000..bceca7b2a
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/.svn/all-wcprops
@@ -0,0 +1,35 @@
+K 25
+svn:wc:ra_dav:version-url
+V 113
+/svnroot/openfoam-extend/!svn/ver/1591/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam
+END
+createFields.H
+K 25
+svn:wc:ra_dav:version-url
+V 127
+/svnroot/openfoam-extend/!svn/ver/915/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/createFields.H
+END
+readControls.H
+K 25
+svn:wc:ra_dav:version-url
+V 127
+/svnroot/openfoam-extend/!svn/ver/915/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/readControls.H
+END
+correctPhi.H
+K 25
+svn:wc:ra_dav:version-url
+V 125
+/svnroot/openfoam-extend/!svn/ver/981/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/correctPhi.H
+END
+icoDyMFoam.C
+K 25
+svn:wc:ra_dav:version-url
+V 126
+/svnroot/openfoam-extend/!svn/ver/1591/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C
+END
+UEqn.H
+K 25
+svn:wc:ra_dav:version-url
+V 119
+/svnroot/openfoam-extend/!svn/ver/912/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/UEqn.H
+END
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/dir-prop-base b/applications/solvers/incompressible/icoDyMFoam/.svn/dir-prop-base
new file mode 100644
index 000000000..3b7a2aa5a
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/.svn/dir-prop-base
@@ -0,0 +1,7 @@
+K 10
+svn:ignore
+V 16
+lnInclude
+*.dep
+
+END
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/entries b/applications/solvers/incompressible/icoDyMFoam/.svn/entries
new file mode 100644
index 000000000..a302fbbf2
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/.svn/entries
@@ -0,0 +1,201 @@
+9
+
+dir
+1730
+https://openfoam-extend.svn.sourceforge.net/svnroot/openfoam-extend/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam
+https://openfoam-extend.svn.sourceforge.net/svnroot/openfoam-extend
+
+
+
+2010-01-20T10:44:39.337304Z
+1591
+hjasak
+has-props
+
+svn:special svn:externals svn:needs-lock
+
+
+
+
+
+
+
+
+
+
+
+e4e07f05-0c2f-0410-a05a-b8ba57e0c909
+
+createFields.H
+file
+
+
+
+
+2009-08-11T21:05:38.000000Z
+00e915519b705cad0056a257a3fd2a92
+2009-01-20T00:30:44.740688Z
+915
+hjasak
+has-props
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+1874
+
+readControls.H
+file
+
+
+
+
+2009-08-11T21:05:38.000000Z
+eea361849937e2d7f8dcb00e72a34c00
+2009-01-20T00:30:44.740688Z
+915
+hjasak
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+362
+
+correctPhi.H
+file
+
+
+
+
+2009-08-11T21:05:38.000000Z
+074a2d4fbd974ca63d4207e2d816104f
+2009-02-06T16:17:03.879174Z
+981
+hjasak
+has-props
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+1048
+
+icoDyMFoam.C
+file
+
+
+
+
+2009-12-13T23:44:56.000000Z
+97fc42064b2254262c779f71eabc4991
+2010-01-20T10:44:39.337304Z
+1591
+hjasak
+has-props
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+4049
+
+Make
+dir
+
+UEqn.H
+file
+
+
+
+
+2009-08-11T21:05:38.000000Z
+84ccfd4a6cf96222b86f76eaacb9f17c
+2008-06-19T13:53:26.834514Z
+684
+hjasak
+has-props
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+190
+
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/format b/applications/solvers/incompressible/icoDyMFoam/.svn/format
new file mode 100644
index 000000000..ec635144f
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/.svn/format
@@ -0,0 +1 @@
+9
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/UEqn.H.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/UEqn.H.svn-base
new file mode 100644
index 000000000..92c8ad747
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/UEqn.H.svn-base
@@ -0,0 +1,5 @@
+K 12
+svn:keywords
+V 2
+Id
+END
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/correctPhi.H.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/correctPhi.H.svn-base
new file mode 100644
index 000000000..92c8ad747
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/correctPhi.H.svn-base
@@ -0,0 +1,5 @@
+K 12
+svn:keywords
+V 2
+Id
+END
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/createFields.H.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/createFields.H.svn-base
new file mode 100644
index 000000000..92c8ad747
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/createFields.H.svn-base
@@ -0,0 +1,5 @@
+K 12
+svn:keywords
+V 2
+Id
+END
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/icoDyMFoam.C.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/icoDyMFoam.C.svn-base
new file mode 100644
index 000000000..92c8ad747
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/icoDyMFoam.C.svn-base
@@ -0,0 +1,5 @@
+K 12
+svn:keywords
+V 2
+Id
+END
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/UEqn.H.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/UEqn.H.svn-base
new file mode 100644
index 000000000..faf508468
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/UEqn.H.svn-base
@@ -0,0 +1,11 @@
+    fvVectorMatrix UEqn
+    (
+        fvm::ddt(U)
+      + fvm::div(phi, U)
+      - fvm::laplacian(nu, U)
+    );
+
+    if (momentumPredictor)
+    {
+        solve(UEqn == -fvc::grad(p));
+    }
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/correctPhi.H.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/correctPhi.H.svn-base
new file mode 100644
index 000000000..69978f0e6
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/correctPhi.H.svn-base
@@ -0,0 +1,47 @@
+{
+    wordList pcorrTypes(p.boundaryField().types());
+
+    for (label i=0; i<p.boundaryField().size(); i++)
+    {
+        if (p.boundaryField()[i].fixesValue())
+        {
+            pcorrTypes[i] = fixedValueFvPatchScalarField::typeName;
+        }
+    }
+
+    volScalarField pcorr
+    (
+        IOobject
+        (
+            "pcorr",
+            runTime.timeName(),
+            mesh,
+            IOobject::NO_READ,
+            IOobject::NO_WRITE
+        ),
+        mesh,
+        dimensionedScalar("pcorr", p.dimensions(), 0.0),
+        pcorrTypes
+    );
+
+    for(int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
+    {
+        fvScalarMatrix pcorrEqn
+        (
+            fvm::laplacian(rAU, pcorr) == fvc::div(phi)
+        );
+
+        pcorrEqn.setReference(pRefCell, pRefValue);
+        pcorrEqn.solve();
+
+        if (nonOrth == nNonOrthCorr)
+        {
+            phi -= pcorrEqn.flux();
+        }
+
+        // Fluxes are corrected to absolute velocity and further corrected
+        // later.  HJ, 6/Feb/2009
+    }
+}
+
+#include "continuityErrs.H"
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/createFields.H.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/createFields.H.svn-base
new file mode 100644
index 000000000..561cd7bf4
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/createFields.H.svn-base
@@ -0,0 +1,90 @@
+    Info<< "Reading transportProperties\n" << endl;
+
+    IOdictionary transportProperties
+    (
+        IOobject
+        (
+            "transportProperties",
+            runTime.constant(),
+            mesh,
+            IOobject::MUST_READ,
+            IOobject::NO_WRITE
+        )
+    );
+
+    dimensionedScalar nu
+    (
+        transportProperties.lookup("nu")
+    );
+
+
+    Info<< "Reading field p\n" << endl;
+    volScalarField p
+    (
+        IOobject
+        (
+            "p",
+            runTime.timeName(),
+            mesh,
+            IOobject::MUST_READ,
+            IOobject::AUTO_WRITE
+        ),
+        mesh
+    );
+
+
+    Info<< "Reading field U\n" << endl;
+    volVectorField U
+    (
+        IOobject
+        (
+            "U",
+            runTime.timeName(),
+            mesh,
+            IOobject::MUST_READ,
+            IOobject::AUTO_WRITE
+        ),
+        mesh
+    );
+
+#   include "createPhi.H"
+
+
+    label pRefCell = 0;
+    scalar pRefValue = 0.0;
+    setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);
+
+
+    Info<< "Reading field rAU if present\n" << endl;
+    volScalarField rAU
+    (
+        IOobject
+        (
+            "rAU",
+            runTime.timeName(),
+            mesh,
+            IOobject::READ_IF_PRESENT,
+            IOobject::AUTO_WRITE
+        ),
+        mesh,
+        runTime.deltaT(),
+        zeroGradientFvPatchScalarField::typeName
+    );
+
+
+    {
+        const dictionary& piso = mesh.solutionDict().subDict("PISO");
+
+        bool momentumPredictor = true;
+        if (piso.found("momentumPredictor"))
+        {
+            momentumPredictor = Switch(piso.lookup("momentumPredictor"));
+        }
+
+        bool momentumPredictorSave = momentumPredictor;
+        momentumPredictor = false;
+#       include "UEqn.H"
+        momentumPredictor = momentumPredictorSave;
+        
+        rAU = 1.0/UEqn.A();
+    }
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/icoDyMFoam.C.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/icoDyMFoam.C.svn-base
new file mode 100644
index 000000000..1c85fdec5
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/icoDyMFoam.C.svn-base
@@ -0,0 +1,144 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2006-2007 Zeljko Tukovic
+     \\/     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 2 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, write to the Free Software Foundation,
+    Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
+
+Application
+    icoDyMFoam
+
+Description
+    Transient solver for incompressible, laminar flow of Newtonian fluids
+    with dynamic mesh.
+
+\*---------------------------------------------------------------------------*/
+
+#include "fvCFD.H"
+#include "dynamicFvMesh.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+int main(int argc, char *argv[])
+{
+#   include "setRootCase.H"
+#   include "createTime.H"
+#   include "createDynamicFvMesh.H"
+#   include "initContinuityErrs.H"
+#   include "initTotalVolume.H"
+#   include "createFields.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+    Info<< "\nStarting time loop\n" << endl;
+
+    while (runTime.run())
+    {
+#       include "readControls.H"
+#       include "checkTotalVolume.H"
+#       include "CourantNo.H"
+
+#       include "setDeltaT.H"
+
+        // Make the fluxes absolute
+        fvc::makeAbsolute(phi, U);
+
+        runTime++;
+
+        Info<< "Time = " << runTime.timeName() << nl << endl;
+
+        bool meshChanged = mesh.update();
+
+#       include "volContinuity.H"
+
+        if (correctPhi && meshChanged)
+        {
+#           include "correctPhi.H"
+        }
+
+        // Make the fluxes relative to the mesh motion
+        fvc::makeRelative(phi, U);
+
+        if (meshChanged)
+        {
+#           include "CourantNo.H"
+        }
+
+#       include "UEqn.H"
+
+        // --- PISO loop
+
+        for (int corr=0; corr<nCorr; corr++)
+        {
+            rAU = 1.0/UEqn.A();
+
+            U = rAU*UEqn.H();
+            phi = (fvc::interpolate(U) & mesh.Sf());
+              //+ fvc::ddtPhiCorr(rAU, U, phi);
+
+            adjustPhi(phi, U, p);
+
+            for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
+            {
+                fvScalarMatrix pEqn
+                (
+                    fvm::laplacian(rAU, p) == fvc::div(phi)
+                );
+
+                pEqn.setReference(pRefCell, pRefValue);
+
+                if (corr == nCorr - 1 && nonOrth == nNonOrthCorr)
+                {
+                    pEqn.solve(mesh.solver(p.name() + "Final"));
+                }
+                else
+                {
+                    pEqn.solve(mesh.solver(p.name()));
+                }
+
+                if (nonOrth == nNonOrthCorr)
+                {
+                    phi -= pEqn.flux();
+                }
+            }
+
+#           include "continuityErrs.H"
+
+            // Make the fluxes relative to the mesh motion
+            fvc::makeRelative(phi, U);
+
+            U -= rAU*fvc::grad(p);
+            U.correctBoundaryConditions();
+        }
+
+        runTime.write();
+
+        Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
+            << "  ClockTime = " << runTime.elapsedClockTime() << " s"
+            << nl << endl;
+    }
+
+    Info<< "End\n" << endl;
+
+    return(0);
+}
+
+
+// ************************************************************************* //
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/readControls.H.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/readControls.H.svn-base
new file mode 100644
index 000000000..9336616a1
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/readControls.H.svn-base
@@ -0,0 +1,14 @@
+#   include "readTimeControls.H"
+#   include "readPISOControls.H"
+
+    bool correctPhi = false;
+    if (piso.found("correctPhi"))
+    {
+        correctPhi = Switch(piso.lookup("correctPhi"));
+    }
+
+    bool checkMeshCourantNo = false;
+    if (piso.found("checkMeshCourantNo"))
+    {
+        checkMeshCourantNo = Switch(piso.lookup("checkMeshCourantNo"));
+    }
diff --git a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/all-wcprops b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/all-wcprops
new file mode 100644
index 000000000..d156d8090
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/all-wcprops
@@ -0,0 +1,17 @@
+K 25
+svn:wc:ra_dav:version-url
+V 117
+/svnroot/openfoam-extend/!svn/ver/912/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/Make
+END
+files
+K 25
+svn:wc:ra_dav:version-url
+V 123
+/svnroot/openfoam-extend/!svn/ver/912/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/Make/files
+END
+options
+K 25
+svn:wc:ra_dav:version-url
+V 125
+/svnroot/openfoam-extend/!svn/ver/912/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/Make/options
+END
diff --git a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/dir-prop-base b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/dir-prop-base
new file mode 100644
index 000000000..91bc687f2
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/dir-prop-base
@@ -0,0 +1,9 @@
+K 10
+svn:ignore
+V 28
+*Opt
+*Debug
+lnInclude
+*.dep
+
+END
diff --git a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/entries b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/entries
new file mode 100644
index 000000000..0c0be6d7a
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/entries
@@ -0,0 +1,96 @@
+9
+
+dir
+1730
+https://openfoam-extend.svn.sourceforge.net/svnroot/openfoam-extend/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/Make
+https://openfoam-extend.svn.sourceforge.net/svnroot/openfoam-extend
+
+
+
+2007-11-08T17:59:31.038014Z
+280
+hjasak
+has-props
+
+svn:special svn:externals svn:needs-lock
+
+
+
+
+
+
+
+
+
+
+
+e4e07f05-0c2f-0410-a05a-b8ba57e0c909
+
+files
+file
+
+
+
+
+2009-08-11T21:05:38.000000Z
+b327ba4ea7390a276aef34f189cc6ba6
+2007-06-05T22:18:47.316404Z
+3
+oforiginal
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+46
+
+options
+file
+
+
+
+
+2009-08-11T21:05:38.000000Z
+290d08439d3bdbe6dae9992a3d3d9f14
+2007-11-08T17:59:31.038014Z
+280
+hjasak
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+302
+
diff --git a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/format b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/format
new file mode 100644
index 000000000..ec635144f
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/format
@@ -0,0 +1 @@
+9
diff --git a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/files.svn-base b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/files.svn-base
new file mode 100644
index 000000000..22873bb94
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/files.svn-base
@@ -0,0 +1,3 @@
+icoDyMFoam.C
+
+EXE = $(FOAM_APPBIN)/icoDyMFoam
diff --git a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/options.svn-base b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/options.svn-base
new file mode 100644
index 000000000..460426a5b
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/options.svn-base
@@ -0,0 +1,13 @@
+EXE_INC = \
+    -I$(LIB_SRC)/dynamicFvMesh/lnInclude \
+    -I$(LIB_SRC)/dynamicMesh/lnInclude \
+    -I$(LIB_SRC)/meshTools/lnInclude \
+    -I$(LIB_SRC)/finiteVolume/lnInclude
+
+EXE_LIBS = \
+    -ldynamicFvMesh \
+    -ldynamicMesh \
+    -lengine \
+    -lmeshTools \
+    -lfiniteVolume \
+    -llduSolvers
diff --git a/applications/solvers/incompressible/icoDyMFoam/Make/files b/applications/solvers/incompressible/icoDyMFoam/Make/files
new file mode 100644
index 000000000..22873bb94
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/Make/files
@@ -0,0 +1,3 @@
+icoDyMFoam.C
+
+EXE = $(FOAM_APPBIN)/icoDyMFoam
diff --git a/applications/solvers/incompressible/icoDyMFoam/Make/options b/applications/solvers/incompressible/icoDyMFoam/Make/options
new file mode 100644
index 000000000..460426a5b
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/Make/options
@@ -0,0 +1,13 @@
+EXE_INC = \
+    -I$(LIB_SRC)/dynamicFvMesh/lnInclude \
+    -I$(LIB_SRC)/dynamicMesh/lnInclude \
+    -I$(LIB_SRC)/meshTools/lnInclude \
+    -I$(LIB_SRC)/finiteVolume/lnInclude
+
+EXE_LIBS = \
+    -ldynamicFvMesh \
+    -ldynamicMesh \
+    -lengine \
+    -lmeshTools \
+    -lfiniteVolume \
+    -llduSolvers
diff --git a/applications/solvers/incompressible/icoDyMFoam/UEqn.H b/applications/solvers/incompressible/icoDyMFoam/UEqn.H
new file mode 100644
index 000000000..faf508468
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/UEqn.H
@@ -0,0 +1,11 @@
+    fvVectorMatrix UEqn
+    (
+        fvm::ddt(U)
+      + fvm::div(phi, U)
+      - fvm::laplacian(nu, U)
+    );
+
+    if (momentumPredictor)
+    {
+        solve(UEqn == -fvc::grad(p));
+    }
diff --git a/applications/solvers/incompressible/icoDyMFoam/correctPhi.H b/applications/solvers/incompressible/icoDyMFoam/correctPhi.H
new file mode 100644
index 000000000..69978f0e6
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/correctPhi.H
@@ -0,0 +1,47 @@
+{
+    wordList pcorrTypes(p.boundaryField().types());
+
+    for (label i=0; i<p.boundaryField().size(); i++)
+    {
+        if (p.boundaryField()[i].fixesValue())
+        {
+            pcorrTypes[i] = fixedValueFvPatchScalarField::typeName;
+        }
+    }
+
+    volScalarField pcorr
+    (
+        IOobject
+        (
+            "pcorr",
+            runTime.timeName(),
+            mesh,
+            IOobject::NO_READ,
+            IOobject::NO_WRITE
+        ),
+        mesh,
+        dimensionedScalar("pcorr", p.dimensions(), 0.0),
+        pcorrTypes
+    );
+
+    for(int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
+    {
+        fvScalarMatrix pcorrEqn
+        (
+            fvm::laplacian(rAU, pcorr) == fvc::div(phi)
+        );
+
+        pcorrEqn.setReference(pRefCell, pRefValue);
+        pcorrEqn.solve();
+
+        if (nonOrth == nNonOrthCorr)
+        {
+            phi -= pcorrEqn.flux();
+        }
+
+        // Fluxes are corrected to absolute velocity and further corrected
+        // later.  HJ, 6/Feb/2009
+    }
+}
+
+#include "continuityErrs.H"
diff --git a/applications/solvers/incompressible/icoDyMFoam/createFields.H b/applications/solvers/incompressible/icoDyMFoam/createFields.H
new file mode 100644
index 000000000..561cd7bf4
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/createFields.H
@@ -0,0 +1,90 @@
+    Info<< "Reading transportProperties\n" << endl;
+
+    IOdictionary transportProperties
+    (
+        IOobject
+        (
+            "transportProperties",
+            runTime.constant(),
+            mesh,
+            IOobject::MUST_READ,
+            IOobject::NO_WRITE
+        )
+    );
+
+    dimensionedScalar nu
+    (
+        transportProperties.lookup("nu")
+    );
+
+
+    Info<< "Reading field p\n" << endl;
+    volScalarField p
+    (
+        IOobject
+        (
+            "p",
+            runTime.timeName(),
+            mesh,
+            IOobject::MUST_READ,
+            IOobject::AUTO_WRITE
+        ),
+        mesh
+    );
+
+
+    Info<< "Reading field U\n" << endl;
+    volVectorField U
+    (
+        IOobject
+        (
+            "U",
+            runTime.timeName(),
+            mesh,
+            IOobject::MUST_READ,
+            IOobject::AUTO_WRITE
+        ),
+        mesh
+    );
+
+#   include "createPhi.H"
+
+
+    label pRefCell = 0;
+    scalar pRefValue = 0.0;
+    setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);
+
+
+    Info<< "Reading field rAU if present\n" << endl;
+    volScalarField rAU
+    (
+        IOobject
+        (
+            "rAU",
+            runTime.timeName(),
+            mesh,
+            IOobject::READ_IF_PRESENT,
+            IOobject::AUTO_WRITE
+        ),
+        mesh,
+        runTime.deltaT(),
+        zeroGradientFvPatchScalarField::typeName
+    );
+
+
+    {
+        const dictionary& piso = mesh.solutionDict().subDict("PISO");
+
+        bool momentumPredictor = true;
+        if (piso.found("momentumPredictor"))
+        {
+            momentumPredictor = Switch(piso.lookup("momentumPredictor"));
+        }
+
+        bool momentumPredictorSave = momentumPredictor;
+        momentumPredictor = false;
+#       include "UEqn.H"
+        momentumPredictor = momentumPredictorSave;
+        
+        rAU = 1.0/UEqn.A();
+    }
diff --git a/applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C b/applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C
new file mode 100644
index 000000000..4a0d628fb
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C
@@ -0,0 +1,149 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright held by original author
+     \\/     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 2 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, write to the Free Software Foundation,
+    Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
+
+Application
+    icoDyMFoam
+
+Description
+    Transient solver for incompressible, laminar flow of Newtonian fluids
+    with dynamic mesh.
+
+Author
+    Hrvoje Jasak, Wikki Ltd.  All rights reserved.
+
+\*---------------------------------------------------------------------------*/
+
+#include "fvCFD.H"
+#include "dynamicFvMesh.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+int main(int argc, char *argv[])
+{
+#   include "setRootCase.H"
+#   include "createTime.H"
+#   include "createDynamicFvMesh.H"
+#   include "initContinuityErrs.H"
+#   include "initTotalVolume.H"
+#   include "createFields.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+    Info<< "\nStarting time loop\n" << endl;
+
+    while (runTime.run())
+    {
+#       include "readControls.H"
+#       include "checkTotalVolume.H"
+#       include "CourantNo.H"
+
+#       include "setDeltaT.H"
+
+        // Make the fluxes absolute
+        fvc::makeAbsolute(phi, U);
+
+        runTime++;
+
+        Info<< "Time = " << runTime.timeName() << nl << endl;
+
+        bool meshChanged = mesh.update();
+
+#       include "volContinuity.H"
+
+        if (correctPhi && meshChanged)
+        {
+            // Fluxes will be corrected to absolute velocity
+            // HJ, 6/Feb/2009
+#           include "correctPhi.H"
+        }
+
+        // Make the fluxes relative to the mesh motion
+        fvc::makeRelative(phi, U);
+
+        if (meshChanged)
+        {
+#           include "CourantNo.H"
+        }
+
+#       include "UEqn.H"
+
+        // --- PISO loop
+
+        for (int corr=0; corr<nCorr; corr++)
+        {
+            rAU = 1.0/UEqn.A();
+
+            U = rAU*UEqn.H();
+            phi = (fvc::interpolate(U) & mesh.Sf());
+              //+ fvc::ddtPhiCorr(rAU, U, phi);
+
+            adjustPhi(phi, U, p);
+
+            for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
+            {
+                fvScalarMatrix pEqn
+                (
+                    fvm::laplacian(rAU, p) == fvc::div(phi)
+                );
+
+                pEqn.setReference(pRefCell, pRefValue);
+
+                if (corr == nCorr - 1 && nonOrth == nNonOrthCorr)
+                {
+                    pEqn.solve(mesh.solver(p.name() + "Final"));
+                }
+                else
+                {
+                    pEqn.solve(mesh.solver(p.name()));
+                }
+
+                if (nonOrth == nNonOrthCorr)
+                {
+                    phi -= pEqn.flux();
+                }
+            }
+
+#           include "continuityErrs.H"
+
+            // Make the fluxes relative to the mesh motion
+            fvc::makeRelative(phi, U);
+
+            U -= rAU*fvc::grad(p);
+            U.correctBoundaryConditions();
+        }
+
+        runTime.write();
+
+        Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
+            << "  ClockTime = " << runTime.elapsedClockTime() << " s"
+            << nl << endl;
+    }
+
+    Info<< "End\n" << endl;
+
+    return(0);
+}
+
+
+// ************************************************************************* //
diff --git a/applications/solvers/incompressible/icoDyMFoam/readControls.H b/applications/solvers/incompressible/icoDyMFoam/readControls.H
new file mode 100644
index 000000000..9336616a1
--- /dev/null
+++ b/applications/solvers/incompressible/icoDyMFoam/readControls.H
@@ -0,0 +1,14 @@
+#   include "readTimeControls.H"
+#   include "readPISOControls.H"
+
+    bool correctPhi = false;
+    if (piso.found("correctPhi"))
+    {
+        correctPhi = Switch(piso.lookup("correctPhi"));
+    }
+
+    bool checkMeshCourantNo = false;
+    if (piso.found("checkMeshCourantNo"))
+    {
+        checkMeshCourantNo = Switch(piso.lookup("checkMeshCourantNo"));
+    }

From 93effde8e98c041f66352502a976c2fa89e08481 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Wed, 13 Oct 2010 12:38:34 +0100
Subject: [PATCH 5/9] Tutorial update

---
 tutorials/basic/sixDOFSolver/Allclean         |  14 --
 tutorials/basic/sixDOFSolver/Allrun           |  13 --
 .../basic/sixDOFSolver/sixDOFmotion/Allclean  |   7 +
 .../basic/sixDOFSolver/sixDOFmotion/Allrun    |   7 +
 .../basic/sixDOFSolver/springDamper/Allclean  |   7 +
 .../basic/sixDOFSolver/springDamper/Allrun    |   7 +
 .../mixer2D/constant/polyMesh/boundary        |  52 -------
 .../icoDyMFoam/mixer2D/system/fvSolution      |  60 +++-----
 .../icoDyMFoam/mixerGgi/Allclean              |   2 +-
 .../icoDyMFoam/mixerGgi/system/controlDict    |   2 +-
 .../icoDyMFoam/mixerGgi/system/fvSolution     |  33 ++---
 .../RBFMotionFunctionObject.C                 |  17 ++-
 .../RBFMotionFunctionObject.H                 |  17 +--
 .../movingBlockRBF/constant/polyMesh/boundary |   8 +-
 .../movingBlockRBF/system/controlDict         |   8 +-
 .../movingBlockRBF/system/fvSolution          |  27 ++--
 .../movingCone/constant/polyMesh/boundary     |   7 +-
 .../icoDyMFoam/movingCone/system/fvSolution   |  12 +-
 .../movingCone/system/tetFemSolution          |  22 ---
 .../icoDyMFoam/movingConeMotion/Allclean      |   8 ++
 .../icoDyMFoam/movingConeMotion/Allrun        |  11 ++
 .../constant/polyMesh/boundary                | 128 ++++++++----------
 .../movingConeMotion/system/fvSolution        |  54 ++------
 .../movingConeMotion/system/tetFemSolution    |   3 +-
 .../movingConeTopo/constant/polyMesh/boundary |   6 +-
 .../movingConeTopo/system/fvSolution          |  32 ++---
 .../icoFoam/cavity/constant/polyMesh/boundary |   7 +-
 .../cavityClipped/constant/polyMesh/boundary  |   7 +-
 .../cavityGrade/constant/polyMesh/boundary    |   7 +-
 .../icoFoam/elbow/constant/polyMesh/boundary  |   7 +-
 .../viscoelasticFluidFoam/Allclean            |  56 --------
 .../viscoelastic/viscoelasticFluidFoam/Allrun |  55 --------
 .../viscoelasticFluidFoam/DCPP/Allclean       |   6 +-
 .../viscoelasticFluidFoam/FENE-CR/Allclean    |   3 +-
 .../viscoelasticFluidFoam/FENE-P/Allclean     |   5 +-
 .../viscoelasticFluidFoam/Feta-PTT/Allclean   |   3 +-
 .../viscoelasticFluidFoam/Giesekus/Allclean   |   5 +-
 .../viscoelasticFluidFoam/Leonov/Allclean     |   4 +-
 .../viscoelasticFluidFoam/Oldroyd-B/Allclean  |   5 +-
 .../PTT-Exponential/Allclean                  |   3 +-
 .../viscoelasticFluidFoam/PTT-Linear/Allclean |   5 +-
 .../viscoelasticFluidFoam/S-MDCPP/Allclean    |   3 +-
 .../viscoelasticFluidFoam/XPP_DE/Allclean     |   3 +-
 .../viscoelasticFluidFoam/XPP_SE/Allclean     |   1 -
 44 files changed, 247 insertions(+), 502 deletions(-)
 delete mode 100755 tutorials/basic/sixDOFSolver/Allclean
 delete mode 100755 tutorials/basic/sixDOFSolver/Allrun
 create mode 100755 tutorials/basic/sixDOFSolver/sixDOFmotion/Allclean
 create mode 100755 tutorials/basic/sixDOFSolver/sixDOFmotion/Allrun
 create mode 100755 tutorials/basic/sixDOFSolver/springDamper/Allclean
 create mode 100755 tutorials/basic/sixDOFSolver/springDamper/Allrun
 delete mode 100644 tutorials/incompressible/icoDyMFoam/mixer2D/constant/polyMesh/boundary
 delete mode 100644 tutorials/incompressible/icoDyMFoam/movingCone/system/tetFemSolution
 create mode 100755 tutorials/incompressible/icoDyMFoam/movingConeMotion/Allclean
 create mode 100755 tutorials/incompressible/icoDyMFoam/movingConeMotion/Allrun
 delete mode 100755 tutorials/viscoelastic/viscoelasticFluidFoam/Allclean
 delete mode 100755 tutorials/viscoelastic/viscoelasticFluidFoam/Allrun

diff --git a/tutorials/basic/sixDOFSolver/Allclean b/tutorials/basic/sixDOFSolver/Allclean
deleted file mode 100755
index a4dd41e90..000000000
--- a/tutorials/basic/sixDOFSolver/Allclean
+++ /dev/null
@@ -1,14 +0,0 @@
-#!/bin/sh
-
-currDir=`pwd`
-application=`basename $currDir`
-cases="springDamper sixDOFmotion"
-
-tutorialPath=`dirname $0`/..
-. $WM_PROJECT_DIR/bin/tools/CleanFunctions
-
-for case in $cases
-do
-    cleanCase $case
-    rm $case/*.dat
-done 
diff --git a/tutorials/basic/sixDOFSolver/Allrun b/tutorials/basic/sixDOFSolver/Allrun
deleted file mode 100755
index d028eaf4a..000000000
--- a/tutorials/basic/sixDOFSolver/Allrun
+++ /dev/null
@@ -1,13 +0,0 @@
-#!/bin/sh
-
-currDir=`pwd`
-application=`basename $currDir`
-cases="springDamper sixDOFmotion"
-
-tutorialPath=`dirname $0`/..
-. $WM_PROJECT_DIR/bin/tools/RunFunctions
-
-for case in $cases
-do
-    (cd $case; runApplication $application)
-done
diff --git a/tutorials/basic/sixDOFSolver/sixDOFmotion/Allclean b/tutorials/basic/sixDOFSolver/sixDOFmotion/Allclean
new file mode 100755
index 000000000..baa1c8da8
--- /dev/null
+++ b/tutorials/basic/sixDOFSolver/sixDOFmotion/Allclean
@@ -0,0 +1,7 @@
+#!/bin/sh
+
+# Source tutorial clean functions
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions
+
+cleanCase
+rm -f *.dat
diff --git a/tutorials/basic/sixDOFSolver/sixDOFmotion/Allrun b/tutorials/basic/sixDOFSolver/sixDOFmotion/Allrun
new file mode 100755
index 000000000..4332c2541
--- /dev/null
+++ b/tutorials/basic/sixDOFSolver/sixDOFmotion/Allrun
@@ -0,0 +1,7 @@
+#!/bin/sh
+# Source tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+application="sixDOFSolver"
+
+runApplication $application
diff --git a/tutorials/basic/sixDOFSolver/springDamper/Allclean b/tutorials/basic/sixDOFSolver/springDamper/Allclean
new file mode 100755
index 000000000..baa1c8da8
--- /dev/null
+++ b/tutorials/basic/sixDOFSolver/springDamper/Allclean
@@ -0,0 +1,7 @@
+#!/bin/sh
+
+# Source tutorial clean functions
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions
+
+cleanCase
+rm -f *.dat
diff --git a/tutorials/basic/sixDOFSolver/springDamper/Allrun b/tutorials/basic/sixDOFSolver/springDamper/Allrun
new file mode 100755
index 000000000..4332c2541
--- /dev/null
+++ b/tutorials/basic/sixDOFSolver/springDamper/Allrun
@@ -0,0 +1,7 @@
+#!/bin/sh
+# Source tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+application="sixDOFSolver"
+
+runApplication $application
diff --git a/tutorials/incompressible/icoDyMFoam/mixer2D/constant/polyMesh/boundary b/tutorials/incompressible/icoDyMFoam/mixer2D/constant/polyMesh/boundary
deleted file mode 100644
index d229dca29..000000000
--- a/tutorials/incompressible/icoDyMFoam/mixer2D/constant/polyMesh/boundary
+++ /dev/null
@@ -1,52 +0,0 @@
-/*--------------------------------*- C++ -*----------------------------------*\
-| =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.5-dev                               |
-|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
-|    \\/     M anipulation  |                                                 |
-\*---------------------------------------------------------------------------*/
-FoamFile
-{
-    version     2.0;
-    format      ascii;
-    class       polyBoundaryMesh;
-    location    "constant/polyMesh";
-    object      boundary;
-}
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-5
-(
-    impellerWall
-    {
-        type            wall;
-        nFaces          68;
-        startFace       1040;
-    }
-    baffleWall
-    {
-        type            wall;
-        nFaces          84;
-        startFace       1108;
-    }
-    insideSlider
-    {
-        type            patch;
-        nFaces          36;
-        startFace       1192;
-    }
-    outsideSlider
-    {
-        type            patch;
-        nFaces          36;
-        startFace       1228;
-    }
-    defaultFaces
-    {
-        type            empty;
-        nFaces          1152;
-        startFace       1264;
-    }
-)
-
-// ************************************************************************* //
diff --git a/tutorials/incompressible/icoDyMFoam/mixer2D/system/fvSolution b/tutorials/incompressible/icoDyMFoam/mixer2D/system/fvSolution
index 6e99f2828..bba0e468b 100644
--- a/tutorials/incompressible/icoDyMFoam/mixer2D/system/fvSolution
+++ b/tutorials/incompressible/icoDyMFoam/mixer2D/system/fvSolution
@@ -10,12 +10,6 @@ FoamFile
 {
     version         2.0;
     format          ascii;
-
-    root            "";
-    case            "";
-    instance        "";
-    local           "";
-
     class           dictionary;
     object          fvSolution;
 }
@@ -24,12 +18,10 @@ FoamFile
 
 solvers
 {
-    pcorr PCG
+    pcorr
     {
-        preconditioner
-        {
-            type         DIC;
-        }
+        solver           BiCGStab;
+        preconditioner   DILU;
 
         minIter          0;
         maxIter          1000;
@@ -37,59 +29,41 @@ solvers
         relTol           0;
     };
 
-    p PCG
+    p
     {
-        preconditioner
-        {
-            type         DIC;
-        }
+        solver           CG;
+        preconditioner   Cholesky;
 
         minIter          0;
         maxIter          1000;
-        tolerance        1e-06;
-        relTol           0.05;
+        tolerance        1e-07;
+        relTol           0.0;
     };
 
-    pFinal PCG
+    pFinal
     {
-        preconditioner
-        {
-            type         DIC;
-        }
+        solver           CG;
+        preconditioner   Cholesky;
 
         minIter          0;
         maxIter          1000;
-        tolerance        1e-06;
+        tolerance        1e-07;
         relTol           0;
     };
 
-    U PBiCG
+    U
     {
-        preconditioner   
-        {
-            type         DILU;
-        }
+        solver            BiCGStab;
+        preconditioner    DILU;
 
         minIter          0;
         maxIter          1000;
-        tolerance        1e-05;
-        relTol           0;
-    };
-
-    cellMotionUx PCG
-    {
-        preconditioner      
-        {
-            type         DIC;
-        }
-
-        minIter          0;
-        maxIter          1000;
-        tolerance        1e-08;
+        tolerance        1e-07;
         relTol           0;
     };
 }
 
+
 PISO
 {
     nCorrectors     4;
diff --git a/tutorials/incompressible/icoDyMFoam/mixerGgi/Allclean b/tutorials/incompressible/icoDyMFoam/mixerGgi/Allclean
index b8e76f777..1892bfc36 100755
--- a/tutorials/incompressible/icoDyMFoam/mixerGgi/Allclean
+++ b/tutorials/incompressible/icoDyMFoam/mixerGgi/Allclean
@@ -5,5 +5,5 @@
 
 cleanCase
 rm -rf VTK
-rm constant/polyMesh/boundary
+rm -f constant/polyMesh/boundary
 rm -rf constant/polyMesh/sets
\ No newline at end of file
diff --git a/tutorials/incompressible/icoDyMFoam/mixerGgi/system/controlDict b/tutorials/incompressible/icoDyMFoam/mixerGgi/system/controlDict
index eeeea40ff..e4bd520e2 100644
--- a/tutorials/incompressible/icoDyMFoam/mixerGgi/system/controlDict
+++ b/tutorials/incompressible/icoDyMFoam/mixerGgi/system/controlDict
@@ -62,7 +62,7 @@ functions
         phi phi;
 
         // Where to load it from (if not already in solver)
-        functionObjectLibs ("libsampling.so");
+        functionObjectLibs ("libcheckFunctionObjects.so");
     }
 );
 
diff --git a/tutorials/incompressible/icoDyMFoam/mixerGgi/system/fvSolution b/tutorials/incompressible/icoDyMFoam/mixerGgi/system/fvSolution
index afa0988fb..8ec6bae19 100644
--- a/tutorials/incompressible/icoDyMFoam/mixerGgi/system/fvSolution
+++ b/tutorials/incompressible/icoDyMFoam/mixerGgi/system/fvSolution
@@ -18,12 +18,10 @@ FoamFile
 
 solvers
 {
-    pcorr BiCGStab
+    pcorr
     {
-        preconditioner
-        {
-            type         DILU;
-        }
+        solver           BiCGStab;
+        preconditioner   DILU;
 
         minIter          0;
         maxIter          1000;
@@ -31,12 +29,10 @@ solvers
         relTol           0;
     };
 
-    p BiCGStab
+    p
     {
-        preconditioner
-        {
-            type         DILU;
-        }
+        solver           CG;
+        preconditioner   Cholesky;
 
         minIter          0;
         maxIter          1000;
@@ -44,12 +40,10 @@ solvers
         relTol           0.0;
     };
 
-    pFinal BiCGStab
+    pFinal
     {
-        preconditioner
-        {
-            type         DILU;
-        }
+        solver           CG;
+        preconditioner   Cholesky;
 
         minIter          0;
         maxIter          1000;
@@ -57,12 +51,10 @@ solvers
         relTol           0;
     };
 
-    U BiCGStab
+    U
     {
-        preconditioner   
-        {
-            type         DILU;
-        }
+        solver            BiCGStab;
+        preconditioner    DILU;
 
         minIter          0;
         maxIter          1000;
@@ -71,6 +63,7 @@ solvers
     };
 }
 
+
 PISO
 {
     nCorrectors     4;
diff --git a/tutorials/incompressible/icoDyMFoam/movingBlockRBF/RBFMotionFunction/RBFMotionFunctionObject.C b/tutorials/incompressible/icoDyMFoam/movingBlockRBF/RBFMotionFunction/RBFMotionFunctionObject.C
index e32b02c78..322546055 100644
--- a/tutorials/incompressible/icoDyMFoam/movingBlockRBF/RBFMotionFunction/RBFMotionFunctionObject.C
+++ b/tutorials/incompressible/icoDyMFoam/movingBlockRBF/RBFMotionFunction/RBFMotionFunctionObject.C
@@ -25,7 +25,7 @@ License
 Author
     Frank Bos, TU Delft.  All rights reserved.
 
-\*----------------------------------------------------------------------------*/
+\*---------------------------------------------------------------------------*/
 
 #include "RBFMotionFunctionObject.H"
 #include "addToRunTimeSelectionTable.H"
@@ -56,15 +56,14 @@ Foam::RBFMotionFunctionObject::RBFMotionFunctionObject
     const dictionary& dict
 )
 :
-    functionObject(),
-    name_(name),
+    functionObject(name),
     time_(t),
     regionName_(polyMesh::defaultRegion),
-	rotationAmplitude_(readScalar(dict.lookup("rotationAmplitude"))),
-	rotationFrequency_(readScalar(dict.lookup("rotationFrequency"))),
-	translationAmplitude_(dict.lookup("translationAmplitude")),
-	translationFrequency_(dict.lookup("translationFrequency")),
-	initialRotationOrigin_(dict.lookup("initialRotationOrigin")),
+    rotationAmplitude_(readScalar(dict.lookup("rotationAmplitude"))),
+    rotationFrequency_(readScalar(dict.lookup("rotationFrequency"))),
+    translationAmplitude_(dict.lookup("translationAmplitude")),
+    translationFrequency_(dict.lookup("translationFrequency")),
+    initialRotationOrigin_(dict.lookup("initialRotationOrigin")),
     statPoints_()
 
 {
@@ -86,7 +85,7 @@ bool Foam::RBFMotionFunctionObject::start()
             mesh.lookupObject<RBFMotionSolver>("dynamicMeshDict")
         );
 
- 	statPoints_ = ms.movingPoints();
+    statPoints_ = ms.movingPoints();
 
     return true;
 }
diff --git a/tutorials/incompressible/icoDyMFoam/movingBlockRBF/RBFMotionFunction/RBFMotionFunctionObject.H b/tutorials/incompressible/icoDyMFoam/movingBlockRBF/RBFMotionFunction/RBFMotionFunctionObject.H
index 2e5446064..8b430483e 100644
--- a/tutorials/incompressible/icoDyMFoam/movingBlockRBF/RBFMotionFunction/RBFMotionFunctionObject.H
+++ b/tutorials/incompressible/icoDyMFoam/movingBlockRBF/RBFMotionFunction/RBFMotionFunctionObject.H
@@ -58,20 +58,17 @@ class RBFMotionFunctionObject
 {
     // Private data
 
-        //- Name
-        const word name_;
-
         //- Reference to main object registry
         const Time& time_;
 
         //- Region name
         word regionName_;
 
-		scalar rotationAmplitude_;		
-		scalar rotationFrequency_;
-		vector translationAmplitude_;
-		vector translationFrequency_;
-		vector initialRotationOrigin_;
+        scalar rotationAmplitude_;
+        scalar rotationFrequency_;
+        vector translationAmplitude_;
+        vector translationFrequency_;
+        vector initialRotationOrigin_;
 
         //- Stationary points
         pointField statPoints_;
@@ -113,8 +110,8 @@ public:
 
         //- Read and set the function object if its data has changed
         virtual bool read(const dictionary& dict);
-        
-		//- Update for changes of mesh
+
+        //- Update for changes of mesh
         virtual void movePoints(const pointField&)
         {}
 };
diff --git a/tutorials/incompressible/icoDyMFoam/movingBlockRBF/constant/polyMesh/boundary b/tutorials/incompressible/icoDyMFoam/movingBlockRBF/constant/polyMesh/boundary
index 42d43b7c0..7dc5fb95e 100644
--- a/tutorials/incompressible/icoDyMFoam/movingBlockRBF/constant/polyMesh/boundary
+++ b/tutorials/incompressible/icoDyMFoam/movingBlockRBF/constant/polyMesh/boundary
@@ -1,9 +1,9 @@
 /*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.5-dev                               |
-|   \\  /    A nd           | Revision: 1708                                  |
-|    \\/     M anipulation  | Web:      http://www.OpenFOAM.org               |
+| \\      /  F ield         | OpenFOAM Extend Project: Open source CFD        |
+|  \\    /   O peration     | Version:  1.6-ext                               |
+|   \\  /    A nd           | Web:      www.extend-project.de                 |
+|    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile
 {
diff --git a/tutorials/incompressible/icoDyMFoam/movingBlockRBF/system/controlDict b/tutorials/incompressible/icoDyMFoam/movingBlockRBF/system/controlDict
index f37de14b1..d0a171bbe 100644
--- a/tutorials/incompressible/icoDyMFoam/movingBlockRBF/system/controlDict
+++ b/tutorials/incompressible/icoDyMFoam/movingBlockRBF/system/controlDict
@@ -58,10 +58,10 @@ maxDeltaT 0.25;
 
 functions
 (
-	RBFMotion
-	{
-		type RBFMotion;
-		functionObjectLibs ("libRBFMotionFunction.so");
+    RBFMotion
+    {
+        type RBFMotion;
+        functionObjectLibs ("libRBFMotionFunction.so");
         rotationAmplitude        0.9;
         rotationFrequency        0.25;
         translationAmplitude    (0.0 3.0 0.0);
diff --git a/tutorials/incompressible/icoDyMFoam/movingBlockRBF/system/fvSolution b/tutorials/incompressible/icoDyMFoam/movingBlockRBF/system/fvSolution
index 2966c2c86..a101d4f5e 100644
--- a/tutorials/incompressible/icoDyMFoam/movingBlockRBF/system/fvSolution
+++ b/tutorials/incompressible/icoDyMFoam/movingBlockRBF/system/fvSolution
@@ -18,32 +18,29 @@ FoamFile
 
 solvers
 {
-    p CG
+    p
     {
-        preconditioner
-        {
-            type         DIC;
-        }
+        solver           CG;
+        preconditioner   DIC;
+
         tolerance        1e-06;
         relTol           0;
     };
 
-    pFinal CG
+    pFinal
     {
-        preconditioner
-        {
-            type         DIC;
-        }
+        solver           CG;
+        preconditioner   DIC;
+
         tolerance        1e-06;
         relTol           0;
     };
 
-    U BiCGStab
+    U
     {
-        preconditioner
-        {
-            type         DILU;
-        }
+        solver           BiCGStab;
+        preconditioner   DILU;
+
         tolerance        1e-7;
         relTol           0;
     };
diff --git a/tutorials/incompressible/icoDyMFoam/movingCone/constant/polyMesh/boundary b/tutorials/incompressible/icoDyMFoam/movingCone/constant/polyMesh/boundary
index f27224369..2d0390ca2 100644
--- a/tutorials/incompressible/icoDyMFoam/movingCone/constant/polyMesh/boundary
+++ b/tutorials/incompressible/icoDyMFoam/movingCone/constant/polyMesh/boundary
@@ -1,8 +1,8 @@
 /*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.5                                   |
-|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+| \\      /  F ield         | OpenFOAM Extend Project: Open source CFD        |
+|  \\    /   O peration     | Version:  1.6-ext                               |
+|   \\  /    A nd           | Web:      www.extend-project.de                 |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile
@@ -10,6 +10,7 @@ FoamFile
     version     2.0;
     format      binary;
     class       polyBoundaryMesh;
+    location    "constant/polyMesh";
     object      boundary;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/tutorials/incompressible/icoDyMFoam/movingCone/system/fvSolution b/tutorials/incompressible/icoDyMFoam/movingCone/system/fvSolution
index a348c8239..c6a864c8b 100644
--- a/tutorials/incompressible/icoDyMFoam/movingCone/system/fvSolution
+++ b/tutorials/incompressible/icoDyMFoam/movingCone/system/fvSolution
@@ -16,29 +16,33 @@ FoamFile
 
 solvers
 {
-    pcorr PCG
+    pcorr
     {
+        solver           PCG;
         preconditioner   DIC;
         tolerance        1e-02;
         relTol           0;
     };
 
-    p PCG
+    p
     {
+        solver           PCG;
         preconditioner   DIC;
         tolerance        1e-06;
         relTol           0.05;
     };
 
-    pFinal PCG
+    pFinal
     {
+        solver           PCG;
         preconditioner   DIC;
         tolerance        1e-06;
         relTol           0;
     };
 
-    U PBiCG
+    U
     {
+        solver           PBiCG;
         preconditioner   DILU;
         tolerance        1e-05;
         relTol           0;
diff --git a/tutorials/incompressible/icoDyMFoam/movingCone/system/tetFemSolution b/tutorials/incompressible/icoDyMFoam/movingCone/system/tetFemSolution
deleted file mode 100644
index 7d2e64868..000000000
--- a/tutorials/incompressible/icoDyMFoam/movingCone/system/tetFemSolution
+++ /dev/null
@@ -1,22 +0,0 @@
-/*--------------------------------*- C++ -*----------------------------------*\
-| =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.5                                   |
-|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
-|    \\/     M anipulation  |                                                 |
-\*---------------------------------------------------------------------------*/
-FoamFile
-{
-    version     2.0;
-    format      ascii;
-    class       dictionary;
-    object      tetFemSolution;
-}
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-solvers
-{
-    motionU         ICCG 1e-06 0;
-}
-
-// ************************************************************************* //
diff --git a/tutorials/incompressible/icoDyMFoam/movingConeMotion/Allclean b/tutorials/incompressible/icoDyMFoam/movingConeMotion/Allclean
new file mode 100755
index 000000000..e81dc6208
--- /dev/null
+++ b/tutorials/incompressible/icoDyMFoam/movingConeMotion/Allclean
@@ -0,0 +1,8 @@
+#!/bin/sh
+
+# Source tutorial clean functions
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions
+
+cleanCase
+rm -f 0/motionGamma
+( cd setMotionMovingCone ; wclean )
\ No newline at end of file
diff --git a/tutorials/incompressible/icoDyMFoam/movingConeMotion/Allrun b/tutorials/incompressible/icoDyMFoam/movingConeMotion/Allrun
new file mode 100755
index 000000000..f1ff0ee65
--- /dev/null
+++ b/tutorials/incompressible/icoDyMFoam/movingConeMotion/Allrun
@@ -0,0 +1,11 @@
+#!/bin/sh
+# Source tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+# Get application name from directory
+application="icoDyMFoam"
+
+compileApplication setMotionMovingCone
+runApplication setMotionMovingCone
+runApplication blockMesh
+runApplication $application
diff --git a/tutorials/incompressible/icoDyMFoam/movingConeMotion/constant/polyMesh/boundary b/tutorials/incompressible/icoDyMFoam/movingConeMotion/constant/polyMesh/boundary
index 7f6658410..635f47626 100644
--- a/tutorials/incompressible/icoDyMFoam/movingConeMotion/constant/polyMesh/boundary
+++ b/tutorials/incompressible/icoDyMFoam/movingConeMotion/constant/polyMesh/boundary
@@ -1,84 +1,70 @@
-/*---------------------------------------------------------------------------*\
+/*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.4.1                                 |
-|   \\  /    A nd           | Web:      http://www.openfoam.org               |
+| \\      /  F ield         | OpenFOAM Extend Project: Open source CFD        |
+|  \\    /   O peration     | Version:  1.6-ext                               |
+|   \\  /    A nd           | Web:      www.extend-project.de                 |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
-
 FoamFile
 {
-    version 2.0;
-    format ascii;
-
-    root "/home/hjasak/OpenFOAM/OpenFOAM-1.4.1/tutorials/icoDyMFoam";
-    case "movingConeMotion";
-    instance "constant";
-    local "polyMesh";
-
-    class polyBoundaryMesh;
-    object boundary;
+    version     2.0;
+    format      ascii;
+    class       polyBoundaryMesh;
+    location    "constant/polyMesh";
+    object      boundary;
 }
-
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 8
 (
-movingWall
-{
-    type wall;
-    nFaces 85;
-    startFace 3665;
-}
-
-farFieldMoving
-{
-    type patch;
-    nFaces 50;
-    startFace 3750;
-}
-
-fixedWall
-{
-    type wall;
-    nFaces 35;
-    startFace 3800;
-}
-
-axis
-{
-    type empty;
-    nFaces 0;
-    startFace 3835;
-}
-
-left
-{
-    type patch;
-    nFaces 30;
-    startFace 3835;
-}
-
-farField
-{
-    type patch;
-    nFaces 35;
-    startFace 3865;
-}
-
-back
-{
-    type wedge;
-    nFaces 1900;
-    startFace 3900;
-}
-
-front
-{
-    type wedge;
-    nFaces 1900;
-    startFace 5800;
-}
+    movingWall
+    {
+        type            wall;
+        nFaces          85;
+        startFace       3665;
+    }
+    farFieldMoving
+    {
+        type            patch;
+        nFaces          50;
+        startFace       3750;
+    }
+    fixedWall
+    {
+        type            wall;
+        nFaces          35;
+        startFace       3800;
+    }
+    axis
+    {
+        type            empty;
+        nFaces          0;
+        startFace       3835;
+    }
+    left
+    {
+        type            patch;
+        nFaces          30;
+        startFace       3835;
+    }
+    farField
+    {
+        type            patch;
+        nFaces          35;
+        startFace       3865;
+    }
+    back
+    {
+        type            wedge;
+        nFaces          1900;
+        startFace       3900;
+    }
+    front
+    {
+        type            wedge;
+        nFaces          1900;
+        startFace       5800;
+    }
 )
 
 // ************************************************************************* //
diff --git a/tutorials/incompressible/icoDyMFoam/movingConeMotion/system/fvSolution b/tutorials/incompressible/icoDyMFoam/movingConeMotion/system/fvSolution
index 9ebe48d0d..88351fb1b 100644
--- a/tutorials/incompressible/icoDyMFoam/movingConeMotion/system/fvSolution
+++ b/tutorials/incompressible/icoDyMFoam/movingConeMotion/system/fvSolution
@@ -1,4 +1,4 @@
-/*---------------------------------------------------------------------------*\
+/*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
 | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
 |  \\    /   O peration     | Version:  1.4                                   |
@@ -10,12 +10,6 @@ FoamFile
 {
     version         2.0;
     format          ascii;
-
-    root            "";
-    case            "";
-    instance        "";
-    local           "";
-
     class           dictionary;
     object          fvSolution;
 }
@@ -24,12 +18,10 @@ FoamFile
 
 solvers
 {
-    pcorr PCG
+    pcorr
     {
-        preconditioner
-        {
-            type         DIC;
-        }
+        solver           PCG;
+        preconditioner   DIC;
 
         minIter          0;
         maxIter          1000;
@@ -37,12 +29,10 @@ solvers
         relTol           0;
     };
 
-    p PCG
+    p
     {
-        preconditioner
-        {
-            type         DIC;
-        }
+        solver           PCG;
+        preconditioner   DIC;
 
         minIter          0;
         maxIter          1000;
@@ -50,12 +40,10 @@ solvers
         relTol           0.05;
     };
 
-    pFinal PCG
+    pFinal
     {
-        preconditioner
-        {
-            type         DIC;
-        }
+        solver           PCG;
+        preconditioner   DIC;
 
         minIter          0;
         maxIter          1000;
@@ -63,33 +51,19 @@ solvers
         relTol           0;
     };
 
-    U PBiCG
+    U
     {
-        preconditioner   
-        {
-            type         DILU;
-        }
+        solver           PBiCG;
+        preconditioner   DILU;
 
         minIter          0;
         maxIter          1000;
         tolerance        1e-05;
         relTol           0;
     };
-
-    cellMotionUx PCG
-    {
-        preconditioner      
-        {
-            type         DIC;
-        }
-
-        minIter          0;
-        maxIter          1000;
-        tolerance        1e-08;
-        relTol           0;
-    };
 }
 
+
 PISO
 {
     nCorrectors     2;
diff --git a/tutorials/incompressible/icoDyMFoam/movingConeMotion/system/tetFemSolution b/tutorials/incompressible/icoDyMFoam/movingConeMotion/system/tetFemSolution
index e3826c1ef..0c354c04e 100644
--- a/tutorials/incompressible/icoDyMFoam/movingConeMotion/system/tetFemSolution
+++ b/tutorials/incompressible/icoDyMFoam/movingConeMotion/system/tetFemSolution
@@ -18,8 +18,9 @@ FoamFile
 
 solvers
 {
-    motionU amgSolver
+    motionU
     {
+        solver          amgSolver;
         cycle           W-cycle;
         policy          AAMG;
         nPreSweeps      0;
diff --git a/tutorials/incompressible/icoDyMFoam/movingConeTopo/constant/polyMesh/boundary b/tutorials/incompressible/icoDyMFoam/movingConeTopo/constant/polyMesh/boundary
index 43b7d9ef5..0afd3439c 100644
--- a/tutorials/incompressible/icoDyMFoam/movingConeTopo/constant/polyMesh/boundary
+++ b/tutorials/incompressible/icoDyMFoam/movingConeTopo/constant/polyMesh/boundary
@@ -1,8 +1,8 @@
 /*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.5-dev                               |
-|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+| \\      /  F ield         | OpenFOAM Extend Project: Open source CFD        |
+|  \\    /   O peration     | Version:  1.6-ext                               |
+|   \\  /    A nd           | Web:      www.extend-project.de                 |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile
diff --git a/tutorials/incompressible/icoDyMFoam/movingConeTopo/system/fvSolution b/tutorials/incompressible/icoDyMFoam/movingConeTopo/system/fvSolution
index 07a7b90fc..1c5f4e3aa 100644
--- a/tutorials/incompressible/icoDyMFoam/movingConeTopo/system/fvSolution
+++ b/tutorials/incompressible/icoDyMFoam/movingConeTopo/system/fvSolution
@@ -18,12 +18,10 @@ FoamFile
 
 solvers
 {
-    pcorr PCG
+    pcorr
     {
-        preconditioner
-        {
-            type         DIC;
-        }
+        solver           PCG;
+        preconditioner   DIC;
 
         minIter          0;
         maxIter          1000;
@@ -31,12 +29,10 @@ solvers
         relTol           0;
     };
 
-    p PCG
+    p
     {
-        preconditioner
-        {
-            type         DIC;
-        }
+        solver           PCG;
+        preconditioner   DIC;
 
         minIter          0;
         maxIter          1000;
@@ -44,12 +40,10 @@ solvers
         relTol           0.0;
     };
 
-    pFinal PCG
+    pFinal
     {
-        preconditioner
-        {
-            type         DIC;
-        }
+        solver           PCG;
+        preconditioner   DIC;
 
         minIter          0;
         maxIter          1000;
@@ -57,12 +51,10 @@ solvers
         relTol           0;
     };
 
-    U PBiCG
+    U
     {
-        preconditioner   
-        {
-            type         DILU;
-        }
+        solver           PBiCG;
+        preconditioner   DILU;
 
         minIter          0;
         maxIter          1000;
diff --git a/tutorials/incompressible/icoFoam/cavity/constant/polyMesh/boundary b/tutorials/incompressible/icoFoam/cavity/constant/polyMesh/boundary
index cc15fe93f..e5f158a50 100644
--- a/tutorials/incompressible/icoFoam/cavity/constant/polyMesh/boundary
+++ b/tutorials/incompressible/icoFoam/cavity/constant/polyMesh/boundary
@@ -1,8 +1,8 @@
 /*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.5                                   |
-|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+| \\      /  F ield         | OpenFOAM Extend Project: Open source CFD        |
+|  \\    /   O peration     | Version:  1.6-ext                               |
+|   \\  /    A nd           | Web:      www.extend-project.de                 |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile
@@ -10,6 +10,7 @@ FoamFile
     version     2.0;
     format      ascii;
     class       polyBoundaryMesh;
+    location    "constant/polyMesh";
     object      boundary;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/tutorials/incompressible/icoFoam/cavityClipped/constant/polyMesh/boundary b/tutorials/incompressible/icoFoam/cavityClipped/constant/polyMesh/boundary
index 63b839343..160211c9a 100644
--- a/tutorials/incompressible/icoFoam/cavityClipped/constant/polyMesh/boundary
+++ b/tutorials/incompressible/icoFoam/cavityClipped/constant/polyMesh/boundary
@@ -1,8 +1,8 @@
 /*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.5                                   |
-|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+| \\      /  F ield         | OpenFOAM Extend Project: Open source CFD        |
+|  \\    /   O peration     | Version:  1.6-ext                               |
+|   \\  /    A nd           | Web:      www.extend-project.de                 |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile
@@ -10,6 +10,7 @@ FoamFile
     version     2.0;
     format      ascii;
     class       polyBoundaryMesh;
+    location    "constant/polyMesh";
     object      boundary;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/tutorials/incompressible/icoFoam/cavityGrade/constant/polyMesh/boundary b/tutorials/incompressible/icoFoam/cavityGrade/constant/polyMesh/boundary
index cc15fe93f..e5f158a50 100644
--- a/tutorials/incompressible/icoFoam/cavityGrade/constant/polyMesh/boundary
+++ b/tutorials/incompressible/icoFoam/cavityGrade/constant/polyMesh/boundary
@@ -1,8 +1,8 @@
 /*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.5                                   |
-|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+| \\      /  F ield         | OpenFOAM Extend Project: Open source CFD        |
+|  \\    /   O peration     | Version:  1.6-ext                               |
+|   \\  /    A nd           | Web:      www.extend-project.de                 |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile
@@ -10,6 +10,7 @@ FoamFile
     version     2.0;
     format      ascii;
     class       polyBoundaryMesh;
+    location    "constant/polyMesh";
     object      boundary;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/tutorials/incompressible/icoFoam/elbow/constant/polyMesh/boundary b/tutorials/incompressible/icoFoam/elbow/constant/polyMesh/boundary
index 0b3c20490..2f02a227f 100644
--- a/tutorials/incompressible/icoFoam/elbow/constant/polyMesh/boundary
+++ b/tutorials/incompressible/icoFoam/elbow/constant/polyMesh/boundary
@@ -1,8 +1,8 @@
 /*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.5                                   |
-|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+| \\      /  F ield         | OpenFOAM Extend Project: Open source CFD        |
+|  \\    /   O peration     | Version:  1.6-ext                               |
+|   \\  /    A nd           | Web:      www.extend-project.de                 |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile
@@ -10,6 +10,7 @@ FoamFile
     version     2.0;
     format      ascii;
     class       polyBoundaryMesh;
+    location    "constant/polyMesh";
     object      boundary;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/tutorials/viscoelastic/viscoelasticFluidFoam/Allclean b/tutorials/viscoelastic/viscoelasticFluidFoam/Allclean
deleted file mode 100755
index 66fe1458f..000000000
--- a/tutorials/viscoelastic/viscoelasticFluidFoam/Allclean
+++ /dev/null
@@ -1,56 +0,0 @@
-#!/bin/sh
-
-cd PTT-Exponential
-   ./Allclean
-
-cd ..
-cd DCPP
-   ./Allclean
-
-cd ..
-cd Giesekus
-   ./Allclean
-
-
-cd ..
-cd FENE-CR
-   ./Allclean
-
-cd ..
-cd FENE-P
-   ./Allclean
-
-cd ..
-cd Feta-PTT
-   ./Allclean
-
-cd ..
-cd Leonov
-   ./Allclean
-
-cd ..
-cd Oldroyd-B
-   ./Allclean
-
-cd ..
-cd S-MDCPP
-   ./Allclean
-
-cd ..
-cd PTT-Linear
-   ./Allclean
-
-cd ..
-cd WhiteMetzner/CarreauYasuda
-   ./Allclean
-
-cd ../..
-cd XPP_DE
-   ./Allclean
-
-cd ..
-cd XPP_SE
-   ./Allclean
-
-cd ..
-
diff --git a/tutorials/viscoelastic/viscoelasticFluidFoam/Allrun b/tutorials/viscoelastic/viscoelasticFluidFoam/Allrun
deleted file mode 100755
index ac2754b9b..000000000
--- a/tutorials/viscoelastic/viscoelasticFluidFoam/Allrun
+++ /dev/null
@@ -1,55 +0,0 @@
-#!/bin/sh
-
-cd PTT-Exponential
-   ./Allrun
-
-cd ..
-cd DCPP
-   ./Allrun
-
-cd ..
-cd Giesekus
-   ./Allrun
-
-
-cd ..
-cd FENE-CR
-   ./Allrun
-
-cd ..
-cd FENE-P
-   ./Allrun
-
-cd ..
-cd Feta-PTT
-   ./Allrun
-
-cd ..
-cd Leonov
-   ./Allrun
-
-cd ..
-cd Oldroyd-B
-   ./Allrun
-
-cd ..
-cd S-MDCPP
-   ./Allrun
-
-cd ..
-cd PTT-Linear
-   ./Allrun
-
-cd ..
-cd White-Metzner/Carreau-Yasuda
-   ./Allrun
-
-cd ../..
-cd XPP_DE
-   ./Allrun
-
-cd ..
-cd XPP_SE
-   ./Allrun
-
-cd ..
diff --git a/tutorials/viscoelastic/viscoelasticFluidFoam/DCPP/Allclean b/tutorials/viscoelastic/viscoelasticFluidFoam/DCPP/Allclean
index 20a51e184..917056ad6 100755
--- a/tutorials/viscoelastic/viscoelasticFluidFoam/DCPP/Allclean
+++ b/tutorials/viscoelastic/viscoelasticFluidFoam/DCPP/Allclean
@@ -3,7 +3,5 @@
 # Source tutorial clean functions
 . $WM_PROJECT_DIR/bin/tools/CleanFunctions
 
- cleanCase
- rm -R -f sets 
- echo "done!"
-
+cleanCase
+rm -rf sets 
diff --git a/tutorials/viscoelastic/viscoelasticFluidFoam/FENE-CR/Allclean b/tutorials/viscoelastic/viscoelasticFluidFoam/FENE-CR/Allclean
index 0d957e05d..e80d74cac 100755
--- a/tutorials/viscoelastic/viscoelasticFluidFoam/FENE-CR/Allclean
+++ b/tutorials/viscoelastic/viscoelasticFluidFoam/FENE-CR/Allclean
@@ -3,6 +3,5 @@
 # Source tutorial clean functions
 . $WM_PROJECT_DIR/bin/tools/CleanFunctions
 
- cleanCase
- echo "done!"
+cleanCase
 
diff --git a/tutorials/viscoelastic/viscoelasticFluidFoam/FENE-P/Allclean b/tutorials/viscoelastic/viscoelasticFluidFoam/FENE-P/Allclean
index 20a51e184..83a0d7beb 100755
--- a/tutorials/viscoelastic/viscoelasticFluidFoam/FENE-P/Allclean
+++ b/tutorials/viscoelastic/viscoelasticFluidFoam/FENE-P/Allclean
@@ -3,7 +3,6 @@
 # Source tutorial clean functions
 . $WM_PROJECT_DIR/bin/tools/CleanFunctions
 
- cleanCase
- rm -R -f sets 
- echo "done!"
+cleanCase
+rm -rf sets 
 
diff --git a/tutorials/viscoelastic/viscoelasticFluidFoam/Feta-PTT/Allclean b/tutorials/viscoelastic/viscoelasticFluidFoam/Feta-PTT/Allclean
index 0d957e05d..e80d74cac 100755
--- a/tutorials/viscoelastic/viscoelasticFluidFoam/Feta-PTT/Allclean
+++ b/tutorials/viscoelastic/viscoelasticFluidFoam/Feta-PTT/Allclean
@@ -3,6 +3,5 @@
 # Source tutorial clean functions
 . $WM_PROJECT_DIR/bin/tools/CleanFunctions
 
- cleanCase
- echo "done!"
+cleanCase
 
diff --git a/tutorials/viscoelastic/viscoelasticFluidFoam/Giesekus/Allclean b/tutorials/viscoelastic/viscoelasticFluidFoam/Giesekus/Allclean
index 20a51e184..83a0d7beb 100755
--- a/tutorials/viscoelastic/viscoelasticFluidFoam/Giesekus/Allclean
+++ b/tutorials/viscoelastic/viscoelasticFluidFoam/Giesekus/Allclean
@@ -3,7 +3,6 @@
 # Source tutorial clean functions
 . $WM_PROJECT_DIR/bin/tools/CleanFunctions
 
- cleanCase
- rm -R -f sets 
- echo "done!"
+cleanCase
+rm -rf sets 
 
diff --git a/tutorials/viscoelastic/viscoelasticFluidFoam/Leonov/Allclean b/tutorials/viscoelastic/viscoelasticFluidFoam/Leonov/Allclean
index 0d957e05d..7d7f43b1c 100755
--- a/tutorials/viscoelastic/viscoelasticFluidFoam/Leonov/Allclean
+++ b/tutorials/viscoelastic/viscoelasticFluidFoam/Leonov/Allclean
@@ -3,6 +3,6 @@
 # Source tutorial clean functions
 . $WM_PROJECT_DIR/bin/tools/CleanFunctions
 
- cleanCase
- echo "done!"
+cleanCase
+echo "done!"
 
diff --git a/tutorials/viscoelastic/viscoelasticFluidFoam/Oldroyd-B/Allclean b/tutorials/viscoelastic/viscoelasticFluidFoam/Oldroyd-B/Allclean
index 20a51e184..304161076 100755
--- a/tutorials/viscoelastic/viscoelasticFluidFoam/Oldroyd-B/Allclean
+++ b/tutorials/viscoelastic/viscoelasticFluidFoam/Oldroyd-B/Allclean
@@ -3,7 +3,6 @@
 # Source tutorial clean functions
 . $WM_PROJECT_DIR/bin/tools/CleanFunctions
 
- cleanCase
- rm -R -f sets 
- echo "done!"
+cleanCase
+rm -R -f sets 
 
diff --git a/tutorials/viscoelastic/viscoelasticFluidFoam/PTT-Exponential/Allclean b/tutorials/viscoelastic/viscoelasticFluidFoam/PTT-Exponential/Allclean
index 0d957e05d..e80d74cac 100755
--- a/tutorials/viscoelastic/viscoelasticFluidFoam/PTT-Exponential/Allclean
+++ b/tutorials/viscoelastic/viscoelasticFluidFoam/PTT-Exponential/Allclean
@@ -3,6 +3,5 @@
 # Source tutorial clean functions
 . $WM_PROJECT_DIR/bin/tools/CleanFunctions
 
- cleanCase
- echo "done!"
+cleanCase
 
diff --git a/tutorials/viscoelastic/viscoelasticFluidFoam/PTT-Linear/Allclean b/tutorials/viscoelastic/viscoelasticFluidFoam/PTT-Linear/Allclean
index 20a51e184..46ce878eb 100755
--- a/tutorials/viscoelastic/viscoelasticFluidFoam/PTT-Linear/Allclean
+++ b/tutorials/viscoelastic/viscoelasticFluidFoam/PTT-Linear/Allclean
@@ -3,7 +3,6 @@
 # Source tutorial clean functions
 . $WM_PROJECT_DIR/bin/tools/CleanFunctions
 
- cleanCase
- rm -R -f sets 
- echo "done!"
+cleanCase
+rm -Rf sets 
 
diff --git a/tutorials/viscoelastic/viscoelasticFluidFoam/S-MDCPP/Allclean b/tutorials/viscoelastic/viscoelasticFluidFoam/S-MDCPP/Allclean
index 0d957e05d..e80d74cac 100755
--- a/tutorials/viscoelastic/viscoelasticFluidFoam/S-MDCPP/Allclean
+++ b/tutorials/viscoelastic/viscoelasticFluidFoam/S-MDCPP/Allclean
@@ -3,6 +3,5 @@
 # Source tutorial clean functions
 . $WM_PROJECT_DIR/bin/tools/CleanFunctions
 
- cleanCase
- echo "done!"
+cleanCase
 
diff --git a/tutorials/viscoelastic/viscoelasticFluidFoam/XPP_DE/Allclean b/tutorials/viscoelastic/viscoelasticFluidFoam/XPP_DE/Allclean
index 0d957e05d..e80d74cac 100755
--- a/tutorials/viscoelastic/viscoelasticFluidFoam/XPP_DE/Allclean
+++ b/tutorials/viscoelastic/viscoelasticFluidFoam/XPP_DE/Allclean
@@ -3,6 +3,5 @@
 # Source tutorial clean functions
 . $WM_PROJECT_DIR/bin/tools/CleanFunctions
 
- cleanCase
- echo "done!"
+cleanCase
 
diff --git a/tutorials/viscoelastic/viscoelasticFluidFoam/XPP_SE/Allclean b/tutorials/viscoelastic/viscoelasticFluidFoam/XPP_SE/Allclean
index 0d957e05d..a690e883f 100755
--- a/tutorials/viscoelastic/viscoelasticFluidFoam/XPP_SE/Allclean
+++ b/tutorials/viscoelastic/viscoelasticFluidFoam/XPP_SE/Allclean
@@ -4,5 +4,4 @@
 . $WM_PROJECT_DIR/bin/tools/CleanFunctions
 
  cleanCase
- echo "done!"
 

From 713c27c5c0a812b0344406dfdda49b757d50ebea Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Wed, 13 Oct 2010 12:39:00 +0100
Subject: [PATCH 6/9] Simple dynamic mesh solver

---
 .../icoDyMSimpleFoam/correctPhi.H             |  8 +------
 .../icoDyMSimpleFoam/createFields.H           | 24 +++++++++++++++++--
 .../icoDyMSimpleFoam/icoDyMSimpleFoam.C       |  1 +
 3 files changed, 24 insertions(+), 9 deletions(-)

diff --git a/applications/solvers/incompressible/icoDyMSimpleFoam/correctPhi.H b/applications/solvers/incompressible/icoDyMSimpleFoam/correctPhi.H
index 98092f501..69978f0e6 100644
--- a/applications/solvers/incompressible/icoDyMSimpleFoam/correctPhi.H
+++ b/applications/solvers/incompressible/icoDyMSimpleFoam/correctPhi.H
@@ -24,12 +24,6 @@
         pcorrTypes
     );
 
-#   include "continuityErrs.H"
-
-    // Flux predictor
-    phi = (fvc::interpolate(U) & mesh.Sf());
-    rAU == runTime.deltaT();
-
     for(int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
     {
         fvScalarMatrix pcorrEqn
@@ -50,4 +44,4 @@
     }
 }
 
-
+#include "continuityErrs.H"
diff --git a/applications/solvers/incompressible/icoDyMSimpleFoam/createFields.H b/applications/solvers/incompressible/icoDyMSimpleFoam/createFields.H
index e9c30de89..561cd7bf4 100644
--- a/applications/solvers/incompressible/icoDyMSimpleFoam/createFields.H
+++ b/applications/solvers/incompressible/icoDyMSimpleFoam/createFields.H
@@ -54,17 +54,37 @@
     scalar pRefValue = 0.0;
     setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);
 
-    scalar totalVolume = sum(mesh.V()).value();
 
+    Info<< "Reading field rAU if present\n" << endl;
     volScalarField rAU
     (
         IOobject
         (
             "rAU",
             runTime.timeName(),
-            mesh
+            mesh,
+            IOobject::READ_IF_PRESENT,
+            IOobject::AUTO_WRITE
         ),
         mesh,
         runTime.deltaT(),
         zeroGradientFvPatchScalarField::typeName
     );
+
+
+    {
+        const dictionary& piso = mesh.solutionDict().subDict("PISO");
+
+        bool momentumPredictor = true;
+        if (piso.found("momentumPredictor"))
+        {
+            momentumPredictor = Switch(piso.lookup("momentumPredictor"));
+        }
+
+        bool momentumPredictorSave = momentumPredictor;
+        momentumPredictor = false;
+#       include "UEqn.H"
+        momentumPredictor = momentumPredictorSave;
+        
+        rAU = 1.0/UEqn.A();
+    }
diff --git a/applications/solvers/incompressible/icoDyMSimpleFoam/icoDyMSimpleFoam.C b/applications/solvers/incompressible/icoDyMSimpleFoam/icoDyMSimpleFoam.C
index ffd88b97d..517f57e2a 100644
--- a/applications/solvers/incompressible/icoDyMSimpleFoam/icoDyMSimpleFoam.C
+++ b/applications/solvers/incompressible/icoDyMSimpleFoam/icoDyMSimpleFoam.C
@@ -46,6 +46,7 @@ int main(int argc, char *argv[])
 #   include "createTime.H"
 #   include "createDynamicFvMesh.H"
 #   include "initContinuityErrs.H"
+#   include "initTotalVolume.H"
 #   include "createFields.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

From 7f631ac3c672d85caab6dd4c5ac132cac281c396 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Wed, 13 Oct 2010 12:39:40 +0100
Subject: [PATCH 7/9] Updating call-back functionality on meshObjects: dynamic
 mesh support

---
 .../GAMGAgglomeration/GAMGAgglomeration.H          |  7 ++++---
 src/OpenFOAM/meshes/pointMesh/pointMesh.H          |  2 +-
 src/OpenFOAM/meshes/polyMesh/polyMesh.C            |  1 -
 src/OpenFOAM/meshes/polyMesh/polyMeshClear.C       |  4 ++++
 src/OpenFOAM/meshes/polyMesh/polyMeshIO.C          | 14 ++++++++++++--
 src/finiteVolume/cfdTools/general/bound/bound.H    |  1 +
 src/finiteVolume/fvMesh/fvMesh.C                   | 10 ++--------
 .../schemes/skewCorrected/skewCorrectionVectors.H  |  2 +-
 .../volPointInterpolation/volPointInterpolation.H  |  2 +-
 9 files changed, 26 insertions(+), 17 deletions(-)

diff --git a/src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/GAMGAgglomerations/GAMGAgglomeration/GAMGAgglomeration.H b/src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/GAMGAgglomerations/GAMGAgglomeration/GAMGAgglomeration.H
index 17a100c3f..05e9da8fc 100644
--- a/src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/GAMGAgglomerations/GAMGAgglomeration/GAMGAgglomeration.H
+++ b/src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/GAMGAgglomerations/GAMGAgglomeration/GAMGAgglomeration.H
@@ -81,10 +81,11 @@ protected:
 
         //- Face restriction addressing array.
         //  Maps from the finer to the coarser level.
-        //  Positive indices map the finer faces which form part of the boundary
-        //  of the coarser cells to the corresponding coarser cell face.
+        //  Positive indices map the finer faces which form part of the
+        //   boundary of the coarser cells to the corresponding coarser
+        //  cell face.
         //  Negative indices map the finer faces which are internal to the
-        //  coarser cells to minus the corresponding coarser cell index minus 1.
+        //  coarser cells to minus the corresponding coarser cell index minus 1
         PtrList<labelList> faceRestrictAddressing_;
 
         //- Hierarchy of mesh addressing
diff --git a/src/OpenFOAM/meshes/pointMesh/pointMesh.H b/src/OpenFOAM/meshes/pointMesh/pointMesh.H
index 2a33965c3..2ab1ecb92 100644
--- a/src/OpenFOAM/meshes/pointMesh/pointMesh.H
+++ b/src/OpenFOAM/meshes/pointMesh/pointMesh.H
@@ -73,7 +73,7 @@ class pointMesh
 public:
 
     // Declare name of the class and its debug switch
-    ClassName("pointMesh");
+    TypeName("pointMesh");
 
     typedef pointMesh Mesh;
     typedef pointBoundaryMesh BoundaryMesh;
diff --git a/src/OpenFOAM/meshes/polyMesh/polyMesh.C b/src/OpenFOAM/meshes/polyMesh/polyMesh.C
index 332f38e2f..375922d43 100644
--- a/src/OpenFOAM/meshes/polyMesh/polyMesh.C
+++ b/src/OpenFOAM/meshes/polyMesh/polyMesh.C
@@ -1226,7 +1226,6 @@ void Foam::polyMesh::setOldPoints
     oldPointsPtr_ = new pointField::subField(oldAllPoints(), nPoints());
     curMotionTimeIndex_ = 0;
     primitiveMesh::clearGeom();
-
 }
 
 
diff --git a/src/OpenFOAM/meshes/polyMesh/polyMeshClear.C b/src/OpenFOAM/meshes/polyMesh/polyMeshClear.C
index 1b6774acd..50f499d9a 100644
--- a/src/OpenFOAM/meshes/polyMesh/polyMeshClear.C
+++ b/src/OpenFOAM/meshes/polyMesh/polyMeshClear.C
@@ -28,6 +28,7 @@ License
 #include "primitiveMesh.H"
 #include "globalMeshData.H"
 #include "demandDrivenData.H"
+#include "meshObjectBase.H"
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
@@ -90,6 +91,9 @@ void Foam::polyMesh::clearGeom()
     // Reset valid directions (could change with rotation)
     geometricD_ = Vector<label>::zero;
     solutionD_ = Vector<label>::zero;
+
+    // Move points all mesh objects.  HJ, 13/Oct/2010
+    meshObjectBase::allMovePoints(*this);
 }
 
 
diff --git a/src/OpenFOAM/meshes/polyMesh/polyMeshIO.C b/src/OpenFOAM/meshes/polyMesh/polyMeshIO.C
index 49ad58aae..8e0ce6e91 100644
--- a/src/OpenFOAM/meshes/polyMesh/polyMeshIO.C
+++ b/src/OpenFOAM/meshes/polyMesh/polyMeshIO.C
@@ -27,6 +27,8 @@ License
 #include "polyMesh.H"
 #include "Time.H"
 #include "cellIOList.H"
+#include "meshObjectBase.H"
+#include "mapPolyMesh.H"
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
@@ -394,6 +396,11 @@ Foam::polyMesh::readUpdateState Foam::polyMesh::readUpdate()
             cellZones_.set(czI, newCellZones[czI].clone(cellZones_));
         }
 
+        // Instantiate a dummy mapPolyMesh
+        autoPtr<mapPolyMesh> mapPtr(new mapPolyMesh(*this));
+
+        // Execute dummy  topo change on all mesh objects
+        meshObjectBase::allUpdateTopology(*this, mapPtr());
 
         if (boundaryChanged)
         {
@@ -432,14 +439,17 @@ Foam::polyMesh::readUpdateState Foam::polyMesh::readUpdate()
 
         // Reset points, mesh is not moved
         points_ = pointField::subField(allPoints_, nPoints());
-        
+
         // Derived info
         bounds_ = boundBox(allPoints_);
 
         // Rotation can cause direction vector to change
         geometricD_ = Vector<label>::zero;
         solutionD_ = Vector<label>::zero;
-        
+
+        // Move points in all mesh objects
+        meshObjectBase::allMovePoints<polyMesh>(*this);
+
         return polyMesh::POINTS_MOVED;
     }
     else
diff --git a/src/finiteVolume/cfdTools/general/bound/bound.H b/src/finiteVolume/cfdTools/general/bound/bound.H
index a49943fa3..29609f7ab 100644
--- a/src/finiteVolume/cfdTools/general/bound/bound.H
+++ b/src/finiteVolume/cfdTools/general/bound/bound.H
@@ -59,6 +59,7 @@ void boundMinMax
     const dimensionedScalar& vsf1
 );
 
+
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 } // End namespace Foam
diff --git a/src/finiteVolume/fvMesh/fvMesh.C b/src/finiteVolume/fvMesh/fvMesh.C
index 587386067..556f93363 100644
--- a/src/finiteVolume/fvMesh/fvMesh.C
+++ b/src/finiteVolume/fvMesh/fvMesh.C
@@ -37,7 +37,6 @@ License
 #include "MapFvFields.H"
 #include "fvMeshMapper.H"
 #include "mapClouds.H"
-#include "meshObjectBase.H"
 
 #include "volPointInterpolation.H"
 #include "extendedLeastSquaresVectors.H"
@@ -78,9 +77,8 @@ void Foam::fvMesh::clearGeom()
     // needs to be saved.
 
     // Geometry dependent object updated through call-back
-    // "Reserve" optional delete.  Reconsider
+    // and handled by polyMesh
     // HJ, 29/Aug/2010
-//     meshObjectBase::allDelete(*this);
 }
 
 
@@ -89,9 +87,8 @@ void Foam::fvMesh::clearAddressing()
     deleteDemandDrivenData(lduPtr_);
 
     // Geometry dependent object updated through call-back
-    // "Reserve" optional delete.  Reconsider
+    // and handled by polyMesh
     // HJ, 29/Aug/2010
-//     meshObjectBase::allDelete(*this);
 }
 
 
@@ -528,9 +525,6 @@ void  Foam::fvMesh::syncUpdateMesh()
     // This is a temporary solution
     surfaceInterpolation::movePoints();
 
-    // Instantiate a dummy mapPolyMesh
-    autoPtr<mapPolyMesh> mapPtr(new mapPolyMesh(*this));
-
     // Function object update moved to polyMesh
     // HJ, 29/Aug/2010
 }
diff --git a/src/finiteVolume/interpolation/surfaceInterpolation/schemes/skewCorrected/skewCorrectionVectors.H b/src/finiteVolume/interpolation/surfaceInterpolation/schemes/skewCorrected/skewCorrectionVectors.H
index 769b60a4f..52df3b587 100644
--- a/src/finiteVolume/interpolation/surfaceInterpolation/schemes/skewCorrected/skewCorrectionVectors.H
+++ b/src/finiteVolume/interpolation/surfaceInterpolation/schemes/skewCorrected/skewCorrectionVectors.H
@@ -57,7 +57,7 @@ class skewCorrectionVectors
     public MeshObject<fvMesh, skewCorrectionVectors>
 {
     // Private data
-    
+
         //- Is mesh skew
         mutable bool skew_;
 
diff --git a/src/finiteVolume/interpolation/volPointInterpolation/volPointInterpolation.H b/src/finiteVolume/interpolation/volPointInterpolation/volPointInterpolation.H
index c9510e2d2..c6fd032ae 100644
--- a/src/finiteVolume/interpolation/volPointInterpolation/volPointInterpolation.H
+++ b/src/finiteVolume/interpolation/volPointInterpolation/volPointInterpolation.H
@@ -81,7 +81,7 @@ class volPointInterpolation
 public:
 
     // Declare name of the class and its debug switch
-    ClassName("volPointInterpolation");
+    TypeName("volPointInterpolation");
 
 
     // Constructors

From ebaaefa484bdc00bf4c6bee6536f776646f54650 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Wed, 13 Oct 2010 12:39:54 +0100
Subject: [PATCH 8/9] Tutorial update

---
 bin/tools/CleanFunctions | 20 +++++++++++++++-----
 bin/tools/RunFunctions   |  6 ++++++
 2 files changed, 21 insertions(+), 5 deletions(-)

diff --git a/bin/tools/CleanFunctions b/bin/tools/CleanFunctions
index 9f85a0f6b..bc1e090f5 100644
--- a/bin/tools/CleanFunctions
+++ b/bin/tools/CleanFunctions
@@ -55,7 +55,7 @@ cleanTimeDirectories ()
         zeros=`printf %0${nZeros}d 0`
         nZeros=$(($nZeros + 1))
     done
-    rm -rf ./{[1-9]*,-[1-9]*,log,log.*,log-*,logSummary.*,.fxLock,*.xml,ParaView*,paraFoam*,*.OpenFOAM} > /dev/null 2>&1
+    rm -rf ./[1-9]* ./-[1-9]* ./log ./log.* ./log-* ./logSummary.* ./.fxLock ./*.xml ./ParaView* ./paraFoam* ./*.OpenFOAM > /dev/null 2>&1
 }
 
 cleanCase ()
@@ -66,10 +66,20 @@ cleanCase ()
     rm -rf forces* > /dev/null 2>&1
 
     rm -rf system/machines \
-           constant/polyMesh/{allOwner*,cell*,face*,meshModifiers*} \
-           constant/polyMesh/{owner*,neighbour*,point*,edge*} \
-           constant/polyMesh/{cellLevel*,pointLevel*,refinementHistory*,surfaceIndex*} \
-           constant/{cellToRegion,cellLevel*,pointLevel*} \
+           constant/polyMesh/allOwner* \
+           constant/polyMesh/cell* \
+           constant/polyMesh/face* \
+           constant/polyMesh/meshModifiers* \
+           constant/polyMesh/owner* \
+           constant/polyMesh/neighbour* \
+           constant/polyMesh/point* \
+           constant/polyMesh/edge* \
+           constant/polyMesh/zoneToPatchName \
+           constant/polyMesh/cellLevel* \
+           constant/polyMesh/pointLevel* \
+           constant/polyMesh/refinementHistory*, \
+           constant/polyMesh/surfaceIndex* \
+           constant/cellToRegion \
            constant/polyMesh/sets/ \
            VTK \
            > /dev/null 2>&1
diff --git a/bin/tools/RunFunctions b/bin/tools/RunFunctions
index 592ef6c4e..0a2b21812 100644
--- a/bin/tools/RunFunctions
+++ b/bin/tools/RunFunctions
@@ -65,6 +65,12 @@ compileApplication ()
     wmake $1
 }
 
+compileLibrary ()
+{
+    echo "Compiling $1 application"
+    wmake libso $1
+}
+
 cloneCase ()
 {
     if [ -d $2 ] ; then

From d1f36140f8c37e93e847172ff80a273d24587380 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Wed, 13 Oct 2010 12:42:42 +0100
Subject: [PATCH 9/9] Cleaning SVN garbage

---
 .../icoDyMFoam/.svn/all-wcprops               |  35 ---
 .../icoDyMFoam/.svn/dir-prop-base             |   7 -
 .../incompressible/icoDyMFoam/.svn/entries    | 201 ------------------
 .../incompressible/icoDyMFoam/.svn/format     |   1 -
 .../icoDyMFoam/.svn/prop-base/UEqn.H.svn-base |   5 -
 .../.svn/prop-base/correctPhi.H.svn-base      |   5 -
 .../.svn/prop-base/createFields.H.svn-base    |   5 -
 .../.svn/prop-base/icoDyMFoam.C.svn-base      |   5 -
 .../icoDyMFoam/.svn/text-base/UEqn.H.svn-base |  11 -
 .../.svn/text-base/correctPhi.H.svn-base      |  47 ----
 .../.svn/text-base/createFields.H.svn-base    |  90 --------
 .../.svn/text-base/icoDyMFoam.C.svn-base      | 144 -------------
 .../.svn/text-base/readControls.H.svn-base    |  14 --
 .../icoDyMFoam/Make/.svn/all-wcprops          |  17 --
 .../icoDyMFoam/Make/.svn/dir-prop-base        |   9 -
 .../icoDyMFoam/Make/.svn/entries              |  96 ---------
 .../icoDyMFoam/Make/.svn/format               |   1 -
 .../Make/.svn/text-base/files.svn-base        |   3 -
 .../Make/.svn/text-base/options.svn-base      |  13 --
 19 files changed, 709 deletions(-)
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/all-wcprops
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/dir-prop-base
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/entries
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/format
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/UEqn.H.svn-base
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/correctPhi.H.svn-base
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/createFields.H.svn-base
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/icoDyMFoam.C.svn-base
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/text-base/UEqn.H.svn-base
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/text-base/correctPhi.H.svn-base
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/text-base/createFields.H.svn-base
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/text-base/icoDyMFoam.C.svn-base
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/.svn/text-base/readControls.H.svn-base
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/Make/.svn/all-wcprops
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/Make/.svn/dir-prop-base
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/Make/.svn/entries
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/Make/.svn/format
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/files.svn-base
 delete mode 100644 applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/options.svn-base

diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/all-wcprops b/applications/solvers/incompressible/icoDyMFoam/.svn/all-wcprops
deleted file mode 100644
index bceca7b2a..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/.svn/all-wcprops
+++ /dev/null
@@ -1,35 +0,0 @@
-K 25
-svn:wc:ra_dav:version-url
-V 113
-/svnroot/openfoam-extend/!svn/ver/1591/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam
-END
-createFields.H
-K 25
-svn:wc:ra_dav:version-url
-V 127
-/svnroot/openfoam-extend/!svn/ver/915/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/createFields.H
-END
-readControls.H
-K 25
-svn:wc:ra_dav:version-url
-V 127
-/svnroot/openfoam-extend/!svn/ver/915/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/readControls.H
-END
-correctPhi.H
-K 25
-svn:wc:ra_dav:version-url
-V 125
-/svnroot/openfoam-extend/!svn/ver/981/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/correctPhi.H
-END
-icoDyMFoam.C
-K 25
-svn:wc:ra_dav:version-url
-V 126
-/svnroot/openfoam-extend/!svn/ver/1591/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C
-END
-UEqn.H
-K 25
-svn:wc:ra_dav:version-url
-V 119
-/svnroot/openfoam-extend/!svn/ver/912/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/UEqn.H
-END
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/dir-prop-base b/applications/solvers/incompressible/icoDyMFoam/.svn/dir-prop-base
deleted file mode 100644
index 3b7a2aa5a..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/.svn/dir-prop-base
+++ /dev/null
@@ -1,7 +0,0 @@
-K 10
-svn:ignore
-V 16
-lnInclude
-*.dep
-
-END
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/entries b/applications/solvers/incompressible/icoDyMFoam/.svn/entries
deleted file mode 100644
index a302fbbf2..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/.svn/entries
+++ /dev/null
@@ -1,201 +0,0 @@
-9
-
-dir
-1730
-https://openfoam-extend.svn.sourceforge.net/svnroot/openfoam-extend/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam
-https://openfoam-extend.svn.sourceforge.net/svnroot/openfoam-extend
-
-
-
-2010-01-20T10:44:39.337304Z
-1591
-hjasak
-has-props
-
-svn:special svn:externals svn:needs-lock
-
-
-
-
-
-
-
-
-
-
-
-e4e07f05-0c2f-0410-a05a-b8ba57e0c909
-
-createFields.H
-file
-
-
-
-
-2009-08-11T21:05:38.000000Z
-00e915519b705cad0056a257a3fd2a92
-2009-01-20T00:30:44.740688Z
-915
-hjasak
-has-props
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-1874
-
-readControls.H
-file
-
-
-
-
-2009-08-11T21:05:38.000000Z
-eea361849937e2d7f8dcb00e72a34c00
-2009-01-20T00:30:44.740688Z
-915
-hjasak
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-362
-
-correctPhi.H
-file
-
-
-
-
-2009-08-11T21:05:38.000000Z
-074a2d4fbd974ca63d4207e2d816104f
-2009-02-06T16:17:03.879174Z
-981
-hjasak
-has-props
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-1048
-
-icoDyMFoam.C
-file
-
-
-
-
-2009-12-13T23:44:56.000000Z
-97fc42064b2254262c779f71eabc4991
-2010-01-20T10:44:39.337304Z
-1591
-hjasak
-has-props
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-4049
-
-Make
-dir
-
-UEqn.H
-file
-
-
-
-
-2009-08-11T21:05:38.000000Z
-84ccfd4a6cf96222b86f76eaacb9f17c
-2008-06-19T13:53:26.834514Z
-684
-hjasak
-has-props
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-190
-
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/format b/applications/solvers/incompressible/icoDyMFoam/.svn/format
deleted file mode 100644
index ec635144f..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/.svn/format
+++ /dev/null
@@ -1 +0,0 @@
-9
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/UEqn.H.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/UEqn.H.svn-base
deleted file mode 100644
index 92c8ad747..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/UEqn.H.svn-base
+++ /dev/null
@@ -1,5 +0,0 @@
-K 12
-svn:keywords
-V 2
-Id
-END
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/correctPhi.H.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/correctPhi.H.svn-base
deleted file mode 100644
index 92c8ad747..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/correctPhi.H.svn-base
+++ /dev/null
@@ -1,5 +0,0 @@
-K 12
-svn:keywords
-V 2
-Id
-END
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/createFields.H.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/createFields.H.svn-base
deleted file mode 100644
index 92c8ad747..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/createFields.H.svn-base
+++ /dev/null
@@ -1,5 +0,0 @@
-K 12
-svn:keywords
-V 2
-Id
-END
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/icoDyMFoam.C.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/icoDyMFoam.C.svn-base
deleted file mode 100644
index 92c8ad747..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/.svn/prop-base/icoDyMFoam.C.svn-base
+++ /dev/null
@@ -1,5 +0,0 @@
-K 12
-svn:keywords
-V 2
-Id
-END
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/UEqn.H.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/UEqn.H.svn-base
deleted file mode 100644
index faf508468..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/UEqn.H.svn-base
+++ /dev/null
@@ -1,11 +0,0 @@
-    fvVectorMatrix UEqn
-    (
-        fvm::ddt(U)
-      + fvm::div(phi, U)
-      - fvm::laplacian(nu, U)
-    );
-
-    if (momentumPredictor)
-    {
-        solve(UEqn == -fvc::grad(p));
-    }
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/correctPhi.H.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/correctPhi.H.svn-base
deleted file mode 100644
index 69978f0e6..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/correctPhi.H.svn-base
+++ /dev/null
@@ -1,47 +0,0 @@
-{
-    wordList pcorrTypes(p.boundaryField().types());
-
-    for (label i=0; i<p.boundaryField().size(); i++)
-    {
-        if (p.boundaryField()[i].fixesValue())
-        {
-            pcorrTypes[i] = fixedValueFvPatchScalarField::typeName;
-        }
-    }
-
-    volScalarField pcorr
-    (
-        IOobject
-        (
-            "pcorr",
-            runTime.timeName(),
-            mesh,
-            IOobject::NO_READ,
-            IOobject::NO_WRITE
-        ),
-        mesh,
-        dimensionedScalar("pcorr", p.dimensions(), 0.0),
-        pcorrTypes
-    );
-
-    for(int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
-    {
-        fvScalarMatrix pcorrEqn
-        (
-            fvm::laplacian(rAU, pcorr) == fvc::div(phi)
-        );
-
-        pcorrEqn.setReference(pRefCell, pRefValue);
-        pcorrEqn.solve();
-
-        if (nonOrth == nNonOrthCorr)
-        {
-            phi -= pcorrEqn.flux();
-        }
-
-        // Fluxes are corrected to absolute velocity and further corrected
-        // later.  HJ, 6/Feb/2009
-    }
-}
-
-#include "continuityErrs.H"
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/createFields.H.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/createFields.H.svn-base
deleted file mode 100644
index 561cd7bf4..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/createFields.H.svn-base
+++ /dev/null
@@ -1,90 +0,0 @@
-    Info<< "Reading transportProperties\n" << endl;
-
-    IOdictionary transportProperties
-    (
-        IOobject
-        (
-            "transportProperties",
-            runTime.constant(),
-            mesh,
-            IOobject::MUST_READ,
-            IOobject::NO_WRITE
-        )
-    );
-
-    dimensionedScalar nu
-    (
-        transportProperties.lookup("nu")
-    );
-
-
-    Info<< "Reading field p\n" << endl;
-    volScalarField p
-    (
-        IOobject
-        (
-            "p",
-            runTime.timeName(),
-            mesh,
-            IOobject::MUST_READ,
-            IOobject::AUTO_WRITE
-        ),
-        mesh
-    );
-
-
-    Info<< "Reading field U\n" << endl;
-    volVectorField U
-    (
-        IOobject
-        (
-            "U",
-            runTime.timeName(),
-            mesh,
-            IOobject::MUST_READ,
-            IOobject::AUTO_WRITE
-        ),
-        mesh
-    );
-
-#   include "createPhi.H"
-
-
-    label pRefCell = 0;
-    scalar pRefValue = 0.0;
-    setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);
-
-
-    Info<< "Reading field rAU if present\n" << endl;
-    volScalarField rAU
-    (
-        IOobject
-        (
-            "rAU",
-            runTime.timeName(),
-            mesh,
-            IOobject::READ_IF_PRESENT,
-            IOobject::AUTO_WRITE
-        ),
-        mesh,
-        runTime.deltaT(),
-        zeroGradientFvPatchScalarField::typeName
-    );
-
-
-    {
-        const dictionary& piso = mesh.solutionDict().subDict("PISO");
-
-        bool momentumPredictor = true;
-        if (piso.found("momentumPredictor"))
-        {
-            momentumPredictor = Switch(piso.lookup("momentumPredictor"));
-        }
-
-        bool momentumPredictorSave = momentumPredictor;
-        momentumPredictor = false;
-#       include "UEqn.H"
-        momentumPredictor = momentumPredictorSave;
-        
-        rAU = 1.0/UEqn.A();
-    }
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/icoDyMFoam.C.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/icoDyMFoam.C.svn-base
deleted file mode 100644
index 1c85fdec5..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/icoDyMFoam.C.svn-base
+++ /dev/null
@@ -1,144 +0,0 @@
-/*---------------------------------------------------------------------------*\
-  =========                 |
-  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
-   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2006-2007 Zeljko Tukovic
-     \\/     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 2 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, write to the Free Software Foundation,
-    Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
-
-Application
-    icoDyMFoam
-
-Description
-    Transient solver for incompressible, laminar flow of Newtonian fluids
-    with dynamic mesh.
-
-\*---------------------------------------------------------------------------*/
-
-#include "fvCFD.H"
-#include "dynamicFvMesh.H"
-
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-int main(int argc, char *argv[])
-{
-#   include "setRootCase.H"
-#   include "createTime.H"
-#   include "createDynamicFvMesh.H"
-#   include "initContinuityErrs.H"
-#   include "initTotalVolume.H"
-#   include "createFields.H"
-
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-    Info<< "\nStarting time loop\n" << endl;
-
-    while (runTime.run())
-    {
-#       include "readControls.H"
-#       include "checkTotalVolume.H"
-#       include "CourantNo.H"
-
-#       include "setDeltaT.H"
-
-        // Make the fluxes absolute
-        fvc::makeAbsolute(phi, U);
-
-        runTime++;
-
-        Info<< "Time = " << runTime.timeName() << nl << endl;
-
-        bool meshChanged = mesh.update();
-
-#       include "volContinuity.H"
-
-        if (correctPhi && meshChanged)
-        {
-#           include "correctPhi.H"
-        }
-
-        // Make the fluxes relative to the mesh motion
-        fvc::makeRelative(phi, U);
-
-        if (meshChanged)
-        {
-#           include "CourantNo.H"
-        }
-
-#       include "UEqn.H"
-
-        // --- PISO loop
-
-        for (int corr=0; corr<nCorr; corr++)
-        {
-            rAU = 1.0/UEqn.A();
-
-            U = rAU*UEqn.H();
-            phi = (fvc::interpolate(U) & mesh.Sf());
-              //+ fvc::ddtPhiCorr(rAU, U, phi);
-
-            adjustPhi(phi, U, p);
-
-            for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
-            {
-                fvScalarMatrix pEqn
-                (
-                    fvm::laplacian(rAU, p) == fvc::div(phi)
-                );
-
-                pEqn.setReference(pRefCell, pRefValue);
-
-                if (corr == nCorr - 1 && nonOrth == nNonOrthCorr)
-                {
-                    pEqn.solve(mesh.solver(p.name() + "Final"));
-                }
-                else
-                {
-                    pEqn.solve(mesh.solver(p.name()));
-                }
-
-                if (nonOrth == nNonOrthCorr)
-                {
-                    phi -= pEqn.flux();
-                }
-            }
-
-#           include "continuityErrs.H"
-
-            // Make the fluxes relative to the mesh motion
-            fvc::makeRelative(phi, U);
-
-            U -= rAU*fvc::grad(p);
-            U.correctBoundaryConditions();
-        }
-
-        runTime.write();
-
-        Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
-            << "  ClockTime = " << runTime.elapsedClockTime() << " s"
-            << nl << endl;
-    }
-
-    Info<< "End\n" << endl;
-
-    return(0);
-}
-
-
-// ************************************************************************* //
diff --git a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/readControls.H.svn-base b/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/readControls.H.svn-base
deleted file mode 100644
index 9336616a1..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/.svn/text-base/readControls.H.svn-base
+++ /dev/null
@@ -1,14 +0,0 @@
-#   include "readTimeControls.H"
-#   include "readPISOControls.H"
-
-    bool correctPhi = false;
-    if (piso.found("correctPhi"))
-    {
-        correctPhi = Switch(piso.lookup("correctPhi"));
-    }
-
-    bool checkMeshCourantNo = false;
-    if (piso.found("checkMeshCourantNo"))
-    {
-        checkMeshCourantNo = Switch(piso.lookup("checkMeshCourantNo"));
-    }
diff --git a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/all-wcprops b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/all-wcprops
deleted file mode 100644
index d156d8090..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/all-wcprops
+++ /dev/null
@@ -1,17 +0,0 @@
-K 25
-svn:wc:ra_dav:version-url
-V 117
-/svnroot/openfoam-extend/!svn/ver/912/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/Make
-END
-files
-K 25
-svn:wc:ra_dav:version-url
-V 123
-/svnroot/openfoam-extend/!svn/ver/912/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/Make/files
-END
-options
-K 25
-svn:wc:ra_dav:version-url
-V 125
-/svnroot/openfoam-extend/!svn/ver/912/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/Make/options
-END
diff --git a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/dir-prop-base b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/dir-prop-base
deleted file mode 100644
index 91bc687f2..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/dir-prop-base
+++ /dev/null
@@ -1,9 +0,0 @@
-K 10
-svn:ignore
-V 28
-*Opt
-*Debug
-lnInclude
-*.dep
-
-END
diff --git a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/entries b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/entries
deleted file mode 100644
index 0c0be6d7a..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/entries
+++ /dev/null
@@ -1,96 +0,0 @@
-9
-
-dir
-1730
-https://openfoam-extend.svn.sourceforge.net/svnroot/openfoam-extend/trunk/Core/OpenFOAM-1.5-dev/applications/solvers/incompressible/icoDyMFoam/Make
-https://openfoam-extend.svn.sourceforge.net/svnroot/openfoam-extend
-
-
-
-2007-11-08T17:59:31.038014Z
-280
-hjasak
-has-props
-
-svn:special svn:externals svn:needs-lock
-
-
-
-
-
-
-
-
-
-
-
-e4e07f05-0c2f-0410-a05a-b8ba57e0c909
-
-files
-file
-
-
-
-
-2009-08-11T21:05:38.000000Z
-b327ba4ea7390a276aef34f189cc6ba6
-2007-06-05T22:18:47.316404Z
-3
-oforiginal
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-46
-
-options
-file
-
-
-
-
-2009-08-11T21:05:38.000000Z
-290d08439d3bdbe6dae9992a3d3d9f14
-2007-11-08T17:59:31.038014Z
-280
-hjasak
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-302
-
diff --git a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/format b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/format
deleted file mode 100644
index ec635144f..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/format
+++ /dev/null
@@ -1 +0,0 @@
-9
diff --git a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/files.svn-base b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/files.svn-base
deleted file mode 100644
index 22873bb94..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/files.svn-base
+++ /dev/null
@@ -1,3 +0,0 @@
-icoDyMFoam.C
-
-EXE = $(FOAM_APPBIN)/icoDyMFoam
diff --git a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/options.svn-base b/applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/options.svn-base
deleted file mode 100644
index 460426a5b..000000000
--- a/applications/solvers/incompressible/icoDyMFoam/Make/.svn/text-base/options.svn-base
+++ /dev/null
@@ -1,13 +0,0 @@
-EXE_INC = \
-    -I$(LIB_SRC)/dynamicFvMesh/lnInclude \
-    -I$(LIB_SRC)/dynamicMesh/lnInclude \
-    -I$(LIB_SRC)/meshTools/lnInclude \
-    -I$(LIB_SRC)/finiteVolume/lnInclude
-
-EXE_LIBS = \
-    -ldynamicFvMesh \
-    -ldynamicMesh \
-    -lengine \
-    -lmeshTools \
-    -lfiniteVolume \
-    -llduSolvers