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
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
|
# Parameter definitions for thorn IDAxiBrillBH
# $Header$
shares: ADMBase
EXTENDS KEYWORD initial_lapse
{
"schwarz" :: "Set lapse to Schwarzschild"
}
EXTENDS KEYWORD initial_data
{
"axibrillbh" :: "Axisymmetry Initial data for Black hole + Brill wave"
}
USES KEYWORD metric_type
private:
############################################################
#
# ***** debugging parameters *****
#
#
# If this parameter is set to true, we output an ASCII data file
# giving psi on the 2D grid. The file format should be directly usable
# by a gnuplot 'splot' command.
#
Boolean output_psi2D \
"should we output the conformal factor psi on the 2D grid?"
{
} false
string output_psi2D_file_name \
"if we output the conformal factor psi on the 2D grid, \
what file name should we use for the output file?"
{
".+" :: "any non-empty string that's a valid file name"
} "psi2D.dat"
#
# this parameter controls the amount of (potentially very detailed)
# debugging information this thorn prints
#
CCTK_INT debug \
"level of debugging information to print \
(0 = none, 2 = a little, 6 = a lot, 10 = huge amounts)"
{
0:* :: "any integer >= 0"
} 0
#
# to keep the output size quasi-finite, some debug printing which is
# logicially "per grid point" on the 2-D (eta,q) grid, is actually only
# done for this single 2-D grid point
#
CCTK_INT debug_ii "i coordinate for per-2D-grid-point debug printing"
{
*:* :: "any integer"
} 14
CCTK_INT debug_jj "j coordinate for per-2D-grid-point debug printing"
{
*:* :: "any integer"
} 15
#
# to keep the output size quasi-finite, some debug printing which is
# logicially "per Cactus grid point" is actually only done for this single
# Cactus grid point
#
CCTK_INT debug_i "i coordinate for per-grid-point debug printing"
{
*:* :: "any integer"
} 14
CCTK_INT debug_j "j coordinate for per-grid-point debug printing"
{
*:* :: "any integer"
} 15
CCTK_INT debug_k "k coordinate for per-grid-point debug printing"
{
*:* :: "any integer"
} 10
############################################################
#
# ***** parameters controlling the Brill wave itself *****
#
REAL amp "Brill wave amplitude"
{
*:* :: "No restriction"
} 0.1
REAL eta0 "Brill wave center (in eta coords)"
{
*:* :: "No restriction"
} 0.0
REAL sigma "Brill wave width (in eta)"
{
*:* :: "No restriction"
} 1.0
REAL etamax "eta value for outer edge of grid"
{
*:* :: "No restriction"
} 5.0
INT n "sin^n theta in Brill wave"
{
*:* :: "No restriction"
} 2
############################################################
#
# ***** parameters for the numerical solution of the *****
# ***** Brill-wave equation on the 2-D (eta,q) grid *****
#
INT ne "eta resolution for solve"
{
*:* :: "No restriction"
} 300
INT nq "theta resolution for solve"
{
*:* :: "No restriction"
} 50
REAL error_tolerance "tolerance parameter for elliptic solver"
{
(0:* :: "any positive real number"
} 1.0e-12
############################################################
#
# ***** interpolation parameters *****
#
#
# This thorn first computes the Brill wave solution on an internal 2-D
# grid, then interpolates this to the 3-D Cactus grid. The following
# parameters control this interpolation.
#
STRING interpolator_name \
"name of CCTK_InterpLocalUniform() interpolation operator"
{
".*" :: "any string"
} "uniform cartesian"
STRING interpolator_pars "parameters for the interpolation operator"
{
".*" :: \
"any nonempty string acceptable to Util_TableSetFromString() \
and to the interpolator, or the empty string to use 'order=n', \
where n is specified by the interpolation_order parameter"
} ""
INT interpolation_order "Order for interpolation" STEERABLE = ALWAYS
{
0:9 :: "any integer accepted by the interpolator"
} 1
############################################################
#
# ***** miscellaneous parameters *****
#
Boolean generate_StaticConformal_metric \
"should we generate a StaticConformal conformal metric (true), \
or a pure ADMBase physical metric (false)"
{
} true
|
