Skip to content

Commit

Permalink
more revert
Browse files Browse the repository at this point in the history
  • Loading branch information
zhichen3 committed Apr 27, 2023
1 parent 047ad09 commit 8074edb
Showing 1 changed file with 7 additions and 6 deletions.
13 changes: 7 additions & 6 deletions pyro/incompressible/incomp_interface.py
Original file line number Diff line number Diff line change
Expand Up @@ -38,6 +38,13 @@ def mac_vels(grid, dt,
ldelta_ux, ldelta_vx,
ldelta_uy, ldelta_vy)

# apply transverse terms on both x- and y- interfaces
u_xl, u_xr, u_yl, u_yr, v_xl, v_xr, v_yl, v_yr = burgers_interface.apply_transverse_corrections(grid, dt,
u_xl, u_xr,
u_yl, u_yr,
v_xl, v_xr,
v_yl, v_yr)

# apply pressure gradient correction terms to the interface state
u_xl, u_xr, u_yl, u_yr, v_xl, v_xr, v_yl, v_yr = apply_gradp_corrections(dt,
u_xl, u_xr,
Expand All @@ -46,12 +53,6 @@ def mac_vels(grid, dt,
v_yl, v_yr,
gradp_x, gradp_y)

# apply transverse terms on both x- and y- interfaces
u_xl, u_xr, u_yl, u_yr, v_xl, v_xr, v_yl, v_yr = burgers_interface.apply_transverse_corrections(grid, dt,
u_xl, u_xr,
u_yl, u_yr,
v_xl, v_xr,
v_yl, v_yr)
# Riemann problem -- this follows Burger's equation. We don't use
# any input velocity for the upwinding. Also, we only care about
# the normal states here (u on x and v on y)
Expand Down

0 comments on commit 8074edb

Please sign in to comment.