From 1b61a46eec35fade3509437fba6f856a23d875af Mon Sep 17 00:00:00 2001
From: laurence <laurence>
Date: Mon, 3 Jun 2013 14:59:05 +0100
Subject: [PATCH] ENH: snappyHexMesh: update outfile location for
 surfaceSimplify

---
 .../generation/snappyHexMesh/snappyHexMesh.C  | 31 +++++++++++++------
 1 file changed, 21 insertions(+), 10 deletions(-)

diff --git a/applications/utilities/mesh/generation/snappyHexMesh/snappyHexMesh.C b/applications/utilities/mesh/generation/snappyHexMesh/snappyHexMesh.C
index fb177e101f..2f6c5b1488 100644
--- a/applications/utilities/mesh/generation/snappyHexMesh/snappyHexMesh.C
+++ b/applications/utilities/mesh/generation/snappyHexMesh/snappyHexMesh.C
@@ -137,8 +137,11 @@ autoPtr<refinementSurfaces> createRefinementSurfaces
                 const scalar surfaceCellSize =
                     readScalar(scsDict.lookup("surfaceCellSizeCoeff"));
 
-                globalMinLevel[surfI] = label(1.0/surfaceCellSize);
-                globalMaxLevel[surfI] = label(1.0/surfaceCellSize);
+                const label calculatedCellLevel =
+                    round(::log(1.0/surfaceCellSize)/::log(2) + 1);
+
+                globalMinLevel[surfI] = calculatedCellLevel;
+                globalMaxLevel[surfI] = calculatedCellLevel;
                 globalLevelIncr[surfI] = shapeDict.lookupOrDefault
                 (
                     "gapLevelIncrement",
@@ -319,8 +322,14 @@ autoPtr<refinementSurfaces> createRefinementSurfaces
                                     scsDict.lookup("surfaceCellSizeCoeff")
                                 );
 
-                            globalMinLevel[surfI] = label(1.0/surfaceCellSize);
-                            globalMaxLevel[surfI] = label(1.0/surfaceCellSize);
+                            const label calculatedCellLevel =
+                                round
+                                (
+                                    ::log(1.0/surfaceCellSize)/::log(2) + 1
+                                );
+
+                            globalMinLevel[surfI] = calculatedCellLevel;
+                            globalMaxLevel[surfI] = calculatedCellLevel;
                             globalLevelIncr[surfI] =
                                 shapeControlRegionDict.lookupOrDefault
                                 (
@@ -1518,12 +1527,14 @@ int main(int argc, char *argv[])
             }
         }
 
-        fileName outFileName("simplifiedSurface.obj");
-
-        if (args.optionFound("outFile"))
-        {
-            outFileName = word(args.optionLookup("outFile")());
-        }
+        fileName outFileName
+        (
+            args.optionLookupOrDefault<fileName>
+            (
+                "outFile",
+                "constant/triSurface/simplifiedSurface.stl"
+            )
+        );
 
         extractSurface
         (