blob: 4162ed6a21126288146163d8868832e4783d6322 (
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
|
# $Header$
# This parameter file sets up Kerr/Kerr-Schild initial data, then
# finds the apparent horizon in it.
# flesh
cactus::cctk_itlast = 2
# 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 = "NoExcision SphericalSurface"
NoExcision::num_regions = 1
NoExcision::radius[0] = 1.01
ActiveThorns = "CarpetIOASCII"
IOASCII::out1D_every = 1
IOASCII::out1D_vars = "ML_BSSN::phi"
IOASCII::out1D_xline_yi = 15
IOASCII::out1D_xline_zi = 3
IOASCII::out1D_yline_xi = 3
IOASCII::out1D_yline_zi = 3
IOASCII::out1D_zline_xi = 3
IOASCII::out1D_zline_yi = 15
IOASCII::out1D_d = no
ActiveThorns = "CarpetIOHDF5"
IOHDF5::checkpoint = yes
IO::checkpoint_dir = $parfile
IO::checkpoint_every = 0
IO::checkpoint_keep = -1
IO::recover = "manual"
IO::recover_dir = "../../../arrangements/EinsteinAnalysis/AHFinderDirect/test/checkpointML-EE"
IO::recover_file = "checkpoint.chkpt.it_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
|
