From ec0f3b05fde492da99d6efe9e5eb71a1b85bee73 Mon Sep 17 00:00:00 2001 From: Daniel_Doehring Date: Thu, 28 Mar 2024 15:16:02 +0100 Subject: [PATCH] Comments --- examples/tree_1d_dgsem/elixir_linearizedeuler_gauss_wall.jl | 3 +++ 1 file changed, 3 insertions(+) diff --git a/examples/tree_1d_dgsem/elixir_linearizedeuler_gauss_wall.jl b/examples/tree_1d_dgsem/elixir_linearizedeuler_gauss_wall.jl index c087989e534..0884249559a 100644 --- a/examples/tree_1d_dgsem/elixir_linearizedeuler_gauss_wall.jl +++ b/examples/tree_1d_dgsem/elixir_linearizedeuler_gauss_wall.jl @@ -18,6 +18,9 @@ mesh = TreeMesh(coordinates_min, coordinates_max, n_cells_max = 100_000, periodicity = false) +# Initialize density and pressure perturbation with a Gaussian bump +# that is advected to left with v - c and to the right with v + c. +# Correspondigly, the bump splits in half. function initial_condition_gauss_wall(x, t, equations::LinearizedEulerEquations1D) v1_prime = 0.0 rho_prime = p_prime = 2 * exp(-(x[1] - 45)^2 / 25)