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;