wrap SphericalAccretion solver into module; recompute bondi_accretion.h5 with less radial points

initialdata/SphericalAccretion.jl

+module SphericalAccretion
+
+
+
 using CairoMakie
 using CairoMakie.FileIO
 using DelimitedFiles
@@ -19,7 +23,7 @@ Alternatively, one could also add an adaptive stepper.
 Note: We implicitly assume θ = π/2, because then sin(θ) = 1, although,
 the sines should cancel out everywhere.
 """
-function solve(; rmin, rmax, rc, ρc, Γ, M=1.0, N_rs=10000,
+function solve(; rmin, rmax, rc, ρc, Γ, M=1.0, N_rs=1000,
                  urc_sign=1)
 
   urc = urc_sign*sqrt(M /(2 * rc))
@@ -336,3 +340,7 @@ function plot_compare_papadopoulus_1999()
 
   return fig
 end
+
+
+
+end # module