blob: bf28d186f6897d8f3d9ec96a610b9bcaa355c85a (
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
|
#------------------------------------------------------------------------------
# cartoon, using BSSN_MoL
#------------------------------------------------------------------------------
ActiveThorns = "admbase admcoupling admmacros coordgauge spacemask staticconformal idanalyticbh time pugh pughreduce pughslab cartgrid3d ioascii ioutil iobasic cartoon2d bssn_mol mol admconstraints boundary nanchecker coordbase symbase localreduce"
IOBasic::outScalar_vars = " admbase::metric admbase::curv admbase::lapse admbase::shift staticconformal::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 admbase::alp "
IOASCII::out1D_vars = " admbase::metric admbase::curv admbase::lapse admbase::shift staticconformal::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
#==============================================================================
admbase::initial_data = "bl_bh"
admbase::metric_type = "static conformal"
idanalyticbh::bl_M_1 = 1.0
#==============================================================================
#generic machine outputdir
#==============================================================================
IO::out_dir=$parfile
#==============================================================================
# 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
adm_bssn::cartoon = "yes"
cartoon2d::cartoon_active = "yes"
admconstraints::cartoon = "yes"
#==============================================================================
# Standard ADM_BSSN evolution
#==============================================================================
admbase::evolution_method = "adm_bssn"
methodoflines::ode_method = "icn"
adm_bssn::ICN_itnum = 3
adm_bssn::stencil_size = 2
adm_bssn::advection = "upwind2"
adm_bssn::bound = "newrad"
boundary::radpower = 3
#==============================================================================
# 1+log
#==============================================================================
admbase::initial_lapse = "one"
admbase::lapse_evolution_method = "1+log"
adm_bssn::lapsesource = "modified"
adm_bssn::harmonic_f = 2.0
#==============================================================================
# Gamma2 shift
#==============================================================================
admbase::initial_shift = "zero"
admbase::shift_evolution_method = "gamma2"
adm_bssn::ShiftGammaCoeff = 0.75
adm_bssn::BetaDriver = 4.0
adm_bssn::ShiftAlpPower = 1
adm_bssn::ShiftPsiPower = 4
#==============================================================================
# 1D output
#==============================================================================
IOBasic::outScalar_every = 1
IOASCII::out1D_every = 1
IOBasic::outInfo_every = 1
IO::out_fileinfo = "none"
|