blob: 1cb441eed11db4cb3759fc8459adc9ac19133abd (
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
|
#==============================================================================
# Cartoon2D, using bitant mode, and ML_BSSN
#==============================================================================
ActiveThorns = "cartoon2d idanalyticbh admbase admcoupling admmacros coordgauge spacemask staticconformal admanalysis ml_bssn ml_bssn_helper newrad genericfd tmunubase loopcontrol mol nanchecker admconstraints pughinterp boundary time pugh pughreduce pughslab cartgrid3d ioascii ioutil iobasic localinterp localreduce coordbase initbase symbase reflectionsymmetry"
## ahfinder
#------------------------------------------------------------------------------
# Run parameters for Cartoon
#------------------------------------------------------------------------------
# Cartoon
cartoon2d::cartoon_active = "yes"
cartoon2d::verbose = "yes"
cartoon2d::stencil = "yes"
cartoon2d::order = 3
# This is for 2 processors, you have to
# change manually by processor numbers.
#driver::processor_topology = "manual"
#driver::processor_topology_3d_x = 1
#driver::processor_topology_3d_y = 1
#driver::processor_topology_3d_z = 2
# Grid
driver::global_nx = 19
driver::global_ny = 5
driver::global_nz = 34
driver::ghost_size_x = 2
driver::ghost_size_y = 2
driver::ghost_size_z = 2
coordbase::boundary_size_x_upper = 2
coordbase::boundary_size_y_upper = 2
coordbase::boundary_size_z_upper = 2
coordbase::boundary_size_x_lower = 2
coordbase::boundary_size_y_lower = 2
coordbase::boundary_size_z_lower = 2
grid::avoid_originy = "no"
grid::type = "byrange"
grid::domain = "full"
grid::xmin = -0.3
grid::xmax = +3.3
grid::ymin = -0.4
grid::ymax = +0.4
grid::zmin = -0.3
grid::zmax = +3.3
coordbase::xmin = -0.3
coordbase::xmax = +3.3
coordbase::ymin = -0.4
coordbase::ymax = +0.4
coordbase::zmin = -0.3
coordbase::zmax = +3.3
reflectionsymmetry::reflection_z = yes
time::dtfac = 0.25
cactus::cctk_itlast = 1
# MODEL
admbase::initial_data = "schwarzschild"
admbase::metric_type = "physical"
#--- Evolution
admbase::evolution_method = "ml_bssn"
methodoflines::ode_method = "icn"
ml_bssn::fdOrder = 2
ml_bssn::my_initial_boundary_condition = "extrapolate-gammas"
ml_bssn::ml_log_confac_bound = "flat"
ml_bssn::ml_metric_bound = "flat"
ml_bssn::ml_gamma_bound = "flat"
ml_bssn::ml_trace_curv_bound = "flat"
ml_bssn::ml_curv_bound = "flat"
ml_bssn::ml_lapse_bound = "flat"
ml_bssn::ml_dtlapse_bound = "flat"
ml_bssn::ml_shift_bound = "flat"
ml_bssn::ml_dtshift_bound = "flat"
admconstraints::cartoon = "yes"
#--- gauge
admbase::lapse_evolution_method = "ML_BSSN"
admbase::dtlapse_evolution_method = "ML_BSSN"
ml_bssn::alphadriver = 1
ml_bssn::lapseadvectioncoeff = 0
ml_bssn::harmonicf = 2.0
## #--- AHFinder.
##
## ahfinder::ahf_active = "yes"
## ahfinder::ahf_cartoon = "yes"
##
## ahfinder::ahf_phi = "no"
## ahfinder::ahf_octant = "high"
##
## ahfinder::ahf_lmax = 0
##
## # Interpolation.
##
## ahfinder::interpolation_order = 1
#------------------------------------------------------------------------------
# Output
#------------------------------------------------------------------------------
IO::out_dir = "test_cartoon_3"
IO::out_fileinfo = "none"
IOBasic::outScalar_every = 1
IOBasic::outScalar_vars = "admconstraints::hamiltonian admconstraints::momentum admbase::alp admanalysis::grr admbase::metric admanalysis::trK ml_bssn::gt11 ml_bssn::phi ml_bssn::trK ml_bssn::ML_Gamma"
IOASCII::out1D_every = 1
IOASCII::out1D_vars = "admconstraints::hamiltonian admconstraints::momentum admbase::alp admanalysis::grr admbase::metric admanalysis::trK ml_bssn::gt11 ml_bssn::phi ml_bssn::trK ml_bssn::Ml_Gamma"
IOBasic::outInfo_every = 1
IOBasic::outInfo_vars = "admbase::alp admbase::dtalp"
#==============================================================================
|