fixes the "findInterfac..." methods, such that they do not overwrite Stoppers

src/GridGenerator/grid/GridBuilder/MultipleGridBuilder.cpp

@@ -343,15 +343,15 @@ void MultipleGridBuilder::buildGrids(LbmOrGks lbmOrGks)
     for (auto grid : grids)
         grid->inital();
 
- //   if (solidObject)
- //   {
- //       grids[grids.size() - 1]->mesh(solidObject);
- //       grids[grids.size() - 1]->findQs(solidObject);
- //   }
+    if (solidObject)
+    {
+        grids[grids.size() - 1]->mesh(solidObject);
+        grids[grids.size() - 1]->findQs(solidObject);
+    }
 
 
- //   for (size_t i = 0; i < grids.size() - 1; i++)
- //       grids[i]->findGridInterface(grids[i + 1], lbmOrGks);
+    for (size_t i = 0; i < grids.size() - 1; i++)
+        grids[i]->findGridInterface(grids[i + 1], lbmOrGks);
 
 	//if (lbmOrGks == LBM) {
 	//	for (size_t i = 0; i < grids.size() - 1; i++)

src/GridGenerator/grid/GridImp.cu

@@ -422,6 +422,23 @@ void GridImp::setCellTo(uint index, char type)
 }
 
 
+void GridImp::setNonStopperOutOfGridCellTo(uint index, char type)
+{
+    real x, y, z;
+    this->transIndexToCoords(index, x, y, z);
+
+    Cell cell(x, y, z, this->delta);
+    for (const auto node : cell)
+    {
+        const uint nodeIndex = transCoordToIndex(node.x, node.y, node.z);
+
+        if( this->getFieldEntry( nodeIndex ) != STOPPER_OUT_OF_GRID && 
+            this->getFieldEntry( nodeIndex ) != STOPPER_OUT_OF_GRID_BOUNDARY )
+            this->field.setFieldEntry(nodeIndex, type);
+    }
+}
+
+
 HOST void GridImp::setPeriodicity(bool periodicityX, bool periodicityY, bool periodicityZ)
 {
     this->periodicityX = periodicityX;

src/GridGenerator/grid/GridImp.h

@@ -81,6 +81,7 @@ public:
     void setPeriodicityZ(bool periodicity) override;
 
     HOSTDEVICE void setCellTo(uint index, char type);
+    HOSTDEVICE void setNonStopperOutOfGridCellTo(uint index, char type);
 
     HOSTDEVICE int transCoordToIndex(const real &x, const real &y, const real &z) const override;
     HOSTDEVICE void transIndexToCoords(int index, real &x, real &y, real &z) const override;

src/GridGenerator/grid/GridInterface.cu

@@ -42,8 +42,8 @@ void GridInterface::findInterfaceCF(const uint& indexOnCoarseGrid, GridImp* coar
 
             cf.numberOfEntries++;
 
-            coarseGrid->setCellTo(indexOnCoarseGrid, FLUID_CFC);
-            fineGrid->setCellTo(indexOnFineGridCF, FLUID_CFF);
+            coarseGrid->setNonStopperOutOfGridCellTo(indexOnCoarseGrid, FLUID_CFC);
+            fineGrid->setNonStopperOutOfGridCellTo(indexOnFineGridCF, FLUID_CFF);
             break;
         }
     }
@@ -101,7 +101,7 @@ void GridInterface::findInterfaceFC(const uint& indexOnCoarseGrid, GridImp* coar
 
             fc.numberOfEntries++;
 
-            fineGrid->setCellTo(indexOnFineGridFC, FLUID_FCF);
+            fineGrid->setNonStopperOutOfGridCellTo(indexOnFineGridFC, FLUID_FCF);
             coarseGrid->getField().setFieldEntry(indexOnCoarseGrid, FLUID_FCC);
             break;
         }