From 66773e4f0378d14e84b084416e90fb1181518515 Mon Sep 17 00:00:00 2001
From: "Todd A. Oliver" <oliver@oden.utexas.edu>
Date: Tue, 9 May 2023 16:21:17 -0500
Subject: [PATCH] Make sure bdrN is initialized

Because of changes added in af1de26, bc types that previously didn't
use bdrN do now.  This mod ensures that we don't return out of
updateMean before it is initialized.
---
 src/inletBC.cpp  | 2 +-
 src/outletBC.cpp | 3 +--
 2 files changed, 2 insertions(+), 3 deletions(-)

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());