add grhd example par files

examples/grhd_bondi/dg.toml

+[EquationOfState]
+eos = "polytrope"
+polytrope_k = 0.02143872868864732
+polytrope_gamma = "$(4/3)"
+
+[GRHD]
+bc = "from_id"
+id = "bondi_accretion"
+id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"bondi_accretion.h5\"))"
+
+[Mesh]
+range = [ 1.8, 10.0 ]
+n = 4
+k = 35
+basis = "lgl"
+periodic = false
+
+[Output]
+every_iteration = 1
+# every_dt  = 0.1
+variables = [ "D", "Sr", "τ", "p", "ρ", "ϵ", "vr", "rhs_D", "rhs_Sr", "rhs_τ", "src_D", "src_Sr", "src_τ" ]
+# aligned_ts      = "$(collect(range(0.01,0.6,step=0.01)))"
+enable1d  = true
+
+[Evolution]
+cfl = 0.4
+tend = 20.0

examples/grhd_bondi/fv.toml

+[EquationOfState]
+eos = "polytrope"
+polytrope_k = 0.02143872868864732
+polytrope_gamma = "$(4/3)"
+
+[GRHD]
+bc = "from_id"
+id = "bondi_accretion"
+id_filename = "bondi_accretion.h5"
+
+[Mesh]
+range = [ 1.8, 10.0 ]
+k = 4048
+scheme = "FV"
+periodic = false
+
+[Output]
+every_iteration = 1
+variables       = [ "D", "Sr", "τ", "p", "ρ", "ϵ", "vr", "rhs_D", "rhs_Sr", "rhs_τ", "src_D", "src_Sr", "src_τ" ]
+# aligned_ts      = "$(collect(range(0.01,0.6,step=0.01)))"
+enable1d        = true
+
+[Evolution]
+cfl = 1.0
+tend = 20.0

examples/grhd_tov/dg.toml

+[EquationOfState]
+eos = "polytrope"
+polytrope_k = 100.0
+polytrope_gamma = 2.0
+
+[GRHD]
+bc = "from_id"
+id = "tov"
+id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))"
+atm_factor = 1e-8
+
+[Mesh]
+range = [ 0.0, 20.0 ]
+n = 6
+k = 125
+basis = "lgl"
+periodic = false
+
+[Output]
+# every_iteration = 1
+every_dt  = 1.0
+variables = [ #"D", "Sr", "τ",
+              #"p", "ρ", "ϵ", "vr",
+              "Sr", "vr",
+              "rD", "rSr", "rτ",
+              "rhs_rD", "rhs_rSr", "rhs_rτ",
+              "src_rD", "src_rSr", "src_rτ",
+              "flr_rD", "flr_rSr", "flr_rτ",
+              #"A", "∂Ar", "B", "∂Br"
+              ]
+# aligned_ts      = "$(collect(range(0.01,0.6,step=0.01)))"
+enable1d  = true
+
+[Evolution]
+cfl = 0.8
+tend = 40.0