blob: 762a51403e45bd88bba5395c07283685d35fede5 (
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
|
# $Header$
# This parameter file sets up Kerr/Kerr-Schild initial data, then
# finds the apparent horizon in it.
# flesh
cactus::cctk_itlast = 1
# create output directory
ActiveThorns = "IOUtil"
IO::out_dir = $parfile
IO::parfile_write = "no"
ActiveThorns = "Carpet CarpetLib"
carpet::domain_from_coordbase = "true"
driver::ghost_size = 3
ActiveThorns = "Boundary CartGrid3D CoordBase RotatingSymmetry180 Slab SymBase ReflectionSymmetry"
grid::domain = "full"
grid::avoid_origin = "false"
grid::type = "coordbase"
CoordBase::domainsize = "minmax"
CoordBase::spacing = "numcells"
CoordBase::xmin = 0.0
CoordBase::ymin = -3.0
CoordBase::zmin = 0.0
CoordBase::xmax = 3.0
CoordBase::ymax = 3.0
CoordBase::zmax = 3.0
CoordBase::ncells_x = 19
CoordBase::ncells_y = 31
CoordBase::ncells_z = 19
ReflectionSymmetry::reflection_z = "yes"
ReflectionSymmetry::avoid_origin_z = "no"
CoordBase::boundary_size_x_lower = 3
CoordBase::boundary_size_y_lower = 3
CoordBase::boundary_size_z_lower = 3
CoordBase::boundary_size_x_upper = 3
CoordBase::boundary_size_y_upper = 3
CoordBase::boundary_size_z_upper = 3
CoordBase::boundary_shiftout_x_lower = 1
CoordBase::boundary_shiftout_z_lower = 1
ActiveThorns = "ADMBase ADMCoupling StaticConformal CoordGauge SpaceMask KerrSchild ADMMacros ML_BSSN ML_BSSN_Helper GenericFD TmunuBase NewRad MoL"
ADMBase::initial_lapse = "KerrSchild"
ADMBase::initial_shift = "KerrSchild"
ADMBase::initial_data = "KerrSchild"
ADMBase::initial_dtlapse = "KerrSchild"
ADMBase::initial_dtshift = "KerrSchild"
ADMBase::lapse_evolution_method = "ML_BSSN"
ADMBase::dtlapse_evolution_method = "ML_BSSN"
ADMBase::shift_evolution_method = "static"
ADMBase::metric_type = "physical"
ADMBase::evolution_method = ML_BSSN
KerrSchild::M = 1.0
KerrSchild::a = 0.6
ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas"
ML_BSSN::my_rhs_boundary_condition = "NewRad"
ML_BSSN::fdOrder = 2
# 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"
ML_BSSN::harmonicN = 1 # 1+log
ML_BSSN::LapseAdvectionCoeff = 1
# the default of 3 fails in Recovery since Carpet tries to time-interpolate the
# third timelevel for which is had allocated no metadata, see Carpet issue #626
ML_BSSN::timelevels = 2 # unigrid run
ActiveThorns = "Time"
Time::dtfac = 0.25
ActiveThorns = "SphericalSurface NoExcision"
NoExcision::num_regions = 1
NoExcision::radius[0] = 1.01
ActiveThorns = "CarpetIOASCII"
IOASCII::out1D_every = 1
IOASCII::out1D_vars = "ML_BSSN::phi"
ActiveThorns = "CarpetIOHDF5"
IOHDF5::checkpoint = yes
IO::checkpoint_dir = $parfile
IO::checkpoint_every = 1
IO::checkpoint_keep = -1
########################################
ActiveThorns = "SphericalSurface"
ActiveThorns = "LocalInterp AEILocalInterp CarpetInterp GSL LoopControl CarpetReduce AHFinderDirect"
AHFinderDirect::print_timing_stats = "true"
AHFinderDirect::h_base_file_name = "h"
AHFinderDirect::BH_diagnostics_base_file_name = "BH_diagnostics"
AHFinderDirect::output_OpenDX_control_files = "no"
AHFinderDirect::Theta_norm_for_convergence = 1.0e-12
AHFinderDirect::geometry_interpolator_name = "Hermite polynomial interpolation"
AHFinderDirect::geometry_interpolator_pars = "order=3"
AHFinderDirect::N_horizons = 1
AHFinderDirect::origin_x[1] = 0.5
AHFinderDirect::origin_y[1] = 0.7
AHFinderDirect::origin_z[1] = 0.0
AHFinderDirect::initial_guess_method[1] = "coordinate sphere"
AHFinderDirect::initial_guess__coord_sphere__x_center[1] = -0.2
AHFinderDirect::initial_guess__coord_sphere__y_center[1] = 0.3
AHFinderDirect::initial_guess__coord_sphere__z_center[1] = 0.0
AHFinderDirect::initial_guess__coord_sphere__radius[1] = 2.0
|
