aboutsummaryrefslogtreecommitdiff
path: root/test/teukolsky.par
diff options
context:
space:
mode:
authorrhaas <rhaas@4f5cb9a8-4dd8-4c2d-9bbd-173fa4467843>2010-11-08 13:08:48 +0000
committerrhaas <rhaas@4f5cb9a8-4dd8-4c2d-9bbd-173fa4467843>2010-11-08 13:08:48 +0000
commit50d42a3f2c07f9f6d8820b32ec58af6c82621c32 (patch)
treef47233ec36cd8f3580ae7f0533a6182e04147bdd /test/teukolsky.par
parente6eeb523c1312cbcc80386571b9577258cc09669 (diff)
fix teukolsky testsuite
needed two ghost zones for upwind derivatives in McLachlan, add initial extrapolatin for Gammas git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/WeylScal4/trunk@66 4f5cb9a8-4dd8-4c2d-9bbd-173fa4467843
Diffstat (limited to 'test/teukolsky.par')
-rw-r--r--test/teukolsky.par52
1 files changed, 33 insertions, 19 deletions
diff --git a/test/teukolsky.par b/test/teukolsky.par
index e3a6190..af713d9 100644
--- a/test/teukolsky.par
+++ b/test/teukolsky.par
@@ -4,10 +4,10 @@
# $Header$
# Required thorns
-ActiveThorns = "CartGrid3D CoordBase SymBase Boundary Time Carpet CarpetLib CarpetSlab CarpetRegrid CarpetReduce InitBase CarpetIOAscii IOUtil IOBasic ADMBase ADMMacros ADMConstraints GenericFD TmunuBase ADMCoupling StaticConformal CoordGauge SpaceMask IDLinearWaves LoopControl GSL WeylScal4 ML_BSSN_O2 ML_BSSN_O2_Helper GenericFD MoL ReflectionSymmetry"
+ActiveThorns = "CartGrid3D CoordBase SymBase Boundary Time Carpet CarpetLib CarpetSlab CarpetRegrid CarpetReduce InitBase CarpetIOAscii IOUtil IOBasic ADMBase ADMMacros ADMConstraints GenericFD TmunuBase ADMCoupling StaticConformal CoordGauge SpaceMask IDLinearWaves LoopControl GSL WeylScal4 ML_BSSN_O2 ML_BSSN_O2_Helper NewRad GenericFD MoL ReflectionSymmetry"
# Grid
-driver::ghost_size = 1
+driver::ghost_size = 2
CartGrid3D::avoid_origin = no
@@ -38,16 +38,16 @@ CoordBase::zmin = 0.0
CoordBase::xmax = 1.9375
CoordBase::ymax = 1.9375
CoordBase::zmax = 2.0
-CoordBase::dx = 0.5
-CoordBase::dy = 0.5
-CoordBase::dz = 0.5
+CoordBase::dx = 0.5
+CoordBase::dy = 0.5
+CoordBase::dz = 0.5
-CoordBase::boundary_size_x_lower = 1
-CoordBase::boundary_size_y_lower = 1
-CoordBase::boundary_size_z_lower = 1
-CoordBase::boundary_size_x_upper = 1
-CoordBase::boundary_size_y_upper = 1
-CoordBase::boundary_size_z_upper = 1
+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
CoordBase::boundary_shiftout_z_lower = 1
CoordBase::boundary_shiftout_z_upper = 0
@@ -58,23 +58,37 @@ ReflectionSymmetry::reflection_z = "yes"
cactus::cctk_initial_time = 0
cactus::cctk_itlast = 4
+# Evolution
+Carpet::use_buffer_zones = "yes"
+MoL::MoL_Intermediate_Steps = 1 # number of iterations for ICN
+
# Initial data
-ADMBase::initial_data = "teukwaves"
-ADMBase::initial_shift = "zero"
-ADMBase::initial_lapse = "one"
+ADMBase::initial_data = "teukwaves"
+ADMBase::initial_shift = "zero"
+ADMBase::initial_lapse = "one"
+ADMBase::initial_dtlapse = "zero"
+ADMBase::initial_dtshift = "zero"
IDLinearWaves::mvalue = 2
# Gauge
-ADMBase::lapse_evolution_method = "ML_BSSN_O2"
-ADMBase::shift_evolution_method = "ML_BSSN_O2"
+ADMBase::lapse_evolution_method = "ML_BSSN_O2"
+ADMBase::shift_evolution_method = "ML_BSSN_O2"
+ADMBase::dtlapse_evolution_method = "ML_BSSN_O2"
ADMBase::dtshift_evolution_method = "ML_BSSN_O2"
# Evolution
ADMBase::evolution_method = "ML_BSSN_O2"
-ML_BSSN_O2::harmonicF = 2
-ML_BSSN_O2::harmonicN = 1 # 1+log
-ML_BSSN_O2::ShiftGammaCoeff = 0
+ML_BSSN_O2::harmonicN = 1 # 1+log
+ML_BSSN_O2::harmonicF = 2.0 # 1+log
+ML_BSSN_O2::ShiftGammaCoeff = 0.75
+ML_BSSN_O2::BetaDriver = 1.0
+ML_BSSN_O2::LapseAdvectionCoeff = 1.0
+ML_BSSN_O2::ShiftAdvectionCoeff = 1.0
+
+ML_BSSN_O2::MinimumLapse = 1.0e-8
+
+ML_BSSN_O2::my_initial_boundary_condition = "extrapolate-gammas"
# These are necessary
ML_BSSN_O2::ML_log_confac_bound = "flat"