blob: 0f1703b1de4693eb1db718b57cacafc2aaddc28f (
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
|
!DESC "distorted BH, no rotation, fixed mesh refinement evolution with carpet"
#
# $Id$
#
# Distorted BH
#
# Frank Herrmann 19.Jun.2003
#
# Edited IH 31.Jul.2003 for Carpet.
#
ActiveThorns = "ahfinderdirect sphericalsurface symbase CoordBase carpetioascii idaxibrillbh spacemask admcoupling admbase staticconformal admmacros coordgauge admanalysis mol bssn_mol admconstraints carpetinterp localinterp aeilocalinterp boundary time carpet carpetlib carpetregrid carpetreduce loopcontrol carpetslab cartgrid3d ioutil carpetiobasic carpetioscalar ellbase nanchecker dissipation localreduce"
#--- grid
time::dtfac = 0.25
grid::type = "coordbase"
grid::avoid_origin = "no"
grid::domain = "bitant"
carpet::domain_from_coordbase = "yes"
coordbase::xmin = -12.45
coordbase::xmax = +12.35
coordbase::ymin = -12.4
coordbase::ymax = +12.4
coordbase::zmin = -0.8
coordbase::zmax = +12.0
coordbase::dx=0.4
coordbase::dy=0.4
coordbase::dz=0.4
coordbase::boundary_size_x_lower = 2 #ghost size
coordbase::boundary_shiftout_x_lower = 1 #due to avoid_origin="no"
coordbase::boundary_size_z_lower = 2
coordbase::boundary_shiftout_z_lower = 1
carpetregrid::smart_outer_boundaries = "yes"
carpet::grid_structure_filename = "grids"
cactus::cctk_itlast = 4
carpet::ghost_size = 2
carpet::max_refinement_levels = 2
carpetregrid::refinement_levels = 2
carpet::prolongation_order_space = 3
carpet::prolongation_order_time = 2
carpet::use_buffer_zones = yes
carpet::additional_buffer_zones = 1
carpet::init_3_timelevels = "yes"
carpetregrid::refined_regions = "manual-coordinate-list"
carpetregrid::coordinates = "[ [ ([-0.2,-0.2,-0.2]:[4.0,4.0,4.0]:[0.2,0.2,0.2]) ] ]"
#--- axisymmetric brill wave+bh
admbase::initial_data = "axibrillbh"
idaxibrillbh::amp = 0.1
idaxibrillbh::eta0 = 0.0
idaxibrillbh::sigma = 1.0
idaxibrillbh::etamax = 6.0
idaxibrillbh::n = 4
idaxibrillbh::ne = 302
idaxibrillbh::nq = 54
admbase::initial_shift = "zero"
admbase::evolution_method = "adm_bssn"
admbase::metric_type = "static conformal"
staticconformal::conformal_storage = "factor+derivs+2nd derivs"
adm_bssn::timelevels = 3
ADMConstraints::constraints_timelevels = 3
ADMConstraints::constraints_persist = yes
#--- numerics
mol::ode_method = "icn"
mol::initial_data_is_crap = "yes"
adm_bssn::bound = "newrad"
adm_bssn::stencil_size = 2
adm_bssn::advection = "upwind2"
boundary::radpower = 3
#--- gauge
admbase::initial_lapse = "one"
admbase::lapse_evolution_method = "1+log"
adm_bssn::lapsesource = "modified"
adm_bssn::harmonic_f = 2.0
adm_bssn::force_lapse_positive = "yes"
adm_bssn::LapsePsiPower = 4
adm_bssn::AlphaDissip = 2.0
admbase::shift_evolution_method = "gamma0"
adm_bssn::ShiftGammaCoeff = 0.75
adm_bssn::BetaDriver = 1.1
adm_bssn::BetaDriverAlpPower = 1
#adm_bssn::BetaDriverAlpPower = 0
adm_bssn::ShiftAlpPower = 1
adm_bssn::ShiftPsiPower = 2
#--- output
io::out_dir = $parfile
io::out_fileinfo = "none"
ioscalar::outScalar_every = 2
ioscalar::outScalar_vars = "admbase::alp admbase::metric admbase::shift admconstraints::momentum admconstraints::ham dissipation::epsdisA"
iobasic::outInfo_every = 2
iobasic::outInfo_vars = "dissipation::epsdisA"
ioascii::out1D_every = 1
ioascii::out1D_vars = "admbase::alp admbase::metric admbase::curv admbase::shift admconstraints::momentum admconstraints::ham dissipation::epsdisA"
## ascii 2D output gets big
#ioascii::out2D_every = 1
#ioascii::out2D_vars = "admbase::alp admbase::metric admbase::curv admbase::shift admconstraints::momentum admconstraints::ham admconstraints::normalized_hamiltonian adm_bssn::adm_bs_phi adm_bssn::adm_bssn_metric adm_bssn::adm_bs_k adm_bssn::adm_bssn_gamma adm_bssn::adm_bssn_curv dissipation::epsdisA"
#--- nanchecker
nanchecker::check_every = 2
nanchecker::check_vars ="admbase::alp admbase::metric admbase::curv admbase::shift"
nanchecker::action_if_found ="terminate"
#--- sphericalsurface
sphericalsurface::nsurfaces = 1
sphericalsurface::maxntheta = 100
sphericalsurface::maxnphi = 100
sphericalsurface::ntheta[0] = 100
sphericalsurface::ntheta[1] = 100
#--- ahfinderdirect
AHFinderDirect::N_horizons = 1
AHFinderDirect::find_every = 1
AHFinderDirect::max_Newton_iterations__initial = 50
AHFinderDirect::max_Newton_iterations__subsequent = 50
AHFinderDirect::initial_guess__coord_sphere__radius [1] = 2.0
ahfinderdirect::which_surface_to_store_info [1] = 0
#--- dissipation
dissipation::ah_slope = 2.0
dissipation::epsdis = 0.2
dissipation::verbose = "no"
dissipation::extra_dissipation_in_horizons = "yes"
dissipation::update_ah_every = 2
dissipation::surface_number[0] = 0
dissipation::horizon_number[0] = 1
dissipation::extra_dissipation_at_outerbound = "yes"
dissipation::outer_bound_npoints = 5
dissipation::ob_slope = 1
Dissipation::vars = "
ADM_BSSN::ADM_BSSN_phi
ADM_BSSN::ADM_BSSN_metric
ADM_BSSN::ADM_BSSN_curv
ADM_BSSN::ADM_BSSN_K
ADM_BSSN::ADM_BSSN_gamma
ADMBase::lapse
ADMBase::shift
"
|
