Work in progress

meshLibrary/utilities/octrees/meshOctree/meshOctreeCreator/meshOctreeCreatorAdjustOctreeToSurface.C

@@ -674,6 +674,10 @@ void meshOctreeCreator::refineBoxesIntersectingEdgeMeshes()
 
         const fileName fName(dict.lookup("edgeFile"));
 
+        Info << "Reading edge mesh " << fName << endl;
+        edgeMesh em(fName);
+        Info << "Finished reading mesh" << endl;
+
         if( meshDictPtr_->found("anisotropicSources") )
         {
             //triSurf origSurf(fName);
@@ -753,6 +757,7 @@ void meshOctreeCreator::refineBoxesIntersectingEdgeMeshes()
         //- select boxes which need to be refined
         forAll(edgeMeshesPtr, emI)
         {
+            Info << "Refining edge mesh " << edgeMeshNames[emI] << endl;
             const edgeMesh& edgeMesh = edgeMeshesPtr[emI];
             const pointField& points = edgeMesh.points();
             const edgeList& edges = edgeMesh.edges();

meshLibrary/utilities/octrees/meshOctree/meshOctreeCreator/meshOctreeCreatorCreateOctreeBoxes.C

@@ -460,9 +460,12 @@ void meshOctreeCreator::createOctreeBoxes()
     Info << "Refining boundary" << endl;
     refineBoundary();
 
-    //- refine parts intersected with surface mesh serving as refinement sources
+    //- refine parts intersected by surface meshes acting as refinement sources
     refineBoxesIntersectingSurfaces();
 
+    //- refine parts intersected by edge meshes acting as refinement sources
+    refineBoxesIntersectingEdgeMeshes();
+
     //- perform automatic octree refinement
     if( !Pstream::parRun() )
     {