spherical accretion:make sure K, Gamma agree with loaded data

src/GRHD/initialdata.jl

@@ -55,8 +55,8 @@ function initialdata_equation!(::Val{:bondi_accretion}, env, P::Project{:spheric
     (; Γ, K, ρc, rc, M, r, ρ, ϵ, p, vr)
   end
 
-  # @show data.K - equation.eos.K
-  # @show data.Γ - equation.eos.Gamma
+  @assert data.K ≈ equation.eos.K
+  @assert data.Γ ≈ equation.eos.Gamma
 
   @unpack D, Sr, τ = get_dynamic_variables(mesh.cache)
   @unpack ρ, p, ϵ, vr, r, max_v,