diff --git a/examples/sod_shock_tube/fv.toml b/examples/sod_shock_tube/fv.toml
new file mode 100644
index 0000000000000000000000000000000000000000..e9880d8a41397700e244a22ad09e683c8a0956b0
--- /dev/null
+++ b/examples/sod_shock_tube/fv.toml
@@ -0,0 +1,23 @@
+[EquationOfState]
+idealgas_gamma = 1.4
+eos = "idealgas"
+
+[EulerEq]
+id = "sod_shock_tube"
+bc = "from_id"
+
+[Evolution]
+cfl = 1.0
+tend = 0.5
+
+[Mesh]
+periodic = [false]
+range = [-1.0, 1.0]
+k = 1536
+scheme = "FV"
+
+[Output]
+every_iteration = 0
+every_dt        = 0.0
+aligned_ts      = "$(collect(range(0.01,0.5,step=0.01)))"
+variables       = [ "rho", "q", "E" ]
diff --git a/examples/sod_shock_tube/fv3.toml b/examples/sod_shock_tube/fv3.toml
new file mode 100644
index 0000000000000000000000000000000000000000..5fa047809259994b6763466b0f9a7cf3848ed95f
--- /dev/null
+++ b/examples/sod_shock_tube/fv3.toml
@@ -0,0 +1,23 @@
+[EquationOfState]
+idealgas_gamma = 1.4
+eos = "idealgas"
+
+[EulerEq]
+id = "sod_shock_tube"
+bc = "from_id"
+
+[Evolution]
+cfl = 1.0
+tend = 0.5
+
+[Mesh]
+periodic = [false]
+range = [-1.0, 1.0]
+k = 768
+scheme = "FV"
+
+[Output]
+every_iteration = 0
+every_dt        = 0.0
+aligned_ts      = "$(collect(range(0.01,0.5,step=0.01)))"
+variables       = [ "rho", "q", "E" ]
diff --git a/examples/sod_shock_tube/fv4.toml b/examples/sod_shock_tube/fv4.toml
new file mode 100644
index 0000000000000000000000000000000000000000..bc12f87c94f53e35ca1a0bce139866727e3daffa
--- /dev/null
+++ b/examples/sod_shock_tube/fv4.toml
@@ -0,0 +1,23 @@
+[EquationOfState]
+idealgas_gamma = 1.4
+eos = "idealgas"
+
+[EulerEq]
+id = "sod_shock_tube"
+bc = "from_id"
+
+[Evolution]
+cfl = 1.0
+tend = 0.5
+
+[Mesh]
+periodic = [false]
+range = [-1.0, 1.0]
+k = 1024
+scheme = "FV"
+
+[Output]
+every_iteration = 0
+every_dt        = 0.0
+aligned_ts      = "$(collect(range(0.01,0.5,step=0.01)))"
+variables       = [ "rho", "q", "E" ]
diff --git a/examples/sod_shock_tube/fv5.toml b/examples/sod_shock_tube/fv5.toml
new file mode 100644
index 0000000000000000000000000000000000000000..f75e9fff714a167396934c1d78c460732aa4072e
--- /dev/null
+++ b/examples/sod_shock_tube/fv5.toml
@@ -0,0 +1,23 @@
+[EquationOfState]
+idealgas_gamma = 1.4
+eos = "idealgas"
+
+[EulerEq]
+id = "sod_shock_tube"
+bc = "from_id"
+
+[Evolution]
+cfl = 1.0
+tend = 0.5
+
+[Mesh]
+periodic = [false]
+range = [-1.0, 1.0]
+k = 1280
+scheme = "FV"
+
+[Output]
+every_iteration = 0
+every_dt        = 0.0
+aligned_ts      = "$(collect(range(0.01,0.5,step=0.01)))"
+variables       = [ "rho", "q", "E" ]
diff --git a/examples/sod_shock_tube/fv6.toml b/examples/sod_shock_tube/fv6.toml
new file mode 100644
index 0000000000000000000000000000000000000000..e9880d8a41397700e244a22ad09e683c8a0956b0
--- /dev/null
+++ b/examples/sod_shock_tube/fv6.toml
@@ -0,0 +1,23 @@
+[EquationOfState]
+idealgas_gamma = 1.4
+eos = "idealgas"
+
+[EulerEq]
+id = "sod_shock_tube"
+bc = "from_id"
+
+[Evolution]
+cfl = 1.0
+tend = 0.5
+
+[Mesh]
+periodic = [false]
+range = [-1.0, 1.0]
+k = 1536
+scheme = "FV"
+
+[Output]
+every_iteration = 0
+every_dt        = 0.0
+aligned_ts      = "$(collect(range(0.01,0.5,step=0.01)))"
+variables       = [ "rho", "q", "E" ]
diff --git a/examples/sod_shock_tube/fv7.toml b/examples/sod_shock_tube/fv7.toml
new file mode 100644
index 0000000000000000000000000000000000000000..0e968a8dc7ad89576e4225d8d13500bc56244bc0
--- /dev/null
+++ b/examples/sod_shock_tube/fv7.toml
@@ -0,0 +1,23 @@
+[EquationOfState]
+idealgas_gamma = 1.4
+eos = "idealgas"
+
+[EulerEq]
+id = "sod_shock_tube"
+bc = "from_id"
+
+[Evolution]
+cfl = 1.0
+tend = 0.5
+
+[Mesh]
+periodic = [false]
+range = [-1.0, 1.0]
+k = 1792
+scheme = "FV"
+
+[Output]
+every_iteration = 0
+every_dt        = 0.0
+aligned_ts      = "$(collect(range(0.01,0.5,step=0.01)))"
+variables       = [ "rho", "q", "E" ]
diff --git a/examples/sod_shock_tube/mono_n3.toml b/examples/sod_shock_tube/mono_n3.toml
new file mode 100644
index 0000000000000000000000000000000000000000..cac7d55ec4643646c6d3d448f270aa8b1cac3d94
--- /dev/null
+++ b/examples/sod_shock_tube/mono_n3.toml
@@ -0,0 +1,33 @@
+[EquationOfState]
+idealgas_gamma = 1.4
+eos = "idealgas"
+
+[EulerEq]
+id = "sod_shock_tube"
+bc = "from_id"
+av_regularization = "mono"
+
+[Mesh]
+periodic = [ false ]
+range = [ -1.0, 1.0 ]
+cfl = 0.75
+n = 3
+k = 256
+basis = "lgl"
+
+[Output]
+every_iteration = 0
+every_dt        = 0.0
+aligned_ts      = "$(collect(range(0.01,0.5,step=0.01)))"
+variables       = [ "rho", "q", "E", "ldg_rho", "ldg_q", "ldg_E", "smoothed_mu" ]
+
+[Evolution]
+tend = 0.5
+
+[Log]
+progress_stdout = true
+
+[HRSC]
+method = "av"
+av_method = "entropy"
+entropy_ce = 5.0
diff --git a/examples/sod_shock_tube/mono_n4.toml b/examples/sod_shock_tube/mono_n4.toml
new file mode 100644
index 0000000000000000000000000000000000000000..756994abfe13df0bd5e030225c3f778311e4c3a3
--- /dev/null
+++ b/examples/sod_shock_tube/mono_n4.toml
@@ -0,0 +1,33 @@
+[EquationOfState]
+idealgas_gamma = 1.4
+eos = "idealgas"
+
+[EulerEq]
+id = "sod_shock_tube"
+bc = "from_id"
+av_regularization = "mono"
+
+[Mesh]
+periodic = [ false ]
+range = [ -1.0, 1.0 ]
+cfl = 0.75
+n = 4
+k = 256
+basis = "lgl"
+
+[Output]
+every_iteration = 0
+every_dt        = 0.0
+aligned_ts      = "$(collect(range(0.01,0.5,step=0.01)))"
+variables       = [ "rho", "q", "E", "ldg_rho", "ldg_q", "ldg_E", "smoothed_mu" ]
+
+[Evolution]
+tend = 0.5
+
+[Log]
+progress_stdout = true
+
+[HRSC]
+method = "av"
+av_method = "entropy"
+entropy_ce = 5.0
diff --git a/examples/sod_shock_tube/mono_n5.toml b/examples/sod_shock_tube/mono_n5.toml
new file mode 100644
index 0000000000000000000000000000000000000000..0a77c7af38e9a1c23d1de2cd5a5a3077de1515a4
--- /dev/null
+++ b/examples/sod_shock_tube/mono_n5.toml
@@ -0,0 +1,33 @@
+[EquationOfState]
+idealgas_gamma = 1.4
+eos = "idealgas"
+
+[EulerEq]
+id = "sod_shock_tube"
+bc = "from_id"
+av_regularization = "mono"
+
+[Mesh]
+periodic = [ false ]
+range = [ -1.0, 1.0 ]
+cfl = 0.75
+n = 5
+k = 256
+basis = "lgl"
+
+[Output]
+every_iteration = 0
+every_dt        = 0.0
+aligned_ts      = "$(collect(range(0.01,0.5,step=0.01)))"
+variables       = [ "rho", "q", "E", "ldg_rho", "ldg_q", "ldg_E", "smoothed_mu" ]
+
+[Evolution]
+tend = 0.5
+
+[Log]
+progress_stdout = true
+
+[HRSC]
+method = "av"
+av_method = "entropy"
+entropy_ce = 5.0
diff --git a/examples/sod_shock_tube/mono_n6.toml b/examples/sod_shock_tube/mono_n6.toml
new file mode 100644
index 0000000000000000000000000000000000000000..e006699f04ef244ccb9b9ac1ce329015a6ac38f2
--- /dev/null
+++ b/examples/sod_shock_tube/mono_n6.toml
@@ -0,0 +1,33 @@
+[EquationOfState]
+idealgas_gamma = 1.4
+eos = "idealgas"
+
+[EulerEq]
+id = "sod_shock_tube"
+bc = "from_id"
+av_regularization = "mono"
+
+[Mesh]
+periodic = [ false ]
+range = [ -1.0, 1.0 ]
+cfl = 0.75
+n = 6
+k = 256
+basis = "lgl"
+
+[Output]
+every_iteration = 0
+every_dt        = 0.0
+aligned_ts      = "$(collect(range(0.01,0.5,step=0.01)))"
+variables       = [ "rho", "q", "E", "ldg_rho", "ldg_q", "ldg_E", "smoothed_mu" ]
+
+[Evolution]
+tend = 0.5
+
+[Log]
+progress_stdout = true
+
+[HRSC]
+method = "av"
+av_method = "entropy"
+entropy_ce = 5.0
diff --git a/examples/sod_shock_tube/mono_n7.toml b/examples/sod_shock_tube/mono_n7.toml
new file mode 100644
index 0000000000000000000000000000000000000000..a3670c3e3696c26c845b3071f18274a30503f593
--- /dev/null
+++ b/examples/sod_shock_tube/mono_n7.toml
@@ -0,0 +1,33 @@
+[EquationOfState]
+idealgas_gamma = 1.4
+eos = "idealgas"
+
+[EulerEq]
+id = "sod_shock_tube"
+bc = "from_id"
+av_regularization = "mono"
+
+[Mesh]
+periodic = [ false ]
+range = [ -1.0, 1.0 ]
+cfl = 0.75
+n = 7
+k = 256
+basis = "lgl"
+
+[Output]
+every_iteration = 0
+every_dt        = 0.0
+aligned_ts      = "$(collect(range(0.01,0.5,step=0.01)))"
+variables       = [ "rho", "q", "E", "ldg_rho", "ldg_q", "ldg_E", "smoothed_mu" ]
+
+[Evolution]
+tend = 0.5
+
+[Log]
+progress_stdout = true
+
+[HRSC]
+method = "av"
+av_method = "entropy"
+entropy_ce = 5.0