fix bug in capillary flux computation

OPM · Aug 29, 2023 · 7394e7f · 7394e7f
1 parent f574140
commit 7394e7f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/models/ColorModel.cpp b/models/ColorModel.cpp
@@ -521,7 +521,7 @@ void ScaLBL_ColorModel::Initialize() {
         double IFT = 6.0 * alpha;
         double CrossSectionalArea =
             (double)(nprocx * (Nx - 2) * nprocy * (Ny - 2));
-        flux = Mask->Porosity() * CrossSectionalArea * (Ny - 2) * IFT *
+        flux = Mask->Porosity() * CrossSectionalArea * IFT *
                capillary_number / MuB;
         if (rank == 0)
             printf("   flux=%f \n", flux);