blob: ada11a066f1c354ed06525e50a2886d35f9ec5e6 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
|
# Initial data and short evolution for gauge wave in Minkowski spacetime
# Author: Mitica Vulcanov <vulcan@aei.mpg.de>
# Author: Roland Haas <roland.haas@physics.gatech.edu>
# $Header$
# Required thorns
ActiveThorns = "CartGrid3D CoordBase SymBase Boundary Time Carpet CarpetLib CarpetSlab CarpetRegrid2 CarpetReduce InitBase CarpetIOAscii IOUtil IOBasic ADMBase ADMMacros ADMConstraints GenericFD TmunuBase ADMCoupling StaticConformal CoordGauge SpaceMask IDLinearWaves LoopControl GSL WeylScal4 ML_BSSN ML_BSSN_Helper NewRad GenericFD MoL ReflectionSymmetry"
# Grid
driver::ghost_size = 2
CartGrid3D::avoid_origin = no
Carpet::max_refinement_levels = 3
Carpet::prolongation_order_space = 3 # 1 is buggy in git version
Carpet::prolongation_order_time = 2
Carpet::domain_from_coordbase = yes
Carpet::poison_new_timelevels = "yes"
Carpet::check_for_poison = "no"
CarpetLib::poison_new_memory = "yes"
CarpetLib::poison_value = 114
InitBase::initial_data_setup_method = init_some_levels
Carpet::init_fill_timelevels = yes
CarpetRegrid2::min_distance = 0
CarpetRegrid2::ensure_proper_nesting = "yes"
CarpetRegrid2::regrid_every = 0
CarpetRegrid2::num_centres = 1
CarpetRegrid2::num_levels_1 = 3
CarpetRegrid2::position_x_1 = -0.0625
CarpetRegrid2::position_y_1 = -0.0625
CarpetRegrid2::position_z_1 = 0
CarpetRegrid2::radius_1[1] = 2.0
CarpetRegrid2::radius_1[2] = 1.0
Grid::type = "coordbase"
CoordBase::xmin = -4.0625
CoordBase::ymin = -4.0625
CoordBase::zmin = 0.0
CoordBase::xmax = 3.9375
CoordBase::ymax = 3.9375
CoordBase::zmax = 4.0
CoordBase::dx = 0.5
CoordBase::dy = 0.5
CoordBase::dz = 0.5
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
ReflectionSymmetry::avoid_origin_z = "no"
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_dtlapse = "zero"
ADMBase::initial_dtshift = "zero"
IDLinearWaves::mvalue = 2
# Gauge
ADMBase::lapse_evolution_method = "ML_BSSN"
ADMBase::shift_evolution_method = "ML_BSSN"
ADMBase::dtlapse_evolution_method = "ML_BSSN"
ADMBase::dtshift_evolution_method = "ML_BSSN"
# Evolution
ADMBase::evolution_method = "ML_BSSN"
ML_BSSN::harmonicN = 1 # 1+log
ML_BSSN::harmonicF = 2.0 # 1+log
ML_BSSN::ShiftGammaCoeff = 0.75
ML_BSSN::BetaDriver = 1.0
ML_BSSN::LapseAdvectionCoeff = 1.0
ML_BSSN::ShiftAdvectionCoeff = 1.0
ML_BSSN::MinimumLapse = 1.0e-8
ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas"
ML_BSSN::fdOrder = 2
# These are necessary
ML_BSSN::ML_log_confac_bound = "flat"
ML_BSSN::ML_metric_bound = "flat"
ML_BSSN::ML_Gamma_bound = "flat"
ML_BSSN::ML_trace_curv_bound = "flat"
ML_BSSN::ML_curv_bound = "flat"
ML_BSSN::ML_lapse_bound = "flat"
ML_BSSN::ML_dtlapse_bound = "flat"
ML_BSSN::ML_shift_bound = "flat"
ML_BSSN::ML_dtshift_bound = "flat"
# Output
IO::out_dir = $parfile
IO::out_fileinfo = "none"
IO::parfile_write = "no"
IOASCII::out1D_every = 1
IOASCII::out1d_yline_x = -0.0625
IOASCII::out1d_xline_y = -0.0625
IOASCII::out1d_zline_x = -0.0625
IOASCII::out1d_zline_y = -0.0625
IOASCII::out1D_vars = "weylscal4::psi4r weylscal4::psi4i"
|
