fix SetForcingBlockVisitor exeption

source/Applications/pChannel/configBombadilpChannel.cfg

@@ -2,9 +2,9 @@
 #Simulation parameters for porous channel
 #
 
-pathOut = d:/temp/ChannelFlow-test
+pathOut = d:/temp/ChannelFlow
 pathGeo = d:/Projects/SFB880/GeometrienPoroeseMedien/isotrop/PA80-110
-numOfThreads = 1
+numOfThreads = 4
 availMem = 14e9
 logToFile = false
 
@@ -83,19 +83,20 @@ Re = 15000
 #real velocity is 54.95 m/s
 u_LB = 0.1
 
-newStart = true
-restartStep = 10
+newStart = false
+restartStep = 230000
 
-cpStep = 10
-cpStart = 10
+cpStep = 100
+cpStart = 100
 
 averaging = false
 averagingReset = false
 timeAvStart = 21000000
 timeAvStop = 2100010000
 
-endTime = 100
-outTime = 10000
+outTime = 100
+endTime = 230000
+
  
-nupsStep = 100 100 10000000 
+nupsStep = 10 10 10000000 
 

source/Applications/pChannel/pChannel.cpp

@@ -505,6 +505,7 @@ void run(string configname)
       calculator->addCoProcessor(nupsCoProcessor);
       calculator->addCoProcessor(AdjForcCoProcessor);
       calculator->addCoProcessor(migCoProcessor);
+      //calculator->addCoProcessor(restartCoProcessor);
       calculator->addCoProcessor(writeMQSelectCoProcessor);
       calculator->addCoProcessor(writeMQCoProcessor);
 

source/VirtualFluidsCore/Visitors/SetForcingBlockVisitor.cpp

@@ -32,12 +32,12 @@ SetForcingBlockVisitor::SetForcingBlockVisitor(const std::string& sForcingX1, co
 //////////////////////////////////////////////////////////////////////////
 void SetForcingBlockVisitor::visit(SPtr<Grid3D> grid, SPtr<Block3D> block)
 {
-    SPtr<LBMKernel> kernel = dynamicPointerCast<LBMKernel>(block->getKernel());
-    if (!kernel)
-        throw std::runtime_error("SetForcingBlockVisitor: Kernel is not a LBMKernel");
-
    if(block->getRank() == grid->getRank())
    {
+      SPtr<LBMKernel> kernel = dynamicPointerCast<LBMKernel>(block->getKernel());
+      if (!kernel)
+         throw UbException(UB_EXARGS, "LBMKernel is not exist");
+
       switch (ftype)
       {
       case 0: