fix some comments

examples/tree_2d_dgsem/elixir_shallowwater_multilayer_dam_break_dry.jl

@@ -29,7 +29,7 @@ function initial_condition_dam_break(x, t, equations::ShallowWaterMultiLayerEqua
     # stays positive. The system would not be stable for h set to a hard 0 due to division by h in
     # the computation of velocity, e.g., (h v) / h. Therefore, a small dry state threshold
     # with a default value of 5*eps() ≈ 1e-15 in double precision, is set in the constructor above
-    # for the ShallowWaterMultiLayerEquations1D and added to the initial condition if h = 0.
+    # for the ShallowWaterMultiLayerEquations2D and added to the initial condition if h = 0.
     # This default value can be changed within the constructor call depending on the simulation setup.
     for i in reverse(eachlayer(equations))
         if i == nlayers(equations)

examples/unstructured_2d_dgsem/elixir_shallowwater_multilayer_dam_break_dry.jl

@@ -124,7 +124,7 @@ function initial_condition_discontinuous_dam_break(x, t, element_id,
     # stays positive. The system would not be stable for h set to a hard 0 due to division by h in
     # the computation of velocity, e.g., (h v) / h. Therefore, a small dry state threshold
     # with a default value of 5*eps() ≈ 1e-15 in double precision, is set in the constructor above
-    # for the ShallowWaterMultiLayerEquations1D and added to the initial condition if h = 0.
+    # for the ShallowWaterMultiLayerEquations2D and added to the initial condition if h = 0.
     # This default value can be changed within the constructor call depending on the simulation setup.
     for i in reverse(eachlayer(equations))
         if i == nlayers(equations)

examples/unstructured_2d_dgsem/elixir_shallowwater_multilayer_well_balanced_wet_dry.jl

@@ -30,7 +30,7 @@ function initial_condition_well_balancedness(x, t,
     # stays positive. The system would not be stable for h set to a hard 0 due to division by h in
     # the computation of velocity, e.g., (h v) / h. Therefore, a small dry state threshold
     # with a default value of 5*eps() ≈ 1e-15 in double precision, is set in the constructor above
-    # for the ShallowWaterMultiLayerEquations1D and added to the initial condition if h = 0.
+    # for the ShallowWaterMultiLayerEquations2D and added to the initial condition if h = 0.
     # This default value can be changed within the constructor call depending on the simulation setup.
     for i in reverse(eachlayer(equations))
         if i == nlayers(equations)

src/equations/numerical_fluxes.jl

@@ -36,13 +36,13 @@ const flux_hll_chen_noelle = FluxHLL(min_max_speed_chen_noelle)
                                                                 normal_direction_ll,
                                                                 normal_direction_average,
                                                                 equations::Trixi.AbstractEquations)
-    @unpack numerical_flux, hydrostatic_reconstruction = numflux
+    @unpack nonconservative_flux, hydrostatic_reconstruction = numflux
 
     # Create the reconstructed left/right solution states in conservative form
     u_ll_star, u_rr_star = hydrostatic_reconstruction(u_ll, u_rr, equations)
 
-    # Use the reconstructed states to compute the numerical surface flux
-    return numerical_flux(u_ll_star, u_rr_star, normal_direction_ll,
+    # Use the reconstructed states to compute the nonconservative surface flux
+    return nonconservative_flux(u_ll_star, u_rr_star, normal_direction_ll,
                           normal_direction_average,
                           equations)
 end