update balsara par files

From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu>

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@146 ac85fae7-cede-4708-beff-ae01c7fa1c26

4 files changed, 1536 insertions, 0 deletions

diff --git a/par/balsara2_2d.par b/par/balsara2_2d.par
new file mode 100644
index 0000000..91068d1
--- /dev/null
+++ b/par/balsara2_2d.par
@@ -0,0 +1,384 @@
+#==============================================================================
+# Balsara magnetized shock tube test #2: GRHydro_InitData
+#==============================================================================
+
+#------------------------------------------------------------------------------
+# Cactus parameters:
+#------------------------------------------------------------------------------
+Cactus::cctk_run_title = "Balsara magnetized shock tube test #2 along z=0 plane square diagonal"
+Cactus::cctk_full_warnings         = "yes"
+
+#Cactus::terminate		= "never"
+Cactus::terminate       = "time"
+Cactus::cctk_final_time = 0.4
+#Cactus::terminate = "iteration"
+#Cactus::cctk_itlast        =  0
+#Cactus::cctk_itlast        =  100
+#Cactus::cctk_itlast        =  256
+#Cactus::cctk_itlast        = 1024
+#Cactus::terminate="runtime"
+#Cactus::max_runtime=28800  #  46 days
+#Cactus::max_runtime=1410  #  23.5 hours
+#Cactus::max_runtime=1380  #  23 hours
+#Cactus::max_runtime=690  #  11.5 hours
+#Cactus::max_runtime=450  #  7.5 hours
+#Cactus::max_runtime=210  #  3.5 hours
+
+#------------------------------------------------------------------------------
+# Activate all necessary thorns:
+#------------------------------------------------------------------------------
+
+ActiveThorns = "Boundary CartGrid3D CoordBase Fortran InitBase IOUtil LocalReduce SymBase Time"
+ActiveThorns = "AEILocalInterp LocalInterp"
+ActiveThorns = "MoL Slab SpaceMask"
+ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetSlab LoopControl"
+#ActiveThorns = "Formaline NaNChecker TerminationTrigger TimerReport"
+ActiveThorns = "NaNChecker TerminationTrigger TimerReport"
+ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal"
+ActiveThorns = "Constants TmunuBase HydroBase "
+#ActiveThorns = "EOS_Base EOS_IdealFluid EOS_Polytrope"
+ActiveThorns = "EOS_Omni"
+ActiveThorns = "GRHydro GRHydro_InitData"
+ActiveThorns = "HDF5"
+ActiveThorns = "SystemStatistics"
+
+#------------------------------------------------------------------------------
+# Diagnostic parameters:
+#------------------------------------------------------------------------------
+Carpet::output_timers_every = 1024
+Carpet::output_timer_tree_every = 256
+Carpet::output_initialise_timer_tree = "yes"
+Carpet::storage_verbose   = "no"
+#Carpet::timers_verbose    = "no"
+#Carpet::verbose           = "yes"  # only for Initial Data
+Carpet::verbose           = "no"
+Carpet::veryverbose       = "no"
+Carpet::grid_structure_filename = "carpet-grid-structure.asc"
+Carpet::grid_coordinates_filename = "carpet-grid-coordinates.asc"
+
+CarpetLib::output_bboxes  = "no"
+CarpetLib::print_memstats_every     = 1024
+CarpetLib::memstat_file    = "carpetlib-memory-statistics.asc"
+#CarpetLib::print_timestats_every    = 1024
+#CarpetLib::timestat_file    = "carpetlib-timing-statistics.asc"
+
+CarpetReduce::verbose    = "no"
+#CarpetRegrid2::verbose    = "yes"
+#CarpetRegrid2::veryverbose    = "yes"
+
+#NaNChecker::verbose         = "all"
+
+TimerReport::out_every    = 1024
+TimerReport::out_filename = "TimerReport"
+#TimerReport::output_all_timers = "yes"
+#TimerReport::output_all_timers_together = "yes"
+#TimerReport::output_all_timers_readable = "yes"
+#TimerReport::before_checkpoint          = "yes"
+TimerReport::n_top_timers = 40
+
+#------------------------------------------------------------------------------
+# Utility parameters:
+#------------------------------------------------------------------------------
+
+#NaNChecker::check_every    =  256
+NaNChecker::check_every   =  1
+NaNChecker::check_vars = "
+            ADMBase::curv 
+            ADMBase::metric 
+            ADMBase::lapse 
+            ADMBase::shift 
+            HydroBase::Bvec
+            HydroBase::rho 
+            HydroBase::eps 
+            HydroBase::press 
+            HydroBase::vel
+            HydroBase::w_lorentz
+"
+NaNChecker::action_if_found   =  "terminate"
+#NaNChecker::action_if_found  =  "abort"
+#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort"
+#NaNChecker::check_after=0
+
+#TerminationTrigger::max_walltime          =  0.0   # hours
+TerminationTrigger::max_walltime          =  24.0   # hours
+TerminationTrigger::on_remaining_walltime = 10.0   # minutes
+TerminationTrigger::termination_from_file   = "yes"
+TerminationTrigger::create_termination_file = "yes"
+TerminationTrigger::termination_file        = "../TERMINATE"
+
+
+#------------------------------------------------------------------------------
+# Run parameters:
+#------------------------------------------------------------------------------
+
+#------
+# Grid:
+#------
+
+#Time::dtfac = 0.8
+Time::dtfac = 0.4
+MoL::ODE_Method             = "rk2"
+MoL::MoL_Intermediate_Steps = 2
+#MoL::ODE_Method             = "rk4"
+#MoL::MoL_Intermediate_Steps = 4
+#MoL::MoL_Num_Scratch_Levels = 1
+
+CartGrid3D::type         = "coordbase"
+CartGrid3D::domain       = "full"
+CartGrid3D::avoid_origin = "no"
+
+CoordBase::domainsize       = "minmax"
+CoordBase::spacing          = "numcells"  # "gridspacing" or "numcells"
+
+#CoordBase::xmin = -0.035355339059327376220
+#CoordBase::ymin = -0.035355339059327376220
+#CoordBase::zmin = -0.0017677669529663688110
+#CoordBase::xmax =  0.035355339059327376220
+#CoordBase::ymax =  0.035355339059327376220
+#CoordBase::zmax =  0.0017677669529663688110
+#CoordBase::ncells_x   = 160  # dx or ncells_x
+#CoordBase::ncells_y   = 160  # dy or ncells_y
+#CoordBase::ncells_z   =   8  # dz or ncells_z
+
+# Approximate memory footprint: 650MB x 8 = 5.2GB
+
+CoordBase::xmin = -0.35355339059327376220
+CoordBase::ymin = -0.35355339059327376220
+CoordBase::zmin = -0.003535533905932737622
+CoordBase::xmax =  0.35355339059327376220
+CoordBase::ymax =  0.35355339059327376220
+CoordBase::zmax =  0.003535533905932737622
+CoordBase::ncells_x   = 800  # dx or ncells_x
+CoordBase::ncells_y   = 800  # dy or ncells_y
+CoordBase::ncells_z   =   8  # dz or ncells_z
+
+Driver::ghost_size                  = 3
+GRHydro::GRHydro_stencil            = 3
+
+# General Carpet parameters:
+Carpet::domain_from_coordbase    = "yes"
+
+CarpetLib::poison_new_memory         = "yes"
+CarpetLib::poison_value              = 114
+
+#------
+# MODEL:
+#------
+
+HydroBase::initial_hydro                = "shocktube"
+GRHydro_InitData::shocktube_type        = "diagshock2d"
+GRHydro_InitData::shock_case            = "Balsara2"
+
+ADMBase::metric_type             = "physical"
+ADMBase::initial_data            = "Cartesian Minkowski"
+ADMBase::initial_lapse           = "one"
+ADMBase::initial_shift           = "zero"
+ADMBase::initial_dtlapse         = "none"
+ADMBase::initial_dtshift         = "none"
+
+EOS_Omni::poly_K     = 100.0
+EOS_Omni::poly_gamma = 1.6666666666666666666
+EOS_Omni::gl_K       = 100.0
+EOS_Omni::gl_gamma   = 1.6666666666666666666
+GRHydro::initial_Gamma = 1.6666666666666666666
+
+#----------
+# Numerics:
+#----------
+#Boundary::radpower                     = 2
+#HydroBase::timelevels            = 3
+HydroBase::timelevels            = 2
+#ADMMacros::spatial_order = 4
+SpaceMask::use_mask      = "yes"
+
+#-----------
+# Evolution:
+#-----------
+
+HydroBase::evolution_method      = "GRHydro"
+HydroBase::Bvec_evolution_method = "GRHydro"
+
+GRHydro::method_type               = "RSA FV"
+GRHydro::recon_method              = "tvd"
+#GRHydro::recon_method              = "ppm"
+#GRHydro::recon_method              = "eno"
+GRHydro::tvd_limiter               = "vanleerMC2"
+#GRHydro::tvd_limiter               = "minmod"
+GRHydro::recon_vars                = "primitive"
+GRHydro::riemann_solver            = "HLLE"
+GRHydro::decouple_normal_Bfield    = "yes"
+#GRHydro::riemann_solver            = "LLF"
+GRHydro::bound                     = "none"
+GRHydro::GRHydro_eos_type          = "General"
+GRHydro::GRHydro_eos_table         = "Ideal_Fluid"
+
+GRHydro::Grhydro_MaxNumConstrainedVars = 13
+GRHydro::GRHydro_MaxNumEvolvedVars     = 10
+
+GRHydro::clean_divergence          = "yes"
+GRHydro::psidcspeed                = "light speed"
+GRHydro::kap_dc                    = 10.0
+GRHydro::track_divB                = "yes"
+GRHydro::calculate_bcom            = "yes" 
+
+#------------------------------------------------------------------------------
+# Output:
+#------------------------------------------------------------------------------
+
+IO::out_dir = $parfile
+IO::out_fileinfo           = "none"
+
+IOBasic::outInfo_every = 1
+IOBasic::outInfo_reductions = "maximum"
+IOBasic::outInfo_vars  = "
+ Carpet::physical_time_per_hour
+ GRHydro::bcom_sq
+ HydroBase::rho
+ HydroBase::press
+ SystemStatistics::maxrss_mb
+"
+
+#IOScalar::outScalar_every      = 1024
+#IOScalar::one_file_per_group   = "yes"
+#IOScalar::outScalar_reductions = "minimum maximum average norm1 norm2"
+#IOScalar::outScalar_vars       = "
+# GRHydro::bcom
+# GRHydro::bcom0
+# GRHydro::bcom_sq
+# GRHydro::Bcons
+# GRHydro::dens
+# GrHydro::divB
+# GRHydro::scon
+# GRHydro::tau
+# HydroBase::Bvec
+# HydroBase::eps
+# HydroBase::press
+# HydroBase::rho
+# HydroBase::vel
+# HydroBase::w_lorentz
+# SystemStatistics::process_memory_mb
+#"
+
+IOASCII::one_file_per_group     = "yes"
+IOASCII::output_symmetry_points = "no"
+IOASCII::output_ghost_points    = "no"  #hg version only!
+IOASCII::output_boundary_points = "no"  #hg version only!
+#IOASCII::out3D_ghosts       = "no"  #hg (DEPRECATED)
+#IOASCII::out3D_outer_ghosts = "no" #hg (DEPRECATED)
+
+IOASCII::out0D_every     = 200
+IOASCII::out0D_vars      = "
+ Carpet::timing
+"
+
+#Set these IOASCII options for initial data only:
+#IOASCII::out1D_every     = 1
+IOASCII::out1D_every     = 1131
+#IOASCII::out1D_dt        = 0.01
+IOASCII::out1D_d         = "yes"
+IOASCII::out1D_vars      = "
+ GRHydro::bcom
+ GRHydro::bcom0
+ GRHydro::bcom_sq
+ GRHydro::Bcons
+ GRHydro::dens
+ GrHydro::divB
+ GrHydro::psidc
+ GRHydro::scon
+ GRHydro::tau
+ HydroBase::Bvec
+ HydroBase::eps
+ HydroBase::press
+ HydroBase::rho
+ HydroBase::vel
+ HydroBase::w_lorentz
+"
+
+IOASCII::out2D_every     = 1131
+IOASCII::out2D_vars      = "
+ GRHydro::bcom
+ GRHydro::bcom0
+ GRHydro::bcom_sq
+ GRHydro::Bcons
+ GRHydro::dens
+ GrHydro::divB
+ GrHydro::psidc
+ GRHydro::scon
+ GRHydro::tau
+ HydroBase::Bvec
+ HydroBase::eps
+ HydroBase::press
+ HydroBase::rho
+ HydroBase::vel
+ HydroBase::w_lorentz
+"
+
+CarpetIOHDF5::one_file_per_group             = "yes"
+CarpetIOHDF5::open_one_input_file_at_a_time  = "yes"
+
+#CarpetIOHDF5::output_symmetry_points         = "no"
+#CarpetIOHDF5::output_ghost_points            = "no"
+#CarpetIOHDF5::output_boundary_points         = "no"
+#CarpetIOHDF5::output_buffer_points           = "no"
+
+#CarpetIOHDF5::out2D_every                    = 1
+CarpetIOHDF5::out2D_every                    = 100
+#CarpetIOHDF5::out2D_dt                       = 0.01
+CarpetIOHDF5::out2D_xy                       = "yes"
+CarpetIOHDF5::out2D_xz                       = "no"
+CarpetIOHDF5::out2D_yz                       = "no"
+CarpetIOHDF5::out2D_xyplane_z                = 0.0
+CarpetIOHDF5::out2D_vars      = "
+  GRHydro::bcom
+  GRHydro::bcom0
+  GRHydro::bcom_sq
+  GRHydro::Bcons
+  GRHydro::dens
+  GRHydro::divB
+  GrHydro::psidc
+  GRHydro::scon
+  GRHydro::tau
+  HydroBase::Bvec
+  HydroBase::eps
+  HydroBase::press
+  HydroBase::rho
+  HydroBase::vel
+  HydroBase::w_lorentz
+"
+
+#CarpetIOHDF5::out3D_every = 1
+CarpetIOHDF5::out3D_every = 100
+#CarpetIOHDF5::out3D_dt = 0.01
+CarpetIOHDF5::out3D_vars  = "
+ GRHydro::bcom
+ GRHydro::bcom0
+ GRHydro::bcom_sq
+ GRHydro::divB
+ HydroBase::Bvec
+ HydroBase::eps
+ HydroBase::press
+ HydroBase::rho
+ HydroBase::vel
+ HydroBase::w_lorentz
+"
+
+#------------------------------------------------------------------------------
+# Analysis:
+#------------------------------------------------------------------------------
+
+
+#------------------------------------------------------------------------------
+# Checkpoint/Recovery:
+#------------------------------------------------------------------------------
+#IOHDF5::checkpoint                  = "yes"
+#IO::checkpoint_dir                  = $parfile
+#IO::checkpoint_ID                   = "yes"
+#IO::checkpoint_every                = 6144
+##IO::checkpoint_every_walltime_hours = 6.0
+#IO::checkpoint_keep=3
+#IO::checkpoint_on_terminate         = "yes"
+#
+#IO::recover     = "autoprobe"
+#IO::recover_dir = $parfile
+
+
diff --git a/par/balsara3_2d.par b/par/balsara3_2d.par
new file mode 100644
index 0000000..3d455bf
--- /dev/null
+++ b/par/balsara3_2d.par
@@ -0,0 +1,384 @@
+#==============================================================================
+# Balsara magnetized shock tube test #3: GRHydro_InitData
+#==============================================================================
+
+#------------------------------------------------------------------------------
+# Cactus parameters:
+#------------------------------------------------------------------------------
+Cactus::cctk_run_title = "Balsara magnetized shock tube test #3 along z=0 plane square diagonal"
+Cactus::cctk_full_warnings         = "yes"
+
+#Cactus::terminate		= "never"
+Cactus::terminate       = "time"
+Cactus::cctk_final_time = 0.4
+#Cactus::terminate = "iteration"
+#Cactus::cctk_itlast        =  0
+#Cactus::cctk_itlast        =  100
+#Cactus::cctk_itlast        =  256
+#Cactus::cctk_itlast        = 1024
+#Cactus::terminate="runtime"
+#Cactus::max_runtime=28800  #  46 days
+#Cactus::max_runtime=1410  #  23.5 hours
+#Cactus::max_runtime=1380  #  23 hours
+#Cactus::max_runtime=690  #  11.5 hours
+#Cactus::max_runtime=450  #  7.5 hours
+#Cactus::max_runtime=210  #  3.5 hours
+
+#------------------------------------------------------------------------------
+# Activate all necessary thorns:
+#------------------------------------------------------------------------------
+
+ActiveThorns = "Boundary CartGrid3D CoordBase Fortran InitBase IOUtil LocalReduce SymBase Time"
+ActiveThorns = "AEILocalInterp LocalInterp"
+ActiveThorns = "MoL Slab SpaceMask"
+ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetSlab LoopControl"
+#ActiveThorns = "Formaline NaNChecker TerminationTrigger TimerReport"
+ActiveThorns = "NaNChecker TerminationTrigger TimerReport"
+ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal"
+ActiveThorns = "Constants TmunuBase HydroBase "
+#ActiveThorns = "EOS_Base EOS_IdealFluid EOS_Polytrope"
+ActiveThorns = "EOS_Omni"
+ActiveThorns = "GRHydro GRHydro_InitData"
+ActiveThorns = "HDF5"
+ActiveThorns = "SystemStatistics"
+
+#------------------------------------------------------------------------------
+# Diagnostic parameters:
+#------------------------------------------------------------------------------
+Carpet::output_timers_every = 1024
+Carpet::output_timer_tree_every = 256
+Carpet::output_initialise_timer_tree = "yes"
+Carpet::storage_verbose   = "no"
+#Carpet::timers_verbose    = "no"
+#Carpet::verbose           = "yes"  # only for Initial Data
+Carpet::verbose           = "no"
+Carpet::veryverbose       = "no"
+Carpet::grid_structure_filename = "carpet-grid-structure.asc"
+Carpet::grid_coordinates_filename = "carpet-grid-coordinates.asc"
+
+CarpetLib::output_bboxes  = "no"
+CarpetLib::print_memstats_every     = 1024
+CarpetLib::memstat_file    = "carpetlib-memory-statistics.asc"
+#CarpetLib::print_timestats_every    = 1024
+#CarpetLib::timestat_file    = "carpetlib-timing-statistics.asc"
+
+CarpetReduce::verbose    = "no"
+#CarpetRegrid2::verbose    = "yes"
+#CarpetRegrid2::veryverbose    = "yes"
+
+#NaNChecker::verbose         = "all"
+
+TimerReport::out_every    = 1024
+TimerReport::out_filename = "TimerReport"
+#TimerReport::output_all_timers = "yes"
+#TimerReport::output_all_timers_together = "yes"
+#TimerReport::output_all_timers_readable = "yes"
+#TimerReport::before_checkpoint          = "yes"
+TimerReport::n_top_timers = 40
+
+#------------------------------------------------------------------------------
+# Utility parameters:
+#------------------------------------------------------------------------------
+
+#NaNChecker::check_every    =  256
+NaNChecker::check_every   =  1
+NaNChecker::check_vars = "
+            ADMBase::curv 
+            ADMBase::metric 
+            ADMBase::lapse 
+            ADMBase::shift 
+            HydroBase::Bvec
+            HydroBase::rho 
+            HydroBase::eps 
+            HydroBase::press 
+            HydroBase::vel
+            HydroBase::w_lorentz
+"
+NaNChecker::action_if_found   =  "terminate"
+#NaNChecker::action_if_found  =  "abort"
+#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort"
+#NaNChecker::check_after=0
+
+#TerminationTrigger::max_walltime          =  0.0   # hours
+TerminationTrigger::max_walltime          =  24.0   # hours
+TerminationTrigger::on_remaining_walltime = 10.0   # minutes
+TerminationTrigger::termination_from_file   = "yes"
+TerminationTrigger::create_termination_file = "yes"
+TerminationTrigger::termination_file        = "../TERMINATE"
+
+
+#------------------------------------------------------------------------------
+# Run parameters:
+#------------------------------------------------------------------------------
+
+#------
+# Grid:
+#------
+
+#Time::dtfac = 0.8
+Time::dtfac = 0.4
+MoL::ODE_Method             = "rk2"
+MoL::MoL_Intermediate_Steps = 2
+#MoL::ODE_Method             = "rk4"
+#MoL::MoL_Intermediate_Steps = 4
+#MoL::MoL_Num_Scratch_Levels = 1
+
+CartGrid3D::type         = "coordbase"
+CartGrid3D::domain       = "full"
+CartGrid3D::avoid_origin = "no"
+
+CoordBase::domainsize       = "minmax"
+CoordBase::spacing          = "numcells"  # "gridspacing" or "numcells"
+
+#CoordBase::xmin = -0.035355339059327376220
+#CoordBase::ymin = -0.035355339059327376220
+#CoordBase::zmin = -0.0017677669529663688110
+#CoordBase::xmax =  0.035355339059327376220
+#CoordBase::ymax =  0.035355339059327376220
+#CoordBase::zmax =  0.0017677669529663688110
+#CoordBase::ncells_x   = 160  # dx or ncells_x
+#CoordBase::ncells_y   = 160  # dy or ncells_y
+#CoordBase::ncells_z   =   8  # dz or ncells_z
+
+# Approximate memory footprint: 650MB x 8 = 5.2GB
+
+CoordBase::xmin = -0.35355339059327376220
+CoordBase::ymin = -0.35355339059327376220
+CoordBase::zmin = -0.003535533905932737622
+CoordBase::xmax =  0.35355339059327376220
+CoordBase::ymax =  0.35355339059327376220
+CoordBase::zmax =  0.003535533905932737622
+CoordBase::ncells_x   = 800  # dx or ncells_x
+CoordBase::ncells_y   = 800  # dy or ncells_y
+CoordBase::ncells_z   =   8  # dz or ncells_z
+
+Driver::ghost_size                  = 3
+GRHydro::GRHydro_stencil            = 3
+
+# General Carpet parameters:
+Carpet::domain_from_coordbase    = "yes"
+
+CarpetLib::poison_new_memory         = "yes"
+CarpetLib::poison_value              = 114
+
+#------
+# MODEL:
+#------
+
+HydroBase::initial_hydro                = "shocktube"
+GRHydro_InitData::shocktube_type        = "diagshock2d"
+GRHydro_InitData::shock_case            = "Balsara3"
+
+ADMBase::metric_type             = "physical"
+ADMBase::initial_data            = "Cartesian Minkowski"
+ADMBase::initial_lapse           = "one"
+ADMBase::initial_shift           = "zero"
+ADMBase::initial_dtlapse         = "none"
+ADMBase::initial_dtshift         = "none"
+
+EOS_Omni::poly_K     = 100.0
+EOS_Omni::poly_gamma = 1.6666666666666666666
+EOS_Omni::gl_K       = 100.0
+EOS_Omni::gl_gamma   = 1.6666666666666666666
+GRHydro::initial_Gamma = 1.6666666666666666666
+
+#----------
+# Numerics:
+#----------
+#Boundary::radpower                     = 2
+#HydroBase::timelevels            = 3
+HydroBase::timelevels            = 2
+#ADMMacros::spatial_order = 4
+SpaceMask::use_mask      = "yes"
+
+#-----------
+# Evolution:
+#-----------
+
+HydroBase::evolution_method      = "GRHydro"
+HydroBase::Bvec_evolution_method = "GRHydro"
+
+GRHydro::method_type               = "RSA FV"
+GRHydro::recon_method              = "tvd"
+#GRHydro::recon_method              = "ppm"
+#GRHydro::recon_method              = "eno"
+GRHydro::tvd_limiter               = "vanleerMC2"
+#GRHydro::tvd_limiter               = "minmod"
+GRHydro::recon_vars                = "primitive"
+GRHydro::riemann_solver            = "HLLE"
+GRHydro::decouple_normal_Bfield    = "yes"
+#GRHydro::riemann_solver            = "LLF"
+GRHydro::bound                     = "none"
+GRHydro::GRHydro_eos_type          = "General"
+GRHydro::GRHydro_eos_table         = "Ideal_Fluid"
+
+GRHydro::Grhydro_MaxNumConstrainedVars = 13
+GRHydro::GRHydro_MaxNumEvolvedVars     = 10
+
+GRHydro::clean_divergence          = "yes"
+GRHydro::psidcspeed                = "light speed"
+GRHydro::kap_dc                    = 10.0
+GRHydro::track_divB                = "yes"
+GRHydro::calculate_bcom            = "yes" 
+
+#------------------------------------------------------------------------------
+# Output:
+#------------------------------------------------------------------------------
+
+IO::out_dir = $parfile
+IO::out_fileinfo           = "none"
+
+IOBasic::outInfo_every = 1
+IOBasic::outInfo_reductions = "maximum"
+IOBasic::outInfo_vars  = "
+ Carpet::physical_time_per_hour
+ GRHydro::bcom_sq
+ HydroBase::rho
+ HydroBase::press
+ SystemStatistics::maxrss_mb
+"
+
+#IOScalar::outScalar_every      = 1024
+#IOScalar::one_file_per_group   = "yes"
+#IOScalar::outScalar_reductions = "minimum maximum average norm1 norm2"
+#IOScalar::outScalar_vars       = "
+# GRHydro::bcom
+# GRHydro::bcom0
+# GRHydro::bcom_sq
+# GRHydro::Bcons
+# GRHydro::dens
+# GrHydro::divB
+# GRHydro::scon
+# GRHydro::tau
+# HydroBase::Bvec
+# HydroBase::eps
+# HydroBase::press
+# HydroBase::rho
+# HydroBase::vel
+# HydroBase::w_lorentz
+# SystemStatistics::process_memory_mb
+#"
+
+IOASCII::one_file_per_group     = "yes"
+IOASCII::output_symmetry_points = "no"
+IOASCII::output_ghost_points    = "no"  #hg version only!
+IOASCII::output_boundary_points = "no"  #hg version only!
+#IOASCII::out3D_ghosts       = "no"  #hg (DEPRECATED)
+#IOASCII::out3D_outer_ghosts = "no" #hg (DEPRECATED)
+
+IOASCII::out0D_every     = 200
+IOASCII::out0D_vars      = "
+ Carpet::timing
+"
+
+#Set these IOASCII options for initial data only:
+#IOASCII::out1D_every     = 1
+IOASCII::out1D_every     = 1131
+#IOASCII::out1D_dt        = 0.01
+IOASCII::out1D_d         = "yes"
+IOASCII::out1D_vars      = "
+ GRHydro::bcom
+ GRHydro::bcom0
+ GRHydro::bcom_sq
+ GRHydro::Bcons
+ GRHydro::dens
+ GrHydro::divB
+ GrHydro::psidc
+ GRHydro::scon
+ GRHydro::tau
+ HydroBase::Bvec
+ HydroBase::eps
+ HydroBase::press
+ HydroBase::rho
+ HydroBase::vel
+ HydroBase::w_lorentz
+"
+
+IOASCII::out2D_every     = 1131
+IOASCII::out2D_vars      = "
+ GRHydro::bcom
+ GRHydro::bcom0
+ GRHydro::bcom_sq
+ GRHydro::Bcons
+ GRHydro::dens
+ GrHydro::divB
+ GrHydro::psidc
+ GRHydro::scon
+ GRHydro::tau
+ HydroBase::Bvec
+ HydroBase::eps
+ HydroBase::press
+ HydroBase::rho
+ HydroBase::vel
+ HydroBase::w_lorentz
+"
+
+CarpetIOHDF5::one_file_per_group             = "yes"
+CarpetIOHDF5::open_one_input_file_at_a_time  = "yes"
+
+#CarpetIOHDF5::output_symmetry_points         = "no"
+#CarpetIOHDF5::output_ghost_points            = "no"
+#CarpetIOHDF5::output_boundary_points         = "no"
+#CarpetIOHDF5::output_buffer_points           = "no"
+
+#CarpetIOHDF5::out2D_every                    = 1
+CarpetIOHDF5::out2D_every                    = 100
+#CarpetIOHDF5::out2D_dt                       = 0.01
+CarpetIOHDF5::out2D_xy                       = "yes"
+CarpetIOHDF5::out2D_xz                       = "no"
+CarpetIOHDF5::out2D_yz                       = "no"
+CarpetIOHDF5::out2D_xyplane_z                = 0.0
+CarpetIOHDF5::out2D_vars      = "
+  GRHydro::bcom
+  GRHydro::bcom0
+  GRHydro::bcom_sq
+  GRHydro::Bcons
+  GRHydro::dens
+  GRHydro::divB
+  GrHydro::psidc
+  GRHydro::scon
+  GRHydro::tau
+  HydroBase::Bvec
+  HydroBase::eps
+  HydroBase::press
+  HydroBase::rho
+  HydroBase::vel
+  HydroBase::w_lorentz
+"
+
+#CarpetIOHDF5::out3D_every = 1
+CarpetIOHDF5::out3D_every = 100
+#CarpetIOHDF5::out3D_dt = 0.01
+CarpetIOHDF5::out3D_vars  = "
+ GRHydro::bcom
+ GRHydro::bcom0
+ GRHydro::bcom_sq
+ GRHydro::divB
+ HydroBase::Bvec
+ HydroBase::eps
+ HydroBase::press
+ HydroBase::rho
+ HydroBase::vel
+ HydroBase::w_lorentz
+"
+
+#------------------------------------------------------------------------------
+# Analysis:
+#------------------------------------------------------------------------------
+
+
+#------------------------------------------------------------------------------
+# Checkpoint/Recovery:
+#------------------------------------------------------------------------------
+#IOHDF5::checkpoint                  = "yes"
+#IO::checkpoint_dir                  = $parfile
+#IO::checkpoint_ID                   = "yes"
+#IO::checkpoint_every                = 6144
+##IO::checkpoint_every_walltime_hours = 6.0
+#IO::checkpoint_keep=3
+#IO::checkpoint_on_terminate         = "yes"
+#
+#IO::recover     = "autoprobe"
+#IO::recover_dir = $parfile
+
+
diff --git a/par/balsara4_2d.par b/par/balsara4_2d.par
new file mode 100644
index 0000000..b7b8c17
--- /dev/null
+++ b/par/balsara4_2d.par
@@ -0,0 +1,384 @@
+#==============================================================================
+# Balsara magnetized shock tube test #4: GRHydro_InitData
+#==============================================================================
+
+#------------------------------------------------------------------------------
+# Cactus parameters:
+#------------------------------------------------------------------------------
+Cactus::cctk_run_title = "Balsara magnetized shock tube test #4 along z=0 plane square diagonal"
+Cactus::cctk_full_warnings         = "yes"
+
+#Cactus::terminate		= "never"
+Cactus::terminate       = "time"
+Cactus::cctk_final_time = 0.4
+#Cactus::terminate = "iteration"
+#Cactus::cctk_itlast        =  0
+#Cactus::cctk_itlast        =  100
+#Cactus::cctk_itlast        =  256
+#Cactus::cctk_itlast        = 1024
+#Cactus::terminate="runtime"
+#Cactus::max_runtime=28800  #  46 days
+#Cactus::max_runtime=1410  #  23.5 hours
+#Cactus::max_runtime=1380  #  23 hours
+#Cactus::max_runtime=690  #  11.5 hours
+#Cactus::max_runtime=450  #  7.5 hours
+#Cactus::max_runtime=210  #  3.5 hours
+
+#------------------------------------------------------------------------------
+# Activate all necessary thorns:
+#------------------------------------------------------------------------------
+
+ActiveThorns = "Boundary CartGrid3D CoordBase Fortran InitBase IOUtil LocalReduce SymBase Time"
+ActiveThorns = "AEILocalInterp LocalInterp"
+ActiveThorns = "MoL Slab SpaceMask"
+ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetSlab LoopControl"
+#ActiveThorns = "Formaline NaNChecker TerminationTrigger TimerReport"
+ActiveThorns = "NaNChecker TerminationTrigger TimerReport"
+ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal"
+ActiveThorns = "Constants TmunuBase HydroBase "
+#ActiveThorns = "EOS_Base EOS_IdealFluid EOS_Polytrope"
+ActiveThorns = "EOS_Omni"
+ActiveThorns = "GRHydro GRHydro_InitData"
+ActiveThorns = "HDF5"
+ActiveThorns = "SystemStatistics"
+
+#------------------------------------------------------------------------------
+# Diagnostic parameters:
+#------------------------------------------------------------------------------
+Carpet::output_timers_every = 1024
+Carpet::output_timer_tree_every = 256
+Carpet::output_initialise_timer_tree = "yes"
+Carpet::storage_verbose   = "no"
+#Carpet::timers_verbose    = "no"
+#Carpet::verbose           = "yes"  # only for Initial Data
+Carpet::verbose           = "no"
+Carpet::veryverbose       = "no"
+Carpet::grid_structure_filename = "carpet-grid-structure.asc"
+Carpet::grid_coordinates_filename = "carpet-grid-coordinates.asc"
+
+CarpetLib::output_bboxes  = "no"
+CarpetLib::print_memstats_every     = 1024
+CarpetLib::memstat_file    = "carpetlib-memory-statistics.asc"
+#CarpetLib::print_timestats_every    = 1024
+#CarpetLib::timestat_file    = "carpetlib-timing-statistics.asc"
+
+CarpetReduce::verbose    = "no"
+#CarpetRegrid2::verbose    = "yes"
+#CarpetRegrid2::veryverbose    = "yes"
+
+#NaNChecker::verbose         = "all"
+
+TimerReport::out_every    = 1024
+TimerReport::out_filename = "TimerReport"
+#TimerReport::output_all_timers = "yes"
+#TimerReport::output_all_timers_together = "yes"
+#TimerReport::output_all_timers_readable = "yes"
+#TimerReport::before_checkpoint          = "yes"
+TimerReport::n_top_timers = 40
+
+#------------------------------------------------------------------------------
+# Utility parameters:
+#------------------------------------------------------------------------------
+
+#NaNChecker::check_every    =  256
+NaNChecker::check_every   =  1
+NaNChecker::check_vars = "
+            ADMBase::curv 
+            ADMBase::metric 
+            ADMBase::lapse 
+            ADMBase::shift 
+            HydroBase::Bvec
+            HydroBase::rho 
+            HydroBase::eps 
+            HydroBase::press 
+            HydroBase::vel
+            HydroBase::w_lorentz
+"
+NaNChecker::action_if_found   =  "terminate"
+#NaNChecker::action_if_found  =  "abort"
+#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort"
+#NaNChecker::check_after=0
+
+#TerminationTrigger::max_walltime          =  0.0   # hours
+TerminationTrigger::max_walltime          =  24.0   # hours
+TerminationTrigger::on_remaining_walltime = 10.0   # minutes
+TerminationTrigger::termination_from_file   = "yes"
+TerminationTrigger::create_termination_file = "yes"
+TerminationTrigger::termination_file        = "../TERMINATE"
+
+
+#------------------------------------------------------------------------------
+# Run parameters:
+#------------------------------------------------------------------------------
+
+#------
+# Grid:
+#------
+
+#Time::dtfac = 0.8
+Time::dtfac = 0.4
+MoL::ODE_Method             = "rk2"
+MoL::MoL_Intermediate_Steps = 2
+#MoL::ODE_Method             = "rk4"
+#MoL::MoL_Intermediate_Steps = 4
+#MoL::MoL_Num_Scratch_Levels = 1
+
+CartGrid3D::type         = "coordbase"
+CartGrid3D::domain       = "full"
+CartGrid3D::avoid_origin = "no"
+
+CoordBase::domainsize       = "minmax"
+CoordBase::spacing          = "numcells"  # "gridspacing" or "numcells"
+
+#CoordBase::xmin = -0.035355339059327376220
+#CoordBase::ymin = -0.035355339059327376220
+#CoordBase::zmin = -0.0017677669529663688110
+#CoordBase::xmax =  0.035355339059327376220
+#CoordBase::ymax =  0.035355339059327376220
+#CoordBase::zmax =  0.0017677669529663688110
+#CoordBase::ncells_x   = 160  # dx or ncells_x
+#CoordBase::ncells_y   = 160  # dy or ncells_y
+#CoordBase::ncells_z   =   8  # dz or ncells_z
+
+# Approximate memory footprint: 650MB x 8 = 5.2GB
+
+CoordBase::xmin = -0.35355339059327376220
+CoordBase::ymin = -0.35355339059327376220
+CoordBase::zmin = -0.003535533905932737622
+CoordBase::xmax =  0.35355339059327376220
+CoordBase::ymax =  0.35355339059327376220
+CoordBase::zmax =  0.003535533905932737622
+CoordBase::ncells_x   = 800  # dx or ncells_x
+CoordBase::ncells_y   = 800  # dy or ncells_y
+CoordBase::ncells_z   =   8  # dz or ncells_z
+
+Driver::ghost_size                  = 3
+GRHydro::GRHydro_stencil            = 3
+
+# General Carpet parameters:
+Carpet::domain_from_coordbase    = "yes"
+
+CarpetLib::poison_new_memory         = "yes"
+CarpetLib::poison_value              = 114
+
+#------
+# MODEL:
+#------
+
+HydroBase::initial_hydro                = "shocktube"
+GRHydro_InitData::shocktube_type        = "diagshock2d"
+GRHydro_InitData::shock_case            = "Balsara4"
+
+ADMBase::metric_type             = "physical"
+ADMBase::initial_data            = "Cartesian Minkowski"
+ADMBase::initial_lapse           = "one"
+ADMBase::initial_shift           = "zero"
+ADMBase::initial_dtlapse         = "none"
+ADMBase::initial_dtshift         = "none"
+
+EOS_Omni::poly_K     = 100.0
+EOS_Omni::poly_gamma = 1.6666666666666666666
+EOS_Omni::gl_K       = 100.0
+EOS_Omni::gl_gamma   = 1.6666666666666666666
+GRHydro::initial_Gamma = 1.6666666666666666666
+
+#----------
+# Numerics:
+#----------
+#Boundary::radpower                     = 2
+#HydroBase::timelevels            = 3
+HydroBase::timelevels            = 2
+#ADMMacros::spatial_order = 4
+SpaceMask::use_mask      = "yes"
+
+#-----------
+# Evolution:
+#-----------
+
+HydroBase::evolution_method      = "GRHydro"
+HydroBase::Bvec_evolution_method = "GRHydro"
+
+GRHydro::method_type               = "RSA FV"
+GRHydro::recon_method              = "tvd"
+#GRHydro::recon_method              = "ppm"
+#GRHydro::recon_method              = "eno"
+GRHydro::tvd_limiter               = "vanleerMC2"
+#GRHydro::tvd_limiter               = "minmod"
+GRHydro::recon_vars                = "primitive"
+GRHydro::riemann_solver            = "HLLE"
+GRHydro::decouple_normal_Bfield    = "yes"
+#GRHydro::riemann_solver            = "LLF"
+GRHydro::bound                     = "none"
+GRHydro::GRHydro_eos_type          = "General"
+GRHydro::GRHydro_eos_table         = "Ideal_Fluid"
+
+GRHydro::Grhydro_MaxNumConstrainedVars = 13
+GRHydro::GRHydro_MaxNumEvolvedVars     = 10
+
+GRHydro::clean_divergence          = "yes"
+GRHydro::psidcspeed                = "light speed"
+GRHydro::kap_dc                    = 10.0
+GRHydro::track_divB                = "yes"
+GRHydro::calculate_bcom            = "yes" 
+
+#------------------------------------------------------------------------------
+# Output:
+#------------------------------------------------------------------------------
+
+IO::out_dir = $parfile
+IO::out_fileinfo           = "none"
+
+IOBasic::outInfo_every = 1
+IOBasic::outInfo_reductions = "maximum"
+IOBasic::outInfo_vars  = "
+ Carpet::physical_time_per_hour
+ GRHydro::bcom_sq
+ HydroBase::rho
+ HydroBase::press
+ SystemStatistics::maxrss_mb
+"
+
+#IOScalar::outScalar_every      = 1024
+#IOScalar::one_file_per_group   = "yes"
+#IOScalar::outScalar_reductions = "minimum maximum average norm1 norm2"
+#IOScalar::outScalar_vars       = "
+# GRHydro::bcom
+# GRHydro::bcom0
+# GRHydro::bcom_sq
+# GRHydro::Bcons
+# GRHydro::dens
+# GrHydro::divB
+# GRHydro::scon
+# GRHydro::tau
+# HydroBase::Bvec
+# HydroBase::eps
+# HydroBase::press
+# HydroBase::rho
+# HydroBase::vel
+# HydroBase::w_lorentz
+# SystemStatistics::process_memory_mb
+#"
+
+IOASCII::one_file_per_group     = "yes"
+IOASCII::output_symmetry_points = "no"
+IOASCII::output_ghost_points    = "no"  #hg version only!
+IOASCII::output_boundary_points = "no"  #hg version only!
+#IOASCII::out3D_ghosts       = "no"  #hg (DEPRECATED)
+#IOASCII::out3D_outer_ghosts = "no" #hg (DEPRECATED)
+
+IOASCII::out0D_every     = 200
+IOASCII::out0D_vars      = "
+ Carpet::timing
+"
+
+#Set these IOASCII options for initial data only:
+#IOASCII::out1D_every     = 1
+IOASCII::out1D_every     = 1131
+#IOASCII::out1D_dt        = 0.01
+IOASCII::out1D_d         = "yes"
+IOASCII::out1D_vars      = "
+ GRHydro::bcom
+ GRHydro::bcom0
+ GRHydro::bcom_sq
+ GRHydro::Bcons
+ GRHydro::dens
+ GrHydro::divB
+ GrHydro::psidc
+ GRHydro::scon
+ GRHydro::tau
+ HydroBase::Bvec
+ HydroBase::eps
+ HydroBase::press
+ HydroBase::rho
+ HydroBase::vel
+ HydroBase::w_lorentz
+"
+
+IOASCII::out2D_every     = 1131
+IOASCII::out2D_vars      = "
+ GRHydro::bcom
+ GRHydro::bcom0
+ GRHydro::bcom_sq
+ GRHydro::Bcons
+ GRHydro::dens
+ GrHydro::divB
+ GrHydro::psidc
+ GRHydro::scon
+ GRHydro::tau
+ HydroBase::Bvec
+ HydroBase::eps
+ HydroBase::press
+ HydroBase::rho
+ HydroBase::vel
+ HydroBase::w_lorentz
+"
+
+CarpetIOHDF5::one_file_per_group             = "yes"
+CarpetIOHDF5::open_one_input_file_at_a_time  = "yes"
+
+#CarpetIOHDF5::output_symmetry_points         = "no"
+#CarpetIOHDF5::output_ghost_points            = "no"
+#CarpetIOHDF5::output_boundary_points         = "no"
+#CarpetIOHDF5::output_buffer_points           = "no"
+
+#CarpetIOHDF5::out2D_every                    = 1
+CarpetIOHDF5::out2D_every                    = 100
+#CarpetIOHDF5::out2D_dt                       = 0.01
+CarpetIOHDF5::out2D_xy                       = "yes"
+CarpetIOHDF5::out2D_xz                       = "no"
+CarpetIOHDF5::out2D_yz                       = "no"
+CarpetIOHDF5::out2D_xyplane_z                = 0.0
+CarpetIOHDF5::out2D_vars      = "
+  GRHydro::bcom
+  GRHydro::bcom0
+  GRHydro::bcom_sq
+  GRHydro::Bcons
+  GRHydro::dens
+  GRHydro::divB
+  GrHydro::psidc
+  GRHydro::scon
+  GRHydro::tau
+  HydroBase::Bvec
+  HydroBase::eps
+  HydroBase::press
+  HydroBase::rho
+  HydroBase::vel
+  HydroBase::w_lorentz
+"
+
+#CarpetIOHDF5::out3D_every = 1
+CarpetIOHDF5::out3D_every = 100
+#CarpetIOHDF5::out3D_dt = 0.01
+CarpetIOHDF5::out3D_vars  = "
+ GRHydro::bcom
+ GRHydro::bcom0
+ GRHydro::bcom_sq
+ GRHydro::divB
+ HydroBase::Bvec
+ HydroBase::eps
+ HydroBase::press
+ HydroBase::rho
+ HydroBase::vel
+ HydroBase::w_lorentz
+"
+
+#------------------------------------------------------------------------------
+# Analysis:
+#------------------------------------------------------------------------------
+
+
+#------------------------------------------------------------------------------
+# Checkpoint/Recovery:
+#------------------------------------------------------------------------------
+#IOHDF5::checkpoint                  = "yes"
+#IO::checkpoint_dir                  = $parfile
+#IO::checkpoint_ID                   = "yes"
+#IO::checkpoint_every                = 6144
+##IO::checkpoint_every_walltime_hours = 6.0
+#IO::checkpoint_keep=3
+#IO::checkpoint_on_terminate         = "yes"
+#
+#IO::recover     = "autoprobe"
+#IO::recover_dir = $parfile
+
+
diff --git a/par/balsara5_2d.par b/par/balsara5_2d.par
new file mode 100644
index 0000000..290f60a
--- /dev/null
+++ b/par/balsara5_2d.par
@@ -0,0 +1,384 @@
+#==============================================================================
+# Balsara magnetized shock tube test #5: GRHydro_InitData
+#==============================================================================
+
+#------------------------------------------------------------------------------
+# Cactus parameters:
+#------------------------------------------------------------------------------
+Cactus::cctk_run_title = "Balsara magnetized shock tube test #5 along z=0 plane square diagonal"
+Cactus::cctk_full_warnings         = "yes"
+
+#Cactus::terminate		= "never"
+Cactus::terminate       = "time"
+Cactus::cctk_final_time = 0.55
+#Cactus::terminate = "iteration"
+#Cactus::cctk_itlast        =  0
+#Cactus::cctk_itlast        =  100
+#Cactus::cctk_itlast        =  256
+#Cactus::cctk_itlast        = 1024
+#Cactus::terminate="runtime"
+#Cactus::max_runtime=28800  #  46 days
+#Cactus::max_runtime=1410  #  23.5 hours
+#Cactus::max_runtime=1380  #  23 hours
+#Cactus::max_runtime=690  #  11.5 hours
+#Cactus::max_runtime=450  #  7.5 hours
+#Cactus::max_runtime=210  #  3.5 hours
+
+#------------------------------------------------------------------------------
+# Activate all necessary thorns:
+#------------------------------------------------------------------------------
+
+ActiveThorns = "Boundary CartGrid3D CoordBase Fortran InitBase IOUtil LocalReduce SymBase Time"
+ActiveThorns = "AEILocalInterp LocalInterp"
+ActiveThorns = "MoL Slab SpaceMask"
+ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetSlab LoopControl"
+#ActiveThorns = "Formaline NaNChecker TerminationTrigger TimerReport"
+ActiveThorns = "NaNChecker TerminationTrigger TimerReport"
+ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal"
+ActiveThorns = "Constants TmunuBase HydroBase "
+#ActiveThorns = "EOS_Base EOS_IdealFluid EOS_Polytrope"
+ActiveThorns = "EOS_Omni"
+ActiveThorns = "GRHydro GRHydro_InitData"
+ActiveThorns = "HDF5"
+ActiveThorns = "SystemStatistics"
+
+#------------------------------------------------------------------------------
+# Diagnostic parameters:
+#------------------------------------------------------------------------------
+Carpet::output_timers_every = 1024
+Carpet::output_timer_tree_every = 256
+Carpet::output_initialise_timer_tree = "yes"
+Carpet::storage_verbose   = "no"
+#Carpet::timers_verbose    = "no"
+#Carpet::verbose           = "yes"  # only for Initial Data
+Carpet::verbose           = "no"
+Carpet::veryverbose       = "no"
+Carpet::grid_structure_filename = "carpet-grid-structure.asc"
+Carpet::grid_coordinates_filename = "carpet-grid-coordinates.asc"
+
+CarpetLib::output_bboxes  = "no"
+CarpetLib::print_memstats_every     = 1024
+CarpetLib::memstat_file    = "carpetlib-memory-statistics.asc"
+#CarpetLib::print_timestats_every    = 1024
+#CarpetLib::timestat_file    = "carpetlib-timing-statistics.asc"
+
+CarpetReduce::verbose    = "no"
+#CarpetRegrid2::verbose    = "yes"
+#CarpetRegrid2::veryverbose    = "yes"
+
+#NaNChecker::verbose         = "all"
+
+TimerReport::out_every    = 1024
+TimerReport::out_filename = "TimerReport"
+#TimerReport::output_all_timers = "yes"
+#TimerReport::output_all_timers_together = "yes"
+#TimerReport::output_all_timers_readable = "yes"
+#TimerReport::before_checkpoint          = "yes"
+TimerReport::n_top_timers = 40
+
+#------------------------------------------------------------------------------
+# Utility parameters:
+#------------------------------------------------------------------------------
+
+#NaNChecker::check_every    =  256
+NaNChecker::check_every   =  1
+NaNChecker::check_vars = "
+            ADMBase::curv 
+            ADMBase::metric 
+            ADMBase::lapse 
+            ADMBase::shift 
+            HydroBase::Bvec
+            HydroBase::rho 
+            HydroBase::eps 
+            HydroBase::press 
+            HydroBase::vel
+            HydroBase::w_lorentz
+"
+NaNChecker::action_if_found   =  "terminate"
+#NaNChecker::action_if_found  =  "abort"
+#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort"
+#NaNChecker::check_after=0
+
+#TerminationTrigger::max_walltime          =  0.0   # hours
+TerminationTrigger::max_walltime          =  24.0   # hours
+TerminationTrigger::on_remaining_walltime = 10.0   # minutes
+TerminationTrigger::termination_from_file   = "yes"
+TerminationTrigger::create_termination_file = "yes"
+TerminationTrigger::termination_file        = "../TERMINATE"
+
+
+#------------------------------------------------------------------------------
+# Run parameters:
+#------------------------------------------------------------------------------
+
+#------
+# Grid:
+#------
+
+#Time::dtfac = 0.8
+Time::dtfac = 0.4
+MoL::ODE_Method             = "rk2"
+MoL::MoL_Intermediate_Steps = 2
+#MoL::ODE_Method             = "rk4"
+#MoL::MoL_Intermediate_Steps = 4
+#MoL::MoL_Num_Scratch_Levels = 1
+
+CartGrid3D::type         = "coordbase"
+CartGrid3D::domain       = "full"
+CartGrid3D::avoid_origin = "no"
+
+CoordBase::domainsize       = "minmax"
+CoordBase::spacing          = "numcells"  # "gridspacing" or "numcells"
+
+#CoordBase::xmin = -0.035355339059327376220
+#CoordBase::ymin = -0.035355339059327376220
+#CoordBase::zmin = -0.0017677669529663688110
+#CoordBase::xmax =  0.035355339059327376220
+#CoordBase::ymax =  0.035355339059327376220
+#CoordBase::zmax =  0.0017677669529663688110
+#CoordBase::ncells_x   = 160  # dx or ncells_x
+#CoordBase::ncells_y   = 160  # dy or ncells_y
+#CoordBase::ncells_z   =   8  # dz or ncells_z
+
+# Approximate memory footprint: 650MB x 8 = 5.2GB
+
+CoordBase::xmin = -0.35355339059327376220
+CoordBase::ymin = -0.35355339059327376220
+CoordBase::zmin = -0.003535533905932737622
+CoordBase::xmax =  0.35355339059327376220
+CoordBase::ymax =  0.35355339059327376220
+CoordBase::zmax =  0.003535533905932737622
+CoordBase::ncells_x   = 800  # dx or ncells_x
+CoordBase::ncells_y   = 800  # dy or ncells_y
+CoordBase::ncells_z   =   8  # dz or ncells_z
+
+Driver::ghost_size                  = 3
+GRHydro::GRHydro_stencil            = 3
+
+# General Carpet parameters:
+Carpet::domain_from_coordbase    = "yes"
+
+CarpetLib::poison_new_memory         = "yes"
+CarpetLib::poison_value              = 114
+
+#------
+# MODEL:
+#------
+
+HydroBase::initial_hydro                = "shocktube"
+GRHydro_InitData::shocktube_type        = "diagshock2d"
+GRHydro_InitData::shock_case            = "Balsara5"
+
+ADMBase::metric_type             = "physical"
+ADMBase::initial_data            = "Cartesian Minkowski"
+ADMBase::initial_lapse           = "one"
+ADMBase::initial_shift           = "zero"
+ADMBase::initial_dtlapse         = "none"
+ADMBase::initial_dtshift         = "none"
+
+EOS_Omni::poly_K     = 100.0
+EOS_Omni::poly_gamma = 1.6666666666666666666
+EOS_Omni::gl_K       = 100.0
+EOS_Omni::gl_gamma   = 1.6666666666666666666
+GRHydro::initial_Gamma = 1.6666666666666666666
+
+#----------
+# Numerics:
+#----------
+#Boundary::radpower                     = 2
+#HydroBase::timelevels            = 3
+HydroBase::timelevels            = 2
+#ADMMacros::spatial_order = 4
+SpaceMask::use_mask      = "yes"
+
+#-----------
+# Evolution:
+#-----------
+
+HydroBase::evolution_method      = "GRHydro"
+HydroBase::Bvec_evolution_method = "GRHydro"
+
+GRHydro::method_type               = "RSA FV"
+GRHydro::recon_method              = "tvd"
+#GRHydro::recon_method              = "ppm"
+#GRHydro::recon_method              = "eno"
+GRHydro::tvd_limiter               = "vanleerMC2"
+#GRHydro::tvd_limiter               = "minmod"
+GRHydro::recon_vars                = "primitive"
+GRHydro::riemann_solver            = "HLLE"
+GRHydro::decouple_normal_Bfield    = "yes"
+#GRHydro::riemann_solver            = "LLF"
+GRHydro::bound                     = "none"
+GRHydro::GRHydro_eos_type          = "General"
+GRHydro::GRHydro_eos_table         = "Ideal_Fluid"
+
+GRHydro::Grhydro_MaxNumConstrainedVars = 13
+GRHydro::GRHydro_MaxNumEvolvedVars     = 10
+
+GRHydro::clean_divergence          = "yes"
+GRHydro::psidcspeed                = "light speed"
+GRHydro::kap_dc                    = 10.0
+GRHydro::track_divB                = "yes"
+GRHydro::calculate_bcom            = "yes" 
+
+#------------------------------------------------------------------------------
+# Output:
+#------------------------------------------------------------------------------
+
+IO::out_dir = $parfile
+IO::out_fileinfo           = "none"
+
+IOBasic::outInfo_every = 1
+IOBasic::outInfo_reductions = "maximum"
+IOBasic::outInfo_vars  = "
+ Carpet::physical_time_per_hour
+ GRHydro::bcom_sq
+ HydroBase::rho
+ HydroBase::press
+ SystemStatistics::maxrss_mb
+"
+
+#IOScalar::outScalar_every      = 1024
+#IOScalar::one_file_per_group   = "yes"
+#IOScalar::outScalar_reductions = "minimum maximum average norm1 norm2"
+#IOScalar::outScalar_vars       = "
+# GRHydro::bcom
+# GRHydro::bcom0
+# GRHydro::bcom_sq
+# GRHydro::Bcons
+# GRHydro::dens
+# GrHydro::divB
+# GRHydro::scon
+# GRHydro::tau
+# HydroBase::Bvec
+# HydroBase::eps
+# HydroBase::press
+# HydroBase::rho
+# HydroBase::vel
+# HydroBase::w_lorentz
+# SystemStatistics::process_memory_mb
+#"
+
+IOASCII::one_file_per_group     = "yes"
+IOASCII::output_symmetry_points = "no"
+IOASCII::output_ghost_points    = "no"  #hg version only!
+IOASCII::output_boundary_points = "no"  #hg version only!
+#IOASCII::out3D_ghosts       = "no"  #hg (DEPRECATED)
+#IOASCII::out3D_outer_ghosts = "no" #hg (DEPRECATED)
+
+IOASCII::out0D_every     = 200
+IOASCII::out0D_vars      = "
+ Carpet::timing
+"
+
+#Set these IOASCII options for initial data only:
+#IOASCII::out1D_every     = 1
+IOASCII::out1D_every     = 1555
+#IOASCII::out1D_dt        = 0.01
+IOASCII::out1D_d         = "yes"
+IOASCII::out1D_vars      = "
+ GRHydro::bcom
+ GRHydro::bcom0
+ GRHydro::bcom_sq
+ GRHydro::Bcons
+ GRHydro::dens
+ GrHydro::divB
+ GrHydro::psidc
+ GRHydro::scon
+ GRHydro::tau
+ HydroBase::Bvec
+ HydroBase::eps
+ HydroBase::press
+ HydroBase::rho
+ HydroBase::vel
+ HydroBase::w_lorentz
+"
+
+IOASCII::out2D_every     = 1555
+IOASCII::out2D_vars      = "
+ GRHydro::bcom
+ GRHydro::bcom0
+ GRHydro::bcom_sq
+ GRHydro::Bcons
+ GRHydro::dens
+ GrHydro::divB
+ GrHydro::psidc
+ GRHydro::scon
+ GRHydro::tau
+ HydroBase::Bvec
+ HydroBase::eps
+ HydroBase::press
+ HydroBase::rho
+ HydroBase::vel
+ HydroBase::w_lorentz
+"
+
+CarpetIOHDF5::one_file_per_group             = "yes"
+CarpetIOHDF5::open_one_input_file_at_a_time  = "yes"
+
+#CarpetIOHDF5::output_symmetry_points         = "no"
+#CarpetIOHDF5::output_ghost_points            = "no"
+#CarpetIOHDF5::output_boundary_points         = "no"
+#CarpetIOHDF5::output_buffer_points           = "no"
+
+#CarpetIOHDF5::out2D_every                    = 1
+CarpetIOHDF5::out2D_every                    = 100
+#CarpetIOHDF5::out2D_dt                       = 0.01
+CarpetIOHDF5::out2D_xy                       = "yes"
+CarpetIOHDF5::out2D_xz                       = "no"
+CarpetIOHDF5::out2D_yz                       = "no"
+CarpetIOHDF5::out2D_xyplane_z                = 0.0
+CarpetIOHDF5::out2D_vars      = "
+  GRHydro::bcom
+  GRHydro::bcom0
+  GRHydro::bcom_sq
+  GRHydro::Bcons
+  GRHydro::dens
+  GRHydro::divB
+  GrHydro::psidc
+  GRHydro::scon
+  GRHydro::tau
+  HydroBase::Bvec
+  HydroBase::eps
+  HydroBase::press
+  HydroBase::rho
+  HydroBase::vel
+  HydroBase::w_lorentz
+"
+
+#CarpetIOHDF5::out3D_every = 1
+CarpetIOHDF5::out3D_every = 100
+#CarpetIOHDF5::out3D_dt = 0.01
+CarpetIOHDF5::out3D_vars  = "
+ GRHydro::bcom
+ GRHydro::bcom0
+ GRHydro::bcom_sq
+ GRHydro::divB
+ HydroBase::Bvec
+ HydroBase::eps
+ HydroBase::press
+ HydroBase::rho
+ HydroBase::vel
+ HydroBase::w_lorentz
+"
+
+#------------------------------------------------------------------------------
+# Analysis:
+#------------------------------------------------------------------------------
+
+
+#------------------------------------------------------------------------------
+# Checkpoint/Recovery:
+#------------------------------------------------------------------------------
+#IOHDF5::checkpoint                  = "yes"
+#IO::checkpoint_dir                  = $parfile
+#IO::checkpoint_ID                   = "yes"
+#IO::checkpoint_every                = 6144
+##IO::checkpoint_every_walltime_hours = 6.0
+#IO::checkpoint_keep=3
+#IO::checkpoint_on_terminate         = "yes"
+#
+#IO::recover     = "autoprobe"
+#IO::recover_dir = $parfile
+
+