From 8d0522633e087ed68d94b2da95191894114688ad Mon Sep 17 00:00:00 2001
From: Mark Olesen <Mark.Olesen@Germany>
Date: Wed, 9 Dec 2009 09:03:49 +0100
Subject: [PATCH 01/16] argList - surround option parameters with <>

- increase usageMin to 20 to compensate for extra <>

- look for _source.html instead of -source.html to match newer doxygen
---
 src/OpenFOAM/global/argList/argList.C | 13 +++++++------
 1 file changed, 7 insertions(+), 6 deletions(-)

diff --git a/src/OpenFOAM/global/argList/argList.C b/src/OpenFOAM/global/argList/argList.C
index 0e3d05b72b..cd2de41fa9 100644
--- a/src/OpenFOAM/global/argList/argList.C
+++ b/src/OpenFOAM/global/argList/argList.C
@@ -43,7 +43,7 @@ Foam::SLList<Foam::string>    Foam::argList::validArgs;
 Foam::HashTable<Foam::string> Foam::argList::validOptions;
 Foam::HashTable<Foam::string> Foam::argList::validParOptions;
 Foam::HashTable<Foam::string> Foam::argList::optionUsage;
-Foam::string::size_type Foam::argList::usageMin = 16;
+Foam::string::size_type Foam::argList::usageMin = 20;
 Foam::string::size_type Foam::argList::usageMax = 80;
 
 
