From 740ae4dc7e479128c270f4e02789fde3e8c4124d Mon Sep 17 00:00:00 2001
From: cgilet <cgilet@gmail.com>
Date: Wed, 22 Oct 2025 12:27:03 -0400
Subject: [PATCH] Ensure compute_divtau sees correct component of velocity.
 Only correct for vel_x perviously.

---
 src/diffusion/DiffusionScalarOp.cpp | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/src/diffusion/DiffusionScalarOp.cpp b/src/diffusion/DiffusionScalarOp.cpp
index b277c1d0..6a7a39df 100644
--- a/src/diffusion/DiffusionScalarOp.cpp
+++ b/src/diffusion/DiffusionScalarOp.cpp
@@ -702,8 +702,6 @@ void DiffusionScalarOp::compute_divtau (Vector<MultiFab*> const& a_divtau,
         m_reg_vel_apply_op->setScalars(0.0, -1.0);
 
         int eta_comp = 0;
-        Vector<MultiFab> divtau_single;
-        Vector<MultiFab>    vel_single;
 
         for (int lev = 0; lev <= finest_level; ++lev)
         {
@@ -719,6 +717,9 @@ void DiffusionScalarOp::compute_divtau (Vector<MultiFab*> const& a_divtau,
             m_reg_vel_apply_op->setDomainBC(m_incflo->get_diffuse_velocity_bc(Orientation::low ,comp),
                                         m_incflo->get_diffuse_velocity_bc(Orientation::high,comp));
 
+            Vector<MultiFab> divtau_single;
+            Vector<MultiFab>    vel_single;
+
             for (int lev = 0; lev <= finest_level; ++lev) {
                 divtau_single.emplace_back(*a_divtau[lev],amrex::make_alias,comp,1);
                    vel_single.emplace_back(      vel[lev],amrex::make_alias,comp,1);