diff --git a/src/inletBC.cpp b/src/inletBC.cpp
index d6ac82f14..1290a1bb0 100644
--- a/src/inletBC.cpp
+++ b/src/inletBC.cpp
@@ -464,8 +464,8 @@ void InletBC::computeBdrPrimitiveStateForGradient(const int i, const Vector &pri
 
 
 void InletBC::updateMean(IntegrationRules *intRules, ParGridFunction *Up) {
-  if (inletType_ == SUB_DENS_VEL) return;
   bdrN = 0;
+  if (inletType_ == SUB_DENS_VEL) return;
 
 #ifdef _GPU_
   DGNonLinearForm::setToZero_gpu(bdrUp, bdrUp.Size());
diff --git a/src/outletBC.cpp b/src/outletBC.cpp
index 7eae260f2..3a904a622 100644
--- a/src/outletBC.cpp
+++ b/src/outletBC.cpp
@@ -472,9 +472,8 @@ void OutletBC::initBoundaryU(ParGridFunction *Up) {
 }
 
 void OutletBC::updateMean(IntegrationRules *intRules, ParGridFunction *Up) {
-  if (outletType_ == SUB_P) return;
-
   bdrN = 0;
+  if (outletType_ == SUB_P) return;
 
 #ifdef _GPU_
   DGNonLinearForm::setToZero_gpu(bdrUp, bdrUp.Size());