@@ -51,7 +51,7 @@ Foam::argList::initValidTables::initValidTables()
 {
     argList::addOption
     (
-        "case", "DIR",
+        "case", "dir",
         "specify alternate case directory, default is the cwd"
     );
     argList::addBoolOption("parallel", "run in parallel");
@@ -769,8 +769,9 @@ void Foam::argList::printUsage() const
 
         if (iter().size())
         {
-            len += iter().size() + 1; // include space between option and param
-            Info<< ' ' << iter().c_str();
+            // length includes space and between option/param and '<>'
+            len += iter().size() + 3;
+            Info<< " <" << iter().c_str() << '>';
         }
 
         HashTable<string>::const_iterator usageIter =
@@ -827,12 +828,12 @@ void Foam::argList::displayDoc(bool source) const
     List<fileName> docDirs(docDict.lookup("doxyDocDirs"));
     List<fileName> docExts(docDict.lookup("doxySourceFileExts"));
 
-    // for source code: change foo_8C.html to foo_8C-source.html
+    // for source code: change foo_8C.html to foo_8C_source.html
     if (source)
     {
         forAll(docExts, extI)
         {
-            docExts[extI].replace(".", "-source.");
+            docExts[extI].replace(".", "_source.");
         }
     }
 

From dd8f4f657a1c11f1007cce12b0b614c1bb52e3f9 Mon Sep 17 00:00:00 2001
From: Mark Olesen <Mark.Olesen@Germany>
Date: Wed, 9 Dec 2009 10:19:53 +0100
Subject: [PATCH 02/16] adjust usage for a few applications to account for new
 style

---
 .../mesh/generation/blockMesh/blockMeshApp.C  |  4 +-
 .../transformPoints/transformPoints.C         | 37 +++++++++++++---
 .../foamToEnsight/foamToEnsight.C             |  6 +--
 .../foamToEnsightParts/foamToEnsightParts.C   |  2 +-
 .../dataConversion/foamToVTK/foamToVTK.C      | 42 +++++++++++++++----
 5 files changed, 72 insertions(+), 19 deletions(-)

diff --git a/applications/utilities/mesh/generation/blockMesh/blockMeshApp.C b/applications/utilities/mesh/generation/blockMesh/blockMeshApp.C
index 93a6536506..9140c090ec 100644
--- a/applications/utilities/mesh/generation/blockMesh/blockMeshApp.C
+++ b/applications/utilities/mesh/generation/blockMesh/blockMeshApp.C
@@ -42,7 +42,7 @@ Usage
     @param -region \<name\> \n
     Specify an alternative mesh region.
 
-    @param -dict \<dictionary\> \n
+    @param -dict \<filename\> \n
     Specify an alternative dictionary for the block mesh description.
 
 \*---------------------------------------------------------------------------*/
@@ -80,7 +80,7 @@ int main(int argc, char *argv[])
     argList::addOption
     (
         "dict",
-        "NAME",
+        "file",
         "specify an alternative dictionary for the blockMesh description"
     );
 
diff --git a/applications/utilities/mesh/manipulation/transformPoints/transformPoints.C b/applications/utilities/mesh/manipulation/transformPoints/transformPoints.C
index 9fdc9681c6..2ea94212ca 100644
--- a/applications/utilities/mesh/manipulation/transformPoints/transformPoints.C
+++ b/applications/utilities/mesh/manipulation/transformPoints/transformPoints.C
@@ -142,12 +142,39 @@ void rotateFields(const Time& runTime, const tensor& T)
 
 int main(int argc, char *argv[])
 {
-    argList::addOption("translate", "vector");
-    argList::addOption("rotate", "(vector vector)");
-    argList::addOption("rollPitchYaw", "(roll pitch yaw)");
-    argList::addOption("yawPitchRoll", "(yaw pitch roll)");
+    argList::addOption
+    (
+        "translate",
+        "vector",
+        "translate by the specified <vector> - eg, '(1 0 0)'"
+    );
+    argList::addOption
+    (
+        "rotate",
+        "(vectorA vectorB)",
+        "transform in terms of a rotation between <vectorA> and <vectorB> "
+        "- eg, '( (1 0 0) (0 0 1) )'"
+    );
+    argList::addOption
+    (
+        "rollPitchYaw",
+        "vector",
+        "transform in terms of '( roll pitch yaw )' in degrees"
+    );
+    argList::addOption
+    (
+        "yawPitchRoll",
+        "vector",
+        "transform in terms of '( yaw pitch roll )' in degrees"
+    );
     argList::addBoolOption("rotateFields");
-    argList::addOption("scale", "vector");
+    argList::addOption
+    (
+        "scale",
+        "vector",
+        "scale by the specified amount - eg, '(0.001 0.001 0.001)' for a "
+        "uniform [mm] to [m] scaling"
+    );
 
 #   include "setRootCase.H"
 #   include "createTime.H"
diff --git a/applications/utilities/postProcessing/dataConversion/foamToEnsight/foamToEnsight.C b/applications/utilities/postProcessing/dataConversion/foamToEnsight/foamToEnsight.C
index 94cf1dfecf..c0ea000959 100644
--- a/applications/utilities/postProcessing/dataConversion/foamToEnsight/foamToEnsight.C
+++ b/applications/utilities/postProcessing/dataConversion/foamToEnsight/foamToEnsight.C
@@ -101,13 +101,13 @@ int main(int argc, char *argv[])
     argList::addBoolOption
     (
         "noPatches",
-        "Suppress writing any patches"
+        "suppress writing any patches"
     );
     argList::addOption
     (
         "patches",
-        "patchList",
-        "Specify particular patches to write. "
+        "wordList",
+        "specify particular patches to write - eg '(inlet outlet)'. "
         "An empty list suppresses writing the internalMesh."
     );
 
diff --git a/applications/utilities/postProcessing/dataConversion/foamToEnsightParts/foamToEnsightParts.C b/applications/utilities/postProcessing/dataConversion/foamToEnsightParts/foamToEnsightParts.C
index 862816a946..22ba937868 100644
--- a/applications/utilities/postProcessing/dataConversion/foamToEnsightParts/foamToEnsightParts.C
+++ b/applications/utilities/postProcessing/dataConversion/foamToEnsightParts/foamToEnsightParts.C
@@ -94,7 +94,7 @@ int main(int argc, char *argv[])
     argList::addBoolOption
     (
         "noMesh",
-        "Suppress writing the geometry. "
+        "suppress writing the geometry. "
         "Can be useful for converting partial results for a static geometry"
     );
 
diff --git a/applications/utilities/postProcessing/dataConversion/foamToVTK/foamToVTK.C b/applications/utilities/postProcessing/dataConversion/foamToVTK/foamToVTK.C
index 9a82419e20..73dc60b762 100644
--- a/applications/utilities/postProcessing/dataConversion/foamToVTK/foamToVTK.C
+++ b/applications/utilities/postProcessing/dataConversion/foamToVTK/foamToVTK.C
@@ -241,20 +241,46 @@ int main(int argc, char *argv[])
 
 #   include "addRegionOption.H"
 
-    argList::addOption("fields", "fields");
-    argList::addOption("cellSet", "cellSet name");
-    argList::addOption("faceSet", "faceSet name");
-    argList::addOption("pointSet", "pointSet name");
-    argList::addBoolOption("ascii");
+    argList::addOption
+    (
+        "fields", "wordList",
+        "only convert the specified fields - eg '(p T U)'"
+    );
+    argList::addOption
+    (
+        "cellSet",
+        "name",
+        "convert a mesh subset corresponding to the specified cellSet"
+    );
+    argList::addOption("faceSet", "name");
+    argList::addOption("pointSet", "name");
+    argList::addBoolOption
+    (
+        "ascii",
+        "write in ASCII format instead of binary"
+    );
     argList::addBoolOption("surfaceFields");
     argList::addBoolOption("nearCellValue");
     argList::addBoolOption("noInternal");
     argList::addBoolOption("noPointValues");
     argList::addBoolOption("allPatches");
-    argList::addOption("excludePatches","patches to exclude");
+    argList::addOption
+    (
+        "excludePatches",
+        "wordReList",
+        "a list of patches to exclude - eg '( inlet \".*Wall\" )' "
+    );
     argList::addBoolOption("noFaceZones");
-    argList::addBoolOption("noLinks");
-    argList::addBoolOption("useTimeName");
+    argList::addBoolOption
+    (
+        "noLinks",
+        "don't link processor VTK files - parallel only"
+    );
+    argList::addBoolOption
+    (
+        "useTimeName",
+        "use the time name instead of the time index when naming the files"
+    );
 
 #   include "setRootCase.H"
 #   include "createTime.H"

From 30eac79d5551483cde101a4bdd3b2192566f9367 Mon Sep 17 00:00:00 2001
From: Mark Olesen <Mark.Olesen@Germany>
Date: Wed, 9 Dec 2009 10:58:32 +0100
Subject: [PATCH 03/16] argList gets addNote() static method

- output any notes in the usage
---
 .../surfaceTransformPoints.C                  | 47 +++++++++++++++----
 src/OpenFOAM/global/argList/argList.C         | 28 ++++++++++-
 src/OpenFOAM/global/argList/argList.H         |  9 +++-
 3 files changed, 74 insertions(+), 10 deletions(-)

diff --git a/applications/utilities/surface/surfaceTransformPoints/surfaceTransformPoints.C b/applications/utilities/surface/surfaceTransformPoints/surfaceTransformPoints.C
index 85e3cc1172..8bb87d3ee9 100644
--- a/applications/utilities/surface/surfaceTransformPoints/surfaceTransformPoints.C
+++ b/applications/utilities/surface/surfaceTransformPoints/surfaceTransformPoints.C
@@ -23,8 +23,8 @@ License
     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
 
 Description
-    Transform (scale/rotate) a surface. Like transformPoints but then for
-    surfaces.
+    Transform (scale/rotate) a surface.
+    Like transformPoints but for surfaces.
 
     The rollPitchYaw option takes three angles (degrees):
     - roll (rotation about x) followed by
@@ -54,16 +54,47 @@ using namespace Foam::constant::mathematical;
 
 int main(int argc, char *argv[])
 {
+    argList::addNote
+    (
+        "Transform (scale/rotate) a surface. "
+        "Like transformPoints but for surfaces."
+    );
     argList::noParallel();
     argList::validArgs.clear();
-
     argList::validArgs.append("surface file");
     argList::validArgs.append("output surface file");
-    argList::addOption("translate", "vector");
-    argList::addOption("rotate", "(vector vector)");
-    argList::addOption("scale", "vector");
-    argList::addOption("rollPitchYaw", "(roll pitch yaw)");
-    argList::addOption("yawPitchRoll", "(yaw pitch roll)");
+    argList::addOption
+    (
+        "translate",
+        "vector",
+        "translate by the specified <vector> - eg, '(1 0 0)'"
+    );
+    argList::addOption
+    (
+        "rotate",
+        "(vectorA vectorB)",
+        "transform in terms of a rotation between <vectorA> and <vectorB> "
+        "- eg, '( (1 0 0) (0 0 1) )'"
+    );
+    argList::addOption
+    (
+        "scale",
+        "vector",
+        "scale by the specified amount - eg, '(0.001 0.001 0.001)' for a "
+        "uniform [mm] to [m] scaling"
+    );
+    argList::addOption
+    (
+        "rollPitchYaw",
+        "vector",
+        "transform in terms of '( roll pitch yaw )' in degrees"
+    );
+    argList::addOption
+    (
+        "yawPitchRoll",
+        "vector",
+        "transform in terms of '( yaw pitch roll )' in degrees"
+    );
     argList args(argc, argv);
 
     fileName surfFileName(args.additionalArgs()[0]);
diff --git a/src/OpenFOAM/global/argList/argList.C b/src/OpenFOAM/global/argList/argList.C
index cd2de41fa9..f574fde3e6 100644
--- a/src/OpenFOAM/global/argList/argList.C
+++ b/src/OpenFOAM/global/argList/argList.C
@@ -43,6 +43,7 @@ Foam::SLList<Foam::string>    Foam::argList::validArgs;
 Foam::HashTable<Foam::string> Foam::argList::validOptions;
 Foam::HashTable<Foam::string> Foam::argList::validParOptions;
 Foam::HashTable<Foam::string> Foam::argList::optionUsage;
+Foam::SLList<Foam::string>    Foam::argList::notes;
 Foam::string::size_type Foam::argList::usageMin = 20;
 Foam::string::size_type Foam::argList::usageMax = 80;
 
@@ -108,6 +109,15 @@ void Foam::argList::addUsage
 }
 
 
+void Foam::argList::addNote(const string& note)
+{
+    if (!note.empty())
+    {
+        notes.append(note);
+    }
+}
+
+
 void Foam::argList::removeOption(const word& opt)
 {
     validOptions.erase(opt);
@@ -765,7 +775,7 @@ void Foam::argList::printUsage() const
 
         HashTable<string>::const_iterator iter = validOptions.find(optionName);
         Info<< "  -" << optionName;
-        label len = optionName.size() + 3;  // include leading "  -"
+        label len = optionName.size() + 3;  // length includes leading '  -'
 
         if (iter().size())
         {
@@ -815,6 +825,22 @@ void Foam::argList::printUsage() const
         "print the usage"
     );
 
+
+    // output notes directly - no automatic text wrapping
+    if (!notes.empty())
+    {
+        Info<< nl;
+        for
+        (
+            SLList<string>::const_iterator iter = notes.begin();
+            iter != notes.end();
+            ++iter
+        )
+        {
+            Info<< iter().c_str() << nl;
+        }
+    }
+
     Info<< nl
         <<"Using OpenFOAM-" << Foam::FOAMversion
         <<" (build: " << Foam::FOAMbuild << ") - see www.OpenFOAM.org"
diff --git a/src/OpenFOAM/global/argList/argList.H b/src/OpenFOAM/global/argList/argList.H
index b1efcb0669..da75c4d066 100644
--- a/src/OpenFOAM/global/argList/argList.H
+++ b/src/OpenFOAM/global/argList/argList.H
@@ -153,7 +153,10 @@ public:
         //- Short usage information for validOptions
         static HashTable<string> optionUsage;
 
-        //- Min offset for displaying usage (default: 16)
+        //- Additional notes for usage
+        static SLList<string> notes;
+
+        //- Min offset for displaying usage (default: 20)
         static string::size_type usageMin;
 
         //- Max screen width for displaying usage (default: 80)
@@ -287,6 +290,10 @@ public:
                 const string& usage
             );
 
+            //- Add extra notes for the usage information
+            //  This string is used "as-is" without additional formatting
+            static void addNote(const string&);
+
             //- Remove option from validOptions and from optionUsage
             static void removeOption(const word& opt);
 

From 0fc464254aa0559ba190d0dc5bf767b224ca0a15 Mon Sep 17 00:00:00 2001
From: Mark Olesen <Mark.Olesen@Germany>
Date: Wed, 9 Dec 2009 13:50:31 +0100
Subject: [PATCH 04/16] Fix minor bug in surfaceFeatureExtract.C

---
 .../surface/surfaceFeatureExtract/surfaceFeatureExtract.C       | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/applications/utilities/surface/surfaceFeatureExtract/surfaceFeatureExtract.C b/applications/utilities/surface/surfaceFeatureExtract/surfaceFeatureExtract.C
index 307edd8568..f38bf13cc5 100644
--- a/applications/utilities/surface/surfaceFeatureExtract/surfaceFeatureExtract.C
+++ b/applications/utilities/surface/surfaceFeatureExtract/surfaceFeatureExtract.C
@@ -199,7 +199,7 @@ int main(int argc, char *argv[])
     }
 
     // Trim away small groups of features
-    if (minLen > 0 || minLen > 0)
+    if (minElem > 0 || minLen > 0)
     {
         set.trimFeatures(minLen, minElem);
         Pout<< endl << "Removed small features" << endl;

From 17846a6bda7d6573a7fe4eb2a18cafffa99b8f6b Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Thu, 10 Dec 2009 15:27:21 +0000
Subject: [PATCH 05/16] Removed unused include file

---
 src/dynamicMesh/fvMeshDistribute/fvMeshDistribute.H | 1 -
 1 file changed, 1 deletion(-)

diff --git a/src/dynamicMesh/fvMeshDistribute/fvMeshDistribute.H b/src/dynamicMesh/fvMeshDistribute/fvMeshDistribute.H
index 0be663b733..595a53de9c 100644
--- a/src/dynamicMesh/fvMeshDistribute/fvMeshDistribute.H
+++ b/src/dynamicMesh/fvMeshDistribute/fvMeshDistribute.H
@@ -54,7 +54,6 @@ SourceFiles
 #define fvMeshDistribute_H
 
 #include "Field.H"
-#include "uLabel.H"
 #include "fvMeshSubset.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

From 35eeedeaa511514b5933ba1d5fa8eb8735c133c5 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Thu, 10 Dec 2009 15:31:04 +0000
Subject: [PATCH 06/16] Added allowFreeStandingZoneFaces flag

This determines when faceZones and cellZones are being meshed any zone faces
can only be on the boundary of the cellZone or also 'free standing'. The
latter option can be used to e.g. create free-standing baffles
(that do not separate two regions)
---
 .../autoHexMeshDriver/autoRefineDriver.C            |  6 +++++-
 .../refinementParameters/refinementParameters.C     |  2 ++
 .../refinementParameters/refinementParameters.H     | 12 ++++++++++++
 .../autoHexMesh/meshRefinement/meshRefinement.H     |  6 +++++-
 .../meshRefinement/meshRefinementBaffles.C          | 13 ++++++++-----
 5 files changed, 32 insertions(+), 7 deletions(-)

diff --git a/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoRefineDriver.C b/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoRefineDriver.C
index a1346b4e8e..cff924d864 100644
--- a/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoRefineDriver.C
+++ b/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoRefineDriver.C
@@ -560,7 +560,11 @@ void Foam::autoRefineDriver::zonify
             const_cast<Time&>(mesh.time())++;
         }
 
-        meshRefiner_.zonify(refineParams.keepPoints()[0]);
+        meshRefiner_.zonify
+        (
+            refineParams.keepPoints()[0],
+            refineParams.allowFreeStandingZoneFaces()
+        );
 
         if (debug)
         {
diff --git a/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/refinementParameters/refinementParameters.C b/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/refinementParameters/refinementParameters.C
index a181662b27..7a8a16cfb7 100644
--- a/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/refinementParameters/refinementParameters.C
+++ b/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/refinementParameters/refinementParameters.C
@@ -44,6 +44,7 @@ Foam::refinementParameters::refinementParameters
     curvature_(readScalar(dict.lookup("curvature"))),
     nBufferLayers_(readLabel(dict.lookup("nBufferLayers"))),
     keepPoints_(dict.lookup("keepPoints")),
+    allowFreeStandingZoneFaces_(dict.lookup("allowFreeStandingZoneFaces")),
     maxLoadUnbalance_(dict.lookupOrDefault<scalar>("maxLoadUnbalance",0))
 {}
 
@@ -55,6 +56,7 @@ Foam::refinementParameters::refinementParameters(const dictionary& dict)
     minRefineCells_(readLabel(dict.lookup("minRefinementCells"))),
     nBufferLayers_(readLabel(dict.lookup("nCellsBetweenLevels"))),
     keepPoints_(pointField(1, dict.lookup("locationInMesh"))),
+    allowFreeStandingZoneFaces_(dict.lookup("allowFreeStandingZoneFaces")),
     maxLoadUnbalance_(dict.lookupOrDefault<scalar>("maxLoadUnbalance",0))
 {
     scalar featAngle(readScalar(dict.lookup("resolveFeatureAngle")));
diff --git a/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/refinementParameters/refinementParameters.H b/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/refinementParameters/refinementParameters.H
index ae09e3fc0e..32d6777276 100644
--- a/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/refinementParameters/refinementParameters.H
+++ b/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/refinementParameters/refinementParameters.H
@@ -38,6 +38,7 @@ SourceFiles
 
 #include "dictionary.H"
 #include "pointField.H"
+#include "Switch.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -73,6 +74,10 @@ class refinementParameters
         //- Areas to keep
         const pointField keepPoints_;
 
+        //- FaceZone faces allowed which have owner and neighbour in same
+        //  cellZone?
+        Switch allowFreeStandingZoneFaces_;
+
         //- Allowed load unbalance
         scalar maxLoadUnbalance_;
 
@@ -137,6 +142,13 @@ public:
                 return keepPoints_;
             }
 
+            //- Are zone faces allowed only inbetween different cell zones
+            //  or also just free standing?
+            bool allowFreeStandingZoneFaces() const
+            {
+                return allowFreeStandingZoneFaces_;
+            }
+
             //- Allowed load unbalance
             scalar maxLoadUnbalance() const
             {
diff --git a/src/mesh/autoMesh/autoHexMesh/meshRefinement/meshRefinement.H b/src/mesh/autoMesh/autoHexMesh/meshRefinement/meshRefinement.H
index 5b603ddbac..6347570cac 100644
--- a/src/mesh/autoMesh/autoHexMesh/meshRefinement/meshRefinement.H
+++ b/src/mesh/autoMesh/autoHexMesh/meshRefinement/meshRefinement.H
@@ -704,7 +704,11 @@ public:
             //- Put faces/cells into zones according to surface specification.
             //  Returns null if no zone surfaces present. Region containing
             //  the keepPoint will not be put into a cellZone.
-            autoPtr<mapPolyMesh> zonify(const point& keepPoint);
+            autoPtr<mapPolyMesh> zonify
+            (
+                const point& keepPoint,
+                const bool allowFreeStandingZoneFaces
+            );
 
 
         // Other topo changes
diff --git a/src/mesh/autoMesh/autoHexMesh/meshRefinement/meshRefinementBaffles.C b/src/mesh/autoMesh/autoHexMesh/meshRefinement/meshRefinementBaffles.C
index 307cb743f8..f7adea091d 100644
--- a/src/mesh/autoMesh/autoHexMesh/meshRefinement/meshRefinementBaffles.C
+++ b/src/mesh/autoMesh/autoHexMesh/meshRefinement/meshRefinementBaffles.C
@@ -2032,7 +2032,8 @@ Foam::autoPtr<Foam::mapPolyMesh> Foam::meshRefinement::dupNonManifoldPoints()
 // Zoning
 Foam::autoPtr<Foam::mapPolyMesh> Foam::meshRefinement::zonify
 (
-    const point& keepPoint
+    const point& keepPoint,
+    const bool allowFreeStandingZoneFaces
 )
 {
     const wordList& cellZoneNames = surfaces_.cellZoneNames();
@@ -2344,9 +2345,11 @@ Foam::autoPtr<Foam::mapPolyMesh> Foam::meshRefinement::zonify
     }
 
 
-    //// Make sure namedSurfaceIndex is unset inbetween same cell cell zones.
-    //makeConsistentFaceIndex(cellToZone, namedSurfaceIndex);
-
+    // Make sure namedSurfaceIndex is unset inbetween same cell cell zones.
+    if (!allowFreeStandingZoneFaces)
+    {
+        makeConsistentFaceIndex(cellToZone, namedSurfaceIndex);
+    }
 
     // Topochange container
     polyTopoChange meshMod(mesh_);
@@ -2440,7 +2443,7 @@ Foam::autoPtr<Foam::mapPolyMesh> Foam::meshRefinement::zonify
                 }
                 else if (ownZone == neiZone)
                 {
-                    // Can only happen for coupled boundaries. Keep master
+                    // Free-standing zone face or coupled boundary. Keep master
                     // face unflipped.
                     flip = !isMasterFace[faceI];
                 }

From 9d817bfbdfb3e512627a66cf7497364f3ca65078 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Fri, 11 Dec 2009 11:55:07 +0000
Subject: [PATCH 07/16] Added -region option

---
 .../redistributeMeshPar/redistributeMeshPar.C | 108 ++++++++++--------
 1 file changed, 63 insertions(+), 45 deletions(-)

diff --git a/applications/utilities/parallelProcessing/redistributeMeshPar/redistributeMeshPar.C b/applications/utilities/parallelProcessing/redistributeMeshPar/redistributeMeshPar.C
index dbcb9e1643..b45843e31c 100644
--- a/applications/utilities/parallelProcessing/redistributeMeshPar/redistributeMeshPar.C
+++ b/applications/utilities/parallelProcessing/redistributeMeshPar/redistributeMeshPar.C
@@ -34,13 +34,19 @@ Description
 
     Can also work like decomposePar:
     @verbatim
+        # Create empty processor directories (have to exist for argList)
         mkdir processor0
+                ..
+        mkdir processorN
+
+        # Copy undecomposed polyMesh
         cp -r constant processor0
+
+        # Distribute
         mpirun -np ddd redistributeMeshPar -parallel
     @endverbatim
 \*---------------------------------------------------------------------------*/
 
-#include "Field.H"
 #include "fvMesh.H"
 #include "decompositionMethod.H"
 #include "PstreamReduceOps.H"
@@ -62,6 +68,7 @@ static const scalar defaultMergeTol = 1E-6;
 autoPtr<fvMesh> createMesh
 (
     const Time& runTime,
+    const word& regionName,
     const fileName& instDir,
     const bool haveMesh
 )
@@ -69,43 +76,33 @@ autoPtr<fvMesh> createMesh
     Pout<< "Create mesh for time = "
         << runTime.timeName() << nl << endl;
 
-    // Create dummy mesh. Only used on procs that don't have mesh.
-    // Note constructed on all processors since does parallel comms.
-    fvMesh dummyMesh
+    IOobject io
     (
-        IOobject
-        (
-            fvMesh::defaultRegion,
-            instDir,
-            runTime,
-            IOobject::MUST_READ
-        ),
-        xferCopy(pointField()),
-        xferCopy(faceList()),
-        xferCopy(labelList()),
-        xferCopy(labelList())
+        regionName,
+        instDir,
+        runTime,
+        IOobject::MUST_READ
     );
 
     if (!haveMesh)
     {
-        Pout<< "Writing dummy mesh to " << runTime.path()/instDir << endl;
+        // Create dummy mesh. Only used on procs that don't have mesh.
+        fvMesh dummyMesh
+        (
+            io,
+            xferCopy(pointField()),
+            xferCopy(faceList()),
+            xferCopy(labelList()),
+            xferCopy(labelList()),
+            false
+        );
+        Pout<< "Writing dummy mesh to " << dummyMesh.polyMesh::objectPath()
+            << endl;
         dummyMesh.write();
     }
 
-    Pout<< "Reading mesh from " << runTime.path()/instDir << endl;
-    autoPtr<fvMesh> meshPtr
-    (
-        new fvMesh
-        (
-            IOobject
-            (
-                fvMesh::defaultRegion,
-                instDir,
-                runTime,
-                IOobject::MUST_READ
-            )
-        )
-    );
+    Pout<< "Reading mesh from " << io.objectPath() << endl;
+    autoPtr<fvMesh> meshPtr(new fvMesh(io));
     fvMesh& mesh = meshPtr();
 
 
@@ -229,8 +226,9 @@ autoPtr<fvMesh> createMesh
     if (!haveMesh)
     {
         // We created a dummy mesh file above. Delete it.
-        Pout<< "Removing dummy mesh in " << runTime.path()/instDir << endl;
-        rmDir(runTime.path()/instDir/polyMesh::meshSubDir);
+        Pout<< "Removing dummy mesh " << io.objectPath()
+            << endl;
+        rmDir(io.objectPath());
     }
 
     // Force recreation of globalMeshData.
@@ -285,7 +283,6 @@ scalar getMergeDistance
 void printMeshData(Ostream& os, const polyMesh& mesh)
 {
     os  << "Number of points:           " << mesh.points().size() << nl
-        << "          edges:            " << mesh.edges().size() << nl
         << "          faces:            " << mesh.faces().size() << nl
         << "          internal faces:   " << mesh.faceNeighbour().size() << nl
         << "          cells:            " << mesh.cells().size() << nl
@@ -506,33 +503,53 @@ void compareFields
 
 int main(int argc, char *argv[])
 {
+#   include "addRegionOption.H"
     argList::addOption("mergeTol", "relative merge distance");
+    // Create argList. This will check for non-existing processor dirs.
 #   include "setRootCase.H"
 
-    // Create processor directory if non-existing
-    if (!Pstream::master() && !isDir(args.path()))
-    {
-        Pout<< "Creating case directory " << args.path() << endl;
-        mkDir(args.path());
-    }
+    //- Not useful anymore. See above.
+    //// Create processor directory if non-existing
+    //if (!Pstream::master() && !isDir(args.path()))
+    //{
+    //    Pout<< "Creating case directory " << args.path() << endl;
+    //    mkDir(args.path());
+    //}
 
 #   include "createTime.H"
 
+    word regionName = polyMesh::defaultRegion;
+    fileName meshSubDir;
+
+    if (args.optionReadIfPresent("region", regionName))
+    {
+        meshSubDir = regionName/polyMesh::meshSubDir;
+    }
+    else
+    {
+        meshSubDir = polyMesh::meshSubDir;
+    }
+    Info<< "Using mesh subdirectory " << meshSubDir << nl << endl;
+
+
     // Get time instance directory. Since not all processors have meshes
     // just use the master one everywhere.
 
     fileName masterInstDir;
     if (Pstream::master())
     {
-        masterInstDir = runTime.findInstance(polyMesh::meshSubDir, "points");
+        masterInstDir = runTime.findInstance(meshSubDir, "points");
     }
     Pstream::scatter(masterInstDir);
 
     // Check who has a mesh
-    const fileName meshDir = runTime.path()/masterInstDir/polyMesh::meshSubDir;
+    const fileName meshPath = runTime.path()/masterInstDir/meshSubDir;
+
+    Info<< "Found points in " << meshPath << nl << endl;
+
 
     boolList haveMesh(Pstream::nProcs(), false);
-    haveMesh[Pstream::myProcNo()] = isDir(meshDir);
+    haveMesh[Pstream::myProcNo()] = isDir(meshPath);
     Pstream::gatherList(haveMesh);
     Pstream::scatterList(haveMesh);
     Info<< "Per processor mesh availability : " << haveMesh << endl;
@@ -542,6 +559,7 @@ int main(int argc, char *argv[])
     autoPtr<fvMesh> meshPtr = createMesh
     (
         runTime,
+        regionName,
         masterInstDir,
         haveMesh[Pstream::myProcNo()]
     );
@@ -799,7 +817,7 @@ int main(int argc, char *argv[])
         << nl
         << "the processor directories with 0 sized meshes in them." << nl
         << "Below is a sample set of commands to do this."
-        << " Take care when issueing these" << nl
+        << " Take care when issuing these" << nl
         << "commands." << nl << endl;
 
     forAll(nFaces, procI)
@@ -812,8 +830,8 @@ int main(int argc, char *argv[])
         }
         else
         {
-            fileName timeDir = procDir/runTime.timeName()/polyMesh::meshSubDir;
-            fileName constDir = procDir/runTime.constant()/polyMesh::meshSubDir;
+            fileName timeDir = procDir/runTime.timeName()/meshSubDir;
+            fileName constDir = procDir/runTime.constant()/meshSubDir;
 
             Info<< "    rm -r " << constDir.c_str() << nl
                 << "    mv " << timeDir.c_str()

From 88e407043e88137c1fbdcb8dd6d4d3f4046a6930 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Fri, 11 Dec 2009 11:55:36 +0000
Subject: [PATCH 08/16] Fixed typo in printing

---
 src/dynamicMesh/fvMeshDistribute/fvMeshDistribute.C | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/dynamicMesh/fvMeshDistribute/fvMeshDistribute.C b/src/dynamicMesh/fvMeshDistribute/fvMeshDistribute.C
index 7dc416b98f..848ec35d21 100644
--- a/src/dynamicMesh/fvMeshDistribute/fvMeshDistribute.C
+++ b/src/dynamicMesh/fvMeshDistribute/fvMeshDistribute.C
@@ -1461,7 +1461,7 @@ Foam::autoPtr<Foam::mapDistributePolyMesh> Foam::fvMeshDistribute::distribute
     getNeighbourData(distribution, sourceFace, sourceProc, sourceNewProc);
 
 
-    // Remove meshPhi. Since this would otherwise dissappear anyway
+    // Remove meshPhi. Since this would otherwise disappear anyway
     // during topo changes and we have to guarantee that all the fields
     // can be sent.
     mesh_.clearOut();

From 834e493812b044e4a500f98184fc98da9ab990cd Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Fri, 11 Dec 2009 11:56:04 +0000
Subject: [PATCH 09/16] Fixed missing arguments in error message

---
 .../polyTopoChange/polyTopoChange/addObject/polyAddFace.H    | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/src/dynamicMesh/polyTopoChange/polyTopoChange/addObject/polyAddFace.H b/src/dynamicMesh/polyTopoChange/polyTopoChange/addObject/polyAddFace.H
index 944e645fe0..8ae310fd36 100644
--- a/src/dynamicMesh/polyTopoChange/polyTopoChange/addObject/polyAddFace.H
+++ b/src/dynamicMesh/polyTopoChange/polyTopoChange/addObject/polyAddFace.H
@@ -264,6 +264,10 @@ public:
                     "    const face& f,\n"
                     "    const label owner,"
                     "    const label neighbour,\n"
+                    "    const label masterPointID,\n"
+                    "    const label masterEdgeID,\n"
+                    "    const label masterFaceID,\n"
+                    "    const bool flipFaceFlux,\n"
                     "    const label patchID,\n"
                     "    const label zoneID"
                     ")"
@@ -292,6 +296,7 @@ public:
                     "    const label masterPointID,\n"
                     "    const label masterEdgeID,\n"
                     "    const label masterFaceID,\n"
+                    "    const bool flipFaceFlux,\n"
                     "    const label patchID,\n"
                     "    const label zoneID,\n"
                     "    const bool zoneFlip\n"

From ce7274639d2761a8b04edb59ec3d4d6b6e096a12 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Fri, 11 Dec 2009 11:58:12 +0000
Subject: [PATCH 10/16] Added missing face-flipping flag

---
 src/dynamicMesh/polyTopoChange/polyTopoChange/removeCells.C | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/dynamicMesh/polyTopoChange/polyTopoChange/removeCells.C b/src/dynamicMesh/polyTopoChange/polyTopoChange/removeCells.C
index 6be85d87c1..770d350e25 100644
--- a/src/dynamicMesh/polyTopoChange/polyTopoChange/removeCells.C
+++ b/src/dynamicMesh/polyTopoChange/polyTopoChange/removeCells.C
@@ -339,7 +339,7 @@ void Foam::removeCells::setRefinement
                         faceI,                  // label of face being modified
                         nei,                    // owner
                         -1,                     // neighbour
-                        false,                  // face flip
+                        true,                   // face flip
                         newPatchID[faceI],      // patch for face
                         false,                  // remove from zone
                         zoneID,                 // zone for face

From 6791d105703f0858d680d1b99ca3bebc7646c135 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Fri, 11 Dec 2009 11:59:28 +0000
Subject: [PATCH 11/16] Changed indentation

---
 src/dynamicMesh/polyTopoChange/polyTopoChange/removeFaces.H | 3 +--
 1 file changed, 1 insertion(+), 2 deletions(-)

diff --git a/src/dynamicMesh/polyTopoChange/polyTopoChange/removeFaces.H b/src/dynamicMesh/polyTopoChange/polyTopoChange/removeFaces.H
index 6d217a3f4f..e691a4f4e3 100644
--- a/src/dynamicMesh/polyTopoChange/polyTopoChange/removeFaces.H
+++ b/src/dynamicMesh/polyTopoChange/polyTopoChange/removeFaces.H
@@ -137,8 +137,7 @@ class removeFaces
             ) const;
 
             //- Return face with all pointsToRemove removed.
-            face filterFace(const labelHashSet& pointsToRemove, const label)
-             const;
+            face filterFace(const labelHashSet&, const label) const;
 
             //- Wrapper for meshMod.modifyFace. Reverses face if own>nei.
             void modFace

From 46a455ccea2d6e6874f8b898f0595de4d24b21cf Mon Sep 17 00:00:00 2001
From: Mark Olesen <Mark.Olesen@Germany>
Date: Fri, 11 Dec 2009 15:04:28 +0100
Subject: [PATCH 12/16] Remove unused Istream methods from surfMesh classes.

Move compilation edgeMesh and surfMesh up in the list
(they are independent of triSurface).
---
 src/Allwmake                                               | 2 +-
 src/surfMesh/UnsortedMeshedSurface/UnsortedMeshedSurface.H | 6 ------
 src/surfMesh/surfMesh/surfMeshIO.C                         | 2 +-
 3 files changed, 2 insertions(+), 8 deletions(-)

diff --git a/src/Allwmake b/src/Allwmake
index 6592f244d7..969687eb83 100755
--- a/src/Allwmake
+++ b/src/Allwmake
@@ -22,9 +22,9 @@ wmake libso OpenFOAM
 
 wmake libso lagrangian/basic
 
-wmake libso triSurface
 wmake libso edgeMesh
 wmake libso surfMesh
+wmake libso triSurface
 
 # Decomposition methods needed by meshTools
 wmake libso parallel/decompositionMethods
diff --git a/src/surfMesh/UnsortedMeshedSurface/UnsortedMeshedSurface.H b/src/surfMesh/UnsortedMeshedSurface/UnsortedMeshedSurface.H
index 7e16cc909d..4b46eb7c40 100644
--- a/src/surfMesh/UnsortedMeshedSurface/UnsortedMeshedSurface.H
+++ b/src/surfMesh/UnsortedMeshedSurface/UnsortedMeshedSurface.H
@@ -110,9 +110,6 @@ private:
         //- Disable setSize with value
         void setSize(const label, const Face&);
 
-        //- Read OpenFOAM Surface format
-        bool read(Istream&);
-
 
 protected:
 
@@ -198,9 +195,6 @@ public:
         //- Construct from file name (uses extension to determine type)
         UnsortedMeshedSurface(const fileName&, const word&);
 
-        //- Construct from Istream
-        UnsortedMeshedSurface(Istream&);
-
         //- Construct from objectRegistry and a named surface
         UnsortedMeshedSurface(const Time&, const word& surfName="");
 
diff --git a/src/surfMesh/surfMesh/surfMeshIO.C b/src/surfMesh/surfMesh/surfMeshIO.C
index ca2d840d07..7c2089f820 100644
--- a/src/surfMesh/surfMesh/surfMeshIO.C
+++ b/src/surfMesh/surfMesh/surfMeshIO.C
@@ -31,7 +31,7 @@ License
 
 void Foam::surfMesh::setInstance(const fileName& inst)
 {
-    if (debug or true)
+    if (debug)
     {
         Info<< "void surfMesh::setInstance(const fileName& inst) : "
             << "Resetting file instance to " << inst << endl;

From d016db1bc6193c609716cd21ed276a74219e7c43 Mon Sep 17 00:00:00 2001
From: Mark Olesen <Mark.Olesen@Germany>
Date: Fri, 11 Dec 2009 16:29:55 +0100
Subject: [PATCH 13/16] Enhance edgeMesh to support more formats.

Read only support:

  .bdf, .nas - NASTRAN format. Handles both CBEAM and CROD as lines.
  CROD is what Hypermesh happens to output.

Write only support:

  .vtk - VTK legacy format in ASCII

Read/write support:
  .eMesh - native format, which is simply a list of points, edges
      with an additional IOobject header that lets them be moved about
      easily to use as a featureEdgeMesh.

  .inp - STAR-CD inp/cel/vrt combination
       IOobject header)

  .obj - Alias waverfront format

Radically simplify surfaceFeatureConvert by using the new edgeMesh
functionality.
---
 applications/utilities/surface/README         |   2 +-
 .../surfaceFeatureConvert.C                   | 293 +++----------
 src/edgeMesh/Make/files                       |  23 +-
 src/edgeMesh/edgeFormats/edgeFormatsCore.C    | 200 +++++++++
 src/edgeMesh/edgeFormats/edgeFormatsCore.H    | 118 ++++++
 .../edgeFormats/emesh/EMESHedgeFormat.C       | 195 +++++++++
 .../edgeFormats/emesh/EMESHedgeFormat.H       | 148 +++++++
 .../emesh/EMESHedgeFormatRunTime.C            |  61 +++
 src/edgeMesh/edgeFormats/nas/NASedgeFormat.C  | 260 ++++++++++++
 src/edgeMesh/edgeFormats/nas/NASedgeFormat.H  | 113 +++++
 .../edgeFormats/nas/NASedgeFormatRunTime.C    |  60 +++
 src/edgeMesh/edgeFormats/obj/OBJedgeFormat.C  | 245 +++++++++++
 src/edgeMesh/edgeFormats/obj/OBJedgeFormat.H  | 121 ++++++
 .../edgeFormats/obj/OBJedgeFormatRunTime.C    |  61 +++
 .../edgeFormats/starcd/STARCDedgeFormat.C     | 387 ++++++++++++++++++
 .../edgeFormats/starcd/STARCDedgeFormat.H     | 159 +++++++
 .../starcd/STARCDedgeFormatRunTime.C          |  61 +++
 src/edgeMesh/edgeFormats/vtk/VTKedgeFormat.C  | 107 +++++
 src/edgeMesh/edgeFormats/vtk/VTKedgeFormat.H  | 116 ++++++
 .../edgeFormats/vtk/VTKedgeFormatRunTime.C    |  51 +++
 src/edgeMesh/edgeMesh.C                       | 184 ++++++++-
 src/edgeMesh/edgeMesh.H                       | 169 +++++++-
 src/edgeMesh/edgeMeshI.H                      |  27 +-
 src/edgeMesh/edgeMeshIO.C                     | 107 ++++-
 src/edgeMesh/edgeMeshNew.C                    |  65 +++
 .../{ => featureEdgeMesh}/featureEdgeMesh.C   |  59 ++-
 .../{ => featureEdgeMesh}/featureEdgeMesh.H   |  19 +-
 27 files changed, 3074 insertions(+), 337 deletions(-)
 create mode 100644 src/edgeMesh/edgeFormats/edgeFormatsCore.C
 create mode 100644 src/edgeMesh/edgeFormats/edgeFormatsCore.H
 create mode 100644 src/edgeMesh/edgeFormats/emesh/EMESHedgeFormat.C
 create mode 100644 src/edgeMesh/edgeFormats/emesh/EMESHedgeFormat.H
 create mode 100644 src/edgeMesh/edgeFormats/emesh/EMESHedgeFormatRunTime.C
 create mode 100644 src/edgeMesh/edgeFormats/nas/NASedgeFormat.C
 create mode 100644 src/edgeMesh/edgeFormats/nas/NASedgeFormat.H
 create mode 100644 src/edgeMesh/edgeFormats/nas/NASedgeFormatRunTime.C
 create mode 100644 src/edgeMesh/edgeFormats/obj/OBJedgeFormat.C
 create mode 100644 src/edgeMesh/edgeFormats/obj/OBJedgeFormat.H
 create mode 100644 src/edgeMesh/edgeFormats/obj/OBJedgeFormatRunTime.C
 create mode 100644 src/edgeMesh/edgeFormats/starcd/STARCDedgeFormat.C
 create mode 100644 src/edgeMesh/edgeFormats/starcd/STARCDedgeFormat.H
 create mode 100644 src/edgeMesh/edgeFormats/starcd/STARCDedgeFormatRunTime.C
 create mode 100644 src/edgeMesh/edgeFormats/vtk/VTKedgeFormat.C
 create mode 100644 src/edgeMesh/edgeFormats/vtk/VTKedgeFormat.H
 create mode 100644 src/edgeMesh/edgeFormats/vtk/VTKedgeFormatRunTime.C
 create mode 100644 src/edgeMesh/edgeMeshNew.C
 rename src/edgeMesh/{ => featureEdgeMesh}/featureEdgeMesh.C (81%)
 rename src/edgeMesh/{ => featureEdgeMesh}/featureEdgeMesh.H (85%)

diff --git a/applications/utilities/surface/README b/applications/utilities/surface/README
index 292d8ca0b5..a66c2463a1 100644
--- a/applications/utilities/surface/README
+++ b/applications/utilities/surface/README
@@ -21,7 +21,7 @@ surfaceFind
 - Finds nearest vertex and face to given point.
 
 surfaceMeshTriangulate
-- Triangulate external facses of mesh and write as surface.
+- Triangulate external faces of mesh and write as surface.
 
 surfacePointMerge
 - Explicit point merge of surface.
diff --git a/applications/utilities/surface/surfaceFeatureConvert/surfaceFeatureConvert.C b/applications/utilities/surface/surfaceFeatureConvert/surfaceFeatureConvert.C
index 046f9cefe8..ee221f992a 100644
--- a/applications/utilities/surface/surfaceFeatureConvert/surfaceFeatureConvert.C
+++ b/applications/utilities/surface/surfaceFeatureConvert/surfaceFeatureConvert.C
@@ -26,266 +26,85 @@ Application
     surfaceFeatureConvert
 
 Description
-    Extracts and writes surface features to file
+    Convert between edgeMesh formats
 
 \*---------------------------------------------------------------------------*/
 
-#include "featureEdgeMesh.H"
 #include "argList.H"
 #include "Time.H"
-#include "IFstream.H"
-#include "IStringStream.H"
-#include "OFstream.H"
-#include "Map.H"
+
+#include "edgeMesh.H"
 
 using namespace Foam;
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-void readNASEdges
-(
-    const fileName& inFileName,
-    pointField& allPoints,
-    edgeList& allEdges
-)
-{
-    IFstream is(inFileName);
-
-    if (!is.good())
-    {
-        FatalErrorIn("readNASEdges")
-            << "Cannot read file " << inFileName
-            << exit(FatalError);
-    }
-
-    // coordinates of point
-    DynamicList<point> points;
-    // Nastran index of point
-    DynamicList<label> pointIndices;
-
-    // beams
-    DynamicList<edge> edges;
-    DynamicList<label> edgeIndices;
-
-
-    while (is.good())
-    {
-        string line;
-        is.getLine(line);
-
-        if (line.empty() || line[0] == '$')
-        {
-            // Skip empty and comment
-            continue;
-        }
-
-        // Check if character 72 is continuation
-        if (line.size() > 72 && line[72] == '+')
-        {
-            line = line.substr(0, 72);
-
-            while (true)
-            {
-                string buf;
-                is.getLine(buf);
-
-                if (buf.size() > 72 && buf[72] == '+')
-                {
-                    line += buf.substr(8, 64);
-                }
-                else
-                {
-                    line += buf.substr(8, buf.size()-8);
-                    break;
-                }
-            }
-        }
-
-        // Read first word
-        IStringStream lineStream(line);
-        word cmd;
-        lineStream >> cmd;
-
-        if (cmd == "GRID")
-        {
-            label index;
-            lineStream >> index;
-            pointIndices.append(index);
-
-            scalar x = readScalar(IStringStream(line.substr(24, 8))());
-            scalar y = readScalar(IStringStream(line.substr(32, 8))());
-            scalar z = readScalar(IStringStream(line.substr(40, 8))());
-            points.append(point(x, y, z));
-        }
-        else if (cmd == "CBEAM")
-        {
-            // Read shell type since gives patchnames.
-            label index, group, v0, v1;
-            lineStream >> index >> group >> v0 >> v1;
-
-            edgeIndices.append(index);
-            edges.append(edge(v0, v1));
-        }
-    }
-
-    points.shrink();
-    pointIndices.shrink();
-    edges.shrink();
-    edgeIndices.shrink();
-
-    Pout<< "Read from " << inFileName
-        << " edges:" << edges.size() << " points:" << points.size()
-        << endl;
-
-    {
-        // Build inverse mapping (index to point)
-        Map<label> indexToPoint(2*pointIndices.size());
-        forAll(pointIndices, i)
-        {
-            indexToPoint.insert(pointIndices[i], i);
-        }
-
-        // Relabel edges
-        forAll(edges, i)
-        {
-            edge& e = edges[i];
-            e[0] = indexToPoint[e[0]];
-            e[1] = indexToPoint[e[1]];
-        }
-    }
-
-    allPoints.transfer(points);
-    allEdges.transfer(edges);
-}
-
-
-
-void write
-(
-    const Time& runTime,
-    const fileName& inFileName,
-    const fileName& outFileName,
-    const edgeMesh& eMesh
-)
-{
-    if (outFileName.ext() == "eMesh")
-    {
-        featureEdgeMesh fem
-        (
-            IOobject
-            (
-                outFileName,        // name
-                runTime.constant(), // instance
-                runTime,            // registry
-                IOobject::NO_READ,
-                IOobject::AUTO_WRITE,
-                false
-            ),
-            eMesh.points(),
-            eMesh.edges()
-        );
-
-        Pout<< "Writing feature edge mesh to " << fem.objectPath()
-            << endl;
-
-        fem.write();
-    }
-    else if (outFileName.ext() == "vtk")
-    {
-        OFstream str(outFileName);
-
-        str << "# vtk DataFile Version 2.0" << nl
-            << "featureEdgeMesh " << inFileName << nl
-            << "ASCII" << nl
-            << "DATASET POLYDATA" << nl;
-
-        str << "POINTS " << eMesh.points().size() << " float" << nl;
-        forAll(eMesh.points(), pointI)
-        {
-            const point& pt = eMesh.points()[pointI];
-
-            str << pt.x() << ' ' << pt.y() << ' ' << pt.z() << nl;
-        }
-
-        str << "LINES " << eMesh.edges().size() << ' '
-            << 3*eMesh.edges().size() << nl;
-        forAll(eMesh.edges(), edgeI)
-        {
-            const edge& e = eMesh.edges()[edgeI];
-
-            str << "2 " << e[0] << ' ' << e[1] << nl;
-        }
-    }
-    else
-    {
-        FatalErrorIn("write")
-            << "Supported output formats: .eMesh, .vtk"
-            << exit(FatalError);
-    }
-}
-
-
 // Main program:
 
 int main(int argc, char *argv[])
 {
+    argList::addNote
+    (
+        "Convert between edgeMesh formats"
+    );
     argList::noParallel();
-    argList::validArgs.append("input file");
-    argList::validArgs.append("output file");
-#   include "setRootCase.H"
-#   include "createTime.H"
+    argList::validArgs.append("inputFile");
+    argList::validArgs.append("outputFile");
+    argList::addOption
+    (
+        "scale",
+        "factor",
+        "specify a scaling factor for the points"
+    );
 
-    const fileName inFileName(args.additionalArgs()[0]);
-    const word outFileName(args.additionalArgs()[1]);
+    argList args(argc, argv);
+    Time runTime(args.rootPath(), args.caseName());
+    const stringList& params = args.additionalArgs();
 
-    Pout<< "Input features file  : " << inFileName << nl
-        << "Output features file : " << outFileName << nl
-        << endl;
+    const fileName importName(params[0]);
+    const fileName exportName(params[1]);
 
-
-    // Read
-    // ~~~~
-
-    if (inFileName.ext() == "nas")
-    {
-        pointField points;
-        edgeList edges;
-        readNASEdges(inFileName, points, edges);
-
-        edgeMesh eMesh(points, edges);
-
-        write(runTime, inFileName, outFileName, eMesh);
-    }
-    else if (inFileName.ext() == "eMesh")
-    {
-        featureEdgeMesh fem
-        (
-            IOobject
-            (
-                inFileName,         // name
-                runTime.constant(), // instance
-                runTime,            // registry
-                IOobject::MUST_READ,
-                IOobject::AUTO_WRITE,
-                false
-            )
-        );
-
-        Pout<< "Read from " << inFileName
-            << " edges:" << fem.edges().size()
-            << " points:" << fem.points().size()
-            << endl;
-
-        write(runTime, inFileName, outFileName, fem);
-    }
-    else
+    // disable inplace editing
+    if (importName == exportName)
     {
         FatalErrorIn(args.executable())
-            << "Can only handle NASTRAN data formats (.nas extension)."
+            << "Output file " << exportName << " would overwrite input file."
             << exit(FatalError);
     }
 
-    Pout<< "End\n" << endl;
+    // check that reading/writing is supported
+    if
+    (
+        !edgeMesh::canReadType(importName.ext(), true)
+     || !edgeMesh::canWriteType(exportName.ext(), true)
+    )
+    {
+        return 1;
+    }
+
+    edgeMesh mesh(importName);
+
+    Info<< "\nRead edgeMesh " << importName << nl;
+    mesh.writeStats(Info);
+    Info<< nl
+        << "\nwriting " << exportName;
+
+    scalar scaleFactor = 0;
+    if (args.optionReadIfPresent("scale", scaleFactor) && scaleFactor > 0)
+    {
+        Info<< " with scaling " << scaleFactor << endl;
+        mesh.scalePoints(scaleFactor);
+    }
+    else
+    {
+        Info<< " without scaling" << endl;
+    }
+
+    mesh.write(exportName);
+    mesh.writeStats(Info);
+    Info<< endl;
+
+    Info<< "\nEnd\n" << endl;
 
     return 0;
 }
diff --git a/src/edgeMesh/Make/files b/src/edgeMesh/Make/files
index 4017257d8b..cdad239999 100644
--- a/src/edgeMesh/Make/files
+++ b/src/edgeMesh/Make/files
@@ -1,5 +1,26 @@
 edgeMesh.C
 edgeMeshIO.C
-featureEdgeMesh.C
+edgeMeshNew.C
+
+edgeFormats = edgeFormats
+$(edgeFormats)/edgeFormatsCore.C
+
+$(edgeFormats)/emesh/EMESHedgeFormat.C
+$(edgeFormats)/emesh/EMESHedgeFormatRunTime.C
+
+$(edgeFormats)/nas/NASedgeFormat.C
+$(edgeFormats)/nas/NASedgeFormatRunTime.C
+
+$(edgeFormats)/obj/OBJedgeFormat.C
+$(edgeFormats)/obj/OBJedgeFormatRunTime.C
+
+$(edgeFormats)/starcd/STARCDedgeFormat.C
+$(edgeFormats)/starcd/STARCDedgeFormatRunTime.C
+
+$(edgeFormats)/vtk/VTKedgeFormat.C
+$(edgeFormats)/vtk/VTKedgeFormatRunTime.C
+
+featureEdgeMesh/featureEdgeMesh.C
+
 
 LIB = $(FOAM_LIBBIN)/libedgeMesh
diff --git a/src/edgeMesh/edgeFormats/edgeFormatsCore.C b/src/edgeMesh/edgeFormats/edgeFormatsCore.C
new file mode 100644
index 0000000000..8af7666cff
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/edgeFormatsCore.C
@@ -0,0 +1,200 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+\*---------------------------------------------------------------------------*/
+
+#include "edgeFormatsCore.H"
+
+#include "Time.H"
+#include "IFstream.H"
+#include "OFstream.H"
+#include "edgeMesh.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+Foam::word Foam::fileFormats::edgeFormatsCore::nativeExt("eMesh");
+
+
+// * * * * * * * * * * * * * Static Member Functions * * * * * * * * * * * * //
+
+Foam::string Foam::fileFormats::edgeFormatsCore::getLineNoComment
+(
+    IFstream& is
+)
+{
+    string line;
+    do
+    {
+        is.getLine(line);
+    }
+    while ((line.empty() || line[0] == '#') && is.good());
+
+    return line;
+}
+
+
+#if 0
+Foam::fileName Foam::fileFormats::edgeFormatsCore::localMeshFileName
+(
+    const word& meshName
+)
+{
+    const word name(meshName.size() ? meshName : surfaceRegistry::defaultName);
+
+    return fileName
+    (
+        surfaceRegistry::prefix/name/surfMesh::meshSubDir
+      / name + "." + nativeExt
+    );
+}
+
+
+Foam::fileName Foam::fileFormats::edgeFormatsCore::findMeshInstance
+(
+    const Time& t,
+    const word& meshName
+)
+{
+    fileName localName = localMeshFileName(meshName);
+
+    // Search back through the time directories list to find the time
+    // closest to and lower than current time
+
+    instantList ts = t.times();
+    label instanceI;
+
+    for (instanceI = ts.size()-1; instanceI >= 0; --instanceI)
+    {
+        if (ts[instanceI].value() <= t.timeOutputValue())
+        {
+            break;
+        }
+    }
+
+    // Noting that the current directory has already been searched
+    // for mesh data, start searching from the previously stored time directory
+
+    if (instanceI >= 0)
+    {
+        for (label i = instanceI; i >= 0; --i)
+        {
+            if (isFile(t.path()/ts[i].name()/localName))
+            {
+                return ts[i].name();
+            }
+        }
+    }
+
+    return "constant";
+}
+
+
+Foam::fileName Foam::fileFormats::edgeFormatsCore::findMeshFile
+(
+    const Time& t,
+    const word& meshName
+)
+{
+    fileName localName = localMeshFileName(meshName);
+
+    // Search back through the time directories list to find the time
+    // closest to and lower than current time
+
+    instantList ts = t.times();
+    label instanceI;
+
+    for (instanceI = ts.size()-1; instanceI >= 0; --instanceI)
+    {
+        if (ts[instanceI].value() <= t.timeOutputValue())
+        {
+            break;
+        }
+    }
+
+    // Noting that the current directory has already been searched
+    // for mesh data, start searching from the previously stored time directory
+
+    if (instanceI >= 0)
+    {
+        for (label i = instanceI; i >= 0; --i)
+        {
+            fileName testName(t.path()/ts[i].name()/localName);
+
+            if (isFile(testName))
+            {
+                return testName;
+            }
+        }
+    }
+
+    // fallback to "constant"
+    return t.path()/"constant"/localName;
+}
+#endif
+
+
+bool Foam::fileFormats::edgeFormatsCore::checkSupport
+(
+    const wordHashSet& available,
+    const word& ext,
+    const bool verbose,
+    const word& functionName
+)
+{
+    if (available.found(ext))
+    {
+        return true;
+    }
+    else if (verbose)
+    {
+        wordList known = available.sortedToc();
+
+        Info<<"Unknown file extension for " << functionName
+            << " : " << ext << nl
+            <<"Valid types: (";
+        // compact output:
+        forAll(known, i)
+        {
+            Info<<" " << known[i];
+        }
+        Info<<" )" << endl;
+    }
+
+    return false;
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::fileFormats::edgeFormatsCore::edgeFormatsCore()
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+Foam::fileFormats::edgeFormatsCore::~edgeFormatsCore()
+{}
+
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/edgeFormatsCore.H b/src/edgeMesh/edgeFormats/edgeFormatsCore.H
new file mode 100644
index 0000000000..a18f4f8a4b
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/edgeFormatsCore.H
@@ -0,0 +1,118 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+Class
+    Foam::fileFormats::edgeFormatsCore
+
+Description
+    A collection of helper functions for reading/writing edge formats.
+
+SourceFiles
+    edgeFormatsCore.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef edgeFormatsCore_H
+#define edgeFormatsCore_H
+
+#include "Map.H"
+#include "HashSet.H"
+#include "labelList.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+// Forward declaration of friend functions and operators
+
+class IFstream;
+class Time;
+
+namespace fileFormats
+{
+
+/*---------------------------------------------------------------------------*\
+                       Class edgeFormatsCore Declaration
+\*---------------------------------------------------------------------------*/
+
+class edgeFormatsCore
+{
+protected:
+
+        //- Read non-comment line
+        static string getLineNoComment(IFstream&);
+
+public:
+    // Static Data
+
+        //- The file extension corresponding to 'native' edge format
+        //  Normally "eMesh" (edge-mesh)
+        static word nativeExt;
+
+    // Static Member Functions
+
+        static bool checkSupport
+        (
+            const wordHashSet& available,
+            const word& ext,
+            const bool verbose,
+            const word& functionName
+        );
+
+//        //- Return the local file name (within time directory)
+//        //  NEEDS FIXING
+//        static fileName localMeshFileName(const word& edgeName="");
+//
+//        //- Find instance with edgeName
+//        //  NEEDS FIXING
+//        static fileName findMeshInstance(const Time&, const word& edgeName="");
+//
+//        //- Find mesh file with edgeName
+//        //  NEEDS FIXING
+//        static fileName findMeshFile(const Time&, const word& edgeName="");
+
+
+    // Constructors
+
+        //- Construct null
+        edgeFormatsCore();
+
+
+    //- Destructor
+    virtual ~edgeFormatsCore();
+
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace fileFormats
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/emesh/EMESHedgeFormat.C b/src/edgeMesh/edgeFormats/emesh/EMESHedgeFormat.C
new file mode 100644
index 0000000000..b085542f1d
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/emesh/EMESHedgeFormat.C
@@ -0,0 +1,195 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+\*---------------------------------------------------------------------------*/
+
+#include "EMESHedgeFormat.H"
+#include "IOobject.H"
+#include "IFstream.H"
+#include "clock.H"
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::fileFormats::EMESHedgeFormat::EMESHedgeFormat
+(
+    const fileName& filename
+)
+{
+    read(filename);
+}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+bool Foam::fileFormats::EMESHedgeFormat::read
+(
+    const fileName& filename
+)
+{
+    clear();
+
+    IFstream is(filename);
+    if (!is.good())
+    {
+        FatalErrorIn
+        (
+            "fileFormats::EMESHedgeFormat::read(const fileName&)"
+        )
+            << "Cannot read file " << filename
+            << exit(FatalError);
+    }
+
+    return read(is, this->storedPoints(), this->storedEdges());
+}
+
+
+bool Foam::fileFormats::EMESHedgeFormat::read
+(
+    Istream& is,
+    pointField& pointLst,
+    edgeList& edgeLst
+)
+{
+    if (!is.good())
+    {
+        FatalErrorIn
+        (
+            "fileFormats::EMESHedgeFormat::read"
+            "(Istream&, pointField&, edgeList&)"
+        )
+            << "read error "
+            << exit(FatalError);
+    }
+
+    token firstToken(is);
+
+    // swallow IOobject header
+    if (!is.good())
+    {
+        FatalIOErrorIn
+        (
+            "fileFormats::EMESHedgeFormat::read"
+            "(Istream&, pointField&, edgeList&)",
+            is
+        )
+            << "First token could not be read" << nl
+            << exit(FatalIOError);
+
+        return false;
+    }
+    else if (firstToken.isWord() && firstToken.wordToken() == "FoamFile")
+    {
+        // read and discard
+        dictionary headerDict(is);
+    }
+    else
+    {
+        is.putBack(firstToken);
+    }
+
+    // read points:
+    is  >> pointLst;
+
+    // read edges:
+    is  >> edgeLst;
+
+    return true;
+}
+
+
+Foam::Ostream&
+Foam::fileFormats::EMESHedgeFormat::write
+(
+    Ostream& os,
+    const pointField& pointLst,
+    const edgeList& edgeLst
+)
+{
+    if (!os.good())
+    {
+        FatalErrorIn
+        (
+            "fileFormats::EMESHedgeFormat::write"
+            "(Ostream&, const fileName&, const edgeMesh&)"
+        )
+            << "bad output stream " << os.name()
+            << exit(FatalError);
+    }
+
+    os  << "\n// points:" << nl << pointLst << nl
+        << "\n// edges:" << nl << edgeLst << nl;
+
+    IOobject::writeDivider(os);
+
+    // Check state of Ostream
+    os.check
+    (
+        "EMESHedgeFormat::write"
+        "(Ostream&, const pointField&, const edgeList&)"
+    );
+
+    return os;
+}
+
+
+void Foam::fileFormats::EMESHedgeFormat::write
+(
+    const fileName& filename,
+    const edgeMesh& mesh
+)
+{
+    OFstream os(filename);
+    if (!os.good())
+    {
+        FatalErrorIn
+        (
+            "fileFormats::EMESHedgeFormat::write"
+            "(const fileName&, const edgeMesh&)"
+        )
+            << "Cannot open file for writing " << filename
+            << exit(FatalError);
+    }
+
+
+    // just emit some information until we get a nicer IOobject
+    IOobject::writeBanner(os)
+        << "FoamFile\n{\n"
+        << "    version     " << os.version() << ";\n"
+        << "    format      " << os.format() << ";\n"
+        << "    class       " << "featureEdgeMesh" << ";\n"
+        << "    note        " << "written " + clock::dateTime() << nl
+        << "    object      " << filename.name() << ";\n"
+        << "}" << nl;
+
+    IOobject::writeDivider(os);
+
+    write(os, mesh.points(), mesh.edges());
+
+    // Check state of Ostream
+    os.check("EMESHedgeFormat::write(Ostream&)");
+}
+
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/emesh/EMESHedgeFormat.H b/src/edgeMesh/edgeFormats/emesh/EMESHedgeFormat.H
new file mode 100644
index 0000000000..370b5e7ae1
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/emesh/EMESHedgeFormat.H
@@ -0,0 +1,148 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+Class
+    Foam::fileFormats::EMESHedgeFormat
+
+Description
+    Provide a means of reading/writing the single-file OpenFOAM edge format.
+
+Note
+   This class provides more methods than the regular edge format interface.
+
+SourceFiles
+    EMESHedgeFormat.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef EMESHedgeFormat_H
+#define EMESHedgeFormat_H
+
+#include "edgeMesh.H"
+#include "Ostream.H"
+#include "OFstream.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace fileFormats
+{
+
+/*---------------------------------------------------------------------------*\
+                      Class EMESHedgeFormat Declaration
+\*---------------------------------------------------------------------------*/
+
+class EMESHedgeFormat
+:
+    public edgeMesh
+{
+    // Private Member Functions
+
+        //- Disallow default bitwise copy construct
+        EMESHedgeFormat(const EMESHedgeFormat&);
+
+        //- Disallow default bitwise assignment
+        void operator=(const EMESHedgeFormat&);
+
+
+protected:
+
+    // Protected Member Functions
+
+        //- Write header information
+        static void writeHeader
+        (
+            Ostream&,
+            const pointField&,
+            const edgeList&
+        );
+
+public:
+
+    // Constructors
+
+        //- Construct from file name
+        EMESHedgeFormat(const fileName&);
+
+
+    // Selectors
+
+        //- Read file and return edgeMesh
+        static autoPtr<edgeMesh> New(const fileName& name)
+        {
+            return autoPtr<edgeMesh>
+            (
+                new EMESHedgeFormat(name)
+            );
+        }
+
+
+    //- Destructor
+    virtual ~EMESHedgeFormat()
+    {}
+
+
+    // Member Functions
+
+        //- Read edgeMesh components from stream
+        static bool read
+        (
+            Istream&,
+            pointField&,
+            edgeList&
+        );
+
+        //- Write edgeMesh components to stream
+        static Ostream& write
+        (
+            Ostream&,
+            const pointField&,
+            const edgeList&
+        );
+
+        //- Write edgeMesh with a mimicked IOobject header
+        static void write(const fileName&, const edgeMesh&);
+
+        //- Read from file
+        virtual bool read(const fileName&);
+
+        //- Write object
+        virtual void write(const fileName& name) const
+        {
+            write(name, *this);
+        }
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace fileFormats
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/emesh/EMESHedgeFormatRunTime.C b/src/edgeMesh/edgeFormats/emesh/EMESHedgeFormatRunTime.C
new file mode 100644
index 0000000000..28e6ddd6d1
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/emesh/EMESHedgeFormatRunTime.C
@@ -0,0 +1,61 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+\*---------------------------------------------------------------------------*/
+
+#include "EMESHedgeFormat.H"
+
+#include "addToRunTimeSelectionTable.H"
+#include "addToMemberFunctionSelectionTable.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace fileFormats
+{
+
+// read edgeMesh
+addNamedToRunTimeSelectionTable
+(
+    edgeMesh,
+    EMESHedgeFormat,
+    fileExtension,
+    eMesh
+);
+
+// write edgeMesh
+addNamedToMemberFunctionSelectionTable
+(
+    edgeMesh,
+    EMESHedgeFormat,
+    write,
+    fileExtension,
+    eMesh
+);
+
+}
+}
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/nas/NASedgeFormat.C b/src/edgeMesh/edgeFormats/nas/NASedgeFormat.C
new file mode 100644
index 0000000000..4b178d1dbd
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/nas/NASedgeFormat.C
@@ -0,0 +1,260 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+\*---------------------------------------------------------------------------*/
+
+#include "NASedgeFormat.H"
+#include "IFstream.H"
+#include "IStringStream.H"
+#include "PackedBoolList.H"
+
+// * * * * * * * * * * * * Protected Member Functions  * * * * * * * * * * * //
+
+// Do weird things to extract a floating point number
+Foam::scalar Foam::fileFormats::NASedgeFormat::parseNASCoord
+(
+    const string& s
+)
+{
+    size_t expSign = s.find_last_of("+-");
+
+    if (expSign != string::npos && expSign > 0 && !isspace(s[expSign-1]))
+    {
+        scalar mantissa = readScalar(IStringStream(s.substr(0, expSign))());
+        scalar exponent = readScalar(IStringStream(s.substr(expSign+1))());
+
+        if (s[expSign] == '-')
+        {
+            exponent = -exponent;
+        }
+        return mantissa * pow(10, exponent);
+    }
+    else
+    {
+        return readScalar(IStringStream(s)());
+    }
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::fileFormats::NASedgeFormat::NASedgeFormat
+(
+    const fileName& filename
+)
+{
+    read(filename);
+}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+bool Foam::fileFormats::NASedgeFormat::read
+(
+    const fileName& filename
+)
+{
+    clear();
+
+    IFstream is(filename);
+    if (!is.good())
+    {
+        FatalErrorIn
+        (
+            "fileFormats::NASedgeFormat::read(const fileName&)"
+        )
+            << "Cannot read file " << filename
+            << exit(FatalError);
+    }
+
+    DynamicList<point>  dynPoints;
+    DynamicList<edge>   dynEdges;
+    DynamicList<label>  pointId;     // Nastran index of points
+
+    while (is.good())
+    {
+        string line;
+        is.getLine(line);
+
+        // Skip empty or comment
+        if (line.empty() || line[0] == '$')
+        {
+            continue;
+        }
+
+        // Check if character 72 is continuation
+        if (line.size() > 72 && line[72] == '+')
+        {
+            line = line.substr(0, 72);
+
+            while (true)
+            {
+                string buf;
+                is.getLine(buf);
+
+                if (buf.size() > 72 && buf[72] == '+')
+                {
+                    line += buf.substr(8, 64);
+                }
+                else
+                {
+                    line += buf.substr(8, buf.size()-8);
+                    break;
+                }
+            }
+        }
+
+
+        // Read first word
+        IStringStream lineStream(line);
+        word cmd;
+        lineStream >> cmd;
+
+        if (cmd == "CBEAM" || cmd == "CROD")
+        {
+            edge e;
+
+            // label groupId = readLabel(IStringStream(line.substr(16,8))());
+            e[0] = readLabel(IStringStream(line.substr(24,8))());
+            e[1] = readLabel(IStringStream(line.substr(32,8))());
+
+            // discard groupID
+            dynEdges.append(e);
+        }
+        else if (cmd == "GRID")
+        {
+            label index = readLabel(IStringStream(line.substr(8,8))());
+            scalar x = parseNASCoord(line.substr(24, 8));
+            scalar y = parseNASCoord(line.substr(32, 8));
+            scalar z = parseNASCoord(line.substr(40, 8));
+
+            pointId.append(index);
+            dynPoints.append(point(x, y, z));
+        }
+        else if (cmd == "GRID*")
+        {
+            // Long format is on two lines with '*' continuation symbol
+            // on start of second line.
+            // Typical line (spaces compacted)
+            // GRID*      126   0 -5.55999875E+02 -5.68730474E+02
+            // *         2.14897901E+02
+
+            label index = readLabel(IStringStream(line.substr(8,16))());
+            scalar x = parseNASCoord(line.substr(40, 16));
+            scalar y = parseNASCoord(line.substr(56, 16));
+
+            is.getLine(line);
+            if (line[0] != '*')
+            {
+                FatalErrorIn
+                (
+                    "fileFormats::NASedgeFormat::read(const fileName&)"
+                )
+                    << "Expected continuation symbol '*' when reading GRID*"
+                    << " (double precision coordinate) format" << nl
+                    << "Read:" << line << nl
+                    << "File:" << is.name() << " line:" << is.lineNumber()
+                    << exit(FatalError);
+            }
+            scalar z = parseNASCoord(line.substr(8, 16));
+
+            pointId.append(index);
+            dynPoints.append(point(x, y, z));
+        }
+    }
+
+    // transfer to normal lists
+    storedPoints().transfer(dynPoints);
+
+    pointId.shrink();
+    dynEdges.shrink();
+
+    // Build inverse mapping (NASTRAN pointId -> index)
+    Map<label> mapPointId(2*pointId.size());
+    forAll(pointId, i)
+    {
+        mapPointId.insert(pointId[i], i);
+    }
+
+    // note which points were really used and which can be culled
+    PackedBoolList usedPoints(points().size());
+
+
+    // Pass1: relabel edges
+    // ~~~~~~~~~~~~~~~~~~~~
+    forAll(dynEdges, i)
+    {
+        edge& e = dynEdges[i];
+        e[0] = mapPointId[e[0]];
+        e[1] = mapPointId[e[1]];
+
+        usedPoints.set(e[0]);
+        usedPoints.set(e[1]);
+    }
+    pointId.clearStorage();
+    mapPointId.clear();
+
+    // not all the points were used, cull them accordingly
+    if (unsigned(points().size()) != usedPoints.count())
+    {
+        label nUsed = 0;
+
+        pointField& pts = storedPoints();
+        forAll(pts, pointI)
+        {
+            if (usedPoints.get(pointI))
+            {
+                if (nUsed != pointI)
+                {
+                    pts[nUsed] = pts[pointI];
+                }
+
+                // map prev -> new id
+                mapPointId[pointI] = nUsed;
+
+                ++nUsed;
+            }
+        }
+
+        pts.setSize(nUsed);
+
+        // renumber edge vertices
+        forAll(dynEdges, edgeI)
+        {
+            edge& e = dynEdges[edgeI];
+
+            e[0] = mapPointId[e[0]];
+            e[1] = mapPointId[e[1]];
+        }
+    }
+
+
+    // transfer to normal lists
+    storedEdges().transfer(dynEdges);
+
+    return true;
+}
+
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/nas/NASedgeFormat.H b/src/edgeMesh/edgeFormats/nas/NASedgeFormat.H
new file mode 100644
index 0000000000..600b880cb7
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/nas/NASedgeFormat.H
@@ -0,0 +1,113 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+Class
+    Foam::fileFormats::NASedgeFormat
+
+Description
+    Nastran edge reader.
+
+SourceFiles
+    NASedgeFormat.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef NASedgeFormat_H
+#define NASedgeFormat_H
+
+#include "edgeMesh.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace fileFormats
+{
+
+/*---------------------------------------------------------------------------*\
+                      Class NASedgeFormat Declaration
+\*---------------------------------------------------------------------------*/
+
+class NASedgeFormat
+:
+    public edgeMesh
+{
+    // Private Member Functions
+
+        //- Disallow default bitwise copy construct
+        NASedgeFormat(const NASedgeFormat&);
+
+        //- Disallow default bitwise assignment
+        void operator=(const NASedgeFormat&);
+
+protected:
+
+    // Protected Member Functions
+
+        //- Do weird things to extract number
+        static scalar parseNASCoord(const string&);
+
+public:
+
+    // Constructors
+
+        //- Construct from file name
+        NASedgeFormat(const fileName&);
+
+
+    // Selectors
+
+        //- Read file and return edge mesh
+        static autoPtr<edgeMesh> New(const fileName& name)
+        {
+            return autoPtr<edgeMesh>
+            (
+                new NASedgeFormat(name)
+            );
+        }
+
+
+    // Destructor
+
+        virtual ~NASedgeFormat()
+        {}
+
+
+    // Member Functions
+
+        //- Read from a file
+        virtual bool read(const fileName&);
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace fileFormats
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/nas/NASedgeFormatRunTime.C b/src/edgeMesh/edgeFormats/nas/NASedgeFormatRunTime.C
new file mode 100644
index 0000000000..ec03e44bbb
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/nas/NASedgeFormatRunTime.C
@@ -0,0 +1,60 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+\*---------------------------------------------------------------------------*/
+
+#include "NASedgeFormat.H"
+
+#include "addToRunTimeSelectionTable.H"
+#include "addToMemberFunctionSelectionTable.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace fileFormats
+{
+
+// read edgeMesh - .bdf (Bulk Data Format)
+addNamedToRunTimeSelectionTable
+(
+    edgeMesh,
+    NASedgeFormat,
+    fileExtension,
+    bdf
+);
+
+// read edgeMesh - .nas (Nastran Data Format)
+addNamedToRunTimeSelectionTable
+(
+    edgeMesh,
+    NASedgeFormat,
+    fileExtension,
+    nas
+);
+
+}
+}
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/obj/OBJedgeFormat.C b/src/edgeMesh/edgeFormats/obj/OBJedgeFormat.C
new file mode 100644
index 0000000000..530d2a2a97
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/obj/OBJedgeFormat.C
@@ -0,0 +1,245 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+\*---------------------------------------------------------------------------*/
+
+#include "OBJedgeFormat.H"
+#include "clock.H"
+#include "IFstream.H"
+#include "IStringStream.H"
+#include "Ostream.H"
+#include "OFstream.H"
+#include "ListOps.H"
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::fileFormats::OBJedgeFormat::OBJedgeFormat
+(
+    const fileName& filename
+)
+{
+    read(filename);
+}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+bool Foam::fileFormats::OBJedgeFormat::read
+(
+    const fileName& filename
+)
+{
+    clear();
+
+    IFstream is(filename);
+    if (!is.good())
+    {
+        FatalErrorIn
+        (
+            "fileFormats::OBJedgeFormat::read(const fileName&)"
+        )
+            << "Cannot read file " << filename
+            << exit(FatalError);
+    }
+
+    DynamicList<point> dynPoints;
+    DynamicList<edge>  dynEdges;
+    DynamicList<label> dynUsedPoints;
+
+    while (is.good())
+    {
+        string line = this->getLineNoComment(is);
+
+        // handle continuations
+        if (line[line.size()-1] == '\\')
+        {
+            line.substr(0, line.size()-1);
+            line += this->getLineNoComment(is);
+        }
+
+        // Read first word
+        IStringStream lineStream(line);
+        word cmd;
+        lineStream >> cmd;
+
+        if (cmd == "v")
+        {
+            scalar x, y, z;
+            lineStream >> x >> y >> z;
+
+            dynPoints.append(point(x, y, z));
+            dynUsedPoints.append(-1);
+        }
+        else if (cmd == "l")
+        {
+            edge edgeRead;
+
+            // Assume 'l' is followed by space.
+            string::size_type endNum = 1;
+
+            int nVerts = 0;
+            for (int count = 0; count < 2; ++count)
+            {
+                string::size_type startNum =
+                    line.find_first_not_of(' ', endNum);
+
+                if (startNum == string::npos)
+                {
+                    break;
+                }
+
+                endNum = line.find(' ', startNum);
+
+                string vertexSpec;
+                if (endNum != string::npos)
+                {
+                    vertexSpec = line.substr(startNum, endNum-startNum);
+                }
+                else
+                {
+                    vertexSpec = line.substr(startNum, line.size() - startNum);
+                }
+
+                string::size_type slashPos = vertexSpec.find('/');
+
+                label vertI = 0;
+                if (slashPos != string::npos)
+                {
+                    IStringStream intStream(vertexSpec.substr(0, slashPos));
+
+                    intStream >> vertI;
+                }
+                else
+                {
+                    IStringStream intStream(vertexSpec);
+
+                    intStream >> vertI;
+                }
+
+                edgeRead[nVerts++] = (vertI - 1); // change to zero-offset
+            }
+
+            if (nVerts >= 2)
+            {
+                dynUsedPoints[edgeRead[0]] = edgeRead[0];
+                dynUsedPoints[edgeRead[1]] = edgeRead[1];
+
+                dynEdges.append(edgeRead);
+            }
+        }
+    }
+
+    // cull unused points
+    label nUsed = 0;
+
+    forAll(dynPoints, pointI)
+    {
+        if (dynUsedPoints[pointI] >= 0)
+        {
+            if (nUsed != pointI)
+            {
+                dynPoints[nUsed] = dynPoints[pointI];
+                dynUsedPoints[pointI] = nUsed;   // new position
+            }
+            ++nUsed;
+        }
+    }
+
+    dynPoints.setSize(nUsed);
+
+    // transfer to normal lists
+    storedPoints().transfer(dynPoints);
+
+    // renumber edge vertices
+    if (nUsed != dynUsedPoints.size())
+    {
+        forAll(dynEdges, edgeI)
+        {
+            edge& e = dynEdges[edgeI];
+
+            e[0] = dynUsedPoints[e[0]];
+            e[1] = dynUsedPoints[e[1]];
+        }
+    }
+    storedEdges().transfer(dynEdges);
+
+    return true;
+}
+
+
+void Foam::fileFormats::OBJedgeFormat::write
+(
+    const fileName& filename,
+    const edgeMesh& mesh
+)
+{
+    const pointField& pointLst = mesh.points();
+    const edgeList&   edgeLst  = mesh.edges();
+
+    OFstream os(filename);
+    if (!os.good())
+    {
+        FatalErrorIn
+        (
+            "fileFormats::OBJedgeFormat::write"
+            "(const fileName&, const edgeMesh&)"
+        )
+            << "Cannot open file for writing " << filename
+            << exit(FatalError);
+    }
+
+
+    os  << "# Wavefront OBJ file written " << clock::dateTime().c_str() << nl
+        << "o " << os.name().lessExt().name() << nl
+        << nl
+        << "# points : " << pointLst.size() << nl
+        << "# lines  : " << edgeLst.size() << nl;
+
+    os  << nl
+        << "# <points count=\"" << pointLst.size() << "\">" << nl;
+
+    // Write vertex coords
+    forAll(pointLst, ptI)
+    {
+        const point& p = pointLst[ptI];
+
+        os  << "v " << p.x() << ' ' << p.y() << ' ' << p.z() << nl;
+    }
+
+    os  << "# </points>" << nl
+        << nl
+        << "# <edges count=\"" << edgeLst.size() << "\">" << endl;
+
+    // Write line connectivity
+    forAll(edgeLst, edgeI)
+    {
+        const edge& e = edgeLst[edgeI];
+
+        os << "l " << (e[0] + 1) << " " << (e[1] + 1) << nl;
+    }
+    os << "# </edges>" << endl;
+}
+
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/obj/OBJedgeFormat.H b/src/edgeMesh/edgeFormats/obj/OBJedgeFormat.H
new file mode 100644
index 0000000000..d6eab70b60
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/obj/OBJedgeFormat.H
@@ -0,0 +1,121 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+Class
+    Foam::fileFormats::OBJedgeFormat
+
+Description
+    Provide a means of reading/writing Alias/Wavefront OBJ format.
+
+    Does not handle negative vertex indices.
+
+SourceFiles
+    OBJedgeFormat.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef OBJedgeFormat_H
+#define OBJedgeFormat_H
+
+#include "edgeMesh.H"
+#include "IFstream.H"
+#include "Ostream.H"
+#include "OFstream.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace fileFormats
+{
+
+/*---------------------------------------------------------------------------*\
+                        Class OBJedgeFormat Declaration
+\*---------------------------------------------------------------------------*/
+
+class OBJedgeFormat
+:
+    public edgeMesh
+{
+    // Private Member Functions
+
+        //- Disallow default bitwise copy construct
+        OBJedgeFormat(const OBJedgeFormat&);
+
+        //- Disallow default bitwise assignment
+        void operator=(const OBJedgeFormat&);
+
+
+public:
+
+    // Constructors
+
+        //- Construct from file name
+        OBJedgeFormat(const fileName&);
+
+
+    // Selectors
+
+        //- Read file and return surface
+        static autoPtr<edgeMesh> New(const fileName& name)
+        {
+            return autoPtr<edgeMesh>
+            (
+                new OBJedgeFormat(name)
+            );
+        }
+
+
+    // Destructor
+
+        virtual ~OBJedgeFormat()
+        {}
+
+
+    // Member Functions
+
+        //- Write surface mesh components by proxy
+        static void write(const fileName&, const edgeMesh&);
+
+        //- Read from file
+        virtual bool read(const fileName&);
+
+        //- Write object file
+        virtual void write(const fileName& name) const
+        {
+            write(name, *this);
+        }
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace fileFormats
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/obj/OBJedgeFormatRunTime.C b/src/edgeMesh/edgeFormats/obj/OBJedgeFormatRunTime.C
new file mode 100644
index 0000000000..b223f44c69
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/obj/OBJedgeFormatRunTime.C
@@ -0,0 +1,61 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+\*---------------------------------------------------------------------------*/
+
+#include "OBJedgeFormat.H"
+
+#include "addToRunTimeSelectionTable.H"
+#include "addToMemberFunctionSelectionTable.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace fileFormats
+{
+
+// read edgeMesh
+addNamedToRunTimeSelectionTable
+(
+    edgeMesh,
+    OBJedgeFormat,
+    fileExtension,
+    obj
+);
+
+// write edgeMesh
+addNamedToMemberFunctionSelectionTable
+(
+    edgeMesh,
+    OBJedgeFormat,
+    write,
+    fileExtension,
+    obj
+);
+
+}
+}
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/starcd/STARCDedgeFormat.C b/src/edgeMesh/edgeFormats/starcd/STARCDedgeFormat.C
new file mode 100644
index 0000000000..b6ecd4d19e
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/starcd/STARCDedgeFormat.C
@@ -0,0 +1,387 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+\*---------------------------------------------------------------------------*/
+
+#include "STARCDedgeFormat.H"
+#include "ListOps.H"
+#include "clock.H"
+#include "PackedBoolList.H"
+#include "IStringStream.H"
+
+// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
+
+inline void Foam::fileFormats::STARCDedgeFormat::writeLines
+(
+    Ostream& os,
+    const edgeList& edges
+)
+{
+    writeHeader(os, "CELL");
+
+    forAll(edges, edgeI)
+    {
+        const edge& e = edges[edgeI];
+        const label cellId = edgeI + 1;
+
+        os  << cellId                    // includes 1 offset
+            << ' ' << starcdLineShape_   // 2(line) shape
+            << ' ' << e.size()
+            << ' ' << 401                // arbitrary value
+            << ' ' << starcdLineType_;   // 5(line)
+
+        os  << nl << "  " << cellId << "  "
+            << (e[0]+1) << "  " << (e[1]+1) << nl;
+    }
+}
+
+
+// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
+
+bool Foam::fileFormats::STARCDedgeFormat::readHeader
+(
+    IFstream& is,
+    const word& signature
+)
+{
+    if (!is.good())
+    {
+        FatalErrorIn
+        (
+            "fileFormats::STARCDedgeFormat::readHeader(...)"
+        )
+            << "cannot read " << signature  << "  " << is.name()
+            << abort(FatalError);
+    }
+
+    word header;
+    label majorVersion;
+
+    string line;
+
+    is.getLine(line);
+    IStringStream(line)() >> header;
+
+    is.getLine(line);
+    IStringStream(line)() >> majorVersion;
+
+    // add other checks ...
+    if (header != signature)
+    {
+        Info<< "header mismatch " << signature << "  " << is.name()
+            << endl;
+    }
+
+    return true;
+}
+
+
+void Foam::fileFormats::STARCDedgeFormat::writeHeader
+(
+    Ostream& os,
+    const char* filetype
+)
+{
+    os  << "PROSTAR_" << filetype << nl
+        << 4000
+        << " " << 0
+        << " " << 0
+        << " " << 0
+        << " " << 0
+        << " " << 0
+        << " " << 0
+        << " " << 0
+        << endl;
+}
+
+
+bool Foam::fileFormats::STARCDedgeFormat::readPoints
+(
+    IFstream& is,
+    pointField& points,
+    labelList& ids
+)
+{
+    //
+    // read .vrt file
+    // ~~~~~~~~~~~~~~
+
+    if (!is.good())
+    {
+        FatalErrorIn
+        (
+            "fileFormats::STARCDedgeFormat::readPoints(...)"
+        )
+            << "Cannot read file " << is.name()
+            << exit(FatalError);
+    }
+
+    readHeader(is, "PROSTAR_VERTEX");
+
+    DynamicList<point> dynPoints;
+    DynamicList<label> dynPointId;    // STAR-CD index of points
+
+    label lineLabel;
+    while ((is >> lineLabel).good())
+    {
+        scalar x, y, z;
+
+        is >> x >> y >> z;
+
+        dynPoints.append(point(x, y, z));
+        dynPointId.append(lineLabel);
+    }
+
+    points.transfer(dynPoints);
+    ids.transfer(dynPointId);
+
+    return true;
+}
+
+
+
+void Foam::fileFormats::STARCDedgeFormat::writePoints
+(
+    Ostream& os,
+    const pointField& pointLst
+)
+{
+    writeHeader(os, "VERTEX");
+
+    // Set the precision of the points data to 10
+    os.precision(10);
+
+    // force decimal point for Fortran input
+    os.setf(std::ios::showpoint);
+
+    forAll(pointLst, ptI)
+    {
+        os
+            << ptI + 1 << " "
+            << pointLst[ptI].x() << " "
+            << pointLst[ptI].y() << " "
+            << pointLst[ptI].z() << nl;
+    }
+    os.flush();
+}
+
+
+void Foam::fileFormats::STARCDedgeFormat::writeCase
+(
+    Ostream& os,
+    const pointField& pointLst,
+    const label nEdges
+)
+{
+    word caseName = os.name().lessExt().name();
+
+    os  << "! STAR-CD file written " << clock::dateTime().c_str() << nl
+        << "! " << pointLst.size() << " points, " << nEdges << " lines" << nl
+        << "! case " << caseName << nl
+        << "! ------------------------------" << nl;
+
+//     forAll(zoneLst, zoneI)
+//     {
+//         os  << "ctable " << zoneI + 1 << " line" << nl
+//             << "ctname " << zoneI + 1 << " "
+//             << zoneLst[zoneI].name() << nl;
+//     }
+
+    os  << "! ------------------------------" << nl
+        << "*set icvo mxv - 1" << nl
+        << "vread " << caseName << ".vrt icvo,,,coded" << nl
+        << "cread " << caseName << ".cel icvo,,,add,coded" << nl
+        << "*set icvo" << nl
+        << "! end" << nl;
+
+    os.flush();
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::fileFormats::STARCDedgeFormat::STARCDedgeFormat
+(
+    const fileName& filename
+)
+{
+    read(filename);
+}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+bool Foam::fileFormats::STARCDedgeFormat::read
+(
+    const fileName& filename
+)
+{
+    clear();
+
+    fileName baseName = filename.lessExt();
+
+    // STAR-CD index of points
+    List<label> pointId;
+
+    // read points from .vrt file
+    readPoints
+    (
+        IFstream(baseName + ".vrt")(),
+        storedPoints(),
+        pointId
+    );
+
+    // Build inverse mapping (STAR-CD pointId -> index)
+    Map<label> mapPointId(2*pointId.size());
+    forAll(pointId, i)
+    {
+        mapPointId.insert(pointId[i], i);
+    }
+    pointId.clear();
+
+    // note which points were really used and which can be culled
+    PackedBoolList usedPoints(points().size());
+
+    //
+    // read .cel file
+    // ~~~~~~~~~~~~~~
+    IFstream is(baseName + ".cel");
+    if (!is.good())
+    {
+        FatalErrorIn
+        (
+            "fileFormats::STARCDedgeFormat::read(const fileName&)"
+        )
+            << "Cannot read file " << is.name()
+            << exit(FatalError);
+    }
+
+    readHeader(is, "PROSTAR_CELL");
+
+    DynamicList<edge>  dynEdges;
+
+    label lineLabel, shapeId, nLabels, cellTableId, typeId;
+    DynamicList<label> vertexLabels(64);
+
+    while ((is >> lineLabel).good())
+    {
+        is >> shapeId >> nLabels >> cellTableId >> typeId;
+
+        vertexLabels.clear();
+        vertexLabels.reserve(nLabels);
+
+        // read indices - max 8 per line
+        for (label i = 0; i < nLabels; ++i)
+        {
+            label vrtId;
+            if ((i % 8) == 0)
+            {
+               is >> lineLabel;
+            }
+            is >> vrtId;
+
+            // convert original vertex id to point label
+            vertexLabels.append(mapPointId[vrtId]);
+        }
+
+        if (typeId == starcdLineType_)
+        {
+            if (vertexLabels.size() >= 2)
+            {
+                dynEdges.append(edge(vertexLabels[0], vertexLabels[1]));
+
+                usedPoints.set(vertexLabels[0]);
+                usedPoints.set(vertexLabels[1]);
+            }
+        }
+    }
+
+    mapPointId.clear();
+
+    // not all the points were used, cull them accordingly
+    if (unsigned(points().size()) != usedPoints.count())
+    {
+        label nUsed = 0;
+
+        pointField& pts = storedPoints();
+        forAll(pts, pointI)
+        {
+            if (usedPoints.get(pointI))
+            {
+                if (nUsed != pointI)
+                {
+                    pts[nUsed] = pts[pointI];
+                }
+
+                // map prev -> new id
+                mapPointId.set(pointI, nUsed);
+
+                ++nUsed;
+            }
+        }
+
+        pts.setSize(nUsed);
+
+        // renumber edge vertices
+        forAll(dynEdges, edgeI)
+        {
+            edge& e = dynEdges[edgeI];
+
+            e[0] = mapPointId[e[0]];
+            e[1] = mapPointId[e[1]];
+        }
+    }
+
+
+    storedEdges().transfer(dynEdges);
+
+    return true;
+}
+
+
+void Foam::fileFormats::STARCDedgeFormat::write
+(
+    const fileName& filename,
+    const edgeMesh& mesh
+)
+{
+    const pointField& pointLst = mesh.points();
+    const edgeList&   edgeLst  = mesh.edges();
+
+    fileName baseName = filename.lessExt();
+
+    writePoints(OFstream(baseName + ".vrt")(), pointLst);
+    writeLines(OFstream(baseName + ".cel")(), edgeLst);
+
+    // write a simple .inp file
+    writeCase
+    (
+        OFstream(baseName + ".inp")(),
+        pointLst,
+        edgeLst.size()
+    );
+}
+
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/starcd/STARCDedgeFormat.H b/src/edgeMesh/edgeFormats/starcd/STARCDedgeFormat.H
new file mode 100644
index 0000000000..31d4a8db5d
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/starcd/STARCDedgeFormat.H
@@ -0,0 +1,159 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+Class
+    Foam::fileFormats::STARCDedgeFormat
+
+Description
+    Read/write the lines from pro-STAR vrt/cel files.
+
+Note
+    Uses the extension @a .inp (input) to denote the format.
+
+See Also
+    Foam::meshReaders::STARCD
+
+SourceFiles
+    STARCDedgeFormat.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef STARCDedgeFormat_H
+#define STARCDedgeFormat_H
+
+#include "edgeMesh.H"
+#include "IFstream.H"
+#include "Ostream.H"
+#include "OFstream.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace fileFormats
+{
+
+/*---------------------------------------------------------------------------*\
+                      Class STARCDedgeFormat Declaration
+\*---------------------------------------------------------------------------*/
+
+class STARCDedgeFormat
+:
+    public edgeMesh
+{
+    // Private Data
+
+        //- STAR-CD identifier for line shapes (1d elements)
+        static const int starcdLineShape_ = 2;
+
+        //- STAR-CD identifier for line type
+        static const int starcdLineType_  = 5;
+
+
+    // Private Member Functions
+
+        static inline void writeLines
+        (
+            Ostream&,
+            const edgeList&
+        );
+
+        //- Disallow default bitwise copy construct
+        STARCDedgeFormat(const STARCDedgeFormat&);
+
+        //- Disallow default bitwise assignment
+        void operator=(const STARCDedgeFormat&);
+
+
+protected:
+
+    // Protected Member Functions
+
+    static bool readHeader(IFstream&, const word&);
+
+    static void writeHeader(Ostream&, const char* filetype);
+
+    static bool readPoints(IFstream&, pointField&, labelList& ids);
+
+    static void writePoints(Ostream&, const pointField&);
+
+    static void writeCase
+    (
+        Ostream&,
+        const pointField&,
+        const label nEdges
+    );
+
+public:
+
+    // Constructors
+
+        //- Construct from file name
+        STARCDedgeFormat(const fileName&);
+
+
+    // Selectors
+
+        //- Read file and return edgeMesh
+        static autoPtr<edgeMesh> New(const fileName& name)
+        {
+            return autoPtr<edgeMesh>
+            (
+                new STARCDedgeFormat(name)
+            );
+        }
+
+
+    //- Destructor
+    virtual ~STARCDedgeFormat()
+    {}
+
+
+    // Member Functions
+
+        //- Write edge mesh
+        static void write(const fileName&, const edgeMesh&);
+
+        //- Read from file
+        virtual bool read(const fileName&);
+
+        //- Write object
+        virtual void write(const fileName& name) const
+        {
+            write(name, *this);
+        }
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace fileFormats
+} // End namespace Foam
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/starcd/STARCDedgeFormatRunTime.C b/src/edgeMesh/edgeFormats/starcd/STARCDedgeFormatRunTime.C
new file mode 100644
index 0000000000..192c0aeaab
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/starcd/STARCDedgeFormatRunTime.C
@@ -0,0 +1,61 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+\*---------------------------------------------------------------------------*/
+
+#include "STARCDedgeFormat.H"
+
+#include "addToRunTimeSelectionTable.H"
+#include "addToMemberFunctionSelectionTable.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace fileFormats
+{
+
+// read edgeMesh
+addNamedToRunTimeSelectionTable
+(
+    edgeMesh,
+    STARCDedgeFormat,
+    fileExtension,
+    inp
+);
+
+// write edgeMesh
+addNamedToMemberFunctionSelectionTable
+(
+    edgeMesh,
+    STARCDedgeFormat,
+    write,
+    fileExtension,
+    inp
+);
+
+}
+}
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/vtk/VTKedgeFormat.C b/src/edgeMesh/edgeFormats/vtk/VTKedgeFormat.C
new file mode 100644
index 0000000000..4d78b20d4e
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/vtk/VTKedgeFormat.C
@@ -0,0 +1,107 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+\*---------------------------------------------------------------------------*/
+
+#include "VTKedgeFormat.H"
+#include "OFstream.H"
+#include "clock.H"
+
+// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
+
+void Foam::fileFormats::VTKedgeFormat::writeHeader
+(
+    Ostream& os,
+    const pointField& pointLst
+)
+{
+    // Write header
+    os  << "# vtk DataFile Version 2.0" << nl
+        << "featureEdgeMesh written " << clock::dateTime().c_str() << nl
+        << "ASCII" << nl
+        << nl
+        << "DATASET POLYDATA" << nl;
+
+    // Write vertex coords
+    os  << "POINTS " << pointLst.size() << " float" << nl;
+    forAll(pointLst, ptI)
+    {
+        const point& pt = pointLst[ptI];
+
+        os  << pt.x() << ' ' << pt.y() << ' ' << pt.z() << nl;
+    }
+}
+
+
+void Foam::fileFormats::VTKedgeFormat::writeEdges
+(
+    Ostream& os,
+    const UList<edge>& edgeLst
+)
+{
+    os  << "LINES " << edgeLst.size() << ' '
+        << 3*edgeLst.size() << nl;
+
+    forAll(edgeLst, edgeI)
+    {
+        const edge& e = edgeLst[edgeI];
+
+        os  << "2 " << e[0] << ' ' << e[1] << nl;
+    }
+}
+
+
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::fileFormats::VTKedgeFormat::VTKedgeFormat()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+void Foam::fileFormats::VTKedgeFormat::write
+(
+    const fileName& filename,
+    const edgeMesh& eMesh
+)
+{
+    OFstream os(filename);
+    if (!os.good())
+    {
+        FatalErrorIn
+        (
+            "fileFormats::VTKedgeFormat::write"
+            "(const fileName&, const edgeMesh&)"
+        )
+            << "Cannot open file for writing " << filename
+            << exit(FatalError);
+    }
+
+    writeHeader(os, eMesh.points());
+    writeEdges(os, eMesh.edges());
+}
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/vtk/VTKedgeFormat.H b/src/edgeMesh/edgeFormats/vtk/VTKedgeFormat.H
new file mode 100644
index 0000000000..af4b4b78ed
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/vtk/VTKedgeFormat.H
@@ -0,0 +1,116 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+Class
+    Foam::fileFormats::VTKedgeFormat
+
+Description
+    Provide a means of writing VTK legacy format.
+
+SourceFiles
+    VTKedgeFormat.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef VTKedgeFormat_H
+#define VTKedgeFormat_H
+
+#include "edgeMesh.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace fileFormats
+{
+
+/*---------------------------------------------------------------------------*\
+                     Class VTKedgeFormat Declaration
+\*---------------------------------------------------------------------------*/
+
+class VTKedgeFormat
+:
+    public edgeMesh
+{
+    // Private Member Functions
+
+        //- Disallow default bitwise copy construct
+        VTKedgeFormat(const VTKedgeFormat&);
+
+        //- Disallow default bitwise assignment
+        void operator=(const VTKedgeFormat&);
+
+
+protected:
+
+    // Protected Member Functions
+
+        //- Write header information with points
+        static void writeHeader
+        (
+            Ostream&,
+            const pointField&
+        );
+
+        //- Write edges
+        static void writeEdges(Ostream&, const UList<edge>&);
+
+public:
+
+    // Constructors
+
+        //- Construct null
+        VTKedgeFormat();
+
+
+    //- Destructor
+    virtual ~VTKedgeFormat()
+    {}
+
+
+    // Member Functions
+
+        // Write
+
+            //- Write edgeMesh
+            static void write(const fileName&, const edgeMesh&);
+
+//            //- Write object
+//            virtual void write(Ostream& os) const
+//            {
+//                write(os, *this);
+//            }
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace fileFormats
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeFormats/vtk/VTKedgeFormatRunTime.C b/src/edgeMesh/edgeFormats/vtk/VTKedgeFormatRunTime.C
new file mode 100644
index 0000000000..7a9f4b859f
--- /dev/null
+++ b/src/edgeMesh/edgeFormats/vtk/VTKedgeFormatRunTime.C
@@ -0,0 +1,51 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+\*---------------------------------------------------------------------------*/
+
+#include "VTKedgeFormat.H"
+
+#include "addToRunTimeSelectionTable.H"
+#include "addToMemberFunctionSelectionTable.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace fileFormats
+{
+
+addNamedToMemberFunctionSelectionTable
+(
+    edgeMesh,
+    VTKedgeFormat,
+    write,
+    fileExtension,
+    vtk
+);
+
+}
+}
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/edgeMesh.C b/src/edgeMesh/edgeMesh.C
index d39810a62d..06502ef282 100644
--- a/src/edgeMesh/edgeMesh.C
+++ b/src/edgeMesh/edgeMesh.C
@@ -26,6 +26,76 @@ License
 
 #include "edgeMesh.H"
 #include "mergePoints.H"
+#include "addToRunTimeSelectionTable.H"
+#include "addToMemberFunctionSelectionTable.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+    defineTypeNameAndDebug(edgeMesh, 0);
+    defineRunTimeSelectionTable(edgeMesh, fileExtension);
+    defineMemberFunctionSelectionTable(edgeMesh,write,fileExtension);
+}
+
+Foam::wordHashSet Foam::edgeMesh::readTypes()
+{
+    return wordHashSet(*fileExtensionConstructorTablePtr_);
+}
+
+Foam::wordHashSet Foam::edgeMesh::writeTypes()
+{
+    return wordHashSet(*writefileExtensionMemberFunctionTablePtr_);
+}
+
+
+// * * * * * * * * * * * * * Static Member Functions * * * * * * * * * * * * //
+
+bool Foam::edgeMesh::canReadType
+(
+    const word& ext,
+    const bool verbose
+)
+{
+    return checkSupport
+    (
+        readTypes(),
+        ext,
+        verbose,
+        "reading"
+   );
+}
+
+bool Foam::edgeMesh::canWriteType
+(
+    const word& ext,
+    const bool verbose
+)
+{
+    return checkSupport
+    (
+        writeTypes(),
+        ext,
+        verbose,
+        "writing"
+    );
+}
+
+
+bool Foam::edgeMesh::canRead
+(
+    const fileName& name,
+    const bool verbose
+)
+{
+    word ext = name.ext();
+    if (ext == "gz")
+    {
+        ext = name.lessExt().ext();
+    }
+    return canReadType(ext, verbose);
+}
+
 
 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
 
@@ -63,10 +133,10 @@ void Foam::edgeMesh::calcPointEdges() const
     forAll(edges_, edgeI)
     {
         const edge& e = edges_[edgeI];
+        const label p0 = e[0];
+        const label p1 = e[1];
 
-        label p0 = e[0];
         pointEdges[p0][nEdgesPerPoint[p0]++] = edgeI;
-        label p1 = e[1];
         pointEdges[p1][nEdgesPerPoint[p1]++] = edgeI;
     }
 }
@@ -74,32 +144,102 @@ void Foam::edgeMesh::calcPointEdges() const
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-// construct from components
-Foam::edgeMesh::edgeMesh(const pointField& points, const edgeList& edges)
+Foam::edgeMesh::edgeMesh()
 :
-    points_(points),
-    edges_(edges)
-{}
-
-
-// construct as copy
-Foam::edgeMesh::edgeMesh(const edgeMesh& em)
-:
-    points_(em.points_),
-    edges_(em.edges_),
+    points_(0),
+    edges_(0),
     pointEdgesPtr_(NULL)
 {}
 
 
+Foam::edgeMesh::edgeMesh
+(
+    const pointField& points,
+    const edgeList& edges
+)
+:
+    points_(points),
+    edges_(edges),
+    pointEdgesPtr_(NULL)
+{}
+
+
+Foam::edgeMesh::edgeMesh
+(
+    const Xfer<pointField>& pointLst,
+    const Xfer<edgeList>& edgeLst
+)
+:
+    points_(0),
+    edges_(0),
+    pointEdgesPtr_(NULL)
+{
+    points_.transfer(pointLst());
+    edges_.transfer(edgeLst());
+}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+Foam::edgeMesh::~edgeMesh()
+{}
+
+
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
+void Foam::edgeMesh::clear()
+{
+    points_.clear();
+    edges_.clear();
+    pointEdgesPtr_.clear();
+}
+
+
+void Foam::edgeMesh::reset
+(
+    const Xfer< pointField >& pointLst,
+    const Xfer< edgeList >& edgeLst
+)
+{
+    // Take over new primitive data.
+    // Optimized to avoid overwriting data at all
+    if (&pointLst)
+    {
+        points_.transfer(pointLst());
+    }
+
+    if (&edgeLst)
+    {
+        edges_.transfer(edgeLst());
+
+        // connectivity likely changed
+        pointEdgesPtr_.clear();
+    }
+
+}
+
+
+void Foam::edgeMesh::transfer(edgeMesh& mesh)
+{
+    points_.transfer(mesh.points_);
+    edges_.transfer(mesh.edges_);
+    pointEdgesPtr_ = mesh.pointEdgesPtr_;
+}
+
+
+Foam::Xfer< Foam::edgeMesh >
+Foam::edgeMesh::xfer()
+{
+    return xferMove(*this);
+}
+
+
 Foam::label Foam::edgeMesh::regions(labelList& edgeRegion) const
 {
     edgeRegion.setSize(edges_.size());
     edgeRegion = -1;
 
     label startEdgeI = 0;
-
     label currentRegion = 0;
 
     while (true)
@@ -159,6 +299,16 @@ Foam::label Foam::edgeMesh::regions(labelList& edgeRegion) const
 }
 
 
+void Foam::edgeMesh::scalePoints(const scalar scaleFactor)
+{
+    // avoid bad scaling
+    if (scaleFactor > 0 && scaleFactor != 1.0)
+    {
+        points_ *= scaleFactor;
+    }
+}
+
+
 void Foam::edgeMesh::mergePoints(const scalar mergeDist)
 {
     pointField newPoints;
@@ -180,7 +330,7 @@ void Foam::edgeMesh::mergePoints(const scalar mergeDist)
 
         points_.transfer(newPoints);
 
-        // Renumber and make sure e[0] < e[1] (not really nessecary)
+        // Renumber and make sure e[0] < e[1] (not really necessary)
         forAll(edges_, edgeI)
         {
             edge& e = edges_[edgeI];
@@ -243,7 +393,7 @@ void Foam::edgeMesh::operator=(const edgeMesh& rhs)
 {
     points_ = rhs.points_;
     edges_ = rhs.edges_;
-    pointEdgesPtr_.reset(NULL);
+    pointEdgesPtr_.clear();
 }
 
 
diff --git a/src/edgeMesh/edgeMesh.H b/src/edgeMesh/edgeMesh.H
index 55934f3824..952a5378ed 100644
--- a/src/edgeMesh/edgeMesh.H
+++ b/src/edgeMesh/edgeMesh.H
@@ -26,7 +26,7 @@ Class
     Foam::edgeMesh
 
 Description
-    points connected by edges.
+    Points connected by edges.
 
 SourceFiles
     edgeMeshI.H
@@ -40,17 +40,33 @@ SourceFiles
 
 #include "pointField.H"
 #include "edgeList.H"
+#include "edgeFormatsCore.H"
+#include "runTimeSelectionTables.H"
+#include "memberFunctionSelectionTables.H"
+#include "HashSet.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 namespace Foam
 {
 
+// Forward declaration of classes
+class Istream;
+class Ostream;
+
+// Forward declaration of friend functions and operators
+class edgeMesh;
+Istream& operator>>(Istream&, edgeMesh&);
+Ostream& operator<<(Ostream&, const edgeMesh&);
+
+
 /*---------------------------------------------------------------------------*\
                            Class edgeMesh Declaration
 \*---------------------------------------------------------------------------*/
 
 class edgeMesh
+:
+    public fileFormats::edgeFormatsCore
 {
     // Private data
 
@@ -68,39 +84,179 @@ class edgeMesh
         //- Calculate point-edge addressing (inverse of edges)
         void calcPointEdges() const;
 
+
+protected:
+
+    // Protected Member Functions
+
+        //- Non-const access to global points
+        inline pointField& storedPoints();
+
+        //- Non-const access to the edges
+        inline edgeList& storedEdges();
+
+
 public:
 
+        //- Runtime type information
+        TypeName("edgeMesh");
+
+    // Static
+
+        //- Can we read this file format?
+        static bool canRead(const fileName&, const bool verbose=false);
+
+        //- Can we read this file format?
+        static bool canReadType(const word& ext, const bool verbose=false);
+
+        //- Can we write this file format type?
+        static bool canWriteType(const word& ext, const bool verbose=false);
+
+        static wordHashSet readTypes();
+        static wordHashSet writeTypes();
+
     // Constructors
 
+        //- Construct null
+        edgeMesh();
+
         //- Construct from components
         edgeMesh(const pointField&, const edgeList&);
 
-        //- Construct from file
-        edgeMesh(const fileName&);
-
-        //- Construct from Istream
-        edgeMesh(Istream&);
+        //- Construct by transferring components (points, edges).
+        edgeMesh
+        (
+            const Xfer< pointField >&,
+            const Xfer< edgeList >&
+        );
 
         //- Construct as copy
         edgeMesh(const edgeMesh&);
 
+        //- Construct from file name (uses extension to determine type)
+        edgeMesh(const fileName&);
+
+        //- Construct from file name (uses extension to determine type)
+        edgeMesh(const fileName&, const word& ext);
+
+        //- Construct from Istream
+        edgeMesh(Istream&);
+
+    // Declare run-time constructor selection table
+
+        declareRunTimeSelectionTable
+        (
+            autoPtr,
+            edgeMesh,
+            fileExtension,
+            (
+                const fileName& name
+            ),
+            (name)
+        );
+
+    // Selectors
+
+        //- Select constructed from filename (explicit extension)
+        static autoPtr<edgeMesh> New
+        (
+            const fileName&,
+            const word& ext
+        );
+
+        //- Select constructed from filename (implicit extension)
+        static autoPtr<edgeMesh> New(const fileName&);
+
+
+    //- Destructor
+    virtual ~edgeMesh();
+
+
+    // Member Function Selectors
+
+        declareMemberFunctionSelectionTable
+        (
+            void,
+            edgeMesh,
+            write,
+            fileExtension,
+            (
+                const fileName& name,
+                const edgeMesh& mesh
+            ),
+            (name, mesh)
+        );
+
+        //- Write to file
+        static void write(const fileName&, const edgeMesh&);
 
 
     // Member Functions
 
+        //- Transfer the contents of the argument and annul the argument
+        void transfer(edgeMesh&);
+
+        //- Transfer contents to the Xfer container
+        Xfer< edgeMesh > xfer();
+
+    // Read
+
+        //- Read from file. Chooses reader based on explicit extension
+        bool read(const fileName&, const word& ext);
+
+        //- Read from file. Chooses reader based on detected extension
+        virtual bool read(const fileName&);
+
+
+    // Access
+
+        //- Return points
         inline const pointField& points() const;
 
+        //- Return edges
         inline const edgeList& edges() const;
 
+        //- Return edges
         inline const labelListList& pointEdges() const;
 
         //- Find connected regions. Set region number per edge.
         //  Returns number of regions.
         label regions(labelList& edgeRegion) const;
 
+
+    // Edit
+
+        //- Clear all storage
+        virtual void clear();
+
+        //- Reset primitive data (points, edges)
+        //  Note, optimized to avoid overwriting data (with Xfer::null)
+        virtual void reset
+        (
+            const Xfer< pointField >& points,
+            const Xfer< edgeList >& edges
+        );
+
+        //- Scale points. A non-positive factor is ignored
+        virtual void scalePoints(const scalar);
+
         //- Merge common points (points within mergeDist)
         void mergePoints(const scalar mergeDist);
 
+
+    // Write
+
+        void writeStats(Ostream&) const;
+
+    // Write
+
+        //- Generic write routine. Chooses writer based on extension.
+        virtual void write(const fileName& name) const
+        {
+            write(name, *this);
+        }
+
+
     // Member Operators
 
         inline void operator=(const edgeMesh&);
@@ -109,6 +265,7 @@ public:
 
         friend Ostream& operator<<(Ostream&, const edgeMesh&);
         friend Istream& operator>>(Istream&, edgeMesh&);
+
 };
 
 
diff --git a/src/edgeMesh/edgeMeshI.H b/src/edgeMesh/edgeMeshI.H
index 85951e8735..e0dbb30a08 100644
--- a/src/edgeMesh/edgeMeshI.H
+++ b/src/edgeMesh/edgeMeshI.H
@@ -24,14 +24,14 @@ License
 
 \*---------------------------------------------------------------------------*/
 
-
-// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
-
-
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-
-// * * * * * * * * * * * * * * * * Destructors * * * * * * * * * * * * * * * //
+inline Foam::edgeMesh::edgeMesh(const edgeMesh& em)
+:
+    points_(em.points_),
+    edges_(em.edges_),
+    pointEdgesPtr_(NULL)
+{}
 
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
@@ -58,15 +58,16 @@ inline const Foam::labelListList& Foam::edgeMesh::pointEdges() const
 }
 
 
-// * * * * * * * * * * * * * * * Friend Operators  * * * * * * * * * * * * * //
+inline Foam::pointField& Foam::edgeMesh::storedPoints()
+{
+    return points_;
+}
 
 
-// * * * * * * * * * * * * * * IOstream Operators  * * * * * * * * * * * * * //
+inline Foam::edgeList& Foam::edgeMesh::storedEdges()
+{
+    return edges_;
+}
 
 
-// * * * * * * * * * * * * * * * Ostream Operator  * * * * * * * * * * * * * //
-
-
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
 // ************************************************************************* //
diff --git a/src/edgeMesh/edgeMeshIO.C b/src/edgeMesh/edgeMeshIO.C
index e2dffb864b..9b81c7b81f 100644
--- a/src/edgeMesh/edgeMeshIO.C
+++ b/src/edgeMesh/edgeMeshIO.C
@@ -25,42 +25,107 @@ License
 \*---------------------------------------------------------------------------*/
 
 #include "edgeMesh.H"
-#include "IFstream.H"
-
+#include "boundBox.H"
+#include "EMESHedgeFormat.H"
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-// construct from file
-Foam::edgeMesh::edgeMesh(const fileName& fname)
+Foam::edgeMesh::edgeMesh
+(
+    const fileName& name,
+    const word& ext
+)
 :
     points_(0),
     edges_(0),
     pointEdgesPtr_(NULL)
 {
-    IFstream is(fname);
+    read(name, ext);
+}
 
-    if (is.good())
+
+Foam::edgeMesh::edgeMesh(const fileName& name)
+:
+    points_(0),
+    edges_(0),
+    pointEdgesPtr_(NULL)
+{
+    read(name);
+}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+bool Foam::edgeMesh::read(const fileName& name)
+{
+    word ext = name.ext();
+    if (ext == "gz")
     {
-        is >> points_ >> edges_;
+        fileName unzipName = name.lessExt();
+        return read(unzipName, unzipName.ext());
     }
     else
     {
-        FatalErrorIn("edgeMesh::edgeMesh(const fileName&)")
-            << "cannot open file " << fname
-            << abort(FatalError);
+        return read(name, ext);
     }
 }
 
 
-// construct from Istream
-Foam::edgeMesh::edgeMesh(Istream& is)
-:
-    points_(is),
-    edges_(is),
-    pointEdgesPtr_(NULL)
+// Read from file in given format
+bool Foam::edgeMesh::read
+(
+    const fileName& name,
+    const word& ext
+)
 {
-    // Check state of Istream
-    is.check("edgeMesh::edgeMesh(Istream&)");
+    // read via selector mechanism
+    transfer(New(name, ext)());
+    return true;
+}
+
+
+void Foam::edgeMesh::write
+(
+    const fileName& name,
+    const edgeMesh& mesh
+)
+{
+    if (debug)
+    {
+        Info<< "edgeMesh::write"
+            "(const fileName&, const edgeMesh&) : "
+            "writing to " << name
+            << endl;
+    }
+
+    const word ext = name.ext();
+
+    writefileExtensionMemberFunctionTable::iterator mfIter =
+        writefileExtensionMemberFunctionTablePtr_->find(ext);
+
+    if (mfIter == writefileExtensionMemberFunctionTablePtr_->end())
+    {
+        FatalErrorIn
+        (
+            "MeshedSurface::write"
+            "(const fileName&, const MeshedSurface&)"
+        )   << "Unknown file extension " << ext << nl << nl
+            << "Valid types are :" << endl
+            << writefileExtensionMemberFunctionTablePtr_->sortedToc()
+            << exit(FatalError);
+    }
+    else
+    {
+        mfIter()(name, mesh);
+    }
+}
+
+
+void Foam::edgeMesh::writeStats(Ostream& os) const
+{
+    os  << "points      : " << points().size() << nl;
+    os  << "edges       : " << edges().size() << nl;
+    os  << "boundingBox : " << boundBox(this->points()) << endl;
 }
 
 
@@ -68,7 +133,7 @@ Foam::edgeMesh::edgeMesh(Istream& is)
 
 Foam::Ostream& Foam::operator<<(Ostream& os, const edgeMesh& em)
 {
-    os  << em.points_ << nl << em.edges_ << endl;
+    fileFormats::EMESHedgeFormat::write(os, em.points_, em.edges_);
 
     // Check state of Ostream
     os.check("Ostream& operator<<(Ostream&, const edgeMesh&)");
@@ -79,7 +144,9 @@ Foam::Ostream& Foam::operator<<(Ostream& os, const edgeMesh& em)
 
 Foam::Istream& Foam::operator>>(Istream& is, edgeMesh& em)
 {
-    is >> em.points_ >> em.edges_;
+    fileFormats::EMESHedgeFormat::read(is, em.points_, em.edges_);
+
+    em.pointEdgesPtr_.clear();
 
     // Check state of Istream
     is.check("Istream& operator>>(Istream&, edgeMesh&)");
diff --git a/src/edgeMesh/edgeMeshNew.C b/src/edgeMesh/edgeMeshNew.C
new file mode 100644
index 0000000000..4d63a61536
--- /dev/null
+++ b/src/edgeMesh/edgeMeshNew.C
@@ -0,0 +1,65 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 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
+
+\*---------------------------------------------------------------------------*/
+
+#include "edgeMesh.H"
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+
+Foam::autoPtr< Foam::edgeMesh >
+Foam::edgeMesh::New(const fileName& name, const word& ext)
+{
+    fileExtensionConstructorTable::iterator cstrIter =
+        fileExtensionConstructorTablePtr_->find(ext);
+
+    if (cstrIter == fileExtensionConstructorTablePtr_->end())
+    {
+        FatalErrorIn
+        (
+            "edgeMesh<Face>::New(const fileName&, const word&) : "
+            "constructing edgeMesh"
+        )   << "Unknown file extension " << ext << nl << nl
+            << "Valid types are :" << nl
+            << fileExtensionConstructorTablePtr_->sortedToc()
+            << exit(FatalError);
+    }
+
+    return autoPtr< edgeMesh >(cstrIter()(name));
+}
+
+
+Foam::autoPtr< Foam::edgeMesh >
+Foam::edgeMesh::New(const fileName& name)
+{
+    word ext = name.ext();
+    if (ext == "gz")
+    {
+        ext = name.lessExt().ext();
+    }
+    return New(name, ext);
+}
+
+// ************************************************************************* //
diff --git a/src/edgeMesh/featureEdgeMesh.C b/src/edgeMesh/featureEdgeMesh/featureEdgeMesh.C
similarity index 81%
rename from src/edgeMesh/featureEdgeMesh.C
rename to src/edgeMesh/featureEdgeMesh/featureEdgeMesh.C
index dab3c2cbb5..a579d4ff14 100644
--- a/src/edgeMesh/featureEdgeMesh.C
+++ b/src/edgeMesh/featureEdgeMesh/featureEdgeMesh.C
@@ -28,21 +28,7 @@ License
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
-namespace Foam
-{
-
-defineTypeNameAndDebug(featureEdgeMesh, 0);
-
-}
-
-
-// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
-
-
-// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
-
-
-// * * * * * * * * * * * * * * * Friend Operators  * * * * * * * * * * * * * //
+defineTypeNameAndDebug(Foam::featureEdgeMesh, 0);
 
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
@@ -73,20 +59,6 @@ Foam::featureEdgeMesh::featureEdgeMesh(const IOobject& io)
 }
 
 
-//- Construct from components
-Foam::featureEdgeMesh::featureEdgeMesh
-(
-    const IOobject& io,
-    const pointField& points,
-    const edgeList& edges
-)
-:
-    regIOobject(io),
-    edgeMesh(points, edges)
-{}
-
-
-// Construct as copy
 Foam::featureEdgeMesh::featureEdgeMesh
 (
     const IOobject& io,
@@ -98,24 +70,43 @@ Foam::featureEdgeMesh::featureEdgeMesh
 {}
 
 
+Foam::featureEdgeMesh::featureEdgeMesh
+(
+    const IOobject& io,
+    const Xfer< pointField >& pointLst,
+    const Xfer< edgeList >& edgeLst
+)
+:
+    regIOobject(io),
+    edgeMesh(pointLst, edgeLst)
+{
+    if
+    (
+        io.readOpt() == IOobject::MUST_READ
+     || (io.readOpt() == IOobject::READ_IF_PRESENT && headerOk())
+    )
+    {
+        readStream(typeName) >> *this;
+        close();
+    }
+}
+
+
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
 bool Foam::featureEdgeMesh::readData(Istream& is)
 {
-    is >> *this;
+    is  >> *this;
     return !is.bad();
 }
 
 
 bool Foam::featureEdgeMesh::writeData(Ostream& os) const
 {
-    os << *this;
+    os  << *this;
 
     return os.good();
 }
 
 
-// * * * * * * * * * * * * * * * Friend Operators  * * * * * * * * * * * * * //
-
-
 // ************************************************************************* //
diff --git a/src/edgeMesh/featureEdgeMesh.H b/src/edgeMesh/featureEdgeMesh/featureEdgeMesh.H
similarity index 85%
rename from src/edgeMesh/featureEdgeMesh.H
rename to src/edgeMesh/featureEdgeMesh/featureEdgeMesh.H
index b22859e799..8fc7b36744 100644
--- a/src/edgeMesh/featureEdgeMesh.H
+++ b/src/edgeMesh/featureEdgeMesh/featureEdgeMesh.H
@@ -45,7 +45,7 @@ namespace Foam
 {
 
 /*---------------------------------------------------------------------------*\
-                           Class featureEdgeMesh Declaration
+                       Class featureEdgeMesh Declaration
 \*---------------------------------------------------------------------------*/
 
 class featureEdgeMesh
@@ -56,25 +56,28 @@ class featureEdgeMesh
 
 public:
 
+    //- Runtime type information
     TypeName("featureEdgeMesh");
 
 
     // Constructors
 
         //- Construct (read) given an IOobject
-        featureEdgeMesh(const IOobject&);
+        explicit featureEdgeMesh(const IOobject&);
 
-        //- Construct from featureEdgeMesh data
+        //- Construct as copy
+        explicit featureEdgeMesh(const IOobject&, const featureEdgeMesh&);
+
+        //- Construct by transferring components (points, edges)
         featureEdgeMesh
         (
             const IOobject&,
-            const pointField&,
-            const edgeList&
+            const Xfer<pointField>&,
+            const Xfer<edgeList>&
         );
 
-        //- Construct as copy
-        featureEdgeMesh(const IOobject&, const featureEdgeMesh&);
-
+        //- Give precedence to the regIOobject write
+        using regIOobject::write;
 
         //- ReadData function required for regIOobject read operation
         virtual bool readData(Istream&);

From bf2865a08765223737381aecb7571c5e65969d91 Mon Sep 17 00:00:00 2001
From: Mark Olesen <Mark.Olesen@Germany>
Date: Fri, 11 Dec 2009 17:31:15 +0100
Subject: [PATCH 14/16] Add missing low-level peek() to ISstream.

---
 applications/test/tokenizeTest/tokenizeTest.C  | 8 +++++++-
 src/OpenFOAM/db/IOstreams/Sstreams/ISstream.H  | 6 ++++++
 src/OpenFOAM/db/IOstreams/Sstreams/ISstreamI.H | 6 ++++++
 3 files changed, 19 insertions(+), 1 deletion(-)

diff --git a/applications/test/tokenizeTest/tokenizeTest.C b/applications/test/tokenizeTest/tokenizeTest.C
index cd80997f61..c8be29cbf3 100644
--- a/applications/test/tokenizeTest/tokenizeTest.C
+++ b/applications/test/tokenizeTest/tokenizeTest.C
@@ -67,9 +67,15 @@ int main(int argc, char *argv[])
             while (is.good())
             {
                 token tok(is);
+                // char ch;
+                // is.get(ch);
+                // is.putback(ch);
+                int lookahead = is.peek();
+
                 if (count == 0)
                 {
-                    Info<< "token: " << tok.info() << endl;
+                    Info<< "token: " << tok.info();
+                    Info<< "  lookahead: '" << char(lookahead) << "'" << endl;
                 }
             }
 
diff --git a/src/OpenFOAM/db/IOstreams/Sstreams/ISstream.H b/src/OpenFOAM/db/IOstreams/Sstreams/ISstream.H
index fe778fa010..0c6bc9d44b 100644
--- a/src/OpenFOAM/db/IOstreams/Sstreams/ISstream.H
+++ b/src/OpenFOAM/db/IOstreams/Sstreams/ISstream.H
@@ -137,6 +137,12 @@ public:
             //- Raw, low-level get character function.
             inline ISstream& get(char&);
 
+            //- Raw, low-level peek function.
+            //  Does not remove the character from the stream.
+            //  Returns the next character in the stream or EOF if the
+            //  end of file is read.
+            inline int peek();
+
             //- Raw, low-level getline into a string function.
             inline ISstream& getLine(string&);
 
diff --git a/src/OpenFOAM/db/IOstreams/Sstreams/ISstreamI.H b/src/OpenFOAM/db/IOstreams/Sstreams/ISstreamI.H
index 1fedf76229..63b1087547 100644
--- a/src/OpenFOAM/db/IOstreams/Sstreams/ISstreamI.H
+++ b/src/OpenFOAM/db/IOstreams/Sstreams/ISstreamI.H
@@ -69,6 +69,12 @@ inline Foam::ISstream& Foam::ISstream::get(char& c)
 }
 
 
+inline int Foam::ISstream::peek()
+{
+    return is_.peek();
+}
+
+
 inline Foam::ISstream& Foam::ISstream::getLine(string& s)
 {
     getline(is_, s);

From 192f8ae542f4efd7d88d088988cd9413a1ba2ee2 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Mon, 14 Dec 2009 11:20:22 +0000
Subject: [PATCH 15/16] Changed tecio build options to use local file

---
 wmake/rules/SiCortex64Gcc/tecio | 3 ---
 wmake/rules/linux64Gcc/tecio    | 3 ---
 wmake/rules/linux64Gcc42/tecio  | 3 ---
 wmake/rules/linux64Gcc43/tecio  | 3 ---
 wmake/rules/linux64Gcc44/tecio  | 3 ---
 wmake/rules/linux64Icc/tecio    | 3 ---
 wmake/rules/linuxGcc/tecio      | 3 ---
 wmake/rules/linuxGcc42/tecio    | 3 ---
 wmake/rules/linuxGcc43/tecio    | 3 ---
 wmake/rules/linuxGcc44/tecio    | 3 ---
 wmake/rules/linuxIA64Gcc/tecio  | 3 ---
 wmake/rules/linuxIA64Icc/tecio  | 3 ---
 wmake/rules/linuxIcc/tecio      | 3 ---
 13 files changed, 39 deletions(-)
 delete mode 100644 wmake/rules/SiCortex64Gcc/tecio
 delete mode 100644 wmake/rules/linux64Gcc/tecio
 delete mode 100644 wmake/rules/linux64Gcc42/tecio
 delete mode 100644 wmake/rules/linux64Gcc43/tecio
 delete mode 100644 wmake/rules/linux64Gcc44/tecio
 delete mode 100644 wmake/rules/linux64Icc/tecio
 delete mode 100644 wmake/rules/linuxGcc/tecio
 delete mode 100644 wmake/rules/linuxGcc42/tecio
 delete mode 100644 wmake/rules/linuxGcc43/tecio
 delete mode 100644 wmake/rules/linuxGcc44/tecio
 delete mode 100644 wmake/rules/linuxIA64Gcc/tecio
 delete mode 100644 wmake/rules/linuxIA64Icc/tecio
 delete mode 100644 wmake/rules/linuxIcc/tecio

diff --git a/wmake/rules/SiCortex64Gcc/tecio b/wmake/rules/SiCortex64Gcc/tecio
deleted file mode 100644
index e0cf0911bc..0000000000
--- a/wmake/rules/SiCortex64Gcc/tecio
+++ /dev/null
@@ -1,3 +0,0 @@
-TECIOFLAGS     =-DMAKEARCHIVE -DLINUX -DLINUX64 -DUSEENUM -DTHREED -U_WIN32
-TECIOINC       =
-TECIOLIBS      =
diff --git a/wmake/rules/linux64Gcc/tecio b/wmake/rules/linux64Gcc/tecio
deleted file mode 100644
index e0cf0911bc..0000000000
--- a/wmake/rules/linux64Gcc/tecio
+++ /dev/null
@@ -1,3 +0,0 @@
-TECIOFLAGS     =-DMAKEARCHIVE -DLINUX -DLINUX64 -DUSEENUM -DTHREED -U_WIN32
-TECIOINC       =
-TECIOLIBS      =
diff --git a/wmake/rules/linux64Gcc42/tecio b/wmake/rules/linux64Gcc42/tecio
deleted file mode 100644
index e0cf0911bc..0000000000
--- a/wmake/rules/linux64Gcc42/tecio
+++ /dev/null
@@ -1,3 +0,0 @@
-TECIOFLAGS     =-DMAKEARCHIVE -DLINUX -DLINUX64 -DUSEENUM -DTHREED -U_WIN32
-TECIOINC       =
-TECIOLIBS      =
diff --git a/wmake/rules/linux64Gcc43/tecio b/wmake/rules/linux64Gcc43/tecio
deleted file mode 100644
index e0cf0911bc..0000000000
--- a/wmake/rules/linux64Gcc43/tecio
+++ /dev/null
@@ -1,3 +0,0 @@
-TECIOFLAGS     =-DMAKEARCHIVE -DLINUX -DLINUX64 -DUSEENUM -DTHREED -U_WIN32
-TECIOINC       =
-TECIOLIBS      =
diff --git a/wmake/rules/linux64Gcc44/tecio b/wmake/rules/linux64Gcc44/tecio
deleted file mode 100644
index e0cf0911bc..0000000000
--- a/wmake/rules/linux64Gcc44/tecio
+++ /dev/null
@@ -1,3 +0,0 @@
-TECIOFLAGS     =-DMAKEARCHIVE -DLINUX -DLINUX64 -DUSEENUM -DTHREED -U_WIN32
-TECIOINC       =
-TECIOLIBS      =
diff --git a/wmake/rules/linux64Icc/tecio b/wmake/rules/linux64Icc/tecio
deleted file mode 100644
index e0cf0911bc..0000000000
--- a/wmake/rules/linux64Icc/tecio
+++ /dev/null
@@ -1,3 +0,0 @@
-TECIOFLAGS     =-DMAKEARCHIVE -DLINUX -DLINUX64 -DUSEENUM -DTHREED -U_WIN32
-TECIOINC       =
-TECIOLIBS      =
diff --git a/wmake/rules/linuxGcc/tecio b/wmake/rules/linuxGcc/tecio
deleted file mode 100644
index ca44623d8f..0000000000
--- a/wmake/rules/linuxGcc/tecio
+++ /dev/null
@@ -1,3 +0,0 @@
-TECIOFLAGS     =-DMAKEARCHIVE -DLINUX -DUSEENUM -DTHREED -U_WIN32
-TECIOINC       =
-TECIOLIBS      =
diff --git a/wmake/rules/linuxGcc42/tecio b/wmake/rules/linuxGcc42/tecio
deleted file mode 100644
index ca44623d8f..0000000000
--- a/wmake/rules/linuxGcc42/tecio
+++ /dev/null
@@ -1,3 +0,0 @@
-TECIOFLAGS     =-DMAKEARCHIVE -DLINUX -DUSEENUM -DTHREED -U_WIN32
-TECIOINC       =
-TECIOLIBS      =
diff --git a/wmake/rules/linuxGcc43/tecio b/wmake/rules/linuxGcc43/tecio
deleted file mode 100644
index ca44623d8f..0000000000
--- a/wmake/rules/linuxGcc43/tecio
+++ /dev/null
@@ -1,3 +0,0 @@
-TECIOFLAGS     =-DMAKEARCHIVE -DLINUX -DUSEENUM -DTHREED -U_WIN32
-TECIOINC       =
-TECIOLIBS      =
diff --git a/wmake/rules/linuxGcc44/tecio b/wmake/rules/linuxGcc44/tecio
deleted file mode 100644
index ca44623d8f..0000000000
--- a/wmake/rules/linuxGcc44/tecio
+++ /dev/null
@@ -1,3 +0,0 @@
-TECIOFLAGS     =-DMAKEARCHIVE -DLINUX -DUSEENUM -DTHREED -U_WIN32
-TECIOINC       =
-TECIOLIBS      =
diff --git a/wmake/rules/linuxIA64Gcc/tecio b/wmake/rules/linuxIA64Gcc/tecio
deleted file mode 100644
index 6306e2098f..0000000000
--- a/wmake/rules/linuxIA64Gcc/tecio
+++ /dev/null
@@ -1,3 +0,0 @@
-TECIOFLAGS     =-DMAKEARCHIVE -DLINUX -DLINUXI64 -DUSEENUM -DTHREED -U_WIN32
-TECIOINC       =
-TECIOLIBS      =
diff --git a/wmake/rules/linuxIA64Icc/tecio b/wmake/rules/linuxIA64Icc/tecio
deleted file mode 100644
index 6306e2098f..0000000000
--- a/wmake/rules/linuxIA64Icc/tecio
+++ /dev/null
@@ -1,3 +0,0 @@
-TECIOFLAGS     =-DMAKEARCHIVE -DLINUX -DLINUXI64 -DUSEENUM -DTHREED -U_WIN32
-TECIOINC       =
-TECIOLIBS      =
diff --git a/wmake/rules/linuxIcc/tecio b/wmake/rules/linuxIcc/tecio
deleted file mode 100644
index ca44623d8f..0000000000
--- a/wmake/rules/linuxIcc/tecio
+++ /dev/null
@@ -1,3 +0,0 @@
-TECIOFLAGS     =-DMAKEARCHIVE -DLINUX -DUSEENUM -DTHREED -U_WIN32
-TECIOINC       =
-TECIOLIBS      =

From 5ca567d7cbfbd23df74e1d2c55a147f1a359508a Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Mon, 14 Dec 2009 11:20:57 +0000
Subject: [PATCH 16/16] Changed indentation

---
 .../autoHexMesh/refinementSurfaces/refinementSurfaces.H         | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/mesh/autoMesh/autoHexMesh/refinementSurfaces/refinementSurfaces.H b/src/mesh/autoMesh/autoHexMesh/refinementSurfaces/refinementSurfaces.H
index cee48bbe1f..86767f5731 100644
--- a/src/mesh/autoMesh/autoHexMesh/refinementSurfaces/refinementSurfaces.H
+++ b/src/mesh/autoMesh/autoHexMesh/refinementSurfaces/refinementSurfaces.H
@@ -291,7 +291,7 @@ public:
             (
                 const labelList& surfacesToTest,
                 const pointField& samples,
-                const  scalarField& nearestDistSqr,
+                const scalarField& nearestDistSqr,
                 labelList& hitSurface,
                 labelList& hitRegion
             ) const;