diff --git a/src/equations/compressible_euler_3d.jl b/src/equations/compressible_euler_3d.jl
index fc65f0e42c1..80de60351df 100644
--- a/src/equations/compressible_euler_3d.jl
+++ b/src/equations/compressible_euler_3d.jl
@@ -154,38 +154,6 @@ function initial_condition_weak_blast_wave(x, t, equations::CompressibleEulerEqu
 end
 
 
-# TODO: ICs, remove this one?
-"""
-    initial_condition_sedov_blast_wave(x, t, equations::CompressibleEulerEquations3D)
-
-The Sedov blast wave setup based on Flash
-- http://flash.uchicago.edu/site/flashcode/user_support/flash_ug_devel/node184.html#SECTION010114000000000000000
-"""
-function initial_condition_sedov_blast_wave(x, t, equations::CompressibleEulerEquations3D)
-  # Calculate radius as distance from origin
-  r = sqrt(x[1]^2 + x[2]^2 + x[3]^2)
-
-  # Setup based on http://flash.uchicago.edu/site/flashcode/user_support/flash4_ug_4p62/node184.html#SECTION010114000000000000000
-  r0 = 0.25 # = 4.0 * smallest dx (for domain length=8 and max-ref=7)
-  E = 1.0
-  p_inner   = (equations.gamma - 1) * E / (4/3 * pi * r0^3)
-  p_ambient = 1e-5 # = true Sedov setup
-
-  # Calculate primitive variables
-  rho = 1.0
-  v1 = 0.0
-  v2 = 0.0
-  v3 = 0.0
-
-  # use a logistic function to tranfer pressure value smoothly
-  k  = -50.0 # sharpness of transfer
-  logistic_function_p = p_inner/(1.0 + exp(-k*(r - r0)))
-  p = max(logistic_function_p, p_ambient)
-
-  return prim2cons(SVector(rho, v1, v2, v3, p), equations)
-end
-
-
 """
     initial_condition_blob(x, t, equations::CompressibleEulerEquations3D)
 
diff --git a/test/test_tree_3d_euler.jl b/test/test_tree_3d_euler.jl
index fdd912c14ee..37761f9c9f8 100644
--- a/test/test_tree_3d_euler.jl
+++ b/test/test_tree_3d_euler.jl
@@ -60,15 +60,6 @@ EXAMPLES_DIR = joinpath(pathof(Trixi) |> dirname |> dirname, "examples", "tree_3
       linf = [0.3954458125573179, 0.26876916180359345, 0.26876916180359345, 0.26933123042178553, 1.3724137121660251])
   end
 
-  # TODO: ICs, remove this one?
-  @trixi_testset "elixir_euler_shockcapturing.jl with initial_condition_sedov_blast_wave" begin
-    # OBS! This setup does not run longer but crashes (also the parameters do not make sense) -> only for testing the IC!
-    @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_euler_shockcapturing.jl"),
-      l2   = [0.03627061709701234, 0.05178777768691398, 0.05178777768691404, 0.051787777686913984, 0.2308467334945545],
-      linf = [0.9542995612960261, 1.2654683412445416, 1.2654683412445418, 1.2654683412445418, 12.805752164787227],
-      initial_condition=Trixi.initial_condition_sedov_blast_wave, cfl=0.7, alpha_max=1.0, tspan=(0.0, 0.1))
-  end
-
   @trixi_testset "elixir_euler_shockcapturing.jl with initial_condition_sedov_self_gravity" begin
     # OBS! This setup does not run longer but crashes (also the parameters do not make sense) -> only for testing the IC!
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_euler_shockcapturing.jl"),