From dd122aab19ceef17e5fb98808f2521b06d558f4b Mon Sep 17 00:00:00 2001
From: "Aaron M. Lattanzi" <103702284+AMLattanzi@users.noreply.github.com>
Date: Mon, 23 Sep 2024 11:05:41 -0700
Subject: [PATCH] Fix WPS regtest. (#1828)

---
 Source/ERF_IndexDefines.H              |  4 ++++
 Source/Initialization/ERF_init_bcs.cpp | 10 ++++++++++
 2 files changed, 14 insertions(+)

diff --git a/Source/ERF_IndexDefines.H b/Source/ERF_IndexDefines.H
index 322a2b884..54d20db69 100644
--- a/Source/ERF_IndexDefines.H
+++ b/Source/ERF_IndexDefines.H
@@ -70,6 +70,10 @@ namespace BCVars {
         RhoScalar_bc_comp,
         RhoQ1_bc_comp,
         RhoQ2_bc_comp,
+        RhoQ3_bc_comp,
+        RhoQ4_bc_comp,
+        RhoQ5_bc_comp,
+        RhoQ6_bc_comp,
         xvel_bc = NBCVAR_max,
         yvel_bc = NBCVAR_max+1,
         zvel_bc = NBCVAR_max+2,
diff --git a/Source/Initialization/ERF_init_bcs.cpp b/Source/Initialization/ERF_init_bcs.cpp
index 569eb5e05..428a72be0 100644
--- a/Source/Initialization/ERF_init_bcs.cpp
+++ b/Source/Initialization/ERF_init_bcs.cpp
@@ -34,6 +34,10 @@ void ERF::init_bcs ()
         m_bc_extdir_vals[BCVars::RhoScalar_bc_comp][ori] = 0.0;
         m_bc_extdir_vals[BCVars::RhoQ1_bc_comp][ori]     = 0.0;
         m_bc_extdir_vals[BCVars::RhoQ2_bc_comp][ori]     = 0.0;
+        m_bc_extdir_vals[BCVars::RhoQ3_bc_comp][ori]     = 0.0;
+        m_bc_extdir_vals[BCVars::RhoQ4_bc_comp][ori]     = 0.0;
+        m_bc_extdir_vals[BCVars::RhoQ5_bc_comp][ori]     = 0.0;
+        m_bc_extdir_vals[BCVars::RhoQ6_bc_comp][ori]     = 0.0;
 
         m_bc_extdir_vals[BCVars::xvel_bc][ori] = 0.0; // default
         m_bc_extdir_vals[BCVars::yvel_bc][ori] = 0.0;
@@ -42,11 +46,17 @@ void ERF::init_bcs ()
         // These are simply defaults for Neumann gradients -- they should be over-written below
         m_bc_neumann_vals[BCVars::Rho_bc_comp][ori]       = 0.0;
         m_bc_neumann_vals[BCVars::RhoTheta_bc_comp][ori]  = 0.0;
+
         m_bc_neumann_vals[BCVars::RhoKE_bc_comp][ori]     = 0.0;
         m_bc_neumann_vals[BCVars::RhoQKE_bc_comp][ori]    = 0.0;
         m_bc_neumann_vals[BCVars::RhoScalar_bc_comp][ori] = 0.0;
         m_bc_neumann_vals[BCVars::RhoQ1_bc_comp][ori]     = 0.0;
         m_bc_neumann_vals[BCVars::RhoQ2_bc_comp][ori]     = 0.0;
+        m_bc_neumann_vals[BCVars::RhoQ3_bc_comp][ori]     = 0.0;
+        m_bc_neumann_vals[BCVars::RhoQ4_bc_comp][ori]     = 0.0;
+        m_bc_neumann_vals[BCVars::RhoQ5_bc_comp][ori]     = 0.0;
+        m_bc_neumann_vals[BCVars::RhoQ6_bc_comp][ori]     = 0.0;
+
         m_bc_neumann_vals[BCVars::xvel_bc][ori] = 0.0;
         m_bc_neumann_vals[BCVars::yvel_bc][ori] = 0.0;
         m_bc_neumann_vals[BCVars::zvel_bc][ori] = 0.0;