fix: set atm density in bondi_accretion initial data

src/GRHD/initialdata.jl

@@ -72,6 +72,10 @@ function initialdata_equation!(::Val{:bondi_accretion}, env, P::Project{:spheric
   ϵ   .= interpolate(r, data.r, data.ϵ)
   vru  = interpolate(r, data.r, data.vr)
 
+  equation.atm_density = maximum(ρ) * equation.atm_factor
+  equation.atm_factor = prms["atm_factor"]
+  equation.atm_threshold = prms["atm_threshold_factor"]
+
   # deps, ideps = newsignature(grhd_sources)
   # deps  = renamevars(deps, Dict( (1,2,3,4) .=> (:r,:θ,:ϕ,:t) ))
   # ideps = renamevars(ideps, Dict( (1,2,3,4) .=> (:r,:θ,:ϕ,:t) ))