From 65aa85db5cb3bbabc767d5dde4b106b7022a0f90 Mon Sep 17 00:00:00 2001
From: Michael B Kuhn <31661049+mbkuhn@users.noreply.github.com>
Date: Mon, 5 Aug 2024 15:30:11 -0600
Subject: [PATCH] Anelastic: multiply/divide all components of velocity (#1184)

---
 .../projection/incflo_apply_nodal_projection.cpp | 16 ++++++++++------
 1 file changed, 10 insertions(+), 6 deletions(-)

diff --git a/amr-wind/projection/incflo_apply_nodal_projection.cpp b/amr-wind/projection/incflo_apply_nodal_projection.cpp
index f4af85f3ea..e8bb944559 100644
--- a/amr-wind/projection/incflo_apply_nodal_projection.cpp
+++ b/amr-wind/projection/incflo_apply_nodal_projection.cpp
@@ -334,9 +334,11 @@ void incflo::ApplyProjection(
 
     if (is_anelastic) {
         for (int lev = 0; lev <= finest_level; ++lev) {
-            amrex::Multiply(
-                velocity(lev), (*ref_density)(lev), 0, 0, density[lev]->nComp(),
-                0);
+            for (int idim = 0; idim < velocity.num_comp(); ++idim) {
+                amrex::Multiply(
+                    velocity(lev), (*ref_density)(lev), 0, idim,
+                    density[lev]->nComp(), 0);
+            }
         }
     }
 
@@ -435,9 +437,11 @@ void incflo::ApplyProjection(
 
     if (is_anelastic) {
         for (int lev = 0; lev <= finest_level; ++lev) {
-            amrex::Divide(
-                velocity(lev), (*ref_density)(lev), 0, 0, density[lev]->nComp(),
-                0);
+            for (int idim = 0; idim < velocity.num_comp(); ++idim) {
+                amrex::Divide(
+                    velocity(lev), (*ref_density)(lev), 0, idim,
+                    density[lev]->nComp(), 0);
+            }
         }
     }