diff --git a/applications/solvers/dfLowMachFoam/YEqn.H b/applications/solvers/dfLowMachFoam/YEqn.H
index 173d0e4d..8e0befd6 100644
--- a/applications/solvers/dfLowMachFoam/YEqn.H
+++ b/applications/solvers/dfLowMachFoam/YEqn.H
@@ -21,7 +21,10 @@ forAll(Y, i)
 }
 const surfaceScalarField phiUc = linearInterpolate(sumYDiffError) & mesh.Sf();
 
-MPI_Barrier(PstreamGlobals::MPI_COMM_FOAM);
+//MPI_Barrier(PstreamGlobals::MPI_COMM_FOAM);
+label flag_mpi_init;
+MPI_Initialized(&flag_mpi_init);
+if(flag_mpi_init) MPI_Barrier(PstreamGlobals::MPI_COMM_FOAM);
 end = std::clock();
 time_monitor_corrDiff += double(end - start) / double(CLOCKS_PER_SEC);
 
@@ -30,8 +33,8 @@ time_monitor_corrDiff += double(end - start) / double(CLOCKS_PER_SEC);
     {
         std::chrono::steady_clock::time_point start = std::chrono::steady_clock::now();
         combustion->correct();
-        label flag_mpi_init;
-        MPI_Initialized(&flag_mpi_init);
+        //label flag_mpi_init;
+        //MPI_Initialized(&flag_mpi_init);
         if(flag_mpi_init) MPI_Barrier(PstreamGlobals::MPI_COMM_FOAM);
         std::chrono::steady_clock::time_point stop = std::chrono::steady_clock::now();
         std::chrono::duration<double> processingTime = std::chrono::duration_cast<std::chrono::duration<double>>(stop - start);