diff --git a/apps/benchmarks/MaterialTransport/utilities/InitializerFunctions.cpp b/apps/benchmarks/MaterialTransport/utilities/InitializerFunctions.cpp
index c210a8d42f011ab0ae4854aef1caf4d71f035e45..a3770e5ab5e0cb8ce998c86c64671ccc590702ad 100644
--- a/apps/benchmarks/MaterialTransport/utilities/InitializerFunctions.cpp
+++ b/apps/benchmarks/MaterialTransport/utilities/InitializerFunctions.cpp
@@ -155,8 +155,8 @@ void initFluidField(const shared_ptr< StructuredBlockStorage >& blocks, BlockDat
       WALBERLA_FOR_ALL_CELLS_INCLUDING_GHOST_LAYER_XYZ(FluidVelocityField, {
          Cell globalCell;
          const auto cellAABB                 = blocks->getBlockLocalCellAABB(block, Cell(x, y, z));
-         FluidVelocityField->get(x, y, z, 0) = uInflow[0];
-         FluidVelocityField->get(x, y, z, 1) = uInflow[1];
+         FluidVelocityField->get(x, y, z, 0) = real_t(0);//uInflow[0];
+         FluidVelocityField->get(x, y, z, 1) = real_t(0);//uInflow[1];
          if(domainSize[2] != 1)
          {
             FluidVelocityField->get(x, y, z, 2) = 0;
diff --git a/src/lbm_mesapd_coupling/partially_saturated_cells_method/codegen/PSMCodegen.py b/src/lbm_mesapd_coupling/partially_saturated_cells_method/codegen/PSMCodegen.py
index a59edeae9d6388fb3deca96cc00ea4a774dee2d6..db7ccd9cb6480016acfeadb3dad7440233d05aa5 100644
--- a/src/lbm_mesapd_coupling/partially_saturated_cells_method/codegen/PSMCodegen.py
+++ b/src/lbm_mesapd_coupling/partially_saturated_cells_method/codegen/PSMCodegen.py
@@ -97,8 +97,8 @@ with CodeGeneration() as ctx:
         stencil=stencil,
         method=methods[config_tokens[0]],
         relaxation_rate=omega,
-        #force=sp.symbols("F_:3"),
-        #force_model=ForceModel.LUO,
+        force=sp.symbols("F_:3"),
+        force_model=ForceModel.LUO,
         compressible=True,
         psm_config=psm_config,
     )