Move geometry print to prepareLattice

This way statistics reflect the cells disabled by refinement.

1 files changed, 6 insertions, 6 deletions

diff --git a/apps/adrian/cylinder2d/cylinder2d.cpp b/apps/adrian/cylinder2d/cylinder2d.cpp
index 0893928..56866c4 100644
--- a/apps/adrian/cylinder2d/cylinder2d.cpp
+++ b/apps/adrian/cylinder2d/cylinder2d.cpp
@@ -93,7 +93,6 @@ void prepareGeometry(Grid2D<T,DESCRIPTOR>& grid)
   sGeometry.clean();
   sGeometry.innerClean();
   sGeometry.checkForErrors();
-  sGeometry.print();
 
   clout << "Prepare Geometry ... OK" << std::endl;
 }
@@ -149,6 +148,7 @@ void prepareLattice(Grid2D<T,DESCRIPTOR>& grid)
   sLattice.initialize();
 
   clout << "Prepare lattice ... OK" << std::endl;
+  sGeometry.print();
 }
 
 void setBoundaryValues(Grid2D<T,DESCRIPTOR>& grid, int iT)
@@ -227,10 +227,12 @@ int main(int argc, char* argv[])
   auto& fineGrid = coarseGrid.refine(fineOrigin, fineExtend);
   prepareGeometry(fineGrid);
 
-  const Vector<T,2> fineExtendB {10*coarseDeltaX, ly};
-  const Vector<T,2> fineOriginB {lx-10*coarseDeltaX, 0};
+  disableRefinedArea(coarseGrid, fineGrid);
+
+  const Vector<T,2> fineOutflowExtend {5*coarseDeltaX, ly};
+  const Vector<T,2> fineOutflowOrigin {lx-5*coarseDeltaX, 0};
 
-  auto& fineOutflowGrid = coarseGrid.refine(fineOriginB, fineExtendB, false);
+  auto& fineOutflowGrid = coarseGrid.refine(fineOutflowOrigin, fineOutflowExtend, false);
   prepareGeometry(fineOutflowGrid);
 
   {
@@ -249,8 +251,6 @@ int main(int argc, char* argv[])
     coarseGrid.getSuperGeometry().rename(4,0,refined);
   }
 
-  disableRefinedArea(coarseGrid, fineGrid);
-
   const Vector<T,2> fineExtend2 {4*cylinderD, 2*cylinderD};
   const Vector<T,2> fineOrigin2 {1*cylinderD, 2*cylinderD-fineExtend2[1]/2};