Don't add rho pert if anelastic. (#2106)

erf-model · Feb 6, 2025 · 93fb408 · 93fb408
1 parent f2f94f3
commit 93fb408
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/Source/Initialization/ERF_InitCustom.cpp b/Source/Initialization/ERF_InitCustom.cpp
@@ -76,7 +76,9 @@ ERF::init_custom (int lev)
     } //mfi
 
     // Add problem-specific perturbation to background flow
-    MultiFab::Add(lev_new[Vars::cons], cons_pert, Rho_comp,      Rho_comp,             1, cons_pert.nGrow());
+    if (!solverChoice.anelastic[lev]) {
+        MultiFab::Add(lev_new[Vars::cons], cons_pert, Rho_comp,      Rho_comp,             1, cons_pert.nGrow());
+    }
     MultiFab::Add(lev_new[Vars::cons], cons_pert, RhoTheta_comp, RhoTheta_comp,        1, cons_pert.nGrow());
     MultiFab::Add(lev_new[Vars::cons], cons_pert, RhoScalar_comp,RhoScalar_comp,NSCALARS, cons_pert.nGrow());