Add Kerr-Schild parameter files

3 files changed, 222 insertions, 2 deletions

diff --git a/par/ks_bssn_1.par b/par/ks_bssn_1.par
new file mode 100644
index 0000000..38c2850
--- /dev/null
+++ b/par/ks_bssn_1.par
@@ -0,0 +1,110 @@
+ActiveThorns = "
+        Boundary CartGrid3D CoordBase IOUtil InitBase MoL SymBase
+        ADMBase ADMCoupling CoordGauge StaticConformal
+        NaNChecker
+        Carpet CarpetIOASCII CarpetIOBasic CarpetIOScalar CarpetLib CarpetReduce LoopControl
+        ReflectionSymmetry
+        Exact Formaline GSL
+        Slab
+        GenericFD
+        ML_BSSN
+"
+
+
+
+Cactus::cctk_itlast = 0
+
+grid::type = "CoordBase"
+
+CoordBase::domainsize = "minmax"
+
+CoordBase::xmax = 4.0
+CoordBase::ymax = 4.0
+CoordBase::zmax = 4.0
+CoordBase::dx   = 0.2
+CoordBase::dy   = 0.2
+CoordBase::dz   = 0.2
+
+ReflectionSymmetry::reflection_x = yes
+ReflectionSymmetry::reflection_y = yes
+ReflectionSymmetry::reflection_z = yes
+
+CoordBase::boundary_shiftout_x_lower = 1
+CoordBase::boundary_shiftout_y_lower = 1
+CoordBase::boundary_shiftout_z_lower = 1
+
+Carpet::domain_from_coordbase = yes
+driver::ghost_size            = 2
+Carpet::use_buffer_zones      = yes
+
+CoordBase::boundary_size_x_lower = 2
+CoordBase::boundary_size_y_lower = 2
+CoordBase::boundary_size_z_lower = 2
+CoordBase::boundary_size_x_upper = 2
+CoordBase::boundary_size_y_upper = 2
+CoordBase::boundary_size_z_upper = 2
+
+
+
+ADMBase::initial_data    = "exact"
+ADMBase::initial_lapse   = "exact"
+ADMBase::initial_shift   = "exact"
+ADMBase::initial_dtlapse = "zero"
+ADMBase::initial_dtshift = "zero"
+
+Exact::exact_model = "Kerr/Kerr-Schild"
+Exact::Kerr_KerrSchild__mass = 1.0
+Exact::Kerr_KerrSchild__spin = 0.0
+
+ML_BSSN::harmonicN       = 0      # 1+log
+ML_BSSN::harmonicF       = 2.0    # 1+log
+ML_BSSN::ShiftGammaCoeff = 0.75
+ML_BSSN::BetaDriver      = 0.5
+
+
+
+IO::out_dir = $parfile
+
+IOBasic::outInfo_every = 1
+IOBasic::outInfo_vars = "
+        ML_BSSN::H
+        ML_BSSN::trKrhs
+"
+
+IOScalar::one_file_per_group = yes
+IOScalar::outScalar_every    = 1
+IOScalar::outScalar_vars     = "
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ML_BSSN::Ham
+        ML_BSSN::mom
+        ML_BSSN::cons_detg
+        ML_BSSN::cons_Gamma
+        ML_BSSN::cons_traceA
+        ML_BSSN::log_confacrhs
+        ML_BSSN::metricrhs
+        ML_BSSN::Gammarhs
+        ML_BSSN::trace_curvrhs
+        ML_BSSN::curvrhs
+"
+
+IOASCII::one_file_per_group = yes
+IOASCII::out1D_every        = 1
+IOASCII::out1D_vars         = "
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ML_BSSN::Ham
+        ML_BSSN::mom
+        ML_BSSN::cons_detg
+        ML_BSSN::cons_Gamma
+        ML_BSSN::cons_traceA
+        ML_BSSN::log_confacrhs
+        ML_BSSN::metricrhs
+        ML_BSSN::Gammarhs
+        ML_BSSN::trace_curvrhs
+        ML_BSSN::curvrhs
+"
diff --git a/par/ks_bssn_2.par b/par/ks_bssn_2.par
new file mode 100644
index 0000000..957c1bd
--- /dev/null
+++ b/par/ks_bssn_2.par
@@ -0,0 +1,110 @@
+ActiveThorns = "
+        Boundary CartGrid3D CoordBase IOUtil InitBase MoL SymBase
+        ADMBase ADMCoupling CoordGauge StaticConformal
+        NaNChecker
+        Carpet CarpetIOASCII CarpetIOBasic CarpetIOScalar CarpetLib CarpetReduce LoopControl
+        ReflectionSymmetry
+        Exact Formaline GSL
+        Slab
+        GenericFD
+        ML_BSSN
+"
+
+
+
+Cactus::cctk_itlast = 0
+
+grid::type = "CoordBase"
+
+CoordBase::domainsize = "minmax"
+
+CoordBase::xmax = 4.0
+CoordBase::ymax = 4.0
+CoordBase::zmax = 4.0
+CoordBase::dx   = 0.1
+CoordBase::dy   = 0.1
+CoordBase::dz   = 0.1
+
+ReflectionSymmetry::reflection_x = yes
+ReflectionSymmetry::reflection_y = yes
+ReflectionSymmetry::reflection_z = yes
+
+CoordBase::boundary_shiftout_x_lower = 1
+CoordBase::boundary_shiftout_y_lower = 1
+CoordBase::boundary_shiftout_z_lower = 1
+
+Carpet::domain_from_coordbase = yes
+driver::ghost_size            = 2
+Carpet::use_buffer_zones      = yes
+
+CoordBase::boundary_size_x_lower = 2
+CoordBase::boundary_size_y_lower = 2
+CoordBase::boundary_size_z_lower = 2
+CoordBase::boundary_size_x_upper = 2
+CoordBase::boundary_size_y_upper = 2
+CoordBase::boundary_size_z_upper = 2
+
+
+
+ADMBase::initial_data    = "exact"
+ADMBase::initial_lapse   = "exact"
+ADMBase::initial_shift   = "exact"
+ADMBase::initial_dtlapse = "zero"
+ADMBase::initial_dtshift = "zero"
+
+Exact::exact_model = "Kerr/Kerr-Schild"
+Exact::Kerr_KerrSchild__mass = 1.0
+Exact::Kerr_KerrSchild__spin = 0.0
+
+ML_BSSN::harmonicN       = 0      # 1+log
+ML_BSSN::harmonicF       = 2.0    # 1+log
+ML_BSSN::ShiftGammaCoeff = 0.75
+ML_BSSN::BetaDriver      = 0.5
+
+
+
+IO::out_dir = $parfile
+
+IOBasic::outInfo_every = 1
+IOBasic::outInfo_vars = "
+        ML_BSSN::H
+        ML_BSSN::trKrhs
+"
+
+IOScalar::one_file_per_group = yes
+IOScalar::outScalar_every    = 1
+IOScalar::outScalar_vars     = "
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ML_BSSN::Ham
+        ML_BSSN::mom
+        ML_BSSN::cons_detg
+        ML_BSSN::cons_Gamma
+        ML_BSSN::cons_traceA
+        ML_BSSN::log_confacrhs
+        ML_BSSN::metricrhs
+        ML_BSSN::Gammarhs
+        ML_BSSN::trace_curvrhs
+        ML_BSSN::curvrhs
+"
+
+IOASCII::one_file_per_group = yes
+IOASCII::out1D_every        = 1
+IOASCII::out1D_vars         = "
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ML_BSSN::Ham
+        ML_BSSN::mom
+        ML_BSSN::cons_detg
+        ML_BSSN::cons_Gamma
+        ML_BSSN::cons_traceA
+        ML_BSSN::log_confacrhs
+        ML_BSSN::metricrhs
+        ML_BSSN::Gammarhs
+        ML_BSSN::trace_curvrhs
+        ML_BSSN::curvrhs
+"
diff --git a/par/m_bssn.par b/par/m_bssn.par
index 5533f4d..61616fc 100644
--- a/par/m_bssn.par
+++ b/par/m_bssn.par
@@ -29,8 +29,8 @@ ADMBase::initial_dtshift = "zero"
 
 ML_BSSN::my_initial_data = "Minkowski"
 
-ML_BSSN::harmonic_n      = 0      # 1+log
-ML_BSSN::harmonic_f      = 2.0    # 1+log
+ML_BSSN::harmonicN       = 0      # 1+log
+ML_BSSN::harmonicF       = 2.0    # 1+log
 ML_BSSN::ShiftGammaCoeff = 0.75
 ML_BSSN::BetaDriver      = 0.5