blob: 3c92e689dec5e7c3d3c7a8516885a6cd96e329d1 (
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
|
#------------------------------------------------------------------------------
# cartoon
#------------------------------------------------------------------------------
ActiveThorns = " einstein idanalyticbh time pugh pughreduce pughslab cartgrid3d ioascii ioutil iobasic cartoon2d adm_bssn admconstraints boundary "
IOBasic::outScalar_vars = " einstein::metric einstein::curv einstein::lapse einstein::shift einstein::psi adm_bssn::adm_bssn_metric adm_bssn::adm_bssn_curv adm_bssn::adm_bssn_k adm_bssn::adm_bssn_gamma adm_bssn::adm_bssn_phi admconstraints::momy admconstraints::ham "
IOBasic::outInfo_vars = " admconstraints::ham einstein::alp "
IOASCII::out1D_vars = " einstein::metric einstein::curv einstein::lapse einstein::shift einstein::psi adm_bssn::adm_bssn_metric adm_bssn::adm_bssn_curv adm_bssn::adm_bssn_k adm_bssn::adm_bssn_gamma adm_bssn::adm_bssn_phi admconstraints::momy admconstraints::ham "
#==============================================================================
# schwarzschild initial data
#==============================================================================
einstein::initial_data = "bl_bh"
idanalyticbh::bl_M_1 = 1.0
#==============================================================================
#generic machine outputdir
#==============================================================================
IO::out_dir="schw-0050"
#==============================================================================
# Cartoon grid
#==============================================================================
driver::global_nx = 51
driver::global_ny = 5
driver::global_nz = 51
grid::type = "byrange"
grid::xmin = -0.15
grid::xmax = 4.85
grid::ymin = -0.2
grid::ymax = 0.2
grid::zmin = -4.85
grid::zmax = 4.85
time::dtfac = 0.25
grid::domain = "bitant"
grid::bitant_plane = "xy"
driver::ghost_size_x = 2
driver::ghost_size_y = 2
driver::ghost_size_z = 2
cactus::cctk_itlast = 10
grid::avoid_originy = "no"
cartoon2d::stencil = "yes"
cartoon2d::order = 3
#cartoon2d::verbose = "yes"
adm_bssn::cartoon = "yes"
cartoon2d::cartoon_active = "yes"
admconstraints::cartoon = "yes"
#==============================================================================
# Standard ADM_BSSN evolution
#==============================================================================
einstein::evolution_system = "adm_bssn"
adm_bssn::method = "icn"
adm_bssn::ICN_itnum = 3
adm_bssn::stencil_size = 2
einstein::advection = "upwind2"
adm_bssn::bound = "newrad"
boundary::radpower = 3
#==============================================================================
# 1+log
#==============================================================================
einstein::initial_lapse = "one"
einstein::slicing = "1+log"
adm_bssn::lapsesource = "modified"
adm_bssn::harmonic_f = 2.0
#==============================================================================
# Gamma2 shift
#==============================================================================
einstein::initial_shift = "zero"
einstein::shift = "gamma2"
adm_bssn::ShiftGammaCoeff = 0.75
adm_bssn::BetaDriver = 4.0
adm_bssn::ShiftAlpPower = 1.0
adm_bssn::ShiftPsiPower = 4.0
#==============================================================================
# 1D output
#==============================================================================
IOBasic::outScalar_every = 1
IOASCII::out1D_every = 1
IOBasic::outInfo_every = 1
IO::out_fileinfo = "none"
|