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
|
#==============================================================================
# --- Test suite: PMR and regridding ---
# --- Activation of new refinement level at arbitrary iteration ---
#
# The basic setup is an octant grid with 3 reflevels. 1 rl is
# activated at a specified iteration. The total number of
# timesteps is 30.
#
# Initial model is a non-rotating n=3 polytrope, which collapses
# with a 5 % pressure reduction.
#
# It is suggested to set carpetregrid::regrid_every to 16 for
# regridding at an iteration involving all reflevels, or to
# 18 for a iteration involving only the finest reflevel.
#==============================================================================
#==============================================================================
# Thorn activation
#==============================================================================
ActiveThorns = "time carpet carpetlib carpetregrid GSL LoopControl carpetreduce carpetslab carpetinterp cartgrid3d coordbase symbase mol boundary admbase staticconformal spacemask admcoupling coordgauge admmacros whisky whisky_init_data eos_base eos_polytrope iobasic carpetioascii ioutil admconstraints bssn_mol whisky_rnsid carpetiohdf5 nanchecker"
#==============================================================================
# Test parameters
# Play around with these to produce different test cases.
#==============================================================================
carpetregrid::regrid_every = 18
#==============================================================================
# Directories and recovery
#==============================================================================
io::out_dir = $parfile
#==============================================================================
# Time
#==============================================================================
cactus::cctk_itlast = 30
time::dtfac = 0.375
#==============================================================================
# Grid setup
#==============================================================================
grid::type = "byrange"
grid::domain = "octant"
grid::avoid_origin = "no"
driver::global_nx = 30
driver::global_ny = 30
driver::global_nz = 30
driver::ghost_size = 3
grid::xmax = 2200.000000
grid::ymax = 2200.000000
grid::zmax = 2200.000000
carpet::buffer_width = 6
carpet::max_refinement_levels = 4
carpetregrid::refinement_levels = 3
carpetregrid::refined_regions = "manual-gridpoint-list"
carpetregrid::gridpoints = "[ [ ([12,12,12]:[128,128,128]:[4,4,4]) ], [ ([18,18,18]:[76,76,76]:[2,2,2]) ], [ ([21,21,21]:[50,50,50]:[1,1,1]) ] ]"
carpetregrid::outerbounds = "[ [ [[1,0],[1,0],[1,0]] ], [ [[1,0],[1,0],[1,0]] ], [ [[1,0],[1,0],[1,0]] ] ]"
#==============================================================================
# Regridding
#==============================================================================
carpetregrid::num_new_levels = 1
carpetregrid::activate_levels_on_regrid = "fixed"
#==============================================================================
# Initial data: Supermassive star at the critical limit
#==============================================================================
admbase::initial_data = "rnsid"
admbase::initial_lapse = "rnsid"
admbase::initial_shift = "rnsid"
whisky::whisky_rho_central = 7.84e-9
rnsid::accuracy = 1e-7
rnsid::axes_ratio = 1.0
rnsid::zero_shift = no
rnsid::rns_gamma = 1.3333333333333333
rnsid::rns_k = 1.0
#==============================================================================
# Equation of state
#==============================================================================
whisky::whisky_eos_type = "Polytype"
whisky::whisky_eos_table = "2D_Polytrope"
whisky::EoS_Change = "yes"
EOS_Polytrope::eos_gamma = 1.3333333333333333
EOS_Polytrope::eos_k = 0.95
EOS_Polytrope::use_cgs = "no"
EOS_Polytrope::gamma_ini = 1.3333333333333333
whisky::whisky_maxnumsandrvars = 0
#==============================================================================
# Carpet parameters
#==============================================================================
carpet::init_3_timelevels = "yes"
Carpet::prolongation_order_space = 3
Carpet::prolongation_order_time = 2
#==============================================================================
# Whisky parameters
#==============================================================================
whisky::riemann_solver = "Marquina"
whisky::bound = "none"
whisky::recon_method = "ppm"
whisky::ppm_detect = "yes"
whisky::whisky_stencil = 3
whisky::rho_rel_min = 1.e-8
whisky::initial_atmosphere_factor = 1.e2
whisky::initial_rho_abs_min = 5e-17
spacemask::use_mask = yes
#==============================================================================
# BSSN parameters
#==============================================================================
admbase::evolution_method = "adm_bssn"
admbase::lapse_evolution_method = "1+log"
admbase::shift_evolution_method = "static"
# admbase::shift_evolution_method = "gamma2"
adm_bssn::lapsesource = "modified"
adm_bssn::bound = "radiative"
boundary::radpower = 2
adm_bssn::harmonic_f = 2.0
adm_bssn::ShiftGammaCoeff = 0.75
adm_bssn::BetaDriver = 4
admconstraints::bound = "static"
#==============================================================================
# MoL parameters
#==============================================================================
mol::ode_method = "icn"
mol::generic_type = "RK"
mol::mol_intermediate_steps = 3
#==============================================================================
# NaNChecker
#==============================================================================
nanchecker::check_vars = "admbase::gxx"
nanchecker::action_if_found = terminate
nanchecker::check_every = 1
#==============================================================================
# Output: General settings
#==============================================================================
iobasic::out_format = ".13e"
#==============================================================================
# Output: Scalar and 1D
#==============================================================================
iobasic::outinfo_every = 1
iobasic::outscalar_every = 1
ioascii::out1d_every = 4 # 1
iobasic::outinfo_vars = "whisky::rho admbase::alp"
iobasic::outscalar_vars = "admconstraints::ham admconstraints::momentum whisky::whisky_prim_scalar_var whisky::whisky_prim_vector_var whisky::whisky_cons_scalar_var whisky::whisky_cons_vector_var admbase::metric admbase::curv admbase::lapse admbase::shift"
#ioascii::out1d_vars = "admconstraints::ham admconstraints::momentum whisky::whisky_prim_scalar_var whisky::whisky_prim_vector_var whisky::whisky_cons_scalar_var whisky::whisky_cons_vector_var admbase::metric admbase::curv admbase::lapse admbase::shift adm_bssn::adm_bssn_curv adm_bssn::adm_bssn_k adm_bssn::adm_bssn_gamma adm_bssn::adm_bssn_phi adm_bssn::adm_bssn_metric"
ioascii::out1d_vars = "admconstraints::ham admconstraints::momentum whisky::whisky_prim_scalar_var whisky::whisky_prim_vector_var whisky::whisky_cons_scalar_var whisky::whisky_cons_vector_var admbase::metric admbase::curv admbase::lapse admbase::shift"
##==============================================================================
## Output: 3D
##==============================================================================
#
#iohdf5::out3d_every = 1
#
#iohdf5::out3d_vars = "admconstraints::ham admconstraints::momentum whisky::whisky_prim_scalar_var whisky::whisky_prim_vector_var whisky::whisky_cons_scalar_var whisky::whisky_cons_vector_var admbase::metric admbase::curv admbase::lapse admbase::shift adm_bssn::adm_bssn_curv adm_bssn::adm_bssn_k adm_bssn::adm_bssn_gamma adm_bssn::adm_bssn_phi adm_bssn::adm_bssn_metric"
|
