Include ghost cells when multiplying or dividing rho0 for anelastic (E…

…xawind#1508)
erf-model · Feb 27, 2025 · 23b314b · 23b314b
1 parent f42c0cf
commit 23b314b
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/amr-wind/projection/incflo_apply_nodal_projection.cpp b/amr-wind/projection/incflo_apply_nodal_projection.cpp
@@ -333,7 +333,7 @@ void incflo::ApplyProjection(
             for (int idim = 0; idim < velocity.num_comp(); ++idim) {
                 amrex::Multiply(
                     velocity(lev), (*ref_density)(lev), 0, idim,
-                    density[lev]->nComp(), 0);
+                    density[lev]->nComp(), 1);
             }
         }
     }
@@ -405,7 +405,7 @@ void incflo::ApplyProjection(
             for (int idim = 0; idim < velocity.num_comp(); ++idim) {
                 amrex::Divide(
                     velocity(lev), (*ref_density)(lev), 0, idim,
-                    density[lev]->nComp(), 0);
+                    density[lev]->nComp(), 1);
             }
         }
     }