diff --git a/examples/elixir_potential_euler_dry_bubble.jl b/examples/elixir_potential_euler_dry_bubble.jl
index 288a797..1e65276 100644
--- a/examples/elixir_potential_euler_dry_bubble.jl
+++ b/examples/elixir_potential_euler_dry_bubble.jl
@@ -71,7 +71,7 @@ mesh = StructuredMesh(cells_per_dimension, coordinates_min, coordinates_max,
 initial_condition = initial_condition_warm_bubble
 
 semi = SemidiscretizationHyperbolic(mesh, equations, initial_condition, solver,
-                                    source_terms = source_terms_bubble,
+                                    source_terms = source_terms_gravity,
                                     boundary_conditions = boundary_conditions)
 
 ###############################################################################
diff --git a/src/equations/compressible_euler_potential_temperature_2d.jl b/src/equations/compressible_euler_potential_temperature_2d.jl
index 8ddea9f..f318072 100644
--- a/src/equations/compressible_euler_potential_temperature_2d.jl
+++ b/src/equations/compressible_euler_potential_temperature_2d.jl
@@ -132,7 +132,7 @@ end
     return boundary_flux
 end
 
-@inline function source_terms_bubble(u, x, t,
+@inline function source_terms_gravity(u, x, t,
                                      equations::CompressibleEulerPotentialTemperatureEquations2D)
     rho, _, _, _ = u
     return SVector(zero(eltype(u)), zero(eltype(u)), -equations.g * rho,