From b3722bb8d28a154011e585dcbb50616a081b5dd5 Mon Sep 17 00:00:00 2001
From: cgilet <cgilet@gmail.com>
Date: Sat, 22 Jul 2023 20:30:02 -0400
Subject: [PATCH] In MLEBABecLap::compGrad, add missing call to
 addInhomogNeumannFlux

For inhomogeneous Neumann BC, we convert the inhomogeneous part
into RHS and treat the BC as homogeneous Neumann during the solve.
Thus, the flux is zero there before inhomogeneous flux is added.
---
 Src/LinearSolvers/MLMG/AMReX_MLEBABecLap.cpp | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/Src/LinearSolvers/MLMG/AMReX_MLEBABecLap.cpp b/Src/LinearSolvers/MLMG/AMReX_MLEBABecLap.cpp
index 7d0b1416d81..170ed49f6d4 100644
--- a/Src/LinearSolvers/MLMG/AMReX_MLEBABecLap.cpp
+++ b/Src/LinearSolvers/MLMG/AMReX_MLEBABecLap.cpp
@@ -837,6 +837,8 @@ MLEBABecLap::compGrad (int amrlev, const Array<MultiFab*,AMREX_SPACEDIM>& grad,
             );
         }
     }
+
+    addInhomogNeumannFlux(amrlev, grad, sol, false);
 }
 
 void