aboutsummaryrefslogtreecommitdiff
path: root/par/CCZ4/rbst_ccz4.par
blob: 3801ba44bc142c2c5ae9dc28c59ff26f1555927c (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
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
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
Cactus::cctk_run_title = "CCZ4 Robust Stability Test"

ActiveThorns = "
  ADMBase
  ADMCoupling
  ADMMacros
  Boundary
  Carpet
  CarpetInterp
  CarpetIOASCII
  CarpetIOBasic
  CarpetIOHDF5
  CarpetIOScalar
  CarpetLib
  CarpetReduce
  CarpetRegrid
  CarpetSlab
  CartGrid3d
  CoordBase
  CoordGauge
  Exact
  GenericFD
  IOUtil
  LoopControl
  ML_BSSN
  ML_BSSN_Helper
  MoL
  NanChecker
  Noise
  Periodic
  Slab
  SpaceMask
  SphericalSurface
  StaticConformal
  SymBase
  Time
  TimerReport
  TMuNuBase
" 

ADMBase::evolution_method         = "ML_BSSN"
ADMBase::lapse_evolution_method   = "ML_BSSN"
ADMBase::shift_evolution_method   = "ML_BSSN"
ADMBase::dtlapse_evolution_method = "ML_BSSN"
ADMBase::dtshift_evolution_method = "ML_BSSN"

ML_BSSN::GammaShift          = 0.0
ML_BSSN::dampk1              = 0.0
ML_BSSN::harmonicN           = 2     
ML_BSSN::harmonicF           = 1.0   
ML_BSSN::ShiftGammaCoeff     = 0.0
ML_BSSN::ShiftBCoeff         = 0.0 # B to zero
ML_BSSN::BetaDriver          = 0.0
ML_BSSN::LapseACoeff         = 0.0 # A to zero
ML_BSSN::LapseAdvectionCoeff = 0.0
ML_BSSN::ShiftAdvectionCoeff = 0.0

ML_BSSN::conformalMethod     = 1 # 1 for W
ML_BSSN::formulation         = 1 # 1 for CCZ4

Boundary::radpower                     = 2

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"
ML_BSSN::ML_Theta_bound    = "none"

#############################################################
# Grid
#############################################################

CoordBase::domainsize                   = minmax

CoordBase::boundary_size_x_lower        = 5
CoordBase::boundary_size_y_lower        = 5
CoordBase::boundary_size_z_lower        = 5
CoordBase::boundary_shiftout_x_lower    = 1
CoordBase::boundary_shiftout_y_lower    = 1
CoordBase::boundary_shiftout_z_lower    = 1

CoordBase::boundary_size_x_upper        = 5
CoordBase::boundary_size_y_upper        = 5
CoordBase::boundary_size_z_upper        = 5
CoordBase::boundary_shiftout_x_upper    = 0
CoordBase::boundary_shiftout_y_upper    = 0
CoordBase::boundary_shiftout_z_upper    = 0

CartGrid3D::type                        = "coordbase"
CartGrid3D::domain                      = "full"
CartGrid3D::avoid_origin                = "no"

Periodic::periodic = "yes"

CoordBase::xmin                         = -0.5
CoordBase::ymin                         = -0.5
CoordBase::zmin                         = -0.5

CoordBase::xmax                         = 0.5
CoordBase::ymax                         = 0.5
CoordBase::zmax                         = 0.5

CoordBase::dx                           = 0.02
CoordBase::dy                           = 0.02
CoordBase::dz                           = 0.02

#############################################################
# CarpetRegrid
#############################################################

CarpetRegrid::refinement_levels           = 1
CarpetRegrid::refined_regions             = "manual-coordinate-list"
CarpetRegrid::smart_outer_boundaries      = yes
CarpetRegrid::keep_same_grid_structure    = "yes"

#############################################################
# Carpet
#############################################################

Carpet::ghost_size                      = 5
Carpet::domain_from_coordbase           = "yes"
Carpet::max_refinement_levels           = 1
Carpet::init_3_timelevels               = "no"
Carpet::prolongation_order_space        = 3
Carpet::prolongation_order_time         = 2
Carpet::use_buffer_zones          = yes
Carpet::num_integrator_substeps         = 4

#############################################################
# CarpetLib
#############################################################

#############################################################
# Time integration
#############################################################

Cactus::terminate                         = "time"
Cactus::cctk_final_time                   = 1000
# Use 0.5 for RK4 and 0.25 for ICN
Time::dtfac                           = 0.5
MethodOfLines::ode_method             = "RK4"
MethodOfLines::MoL_Intermediate_Steps = 4
MethodOfLines::MoL_Num_Scratch_Levels = 1

MethodOfLines::MoL_NaN_Check          = "yes"
MethodOfLines::initial_data_is_crap   = "no"

#############################################################
# Initial data
#############################################################

ADMBase::initial_data = "exact"
ADMBase::initial_lapse = "exact"
ADMBase::initial_shift = "exact"

Exact::exact_model = "Minkowski"

Noise::apply_id_noise = yes
Noise::id_vars = "
  ML_BSSN::ML_log_confac
  ML_BSSN::ML_metric
  ML_BSSN::ML_trace_curv
  ML_BSSN::ML_curv
  ML_BSSN::ML_Gamma
  ML_BSSN::ML_lapse
  ML_BSSN::ML_Theta
"
Noise::Amplitude = 1e-10

#############################################################
# Output
#############################################################

IO::out_dir                   = $parfile
IO::out_fileinfo              = "none"

CarpetIOScalar::outScalar_every      = 100
CarpetIOScalar::outScalar_reductions = "norm_inf norm2"
CarpetIOScalar::outScalar_vars       = "
  ADMBase::alp
  ADMBase::shift
  ADMBase::metric
  ADMBase::curv
  ML_BSSN::ML_log_confac
  ML_BSSN::ML_metric
  ML_BSSN::ML_trace_curv
  ML_BSSN::ML_curv
  ML_BSSN::ML_Gamma
  ML_BSSN::ML_lapse
  ML_BSSN::ML_dtlapse
  ML_BSSN::ML_dtshift
  ML_BSSN::ML_shift
  ML_BSSN::ML_Theta
  ML_BSSN::ML_Ham
  ML_BSSN::ML_mom
  ML_BSSN::ML_cons_detg
  ML_BSSN::ML_cons_traceA
  ML_BSSN::ML_cons_Gamma
"

CarpetIOBasic::outInfo_every        = 100
CarpetIOBasic::outInfo_vars         = "ADMBase::alp"

IOASCII::out1D_every          = 1000
IOASCII::out1D_x              = "yes"
IOASCII::out1D_y              = "yes"
IOASCII::out1D_z              = "yes"
IOASCII::out1D_d              = "no"
IOASCII::out1D_vars           = "
  ADMBase::metric
  ADMBase::curv
  ADMBase::alp
  ADMBase::shift
  ML_BSSN::ML_log_confac
  ML_BSSN::ML_metric
  ML_BSSN::ML_trace_curv
  ML_BSSN::ML_curv
  ML_BSSN::ML_Gamma
  ML_BSSN::ML_lapse
  ML_BSSN::ML_dtlapse
  ML_BSSN::ML_Theta
"

CarpetIOASCII::out_precision = 19