diff --git a/src/EulerEq/rhs.jl b/src/EulerEq/rhs.jl index 13ba6d1f2f6d9a6a5669b9a955989084b041bfc9..d2c5808b311d3ea951a53a67fc25fa412ce51cac 100644 --- a/src/EulerEq/rhs.jl +++ b/src/EulerEq/rhs.jl @@ -235,6 +235,7 @@ function rhs!(env, P::Project2d, hrsc::HRSC.AbstractArtificialViscosity, mesh::M ## solve auxiliary equation: q + ∂x u = 0 broadcast_faces_2!(ldg_nflux, equation, mesh) + # TODO Need to revise Bdry conditions # broadcast_bdry_2!(ldg_bdryllf, equation, P.bdrycond, mesh) compute_rhs_weak_form!(ldg_rho_x, rho, 0, nflx_rho_x, mesh) compute_rhs_weak_form!(ldg_rho_y, 0, rho, nflx_rho_y, mesh) @@ -248,6 +249,7 @@ function rhs!(env, P::Project2d, hrsc::HRSC.AbstractArtificialViscosity, mesh::M ## compute rhs of regularized equation: ∂t u + ∂x f + ∂x μ q = 0 broadcast_volume_2!(flux, equation, mesh) broadcast_faces_2!(llf, equation, mesh) + # TODO Need to revise Bdry conditions # broadcast_bdry_2!(bdryllf, equation, mesh) broadcast_volume_2!(av_flux_2d, equation, mesh)