blob: c3b95a486face39ff785efebddd07c2f6493a2c7 (
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
|
ActiveThorns = "CoordBase SymBase Boundary NanChecker CartGrid3d Time MoL CarpetIOBasic CarpetIOScalar IOUtil Carpet CarpetLib CarpetReduce CarpetInterp CarpetSlab CarpetIOASCII ADMBase StaticConformal SpaceMask Slab GenericFD CoordGauge ADMCoupling LoopControl ML_BSSN ML_BSSN_Helper GenericFD TMuNuBase SphericalSurface ADMMacros Coordinates CarpetInterp2 Interpolate2 TimerReport IDAnalyticBH NewRad"
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"
# We use the BBH gauge to ensure it is regression-tested. Also, there
# is no harmonic shift in McLachlan
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::conformalMethod = 1 # 1 for W, 0 for phi
ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas"
ML_BSSN::my_rhs_boundary_condition = "static"
ML_BSSN::apply_dissipation = "always"
ML_BSSN::EpsDiss = 0.1
ML_BSSN::fdOrder = 8
GenericFD::jacobian_group = "Coordinates::jacobian"
GenericFD::jacobian_derivative_group = "Coordinates::jacobian2"
GenericFD::jacobian_identity_map = 0
Boundary::radpower = 2
#############################################################
# Grid
#############################################################
Carpet::domain_from_multipatch = yes
CartGrid3D::type = "multipatch"
CartGrid3D::set_coordinate_ranges_on = "all maps"
Coordinates::coordinate_system = "Thornburg04"
Coordinates::h_cartesian = 0.06
Coordinates::h_radial = 0.15
Coordinates::sphere_inner_radius = 0.6
Coordinates::sphere_outer_radius = 1.8
Coordinates::n_angular = 11
Coordinates::patch_boundary_size = 5
Coordinates::additional_overlap_size = 2
Coordinates::outer_boundary_size = 5
Time::timestep_method = "given"
Time::timestep = 0.03
Interpolate::interpolator_order = 5 # Interpatch interpolation
#############################################################
# Carpet
#############################################################
Carpet::ghost_size = 5
Carpet::init_3_timelevels = "no"
#CarpetLib::max_memory_size_MB = 500
Carpet::poison_new_timelevels = yes
CarpetLib::poison_value = 0 # This is needed because the intersection of the outer and interpatch boundary is not currently initialised by anything
NaNChecker::check_every = 1
NaNChecker::verbose = "all"
NaNChecker::action_if_found = "just warn"
NaNChecker::check_vars = "ADMBase::metric ADMBase::curv ADMBase::lapse ADMBase::shift ML_BSSN::ML_curv ML_BSSN::ML_dtlapse ML_BSSN::ML_dtshift ML_BSSN::ML_Gamma ML_BSSN::ML_lapse ML_BSSN::ML_log_confac ML_BSSN::ML_metric ML_BSSN::ML_shift ML_BSSN::ML_trace_curv ML_BSSN::ML_cons_detg ML_BSSN::ML_cons_Gamma ML_BSSN::ML_cons_traceA"
NaNChecker::report_max = 10
#############################################################
# Timers
#############################################################
# Do not commit these as testsuite output as they will change on
# different machines. They are useful for performance monitoring in
# automated build and test systems.
TimerReport::output_all_timers_readable = yes
TimerReport::n_top_timers = 40
TimerReport::output_schedule_timers = no
#############################################################
# Time integration
#############################################################
Cactus::terminate = "iteration"
Cactus::cctk_itlast = 2
# Use 0.5 for RK4 and 0.25 for ICN
#Time::dtfac = 0.5
MethodOfLines::ode_method = "RK4"
MethodOfLines::MoL_Intermediate_Steps = 4
# MethodOfLines::ode_method = "generic"
# MethodOfLines::generic_type = "RK"
# MethodOfLines::MoL_Intermediate_Steps = 1
MethodOfLines::MoL_Num_Scratch_Levels = 1
MethodOfLines::MoL_NaN_Check = "no"
MethodOfLines::initial_data_is_crap = "no"
#############################################################
# Initial data
#############################################################
ADMBase::initial_data = "schwarzschild"
ADMBase::initial_lapse = "one"
ADMBase::initial_shift = "zero"
ADMBase::initial_dtlapse = "zero"
ADMBase::initial_dtshift = "zero"
ADMBase::metric_type = "physical"
IDAnalyticBH::mass = 1
IDAnalyticBH::epsilon = 1e-3
#############################################################
# Output
#############################################################
IO::out_dir = $parfile
IO::out_fileinfo = "none"
IO::parfile_write = "no"
CarpetIOBasic::outInfo_every = 1
CarpetIOBasic::outInfo_vars = "ML_BSSN::ML_log_confac"
IOASCII::out1D_every = 1
IOASCII::out_precision = 19
IOASCII::out1D_x = "no"
IOASCII::out1D_y = "no"
IOASCII::out1D_z = "yes"
IOASCII::out1D_d = "no"
IOASCII::out1D_vars = "ADMBase::gxx ADMBase::kxx ML_BSSN::ML_cons_detg ML_BSSN::ML_cons_Gamma ML_BSSN::ML_cons_traceA ML_BSSN::ML_curv ML_BSSN::ML_dtlapse ML_BSSN::ML_dtshift ML_BSSN::ML_Gamma ML_BSSN::ML_lapse ML_BSSN::ML_log_confac ML_BSSN::ML_metric ML_BSSN::ML_shift ML_BSSN::ML_trace_curv ML_BSSN::At11rhs"
|
