From e40a8a544cff99ace85f5ab4a55765a4e9576c8b Mon Sep 17 00:00:00 2001 From: Zhi Date: Wed, 24 Jul 2024 21:44:01 -0400 Subject: [PATCH] update BC with dens_src and xmom_src --- pyro/compressible/BC.py | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) diff --git a/pyro/compressible/BC.py b/pyro/compressible/BC.py index 812ae2db2..8b841d9b5 100644 --- a/pyro/compressible/BC.py +++ b/pyro/compressible/BC.py @@ -53,7 +53,9 @@ def user(bc_name, bc_edge, variable, ccdata): # we will take the density to be constant, the velocity to # be outflow, and the pressure to be in HSE - if variable in ["density", "x-momentum", "y-momentum", "ymom_src", "E_src", "fuel", "ash"]: + if variable in ["density", "x-momentum", "y-momentum", + "dens_src", "xmom_src", "ymom_src", "E_src", + "fuel", "ash"]: v = ccdata.get_var(variable) j = myg.jlo-1 while j >= 0: @@ -99,7 +101,9 @@ def user(bc_name, bc_edge, variable, ccdata): # we will take the density to be constant, the velocity to # be outflow, and the pressure to be in HSE - if variable in ["density", "x-momentum", "y-momentum", "ymom_src", "E_src", "fuel", "ash"]: + if variable in ["density", "x-momentum", "y-momentum", + "dens_src", "xmom_src", "ymom_src", "E_src", + "fuel", "ash"]: v = ccdata.get_var(variable) for j in range(myg.jhi+1, myg.jhi+myg.ng+1): v[:, j] = v[:, myg.jhi]