pass boundary_conditions

src/equations/serre_green_naghdi_1d.jl

@@ -314,19 +314,20 @@ end
 function rhs!(dq, q, t, mesh,
               equations::SerreGreenNaghdiEquations1D,
               initial_condition,
-              ::BoundaryConditionPeriodic,
+              boundary_conditions::BoundaryConditionPeriodic,
               source_terms::Nothing,
               solver, cache)
     if cache.D1 isa PeriodicUpwindOperators
-        rhs_sgn_upwind!(dq, q, equations, source_terms, cache, equations.bathymetry_type)
+        rhs_sgn_upwind!(dq, q, equations, source_terms, cache, equations.bathymetry_type, boundary_conditions)
     else
-        rhs_sgn_central!(dq, q, equations, source_terms, cache, equations.bathymetry_type)
+        rhs_sgn_central!(dq, q, equations, source_terms, cache, equations.bathymetry_type, boundary_conditions)
     end
 
     return nothing
 end
 
-function rhs_sgn_central!(dq, q, equations, source_terms, cache, ::BathymetryFlat)
+function rhs_sgn_central!(dq, q, equations, source_terms, cache, ::BathymetryFlat,
+                          boundary_conditions::BoundaryConditionPeriodic)
     # Unpack physical parameters and SBP operator `D1` as well as the
     # SBP operator in sparse matrix form `D1mat`
     g = gravity_constant(equations)
@@ -406,7 +407,8 @@ function rhs_sgn_central!(dq, q, equations, source_terms, cache, ::BathymetryFla
     return nothing
 end
 
-function rhs_sgn_upwind!(dq, q, equations, source_terms, cache, ::BathymetryFlat)
+function rhs_sgn_upwind!(dq, q, equations, source_terms, cache, ::BathymetryFlat,
+                         boundary_conditions::BoundaryConditionPeriodic)
     # Unpack physical parameters and SBP operator `D1` as well as the
     # SBP upwind operator in sparse matrix form `D1mat_minus`
     g = gravity_constant(equations)
@@ -492,7 +494,8 @@ function rhs_sgn_upwind!(dq, q, equations, source_terms, cache, ::BathymetryFlat
 end
 
 function rhs_sgn_central!(dq, q, equations, source_terms, cache,
-                          ::Union{BathymetryMildSlope, BathymetryVariable})
+                          ::Union{BathymetryMildSlope, BathymetryVariable},
+                          boundary_conditions::BoundaryConditionPeriodic)
     # Unpack physical parameters and SBP operator `D1` as well as the
     # SBP operator in sparse matrix form `D1mat`
     g = gravity_constant(equations)
@@ -597,7 +600,8 @@ function rhs_sgn_central!(dq, q, equations, source_terms, cache,
 end
 
 function rhs_sgn_upwind!(dq, q, equations, source_terms, cache,
-                         ::Union{BathymetryMildSlope, BathymetryVariable})
+                         ::Union{BathymetryMildSlope, BathymetryVariable},
+                         boundary_conditions::BoundaryConditionPeriodic)
     # Unpack physical parameters and SBP operator `D1` as well as the
     # SBP operator in sparse matrix form `D1mat`
     g = gravity_constant(equations)