diff --git a/Source/hydro/Castro_ctu_hydro.cpp b/Source/hydro/Castro_ctu_hydro.cpp index b7db069e41..68c93bf857 100644 --- a/Source/hydro/Castro_ctu_hydro.cpp +++ b/Source/hydro/Castro_ctu_hydro.cpp @@ -259,16 +259,13 @@ Castro::construct_ctu_hydro_source(Real time, Real dt) // NOLINT(readability-co shk.resize(obx, 1); fab_size += shk.nBytes(); - Array4 const shk_arr = shk.array(); - // Multidimensional shock detection - // Used for the hybrid Riemann solver -#ifdef SHOCK_VAR - bool compute_shock = true; -#else - bool compute_shock = false; -#endif + // this is a local shock variable used only for the Riemann + // solver -- this will never be used to update the value in the + // conserved state. + + Array4 const shk_arr = shk.array(); // get the primitive variable hydro sources @@ -285,7 +282,7 @@ Castro::construct_ctu_hydro_source(Real time, Real dt) // NOLINT(readability-co hydro::src_to_prim(i, j, k, dt, U_old_arr, q_arr, old_src_arr, src_corr_arr, src_q_arr); }); - if (hybrid_riemann == 1 || compute_shock) { + if (hybrid_riemann == 1) { shock(obx, q_arr, old_src_arr, shk_arr); } else {