Merge branch 'master' of /home/dm4/OpenFOAM/OpenFOAM-dev

applications/solvers/compressible/rhoCentralFoam/rhoCentralFoam.C

@@ -157,7 +157,7 @@ int main(int argc, char *argv[])
 
         Info<< "Time = " << runTime.timeName() << nl << endl;
 
-        surfaceScalarField phi("phi", aphiv_pos*rho_pos + aphiv_neg*rho_neg);
+        phi = aphiv_pos*rho_pos + aphiv_neg*rho_neg;
 
         surfaceVectorField phiUp
         (

applications/utilities/mesh/generation/snappyHexMesh/snappyHexMeshDict

@@ -212,7 +212,7 @@ snapControls
     //- Maximum relative distance for points to be attracted by surface.
     //  True distance is this factor times local maximum edge length.
     // Note: changed(corrected) w.r.t 17x! (17x used 2* tolerance)
-    tolerance 1.0;
+    tolerance 2.0;
 
     //- Number of mesh displacement relaxation iterations.
     nSolveIter 30;

applications/utilities/mesh/manipulation/setSet/setSet.C

@@ -1018,7 +1018,7 @@ int main(int argc, char *argv[])
                     is
                 );
 
-                if (!ok)
+                if (!ok && batch)
                 {
                     // Exit with error.
                     quit = true;

applications/utilities/mesh/manipulation/topoSet/topoSet.C

@@ -35,6 +35,10 @@ Description
 #include "pointSet.H"
 #include "globalMeshData.H"
 #include "timeSelector.H"
+#include "IOobjectList.H"
+#include "cellZoneSet.H"
+#include "faceZoneSet.H"
+#include "pointZoneSet.H"
 
 using namespace Foam;
 
@@ -51,6 +55,96 @@ void printMesh(const Time& runTime, const polyMesh& mesh)
 }
 
 
+template<class ZoneType>
+void removeZone
+(
+    ZoneMesh<ZoneType, polyMesh>& zones,
+    const word& setName
+)
+{
+    label zoneID = zones.findZoneID(setName);
+
+    if (zoneID != -1)
+    {
+        Info<< "Removing zone " << setName << " at index " << zoneID << endl;
+        // Shuffle to last position
+        labelList oldToNew(zones.size());
+        label newI = 0;
+        forAll(oldToNew, i)
+        {
+            if (i != zoneID)
+            {
+                oldToNew[i] = newI++;
+            }
+        }
+        oldToNew[zoneID] = newI;
+        zones.reorder(oldToNew);
+        // Remove last element
+        zones.setSize(zones.size()-1);
+        zones.clearAddressing();
+        zones.write();
+    }
+}
+
+
+// Physically remove a set
+void removeSet
+(
+    const polyMesh& mesh,
+    const word& setType,
+    const word& setName
+)
+{
+    // Remove the file
+    IOobjectList objects
+    (
+        mesh,
+        mesh.time().findInstance
+        (
+            polyMesh::meshSubDir/"sets",
+            word::null,
+            IOobject::READ_IF_PRESENT,
+            mesh.facesInstance()
+        ),
+        polyMesh::meshSubDir/"sets"
+    );
+
+    if (objects.found(setName))
+    {
+        // Remove file
+        fileName object = objects[setName]->objectPath();
+        Info<< "Removing file " << object << endl;
+        rm(object);
+    }
+
+    // See if zone
+    if (setType == cellZoneSet::typeName)
+    {
+        removeZone
+        (
+            const_cast<cellZoneMesh&>(mesh.cellZones()),
+            setName
+        );
+    }
+    else if (setType == faceZoneSet::typeName)
+    {
+        removeZone
+        (
+            const_cast<faceZoneMesh&>(mesh.faceZones()),
+            setName
+        );
+    }
+    else if (setType == pointZoneSet::typeName)
+    {
+        removeZone
+        (
+            const_cast<pointZoneMesh&>(mesh.pointZones()),
+            setName
+        );
+    }
+}
+
+
 polyMesh::readUpdateState meshReadUpdate(polyMesh& mesh)
 {
     polyMesh::readUpdateState stat = mesh.readUpdate();
@@ -284,6 +378,12 @@ int main(int argc, char *argv[])
                     currentSet().write();
                 break;
 
+                case topoSetSource::REMOVE:
+                    Info<< "    Removing set" << endl;
+                    removeSet(mesh, setType, setName);
+                break;
+
+
                 default:
                     WarningIn(args.executable())
                         << "Unhandled action " << action << endl;

applications/utilities/mesh/manipulation/topoSet/topoSetDict

@@ -22,8 +22,15 @@ FoamFile
 //     type    cellSet;
 //
 //     // action to perform on set. Two types:
-//     // - require no source : clear/invert
+//     // - require no source : clear/invert/remove
+//     //       clear  : clears set or zone
+//     //       invert : select all currently non-selected elements
+//     //       remove : removes set or zone
 //     // - require source    : new/add/delete/subset
+//     //       new    : create new set or zone from source
+//     //       add    : add source to contents
+//     //       delete : deletes source from contents
+//     //       subset : keeps elements both in contents and source
 //     action  new;
 //
 // The source entry varies according to the type of set:
@@ -329,6 +336,7 @@ FoamFile
 //        faceSet f0;       // name of faceSet
 //        cellSet c0;       // name of cellSet of slave side
 //    }
+//
 
 actions
 (