Add several parameter files

12 files changed, 1774 insertions, 0 deletions

diff --git a/Examples/advect_shock_2.par b/Examples/advect_shock_2.par
new file mode 100644
index 0000000..f4ad696
--- /dev/null
+++ b/Examples/advect_shock_2.par
@@ -0,0 +1,129 @@
+#!/usr/bin/perl -W
+
+ActiveThorns = "
+Boundary
+Carpet
+CarpetIOASCII
+CarpetIOBasic
+CarpetIOScalar
+CarpetLib
+CarpetReduce
+CarpetSlab
+CartGrid3d
+CoordBase
+Dissipation
+GenericFD
+IOUtil
+LoopControl
+MoL
+NanChecker
+Periodic
+Advect
+Slab
+SpaceMask
+SphericalSurface
+SymBase
+Time
+"
+
+#############################################################
+# Grid
+#############################################################
+
+CoordBase::domainsize                   = minmax
+
+CoordBase::xmin                         = 0
+CoordBase::ymin                         = 0
+CoordBase::zmin                         = 0
+
+CoordBase::xmax                         = 1
+CoordBase::ymax                         = 0.1
+CoordBase::zmax                         = 0.1
+
+CoordBase::dx                           = 0.0025
+CoordBase::dy                           = 0.1
+CoordBase::dz                           = 0.1
+
+CoordBase::boundary_size_x_lower        = 2
+CoordBase::boundary_size_y_lower        = 2
+CoordBase::boundary_size_z_lower        = 2
+CoordBase::boundary_shiftout_x_lower    = 1
+CoordBase::boundary_shiftout_y_lower    = 1
+CoordBase::boundary_shiftout_z_lower    = 1
+
+CoordBase::boundary_size_x_upper        = 2
+CoordBase::boundary_size_y_upper        = 2
+CoordBase::boundary_size_z_upper        = 2
+CoordBase::boundary_shiftout_x_upper    = 0
+CoordBase::boundary_shiftout_y_upper    = 0
+CoordBase::boundary_shiftout_z_upper    = 0
+
+CartGrid3D::type                        = "coordbase"
+CartGrid3D::domain                      = "full"
+CartGrid3D::avoid_origin                = "no"
+
+Periodic::periodic = "yes"
+
+#############################################################
+# Carpet
+#############################################################
+
+Carpet::ghost_size                      = 2
+Carpet::domain_from_coordbase           = "yes"
+Carpet::max_refinement_levels           = 1
+#Carpet::init_each_timelevel             = "yes"
+Carpet::num_integrator_substeps         = 4
+
+#############################################################
+# Time integration
+#############################################################
+
+Cactus::terminate                     = "time"
+Cactus::cctk_final_time               = 0.25
+
+Time::dtfac                           = 0.5
+MethodOfLines::ode_method             = "RK4"
+MethodOfLines::MoL_Intermediate_Steps = 4
+MethodOfLines::MoL_Num_Scratch_Levels = 1
+MethodOfLines::MoL_NaN_Check          = "yes"
+
+Dissipation::vars = "Advect::rho"
+Dissipation::order = 1
+Dissipation::epsdis = 0.5
+
+Advect::amp = 1
+Advect::v0 = 1
+Advect::initial_data = "shock"
+
+#############################################################
+# Boundary conditions
+#############################################################
+
+Advect::rho_group_bound   = "none"
+
+Carpet::poison_new_timelevels           = yes
+Carpet::poison_value                    = 100
+CarpetLib::poison_new_memory            = yes
+CarpetLib::poison_value                 = 200
+Carpet::check_for_poison                = no
+
+#############################################################
+# Output
+#############################################################
+
+IO::out_dir                   = $parfile
+IO::out_fileinfo              = "all"
+
+CarpetIOBasic::outInfo_every        = 1
+CarpetIOBasic::outInfo_vars         = "Advect::rho"
+
+CarpetIOScalar::outScalar_every = 1
+CarpetIOScalar::outScalar_vars = ""
+CarpetIOScalar::outScalar_reductions = "minimum maximum norm2"
+
+IOASCII::out1D_every          = 2
+IOASCII::out1D_x              = "yes"
+IOASCII::out1D_vars           = "Advect::rho Advect::v_group"
+
+CarpetIOASCII::out_precision = 19
+CarpetIOASCII::out3D_ghosts   = "yes"
diff --git a/Examples/advect_test_2.par b/Examples/advect_test_2.par
new file mode 100644
index 0000000..3db6849
--- /dev/null
+++ b/Examples/advect_test_2.par
@@ -0,0 +1,128 @@
+#!/usr/bin/perl -W
+
+ActiveThorns = "
+Boundary
+Carpet
+CarpetIOASCII
+CarpetIOBasic
+CarpetIOScalar
+CarpetLib
+CarpetReduce
+CarpetSlab
+CartGrid3d
+CoordBase
+Dissipation
+GenericFD
+IOUtil
+LoopControl
+MoL
+NanChecker
+Periodic
+Advect
+Slab
+SpaceMask
+SphericalSurface
+SymBase
+Time
+"
+
+#############################################################
+# Grid
+#############################################################
+
+CoordBase::domainsize                   = minmax
+
+CoordBase::xmin                         = 0
+CoordBase::ymin                         = 0
+CoordBase::zmin                         = 0
+
+CoordBase::xmax                         = 1
+CoordBase::ymax                         = 0.1
+CoordBase::zmax                         = 0.1
+
+CoordBase::dx                           = 0.025
+CoordBase::dy                           = 0.1
+CoordBase::dz                           = 0.1
+
+CoordBase::boundary_size_x_lower        = 2
+CoordBase::boundary_size_y_lower        = 2
+CoordBase::boundary_size_z_lower        = 2
+CoordBase::boundary_shiftout_x_lower    = 1
+CoordBase::boundary_shiftout_y_lower    = 1
+CoordBase::boundary_shiftout_z_lower    = 1
+
+CoordBase::boundary_size_x_upper        = 2
+CoordBase::boundary_size_y_upper        = 2
+CoordBase::boundary_size_z_upper        = 2
+CoordBase::boundary_shiftout_x_upper    = 0
+CoordBase::boundary_shiftout_y_upper    = 0
+CoordBase::boundary_shiftout_z_upper    = 0
+
+CartGrid3D::type                        = "coordbase"
+CartGrid3D::domain                      = "full"
+CartGrid3D::avoid_origin                = "no"
+
+Periodic::periodic = "yes"
+
+#############################################################
+# Carpet
+#############################################################
+
+Carpet::ghost_size                      = 2
+Carpet::domain_from_coordbase           = "yes"
+Carpet::max_refinement_levels           = 1
+#Carpet::init_each_timelevel             = "yes"
+Carpet::num_integrator_substeps         = 4
+
+#############################################################
+# Time integration
+#############################################################
+
+Cactus::terminate                     = "time"
+Cactus::cctk_final_time               = 1
+
+Time::dtfac                           = 0.5
+MethodOfLines::ode_method             = "RK4"
+MethodOfLines::MoL_Intermediate_Steps = 4
+MethodOfLines::MoL_Num_Scratch_Levels = 1
+MethodOfLines::MoL_NaN_Check          = "yes"
+
+Dissipation::vars = "Advect::rho"
+Dissipation::order = 3
+Dissipation::epsdis = 0.1
+
+Advect::amp = 1e-8
+Advect::v0 = 1
+
+#############################################################
+# Boundary conditions
+#############################################################
+
+Advect::rho_group_bound   = "none"
+
+Carpet::poison_new_timelevels           = yes
+Carpet::poison_value                    = 100
+CarpetLib::poison_new_memory            = yes
+CarpetLib::poison_value                 = 200
+Carpet::check_for_poison                = no
+
+#############################################################
+# Output
+#############################################################
+
+IO::out_dir                   = $parfile
+IO::out_fileinfo              = "all"
+
+CarpetIOBasic::outInfo_every        = 1
+CarpetIOBasic::outInfo_vars         = "Advect::rho"
+
+CarpetIOScalar::outScalar_every = 1
+CarpetIOScalar::outScalar_vars = ""
+CarpetIOScalar::outScalar_reductions = "minimum maximum norm2"
+
+IOASCII::out1D_every          = 1
+IOASCII::out1D_x              = "yes"
+IOASCII::out1D_vars           = "Advect::rho Advect::v_group"
+
+CarpetIOASCII::out_precision = 19
+CarpetIOASCII::out3D_ghosts   = "yes"
diff --git a/Examples/burgers_shock_2.par b/Examples/burgers_shock_2.par
new file mode 100644
index 0000000..ef4fdcb
--- /dev/null
+++ b/Examples/burgers_shock_2.par
@@ -0,0 +1,134 @@
+#!/usr/bin/perl -W
+
+ActiveThorns = "
+Boundary
+Carpet
+CarpetIOASCII
+CarpetIOBasic
+CarpetIOScalar
+CarpetLib
+CarpetReduce
+CarpetSlab
+CartGrid3d
+CoordBase
+Dissipation
+GenericFD
+IOUtil
+LoopControl
+MoL
+NanChecker
+Periodic
+Burgers
+Slab
+SpaceMask
+SphericalSurface
+SymBase
+Time
+"
+
+#############################################################
+# Grid
+#############################################################
+
+CoordBase::domainsize                   = minmax
+
+CoordBase::xmin                         = -1
+CoordBase::ymin                         = 0
+CoordBase::zmin                         = 0
+
+CoordBase::xmax                         = 2
+CoordBase::ymax                         = 0.1
+CoordBase::zmax                         = 0.1
+
+CoordBase::dx                           = 0.02
+CoordBase::dy                           = 0.1
+CoordBase::dz                           = 0.1
+
+CoordBase::boundary_size_x_lower        = 2
+CoordBase::boundary_size_y_lower        = 2
+CoordBase::boundary_size_z_lower        = 2
+CoordBase::boundary_shiftout_x_lower    = 1
+CoordBase::boundary_shiftout_y_lower    = 1
+CoordBase::boundary_shiftout_z_lower    = 1
+
+CoordBase::boundary_size_x_upper        = 2
+CoordBase::boundary_size_y_upper        = 2
+CoordBase::boundary_size_z_upper        = 2
+CoordBase::boundary_shiftout_x_upper    = 0
+CoordBase::boundary_shiftout_y_upper    = 0
+CoordBase::boundary_shiftout_z_upper    = 0
+
+CartGrid3D::type                        = "coordbase"
+CartGrid3D::domain                      = "full"
+CartGrid3D::avoid_origin                = "no"
+
+Periodic::periodic = "yes"
+
+#############################################################
+# Carpet
+#############################################################
+
+Carpet::ghost_size                      = 2
+Carpet::domain_from_coordbase           = "yes"
+Carpet::max_refinement_levels           = 1
+#Carpet::init_each_timelevel             = "yes"
+Carpet::num_integrator_substeps         = 4
+
+#############################################################
+# Time integration
+#############################################################
+
+Cactus::terminate                     = "time"
+Cactus::cctk_final_time               = 2
+
+Time::dtfac                           = 0.25
+MethodOfLines::ode_method             = "RK4"
+MethodOfLines::MoL_Intermediate_Steps = 4
+MethodOfLines::MoL_Num_Scratch_Levels = 1
+MethodOfLines::MoL_NaN_Check          = "yes"
+
+# MethodOfLines::ode_method             = "generic"
+# MethodOfLines::generic_type             = "RK"
+# MethodOfLines::MoL_Intermediate_Steps = 1
+
+Dissipation::vars = ""
+Dissipation::order = 1
+Dissipation::epsdis = 0.5
+
+Burgers::uL0 = 1
+Burgers::uR0 = 0
+Burgers::alpha = 1
+Burgers::initial_data = "shock"
+
+#############################################################
+# Boundary conditions
+#############################################################
+
+Burgers::u_group_bound   = "none"
+
+Carpet::poison_new_timelevels           = yes
+Carpet::poison_value                    = 100
+CarpetLib::poison_new_memory            = yes
+CarpetLib::poison_value                 = 200
+Carpet::check_for_poison                = no
+
+#############################################################
+# Output
+#############################################################
+
+IO::out_dir                   = $parfile
+IO::out_fileinfo              = "all"
+
+CarpetIOBasic::outInfo_every        = 1
+CarpetIOBasic::outInfo_vars         = "Burgers::u Burgers::uF Burgers::uR Burgers::uLeft Burgers::urhs"
+
+CarpetIOScalar::outScalar_every = 1
+CarpetIOScalar::outScalar_vars = ""
+CarpetIOScalar::outScalar_reductions = "minimum maximum norm2"
+
+IOASCII::out1D_every          = 2
+IOASCII::out1D_x              = "yes"
+IOASCII::out1D_vars           = "Burgers::u Burgers::uF Burgers::uR Burgers::uLeft Burgers::urhs"
+
+CarpetIOASCII::out_precision = 19
+CarpetIOASCII::out3D_ghosts   = "yes"
diff --git a/Examples/burgers_shock_3.par b/Examples/burgers_shock_3.par
new file mode 100644
index 0000000..90c1c02
--- /dev/null
+++ b/Examples/burgers_shock_3.par
@@ -0,0 +1,134 @@
+#!/usr/bin/perl -W
+
+ActiveThorns = "
+Boundary
+Carpet
+CarpetIOASCII
+CarpetIOBasic
+CarpetIOScalar
+CarpetLib
+CarpetReduce
+CarpetSlab
+CartGrid3d
+CoordBase
+Dissipation
+GenericFD
+IOUtil
+LoopControl
+MoL
+NanChecker
+Periodic
+Burgers
+Slab
+SpaceMask
+SphericalSurface
+SymBase
+Time
+"
+
+#############################################################
+# Grid
+#############################################################
+
+CoordBase::domainsize                   = minmax
+
+CoordBase::xmin                         = -1
+CoordBase::ymin                         = 0
+CoordBase::zmin                         = 0
+
+CoordBase::xmax                         = 2
+CoordBase::ymax                         = 0.1
+CoordBase::zmax                         = 0.1
+
+CoordBase::dx                           = 0.01
+CoordBase::dy                           = 0.1
+CoordBase::dz                           = 0.1
+
+CoordBase::boundary_size_x_lower        = 2
+CoordBase::boundary_size_y_lower        = 2
+CoordBase::boundary_size_z_lower        = 2
+CoordBase::boundary_shiftout_x_lower    = 1
+CoordBase::boundary_shiftout_y_lower    = 1
+CoordBase::boundary_shiftout_z_lower    = 1
+
+CoordBase::boundary_size_x_upper        = 2
+CoordBase::boundary_size_y_upper        = 2
+CoordBase::boundary_size_z_upper        = 2
+CoordBase::boundary_shiftout_x_upper    = 0
+CoordBase::boundary_shiftout_y_upper    = 0
+CoordBase::boundary_shiftout_z_upper    = 0
+
+CartGrid3D::type                        = "coordbase"
+CartGrid3D::domain                      = "full"
+CartGrid3D::avoid_origin                = "no"
+
+Periodic::periodic = "yes"
+
+#############################################################
+# Carpet
+#############################################################
+
+Carpet::ghost_size                      = 2
+Carpet::domain_from_coordbase           = "yes"
+Carpet::max_refinement_levels           = 1
+#Carpet::init_each_timelevel             = "yes"
+Carpet::num_integrator_substeps         = 4
+
+#############################################################
+# Time integration
+#############################################################
+
+Cactus::terminate                     = "time"
+Cactus::cctk_final_time               = 2
+
+Time::dtfac                           = 0.25
+MethodOfLines::ode_method             = "RK4"
+MethodOfLines::MoL_Intermediate_Steps = 4
+MethodOfLines::MoL_Num_Scratch_Levels = 1
+MethodOfLines::MoL_NaN_Check          = "yes"
+
+# MethodOfLines::ode_method             = "generic"
+# MethodOfLines::generic_type             = "RK"
+# MethodOfLines::MoL_Intermediate_Steps = 1
+
+Dissipation::vars = ""
+Dissipation::order = 1
+Dissipation::epsdis = 0.5
+
+Burgers::uL0 = 1
+Burgers::uR0 = 0
+Burgers::alpha = 1
+Burgers::initial_data = "shock"
+
+#############################################################
+# Boundary conditions
+#############################################################
+
+Burgers::u_group_bound   = "none"
+
+Carpet::poison_new_timelevels           = yes
+Carpet::poison_value                    = 100
+CarpetLib::poison_new_memory            = yes
+CarpetLib::poison_value                 = 200
+Carpet::check_for_poison                = no
+
+#############################################################
+# Output
+#############################################################
+
+IO::out_dir                   = $parfile
+IO::out_fileinfo              = "all"
+
+CarpetIOBasic::outInfo_every        = 1
+CarpetIOBasic::outInfo_vars         = "Burgers::u Burgers::uF Burgers::uR Burgers::uLeft Burgers::urhs"
+
+CarpetIOScalar::outScalar_every = 1
+CarpetIOScalar::outScalar_vars = ""
+CarpetIOScalar::outScalar_reductions = "minimum maximum norm2"
+
+IOASCII::out1D_every          = 4
+IOASCII::out1D_x              = "yes"
+IOASCII::out1D_vars           = "Burgers::u Burgers::uF Burgers::uR Burgers::uLeft Burgers::urhs"
+
+CarpetIOASCII::out_precision = 19
+CarpetIOASCII::out3D_ghosts   = "yes"
diff --git a/Examples/burgers_sine_2.par b/Examples/burgers_sine_2.par
new file mode 100644
index 0000000..d9a65e2
--- /dev/null
+++ b/Examples/burgers_sine_2.par
@@ -0,0 +1,133 @@
+#!/usr/bin/perl -W
+
+ActiveThorns = "
+Boundary
+Carpet
+CarpetIOASCII
+CarpetIOBasic
+CarpetIOScalar
+CarpetLib
+CarpetReduce
+CarpetSlab
+CartGrid3d
+CoordBase
+Dissipation
+GenericFD
+IOUtil
+LoopControl
+MoL
+NanChecker
+Periodic
+Burgers
+Slab
+SpaceMask
+SphericalSurface
+SymBase
+Time
+"
+
+#############################################################
+# Grid
+#############################################################
+
+CoordBase::domainsize                   = minmax
+
+CoordBase::xmin                         = 0
+CoordBase::ymin                         = 0
+CoordBase::zmin                         = 0
+
+CoordBase::xmax                         = 1
+CoordBase::ymax                         = 0.1
+CoordBase::zmax                         = 0.1
+
+CoordBase::dx                           = 0.01
+CoordBase::dy                           = 0.1
+CoordBase::dz                           = 0.1
+
+CoordBase::boundary_size_x_lower        = 2
+CoordBase::boundary_size_y_lower        = 2
+CoordBase::boundary_size_z_lower        = 2
+CoordBase::boundary_shiftout_x_lower    = 1
+CoordBase::boundary_shiftout_y_lower    = 1
+CoordBase::boundary_shiftout_z_lower    = 1
+
+CoordBase::boundary_size_x_upper        = 2
+CoordBase::boundary_size_y_upper        = 2
+CoordBase::boundary_size_z_upper        = 2
+CoordBase::boundary_shiftout_x_upper    = 0
+CoordBase::boundary_shiftout_y_upper    = 0
+CoordBase::boundary_shiftout_z_upper    = 0
+
+CartGrid3D::type                        = "coordbase"
+CartGrid3D::domain                      = "full"
+CartGrid3D::avoid_origin                = "no"
+
+Periodic::periodic = "yes"
+
+#############################################################
+# Carpet
+#############################################################
+
+Carpet::ghost_size                      = 2
+Carpet::domain_from_coordbase           = "yes"
+Carpet::max_refinement_levels           = 1
+#Carpet::init_each_timelevel             = "yes"
+Carpet::num_integrator_substeps         = 4
+
+#############################################################
+# Time integration
+#############################################################
+
+Cactus::terminate                     = "time"
+Cactus::cctk_final_time               = 2
+
+Time::dtfac                           = 0.25
+MethodOfLines::ode_method             = "RK4"
+MethodOfLines::MoL_Intermediate_Steps = 4
+MethodOfLines::MoL_Num_Scratch_Levels = 1
+MethodOfLines::MoL_NaN_Check          = "yes"
+
+# MethodOfLines::ode_method             = "generic"
+# MethodOfLines::generic_type             = "RK"
+# MethodOfLines::MoL_Intermediate_Steps = 1
+
+Dissipation::vars = ""
+Dissipation::order = 1
+Dissipation::epsdis = 0.5
+
+Burgers::amp = 0.01
+Burgers::alpha = 1
+Burgers::initial_data = "sine"
+
+#############################################################
+# Boundary conditions
+#############################################################
+
+Burgers::u_group_bound   = "none"
+
+Carpet::poison_new_timelevels           = yes
+Carpet::poison_value                    = 100
+CarpetLib::poison_new_memory            = yes
+CarpetLib::poison_value                 = 200
+Carpet::check_for_poison                = no
+
+#############################################################
+# Output
+#############################################################
+
+IO::out_dir                   = $parfile
+IO::out_fileinfo              = "all"
+
+CarpetIOBasic::outInfo_every        = 1
+CarpetIOBasic::outInfo_vars         = "Burgers::u Burgers::uF Burgers::uR Burgers::uLeft Burgers::urhs"
+
+CarpetIOScalar::outScalar_every = 1
+CarpetIOScalar::outScalar_vars = ""
+CarpetIOScalar::outScalar_reductions = "minimum maximum norm2"
+
+IOASCII::out1D_every          = 2
+IOASCII::out1D_x              = "yes"
+IOASCII::out1D_vars           = "Burgers::u Burgers::uF Burgers::uR Burgers::uLeft Burgers::urhs"
+
+CarpetIOASCII::out_precision = 19
+CarpetIOASCII::out3D_ghosts   = "yes"
diff --git a/Examples/burgers_sine_3.par b/Examples/burgers_sine_3.par
new file mode 100644
index 0000000..279ee00
--- /dev/null
+++ b/Examples/burgers_sine_3.par
@@ -0,0 +1,133 @@
+#!/usr/bin/perl -W
+
+ActiveThorns = "
+Boundary
+Carpet
+CarpetIOASCII
+CarpetIOBasic
+CarpetIOScalar
+CarpetLib
+CarpetReduce
+CarpetSlab
+CartGrid3d
+CoordBase
+Dissipation
+GenericFD
+IOUtil
+LoopControl
+MoL
+NanChecker
+Periodic
+Burgers
+Slab
+SpaceMask
+SphericalSurface
+SymBase
+Time
+"
+
+#############################################################
+# Grid
+#############################################################
+
+CoordBase::domainsize                   = minmax
+
+CoordBase::xmin                         = 0
+CoordBase::ymin                         = 0
+CoordBase::zmin                         = 0
+
+CoordBase::xmax                         = 1
+CoordBase::ymax                         = 0.1
+CoordBase::zmax                         = 0.1
+
+CoordBase::dx                           = 0.005
+CoordBase::dy                           = 0.1
+CoordBase::dz                           = 0.1
+
+CoordBase::boundary_size_x_lower        = 2
+CoordBase::boundary_size_y_lower        = 2
+CoordBase::boundary_size_z_lower        = 2
+CoordBase::boundary_shiftout_x_lower    = 1
+CoordBase::boundary_shiftout_y_lower    = 1
+CoordBase::boundary_shiftout_z_lower    = 1
+
+CoordBase::boundary_size_x_upper        = 2
+CoordBase::boundary_size_y_upper        = 2
+CoordBase::boundary_size_z_upper        = 2
+CoordBase::boundary_shiftout_x_upper    = 0
+CoordBase::boundary_shiftout_y_upper    = 0
+CoordBase::boundary_shiftout_z_upper    = 0
+
+CartGrid3D::type                        = "coordbase"
+CartGrid3D::domain                      = "full"
+CartGrid3D::avoid_origin                = "no"
+
+Periodic::periodic = "yes"
+
+#############################################################
+# Carpet
+#############################################################
+
+Carpet::ghost_size                      = 2
+Carpet::domain_from_coordbase           = "yes"
+Carpet::max_refinement_levels           = 1
+#Carpet::init_each_timelevel             = "yes"
+Carpet::num_integrator_substeps         = 4
+
+#############################################################
+# Time integration
+#############################################################
+
+Cactus::terminate                     = "time"
+Cactus::cctk_final_time               = 2
+
+Time::dtfac                           = 0.25
+MethodOfLines::ode_method             = "RK4"
+MethodOfLines::MoL_Intermediate_Steps = 4
+MethodOfLines::MoL_Num_Scratch_Levels = 1
+MethodOfLines::MoL_NaN_Check          = "yes"
+
+# MethodOfLines::ode_method             = "generic"
+# MethodOfLines::generic_type             = "RK"
+# MethodOfLines::MoL_Intermediate_Steps = 1
+
+Dissipation::vars = ""
+Dissipation::order = 1
+Dissipation::epsdis = 0.5
+
+Burgers::amp = 0.01
+Burgers::alpha = 1
+Burgers::initial_data = "sine"
+
+#############################################################
+# Boundary conditions
+#############################################################
+
+Burgers::u_group_bound   = "none"
+
+Carpet::poison_new_timelevels           = yes
+Carpet::poison_value                    = 100
+CarpetLib::poison_new_memory            = yes
+CarpetLib::poison_value                 = 200
+Carpet::check_for_poison                = no
+
+#############################################################
+# Output
+#############################################################
+
+IO::out_dir                   = $parfile
+IO::out_fileinfo              = "all"
+
+CarpetIOBasic::outInfo_every        = 1
+CarpetIOBasic::outInfo_vars         = "Burgers::u Burgers::uF Burgers::uR Burgers::uLeft Burgers::urhs"
+
+CarpetIOScalar::outScalar_every = 1
+CarpetIOScalar::outScalar_vars = ""
+CarpetIOScalar::outScalar_reductions = "minimum maximum norm2"
+
+IOASCII::out1D_every          = 4
+IOASCII::out1D_x              = "yes"
+IOASCII::out1D_vars           = "Burgers::u Burgers::uF Burgers::uR Burgers::uLeft Burgers::urhs"
+
+CarpetIOASCII::out_precision = 19
+CarpetIOASCII::out3D_ghosts   = "yes"
diff --git a/Examples/euler_shock_2.par b/Examples/euler_shock_2.par
new file mode 100644
index 0000000..c02e11d
--- /dev/null
+++ b/Examples/euler_shock_2.par
@@ -0,0 +1,141 @@
+#!/usr/bin/perl -W
+
+ActiveThorns = "
+Boundary
+Carpet
+CarpetIOASCII
+CarpetIOBasic
+CarpetIOScalar
+CarpetLib
+CarpetReduce
+CarpetSlab
+CartGrid3d
+CoordBase
+Dissipation
+GenericFD
+IOUtil
+LoopControl
+MoL
+NanChecker
+Periodic
+Euler
+Slab
+SpaceMask
+SphericalSurface
+SymBase
+Time
+"
+
+#############################################################
+# Grid
+#############################################################
+
+CoordBase::domainsize                   = minmax
+
+CoordBase::xmin                         = -1
+CoordBase::ymin                         = 0
+CoordBase::zmin                         = 0
+
+CoordBase::xmax                         = 2
+CoordBase::ymax                         = 0.1
+CoordBase::zmax                         = 0.1
+
+CoordBase::dx                           = 0.02
+CoordBase::dy                           = 0.1
+CoordBase::dz                           = 0.1
+
+CoordBase::boundary_size_x_lower        = 2
+CoordBase::boundary_size_y_lower        = 2
+CoordBase::boundary_size_z_lower        = 2
+CoordBase::boundary_shiftout_x_lower    = 1
+CoordBase::boundary_shiftout_y_lower    = 1
+CoordBase::boundary_shiftout_z_lower    = 1
+
+CoordBase::boundary_size_x_upper        = 2
+CoordBase::boundary_size_y_upper        = 2
+CoordBase::boundary_size_z_upper        = 2
+CoordBase::boundary_shiftout_x_upper    = 0
+CoordBase::boundary_shiftout_y_upper    = 0
+CoordBase::boundary_shiftout_z_upper    = 0
+
+CartGrid3D::type                        = "coordbase"
+CartGrid3D::domain                      = "full"
+CartGrid3D::avoid_origin                = "no"
+
+Periodic::periodic = "yes"
+
+#############################################################
+# Carpet
+#############################################################
+
+Carpet::ghost_size                      = 2
+Carpet::domain_from_coordbase           = "yes"
+Carpet::max_refinement_levels           = 1
+#Carpet::init_each_timelevel             = "yes"
+Carpet::num_integrator_substeps         = 4
+
+#############################################################
+# Time integration
+#############################################################
+
+Cactus::terminate                     = "time"
+Cactus::cctk_final_time               = 0.2
+
+Time::dtfac                           = 0.25
+MethodOfLines::ode_method             = "RK4"
+MethodOfLines::MoL_Intermediate_Steps = 4
+MethodOfLines::MoL_Num_Scratch_Levels = 1
+MethodOfLines::MoL_NaN_Check          = "yes"
+
+# MethodOfLines::ode_method             = "generic"
+# MethodOfLines::generic_type             = "RK"
+# MethodOfLines::MoL_Intermediate_Steps = 1
+
+Dissipation::vars = ""
+Dissipation::order = 1
+Dissipation::epsdis = 0.5
+
+Euler::rhoL0 = 1
+Euler::rhoR0 = 0.125
+Euler::vL0 = 0
+Euler::vR0 = 0
+Euler::pL0 = 1
+Euler::pR0 = 0.1
+Euler::alpha = 2
+Euler::gamma = 1.4
+Euler::initial_data = "shock"
+
+#############################################################
+# Boundary conditions
+#############################################################
+
+Euler::Den_group_bound   = "none"
+Euler::S_group_bound   = "none"
+Euler::En_group_bound   = "none"
+
+Carpet::poison_new_timelevels           = yes
+Carpet::poison_value                    = 100
+CarpetLib::poison_new_memory            = yes
+CarpetLib::poison_value                 = 200
+Carpet::check_for_poison                = no
+
+#############################################################
+# Output
+#############################################################
+
+IO::out_dir                   = $parfile
+IO::out_fileinfo              = "all"
+
+CarpetIOBasic::outInfo_every        = 1
+CarpetIOBasic::outInfo_vars         = "Euler::rho Euler::DenF Euler::v1 Euler::p Euler::Den Euler::S1 Euler::En"
+
+CarpetIOScalar::outScalar_every = 1
+CarpetIOScalar::outScalar_vars = ""
+CarpetIOScalar::outScalar_reductions = "minimum maximum norm2"
+
+IOASCII::out1D_every          = 2
+IOASCII::out1D_x              = "yes"
+IOASCII::out1D_vars           = "Euler::rho Euler::v1 Euler::p Euler::Den Euler::S1 Euler::En"
+
+CarpetIOASCII::out_precision = 19
+CarpetIOASCII::out3D_ghosts   = "yes"
diff --git a/Examples/euler_shock_3.par b/Examples/euler_shock_3.par
new file mode 100644
index 0000000..78a6f3a
--- /dev/null
+++ b/Examples/euler_shock_3.par
@@ -0,0 +1,141 @@
+#!/usr/bin/perl -W
+
+ActiveThorns = "
+Boundary
+Carpet
+CarpetIOASCII
+CarpetIOBasic
+CarpetIOScalar
+CarpetLib
+CarpetReduce
+CarpetSlab
+CartGrid3d
+CoordBase
+Dissipation
+GenericFD
+IOUtil
+LoopControl
+MoL
+NanChecker
+Periodic
+Euler
+Slab
+SpaceMask
+SphericalSurface
+SymBase
+Time
+"
+
+#############################################################
+# Grid
+#############################################################
+
+CoordBase::domainsize                   = minmax
+
+CoordBase::xmin                         = -1
+CoordBase::ymin                         = 0
+CoordBase::zmin                         = 0
+
+CoordBase::xmax                         = 2
+CoordBase::ymax                         = 0.1
+CoordBase::zmax                         = 0.1
+
+CoordBase::dx                           = 0.0025
+CoordBase::dy                           = 0.1
+CoordBase::dz                           = 0.1
+
+CoordBase::boundary_size_x_lower        = 2
+CoordBase::boundary_size_y_lower        = 2
+CoordBase::boundary_size_z_lower        = 2
+CoordBase::boundary_shiftout_x_lower    = 1
+CoordBase::boundary_shiftout_y_lower    = 1
+CoordBase::boundary_shiftout_z_lower    = 1
+
+CoordBase::boundary_size_x_upper        = 2
+CoordBase::boundary_size_y_upper        = 2
+CoordBase::boundary_size_z_upper        = 2
+CoordBase::boundary_shiftout_x_upper    = 0
+CoordBase::boundary_shiftout_y_upper    = 0
+CoordBase::boundary_shiftout_z_upper    = 0
+
+CartGrid3D::type                        = "coordbase"
+CartGrid3D::domain                      = "full"
+CartGrid3D::avoid_origin                = "no"
+
+Periodic::periodic = "yes"
+
+#############################################################
+# Carpet
+#############################################################
+
+Carpet::ghost_size                      = 2
+Carpet::domain_from_coordbase           = "yes"
+Carpet::max_refinement_levels           = 1
+#Carpet::init_each_timelevel             = "yes"
+Carpet::num_integrator_substeps         = 4
+
+#############################################################
+# Time integration
+#############################################################
+
+Cactus::terminate                     = "time"
+Cactus::cctk_final_time               = 0.2
+
+Time::dtfac                           = 0.25
+MethodOfLines::ode_method             = "RK4"
+MethodOfLines::MoL_Intermediate_Steps = 4
+MethodOfLines::MoL_Num_Scratch_Levels = 1
+MethodOfLines::MoL_NaN_Check          = "yes"
+
+# MethodOfLines::ode_method             = "generic"
+# MethodOfLines::generic_type             = "RK"
+# MethodOfLines::MoL_Intermediate_Steps = 1
+
+Dissipation::vars = ""
+Dissipation::order = 1
+Dissipation::epsdis = 0.5
+
+Euler::rhoL0 = 1
+Euler::rhoR0 = 0.125
+Euler::vL0 = 0
+Euler::vR0 = 0
+Euler::pL0 = 1
+Euler::pR0 = 0.1
+Euler::alpha = 3
+Euler::gamma = 1.4
+Euler::initial_data = "shock"
+
+#############################################################
+# Boundary conditions
+#############################################################
+
+Euler::Den_group_bound   = "none"
+Euler::S_group_bound   = "none"
+Euler::En_group_bound   = "none"
+
+Carpet::poison_new_timelevels           = yes
+Carpet::poison_value                    = 100
+CarpetLib::poison_new_memory            = yes
+CarpetLib::poison_value                 = 200
+Carpet::check_for_poison                = no
+
+#############################################################
+# Output
+#############################################################
+
+IO::out_dir                   = $parfile
+IO::out_fileinfo              = "all"
+
+CarpetIOBasic::outInfo_every        = 1
+CarpetIOBasic::outInfo_vars         = "Euler::rho Euler::DenF Euler::v1 Euler::p Euler::Den Euler::S1 Euler::En"
+
+CarpetIOScalar::outScalar_every = 1
+CarpetIOScalar::outScalar_vars = ""
+CarpetIOScalar::outScalar_reductions = "minimum maximum norm2"
+
+IOASCII::out1D_every          = 16
+IOASCII::out1D_x              = "yes"
+IOASCII::out1D_vars           = "Euler::rho Euler::v1 Euler::p Euler::Den Euler::S1 Euler::En"
+
+CarpetIOASCII::out_precision = 19
+CarpetIOASCII::out3D_ghosts   = "yes"
diff --git a/Examples/euler_test_2.par b/Examples/euler_test_2.par
new file mode 100644
index 0000000..efb6d4b
--- /dev/null
+++ b/Examples/euler_test_2.par
@@ -0,0 +1,130 @@
+#!/usr/bin/perl -W
+
+ActiveThorns = "
+Boundary
+Carpet
+CarpetIOASCII
+CarpetIOBasic
+CarpetIOScalar
+CarpetLib
+CarpetReduce
+CarpetSlab
+CartGrid3d
+CoordBase
+Dissipation
+GenericFD
+IOUtil
+LoopControl
+MoL
+NanChecker
+Periodic
+Euler
+Slab
+SpaceMask
+SphericalSurface
+SymBase
+Time
+"
+
+#############################################################
+# Grid
+#############################################################
+
+CoordBase::domainsize                   = minmax
+
+CoordBase::xmin                         = 0
+CoordBase::ymin                         = 0
+CoordBase::zmin                         = 0
+
+CoordBase::xmax                         = 1
+CoordBase::ymax                         = 0.1
+CoordBase::zmax                         = 0.1
+
+CoordBase::dx                           = 0.025
+CoordBase::dy                           = 0.1
+CoordBase::dz                           = 0.1
+
+CoordBase::boundary_size_x_lower        = 2
+CoordBase::boundary_size_y_lower        = 2
+CoordBase::boundary_size_z_lower        = 2
+CoordBase::boundary_shiftout_x_lower    = 1
+CoordBase::boundary_shiftout_y_lower    = 1
+CoordBase::boundary_shiftout_z_lower    = 1
+
+CoordBase::boundary_size_x_upper        = 2
+CoordBase::boundary_size_y_upper        = 2
+CoordBase::boundary_size_z_upper        = 2
+CoordBase::boundary_shiftout_x_upper    = 0
+CoordBase::boundary_shiftout_y_upper    = 0
+CoordBase::boundary_shiftout_z_upper    = 0
+
+CartGrid3D::type                        = "coordbase"
+CartGrid3D::domain                      = "full"
+CartGrid3D::avoid_origin                = "no"
+
+Periodic::periodic = "yes"
+
+#############################################################
+# Carpet
+#############################################################
+
+Carpet::ghost_size                      = 2
+Carpet::domain_from_coordbase           = "yes"
+Carpet::max_refinement_levels           = 1
+#Carpet::init_each_timelevel             = "yes"
+Carpet::num_integrator_substeps         = 4
+
+#############################################################
+# Time integration
+#############################################################
+
+Cactus::terminate                     = "time"
+Cactus::cctk_final_time               = 1
+
+Time::dtfac                           = 0.5
+MethodOfLines::ode_method             = "RK4"
+MethodOfLines::MoL_Intermediate_Steps = 4
+MethodOfLines::MoL_Num_Scratch_Levels = 1
+MethodOfLines::MoL_NaN_Check          = "yes"
+
+Dissipation::vars = "Euler::rho Euler::w_group Euler::En"
+Dissipation::order = 3
+Dissipation::epsdis = 0.1
+
+Euler::amp = 1e-8
+Euler::v0 = 1
+
+#############################################################
+# Boundary conditions
+#############################################################
+
+Euler::rho_group_bound   = "none"
+Euler::w_group_bound   = "none"
+Euler::En_group_bound   = "none"
+
+Carpet::poison_new_timelevels           = yes
+Carpet::poison_value                    = 100
+CarpetLib::poison_new_memory            = yes
+CarpetLib::poison_value                 = 200
+Carpet::check_for_poison                = no
+
+#############################################################
+# Output
+#############################################################
+
+IO::out_dir                   = $parfile
+IO::out_fileinfo              = "all"
+
+CarpetIOBasic::outInfo_every        = 1
+CarpetIOBasic::outInfo_vars         = "Euler::rho"
+
+CarpetIOScalar::outScalar_every = 1
+CarpetIOScalar::outScalar_vars = ""
+CarpetIOScalar::outScalar_reductions = "minimum maximum norm2"
+
+IOASCII::out1D_every          = 1
+IOASCII::out1D_x              = "yes"
+IOASCII::out1D_vars           = "Euler::rho Euler::w_group Euler::En Euler::v_group"
+
+CarpetIOASCII::out_precision = 19
+CarpetIOASCII::out3D_ghosts   = "yes"
diff --git a/Examples/eulerauto_shock.par b/Examples/eulerauto_shock.par
new file mode 100644
index 0000000..42d1cfa
--- /dev/null
+++ b/Examples/eulerauto_shock.par
@@ -0,0 +1,170 @@
+#!/usr/bin/perl -W
+
+ActiveThorns = "
+Boundary
+Carpet
+CarpetIOASCII
+CarpetIOBasic
+CarpetIOScalar
+CarpetLib
+CarpetReduce
+CarpetSlab
+CartGrid3d
+CoordBase
+Dissipation
+GenericFD
+IOUtil
+LoopControl
+MoL
+NanChecker
+Periodic
+EulerAuto
+Slab
+SpaceMask
+SphericalSurface
+SymBase
+Time
+"
+
+#############################################################
+# Grid
+#############################################################
+
+CoordBase::domainsize                   = minmax
+
+CoordBase::xmin                         = -1
+CoordBase::ymin                         = 0
+CoordBase::zmin                         = 0
+
+CoordBase::xmax                         = 2
+CoordBase::ymax                         = 0.1
+CoordBase::zmax                         = 0.1
+
+CoordBase::dx                           = 0.0025
+CoordBase::dy                           = 0.1
+CoordBase::dz                           = 0.1
+
+CoordBase::boundary_size_x_lower        = 2
+CoordBase::boundary_size_y_lower        = 2
+CoordBase::boundary_size_z_lower        = 2
+CoordBase::boundary_shiftout_x_lower    = 1
+CoordBase::boundary_shiftout_y_lower    = 1
+CoordBase::boundary_shiftout_z_lower    = 1
+
+CoordBase::boundary_size_x_upper        = 2
+CoordBase::boundary_size_y_upper        = 2
+CoordBase::boundary_size_z_upper        = 2
+CoordBase::boundary_shiftout_x_upper    = 0
+CoordBase::boundary_shiftout_y_upper    = 0
+CoordBase::boundary_shiftout_z_upper    = 0
+
+CartGrid3D::type                        = "coordbase"
+CartGrid3D::domain                      = "full"
+CartGrid3D::avoid_origin                = "no"
+
+Periodic::periodic = "yes"
+
+#############################################################
+# Carpet
+#############################################################
+
+Carpet::ghost_size                      = 2
+Carpet::domain_from_coordbase           = "yes"
+Carpet::max_refinement_levels           = 1
+#Carpet::init_each_timelevel             = "yes"
+Carpet::num_integrator_substeps         = 4
+
+#############################################################
+# Time integration
+#############################################################
+
+Cactus::terminate                     = "time"
+Cactus::cctk_final_time               = 0.2
+Cactus::cctk_itlast               = 2
+
+Time::dtfac                           = 0.25
+MethodOfLines::ode_method             = "RK4"
+MethodOfLines::MoL_Intermediate_Steps = 4
+MethodOfLines::MoL_Num_Scratch_Levels = 1
+MethodOfLines::MoL_NaN_Check          = "yes"
+
+# MethodOfLines::ode_method             = "generic"
+# MethodOfLines::generic_type             = "RK"
+# MethodOfLines::MoL_Intermediate_Steps = 1
+
+Dissipation::vars = ""
+Dissipation::order = 1
+Dissipation::epsdis = 0.5
+
+EulerAuto::rhoL0 = 1
+EulerAuto::rhoR0 = 0.125
+EulerAuto::vL0 = 0
+EulerAuto::vR0 = 0
+EulerAuto::pL0 = 1
+EulerAuto::pR0 = 0.1
+EulerAuto::hlleAlpha = 3
+EulerAuto::gamma = 1.4
+EulerAuto::initial_data = "shock"
+
+#############################################################
+# Boundary conditions
+#############################################################
+
+EulerAuto::Den_group_bound   = "none"
+EulerAuto::S_group_bound   = "none"
+EulerAuto::En_group_bound   = "none"
+
+Carpet::poison_new_timelevels           = yes
+Carpet::poison_value                    = 100
+CarpetLib::poison_new_memory            = yes
+CarpetLib::poison_value                 = 200
+Carpet::check_for_poison                = no
+
+#############################################################
+# Output
+#############################################################
+
+IO::out_dir                   = $parfile
+IO::out_fileinfo              = "none"
+
+CarpetIOBasic::outInfo_every        = 1
+CarpetIOBasic::outInfo_vars         = "EulerAuto::rho EulerAuto::v1 EulerAuto::p EulerAuto::Den EulerAuto::S1 EulerAuto::En"
+
+CarpetIOScalar::outScalar_every = 1
+CarpetIOScalar::outScalar_vars = ""
+CarpetIOScalar::outScalar_reductions = "minimum maximum norm2"
+
+IOASCII::out1D_every          = 16
+IOASCII::out1D_x              = "yes"
+IOASCII::out1D_y              = "no"
+IOASCII::out1D_z              = "no"
+IOASCII::out1D_d              = "no"
+IOASCII::out1D_vars           = "
+EulerAuto::Den
+EulerAuto::S1
+EulerAuto::En
+EulerAuto::rho
+EulerAuto::v1
+EulerAuto::p
+
+EulerAuto::DenLeft
+EulerAuto::S1Left
+EulerAuto::EnLeft
+EulerAuto::rhoLeft
+EulerAuto::v1Left
+EulerAuto::pLeft
+
+EulerAuto::DenRight
+EulerAuto::S1Right
+EulerAuto::EnRight
+EulerAuto::rhoRight
+EulerAuto::v1Right
+EulerAuto::pRight
+
+EulerAuto::DenFlux
+EulerAuto::S1Flux
+EulerAuto::EnFlux
+"
+
+CarpetIOASCII::out_precision = 19
+CarpetIOASCII::out3D_ghosts   = "yes"
diff --git a/Examples/eulersr_shock.par b/Examples/eulersr_shock.par
new file mode 100644
index 0000000..7d6d652
--- /dev/null
+++ b/Examples/eulersr_shock.par
@@ -0,0 +1,209 @@
+#!/usr/bin/perl -W
+
+ActiveThorns = "
+Boundary
+Carpet
+CarpetIOASCII
+CarpetIOBasic
+CarpetIOScalar
+CarpetLib
+CarpetReduce
+CarpetSlab
+CarpetIOHDF5
+CartGrid3d
+CoordBase
+Dissipation
+GenericFD
+IOUtil
+LoopControl
+MoL
+NanChecker
+EulerSR
+Slab
+SpaceMask
+SphericalSurface
+SymBase
+Time
+TimerReport
+"
+#Periodic
+
+#############################################################
+# Grid
+#############################################################
+
+CoordBase::domainsize                   = minmax
+
+CoordBase::xmin                         = -1
+CoordBase::ymin                         = 0
+CoordBase::zmin                         = 0
+
+CoordBase::xmax                         = 2
+CoordBase::ymax                         = 0.1
+CoordBase::zmax                         = 0.1
+
+CoordBase::dx                           = 0.005
+CoordBase::dy                           = 0.1
+CoordBase::dz                           = 0.1
+
+CoordBase::boundary_size_x_lower        = 2
+CoordBase::boundary_size_y_lower        = 2
+CoordBase::boundary_size_z_lower        = 2
+CoordBase::boundary_shiftout_x_lower    = 1
+CoordBase::boundary_shiftout_y_lower    = 1
+CoordBase::boundary_shiftout_z_lower    = 1
+
+CoordBase::boundary_size_x_upper        = 2
+CoordBase::boundary_size_y_upper        = 2
+CoordBase::boundary_size_z_upper        = 2
+CoordBase::boundary_shiftout_x_upper    = 0
+CoordBase::boundary_shiftout_y_upper    = 0
+CoordBase::boundary_shiftout_z_upper    = 0
+
+CartGrid3D::type                        = "coordbase"
+CartGrid3D::domain                      = "full"
+CartGrid3D::avoid_origin                = "no"
+
+#Periodic::periodic = "yes"
+
+#############################################################
+# Carpet
+#############################################################
+
+Carpet::ghost_size                      = 2
+Carpet::domain_from_coordbase           = "yes"
+Carpet::max_refinement_levels           = 1
+#Carpet::init_each_timelevel             = "yes"
+Carpet::num_integrator_substeps         = 4
+
+#############################################################
+# Time integration
+#############################################################
+
+Cactus::terminate                     = "time"
+Cactus::cctk_final_time               = 0.5
+#Cactus::cctk_itlast = 40
+
+Time::dtfac                           = 0.5
+MethodOfLines::ode_method             = "RK4"
+MethodOfLines::MoL_Intermediate_Steps = 4
+MethodOfLines::MoL_Num_Scratch_Levels = 1
+MethodOfLines::MoL_NaN_Check          = "yes"
+
+# MethodOfLines::ode_method             = "generic"
+# MethodOfLines::generic_type             = "RK"
+# MethodOfLines::MoL_Intermediate_Steps = 1
+
+Dissipation::vars = ""
+Dissipation::order = 1
+Dissipation::epsdis = 0.5
+
+EulerSR::rhoL0 = 1
+EulerSR::rhoR0 = 0.1
+EulerSR::vL0 = 0
+EulerSR::vR0 = 0
+EulerSR::epsiL0 = 2.5
+EulerSR::epsiR0 = 3.125
+EulerSR::hlleAlpha = 1
+EulerSR::gamma = 1.4
+EulerSR::initial_data = "shock"
+
+#############################################################
+# Boundary conditions
+#############################################################
+
+EulerSR::Den_group_bound   = "none"
+EulerSR::S_group_bound   = "none"
+EulerSR::tau_group_bound   = "none"
+
+Carpet::poison_new_timelevels           = yes
+Carpet::poison_value                    = 0
+CarpetLib::poison_new_memory            = yes
+CarpetLib::poison_value                 = 0
+Carpet::check_for_poison                = no
+
+#############################################################
+# Output
+#############################################################
+
+IO::out_dir                   = $parfile
+IO::out_fileinfo              = "none"
+
+CarpetIOBasic::outInfo_every        = 1
+CarpetIOBasic::outInfo_vars         = "EulerSR::rho EulerSR::v1 EulerSR::epsi EulerSR::Den EulerSR::S1 EulerSR::tau"
+
+CarpetIOScalar::outScalar_every = 1
+CarpetIOScalar::outScalar_vars = ""
+CarpetIOScalar::outScalar_reductions = "minimum maximum norm2"
+
+IOHDF5::out1D_every                    = 40
+IOHDF5::out1D_d                        = no
+IOHDF5::out1D_x                        = yes
+IOHDF5::out1D_y                        = no
+IOHDF5::out1D_z                        = no
+IOHDF5::out1D_vars                     = "
+EulerSR::Den
+EulerSR::S1
+EulerSR::tau
+EulerSR::rho
+EulerSR::v1
+EulerSR::epsi
+EulerSR::p
+
+EulerSR::DenLeft
+EulerSR::S1Left
+EulerSR::tauLeft
+EulerSR::rhoLeft
+EulerSR::v1Left
+EulerSR::epsiLeft
+
+EulerSR::DenRight
+EulerSR::S1Right
+EulerSR::tauRight
+EulerSR::rhoRight
+EulerSR::v1Right
+EulerSR::epsiRight
+
+EulerSR::DenFlux
+EulerSR::S1Flux
+EulerSR::tauFlux
+"
+
+# IOASCII::out1D_every          = 1
+# IOASCII::out1D_x              = "yes"
+# IOASCII::out1D_y              = "no"
+# IOASCII::out1D_z              = "no"
+# IOASCII::out1D_d              = "no"
+# IOASCII::out1D_vars           = "
+# EulerSR::Den
+# EulerSR::S1
+# EulerSR::tau
+# EulerSR::rho
+# EulerSR::v1
+# EulerSR::epsi
+# EulerSR::p
+
+# EulerSR::DenLeft
+# EulerSR::S1Left
+# EulerSR::tauLeft
+# EulerSR::rhoLeft
+# EulerSR::v1Left
+# EulerSR::epsiLeft
+
+# EulerSR::DenRight
+# EulerSR::S1Right
+# EulerSR::tauRight
+# EulerSR::rhoRight
+# EulerSR::v1Right
+# EulerSR::epsiRight
+
+# EulerSR::DenFlux
+# EulerSR::S1Flux
+# EulerSR::tauFlux
+# "
+
+CarpetIOASCII::out_precision = 19
+CarpetIOASCII::out3D_ghosts   = "yes"
+
+TimerReport::n_top_timers = 40
+TimerReport::out_every = 200
diff --git a/Examples/eulersr_shock_diag.par b/Examples/eulersr_shock_diag.par
new file mode 100644
index 0000000..260b52d
--- /dev/null
+++ b/Examples/eulersr_shock_diag.par
@@ -0,0 +1,192 @@
+#!/usr/bin/perl -W
+
+ActiveThorns = "
+Boundary
+Carpet
+CarpetIOASCII
+CarpetIOBasic
+CarpetIOScalar
+CarpetLib
+CarpetReduce
+CarpetSlab
+CarpetIOHDF5
+CartGrid3d
+CoordBase
+Dissipation
+GenericFD
+IOUtil
+LoopControl
+MoL
+NanChecker
+EulerSR
+Slab
+SpaceMask
+SphericalSurface
+SymBase
+Time
+TimerReport
+Periodic
+"
+
+#############################################################
+# Grid
+#############################################################
+
+CoordBase::domainsize                   = minmax
+
+CoordBase::xmin                         = -1
+CoordBase::ymin                         = -1
+CoordBase::zmin                         = -1
+
+CoordBase::xmax                         = 1
+CoordBase::ymax                         = 1
+CoordBase::zmax                         = 1
+
+CoordBase::dx                           = 0.05
+CoordBase::dy                           = 0.05
+CoordBase::dz                           = 0.05
+
+CoordBase::boundary_size_x_lower        = 2
+CoordBase::boundary_size_y_lower        = 2
+CoordBase::boundary_size_z_lower        = 2
+CoordBase::boundary_shiftout_x_lower    = 1
+CoordBase::boundary_shiftout_y_lower    = 1
+CoordBase::boundary_shiftout_z_lower    = 1
+
+CoordBase::boundary_size_x_upper        = 2
+CoordBase::boundary_size_y_upper        = 2
+CoordBase::boundary_size_z_upper        = 2
+CoordBase::boundary_shiftout_x_upper    = 0
+CoordBase::boundary_shiftout_y_upper    = 0
+CoordBase::boundary_shiftout_z_upper    = 0
+
+CartGrid3D::type                        = "coordbase"
+CartGrid3D::domain                      = "full"
+CartGrid3D::avoid_origin                = "no"
+
+Periodic::periodic = "no"
+
+#############################################################
+# Carpet
+#############################################################
+
+Carpet::ghost_size                      = 2
+Carpet::domain_from_coordbase           = "yes"
+Carpet::max_refinement_levels           = 1
+#Carpet::init_each_timelevel             = "yes"
+Carpet::num_integrator_substeps         = 4
+
+#############################################################
+# Time integration
+#############################################################
+
+Cactus::terminate                     = "time"
+Cactus::cctk_final_time               = 2
+#Cactus::cctk_itlast = 40
+
+Time::dtfac                           = 0.25
+MethodOfLines::ode_method             = "RK4"
+MethodOfLines::MoL_Intermediate_Steps = 4
+MethodOfLines::MoL_Num_Scratch_Levels = 1
+MethodOfLines::MoL_NaN_Check          = "yes"
+
+# MethodOfLines::ode_method             = "generic"
+# MethodOfLines::generic_type             = "RK"
+# MethodOfLines::MoL_Intermediate_Steps = 1
+
+Dissipation::vars = ""
+Dissipation::order = 1
+Dissipation::epsdis = 0.5
+
+EulerSR::rhoL0 = 1
+EulerSR::rhoR0 = 0.1
+EulerSR::vL0 = 0
+EulerSR::vR0 = 0
+EulerSR::epsiL0 = 2.5
+EulerSR::epsiR0 = 3.125
+EulerSR::hlleAlpha = 1
+EulerSR::gamma = 1.4
+EulerSR::initial_data = "shock"
+
+#############################################################
+# Boundary conditions
+#############################################################
+
+EulerSR::Den_group_bound   = "none"
+EulerSR::S_group_bound   = "none"
+EulerSR::tau_group_bound   = "none"
+
+Carpet::poison_new_timelevels           = yes
+Carpet::poison_value                    = 0
+CarpetLib::poison_new_memory            = yes
+CarpetLib::poison_value                 = 0
+Carpet::check_for_poison                = no
+
+#############################################################
+# Output
+#############################################################
+
+IO::out_dir                   = $parfile
+IO::out_fileinfo              = "none"
+
+CarpetIOBasic::outInfo_every        = 1
+CarpetIOBasic::outInfo_vars         = "EulerSR::rho EulerSR::v1 EulerSR::epsi EulerSR::Den EulerSR::S1 EulerSR::tau"
+
+CarpetIOScalar::outScalar_every = 1
+CarpetIOScalar::outScalar_vars = ""
+CarpetIOScalar::outScalar_reductions = "minimum maximum norm2"
+
+IOHDF5::out1D_every                    = 4
+IOHDF5::out1D_d                        = no
+IOHDF5::out1D_x                        = yes
+IOHDF5::out1D_y                        = yes
+IOHDF5::out1D_z                        = yes
+IOHDF5::out1D_vars                     = "
+EulerSR::Den
+EulerSR::S1
+EulerSR::tau
+EulerSR::rho
+EulerSR::v1
+EulerSR::v2
+EulerSR::v3
+EulerSR::epsi
+EulerSR::p
+
+EulerSR::DenLeft
+EulerSR::S1Left
+EulerSR::tauLeft
+EulerSR::rhoLeft
+EulerSR::v1Left
+EulerSR::epsiLeft
+
+EulerSR::DenRight
+EulerSR::S1Right
+EulerSR::tauRight
+EulerSR::rhoRight
+EulerSR::v1Right
+EulerSR::epsiRight
+
+EulerSR::DenFlux
+EulerSR::S1Flux
+EulerSR::tauFlux
+"
+
+IOHDF5::out2D_every                    = 1
+IOHDF5::out2D_xy                       = yes
+IOHDF5::out2D_vars                     = "
+EulerSR::Den
+EulerSR::S1
+EulerSR::tau
+EulerSR::rho
+EulerSR::v1
+EulerSR::v2
+EulerSR::v3
+EulerSR::epsi
+EulerSR::p
+"
+
+CarpetIOASCII::out_precision = 19
+CarpetIOASCII::out3D_ghosts   = "yes"
+
+TimerReport::n_top_timers = 40
+TimerReport::out_every = 200