add formulation parameter

src/GRHD/GRHD.jl

@@ -76,6 +76,14 @@ dg1d.@parameters GRHD begin
   id_atmosphere_density = 1e-11
   @check id_atmosphere_density > 0.0
 
+  """
+  Choice of formulation of evolution equations. Available options:
+  - `spherical1d` ... singular at `r=0`
+  - `rescaled_spherical1d` ... regular at `r=0`
+  """
+  formulation = "spherical1d"
+  @check formulation in [ "spherical1d", "rescaled_spherical1d" ]
+
 end
 
 

src/GRHD/setup.jl

@@ -7,14 +7,8 @@ function Project(env::Environment, prms)
   hrsc = HRSC.make_HRSC(env.mesh, prms["HRSC"])
 
   fixedprms = (;)
-  formulation = if prms["GRHD"]["id"] == "bondi_accretion"
-    :spherical1d
-  elseif prms["GRHD"]["id"] == "tov"
-    :rescaled_spherical1d
-  else
-    TODO()
-  end
-  P = Project(formulation, equation, hrsc, tci, fixedprms)
+  formulation = prms["GRHD"]["formulation"]
+  P = Project(Symbol(formulation), equation, hrsc, tci, fixedprms)
 
   # register variables
   @unpack cache = env