diff options
| -rw-r--r-- | apps/adrian/cylinder2d/cylinder2d.cpp | 12 |
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}; |