blob: ae5c73c516e397539771de958a9e01a9aa38753a (
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
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
|
Cactus::cctk_run_title = "Kerr with a 180 degree rotational symmetry"
Cactus::cctk_full_warnings = yes
Cactus::cctk_itlast = 3
ActiveThorns = "InitBase"
ActiveThorns = "Slab"
ActiveThorns = "IOUtil"
IO::out_dir = $parfile
IO::parfile_write = no
ActiveThorns = "PUGH PUGHReduce PUGHSlab"
driver::global_nx = 32
driver::global_ny = 32
driver::global_nz = 18
driver::ghost_size = 2
ActiveThorns = "NaNChecker"
ActiveThorns = "Boundary CartGrid3D CoordBase SymBase ReflectionSymmetry"
CartGrid3D::type = coordbase
CartGrid3D::domain = full
CartGrid3D::avoid_origin = yes
CoordBase::xmin = -4.0+0.5
CoordBase::ymin = -4.0+0.5
CoordBase::zmin = 0.0
CoordBase::xmax = 4.0-0.5
CoordBase::ymax = 4.0-0.5
CoordBase::zmax = 4.0-0.5
CoordBase::dx = 0.25
CoordBase::dy = 0.25
CoordBase::dz = 0.25
CoordBase::boundary_size_x_upper = 2
CoordBase::boundary_size_y_upper = 2
CoordBase::boundary_size_z_upper = 2
CoordBase::boundary_size_x_lower = 2
CoordBase::boundary_size_y_lower = 2
CoordBase::boundary_size_z_lower = 2
CoordBase::boundary_staggered_x_upper = yes
CoordBase::boundary_staggered_y_upper = yes
CoordBase::boundary_staggered_z_upper = yes
CoordBase::boundary_staggered_x_lower = yes
CoordBase::boundary_staggered_y_lower = yes
CoordBase::boundary_staggered_z_lower = yes
ReflectionSymmetry::avoid_origin_x = yes
ReflectionSymmetry::avoid_origin_y = yes
ReflectionSymmetry::avoid_origin_z = yes
ReflectionSymmetry::reflection_z = yes
ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal"
SpaceMask::use_mask = yes
ActiveThorns = "Exact"
ADMBase::initial_data = exact
ADMBase::initial_lapse = exact
ADMBase::initial_dtlapse = exact
ADMBase::initial_shift = exact
ADMBase::initial_dtshift = exact
Exact::exact_model = "Kerr/Kerr-Schild"
Exact::Kerr_KerrSchild__mass = 1.0
Exact::Kerr_KerrSchild__spin = 0.6
Exact::exact_eps = 1e-3
ActiveThorns = "NoExcision SphericalSurface LoopControl"
NoExcision::num_regions = 1
NoExcision::radius[0] = 0.6
NoExcision::Minkowski_scale[0] = 2.6 # value found at r=0.75
NoExcision::reduce_rhs[0] = yes
ActiveThorns = "Time"
Time::dtfac = 0.25
ActiveThorns = "MoL"
MoL::ODE_Method = RK3
ActiveThorns = "ML_BSSN ML_BSSN_Helper GenericFD NewRad TmunuBase LoopControl"
ADMBase::evolution_method = ML_BSSN
ML_BSSN::fdOrder = 2
ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas"
ML_BSSN::my_rhs_boundary_condition = "NewRad"
# These are necessary
ML_BSSN::ML_log_confac_bound = "none"
ML_BSSN::ML_metric_bound = "none"
ML_BSSN::ML_Gamma_bound = "none"
ML_BSSN::ML_trace_curv_bound = "none"
ML_BSSN::ML_curv_bound = "none"
ML_BSSN::ML_lapse_bound = "none"
ML_BSSN::ML_dtlapse_bound = "none"
ML_BSSN::ML_shift_bound = "none"
ML_BSSN::ML_dtshift_bound = "none"
ADMBase::lapse_evolution_method = ML_BSSN
ADMBase::dtlapse_evolution_method = "ML_BSSN"
ML_BSSN::harmonicN = 1 # 1+log
ML_BSSN::harmonicF = 2.0
ML_BSSN::AlphaDriver = 1
ML_BSSN::LapseAdvectionCoeff = 0
ADMBase::shift_evolution_method = ML_BSSN
ADMBase::dtshift_evolution_method = ML_BSSN
ML_BSSN::ShiftGammaCoeff = 0
ActiveThorns = "ADMConstraints"
ActiveThorns = "IOBasic LocalReduce"
IOBasic::outInfo_every = 1
IOBasic::outInfo_reductions = "norm2"
IOBasic::outInfo_vars = "
ADMConstraints::ham
ADMConstraints::momx
"
IOBasic::outScalar_every = 1
IOBasic::outScalar_reductions = "minimum maximum average norm1 norm2 norm_inf"
IOBasic::outScalar_style = gnuplot
IOBasic::outScalar_vars = "
ADMBase::metric
ADMBase::lapse
ADMBase::shift
ADMConstraints::hamiltonian
ADMConstraints::momentum
"
ActiveThorns = "IOASCII"
IOASCII::out1D_every = 1
IOASCII::out1D_style = "gnuplot f(x)"
IOASCII::out1D_vars = "
ADMBase::metric
ADMBase::lapse
ADMBase::shift
ADMConstraints::hamiltonian
ADMConstraints::momentum
SpaceMask::mask
"
|