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
|
# Initial data and short evolution for gauge wave in Minkowski spacetime
# Author: Mitica Vulcanov <vulcan@aei.mpg.de>
# Author: Roland Haas <roland.haas@physics.gatech.edu>
# $Header$
# Required thorns
ActiveThorns = "CartGrid3D CoordBase SymBase Boundary Time Carpet CarpetLib CarpetSlab CarpetRegrid CarpetReduce InitBase CarpetIOAscii IOUtil IOBasic ADMBase ADMMacros ADMConstraints GenericFD TmunuBase ADMCoupling StaticConformal CoordGauge SpaceMask LocalReduce IDLinearWaves LoopControl GSL WeylScal4 ML_BSSN_O2 ML_BSSN_O2_Helper GenericFD MoL ReflectionSymmetry"
# Grid
driver::ghost_size = 1
CartGrid3D::avoid_origin = no
Carpet::max_refinement_levels = 3
Carpet::prolongation_order_space = 1
Carpet::prolongation_order_time = 2
Carpet::domain_from_coordbase = yes
Carpet::poison_new_timelevels = "yes"
Carpet::check_for_poison = "no"
Carpet::poison_value = 113
CarpetLib::poison_new_memory = "yes"
CarpetLib::poison_value = 114
InitBase::initial_data_setup_method = init_some_levels
Carpet::init_fill_timelevels = yes
CarpetRegrid::refinement_levels = 3
CarpetRegrid::regrid_every = 0
CarpetRegrid::refined_regions = "manual-coordinate-list"
CarpetRegrid::coordinates = "[[([-1.0625,-1.0625,0.0]:[0.9375,0.9375,1.0]:[0.25,0.25,0.25])],[([-0.5625,-0.5625,0.0]:[0.4375,0.4375,0.5]:[0.125,0.125,0.125])]]"
CarpetRegrid::smart_outer_boundaries = "yes"
Grid::type = "coordbase"
CoordBase::xmin = -2.0625
CoordBase::ymin = -2.0625
CoordBase::zmin = 0.0
CoordBase::xmax = 1.9375
CoordBase::ymax = 1.9375
CoordBase::zmax = 2.0
CoordBase::dx = 0.5
CoordBase::dy = 0.5
CoordBase::dz = 0.5
CoordBase::boundary_size_x_lower = 1
CoordBase::boundary_size_y_lower = 1
CoordBase::boundary_size_z_lower = 1
CoordBase::boundary_size_x_upper = 1
CoordBase::boundary_size_y_upper = 1
CoordBase::boundary_size_z_upper = 1
CoordBase::boundary_shiftout_z_lower = 1
CoordBase::boundary_shiftout_z_upper = 0
ReflectionSymmetry::avoid_origin_z = "no"
ReflectionSymmetry::reflection_z = "yes"
cactus::cctk_initial_time = 0
cactus::cctk_itlast = 4
# Initial data
ADMBase::initial_data = "teukwaves"
ADMBase::initial_shift = "zero"
ADMBase::initial_lapse = "one"
IDLinearWaves::mvalue = 2
# Gauge
ADMBase::lapse_evolution_method = "ML_BSSN_O2"
ADMBase::shift_evolution_method = "ML_BSSN_O2"
ADMBase::dtshift_evolution_method = "ML_BSSN_O2"
# Evolution
ADMBase::evolution_method = "ML_BSSN_O2"
ML_BSSN_O2::harmonicF = 2
ML_BSSN_O2::harmonicN = 1 # 1+log
ML_BSSN_O2::ShiftGammaCoeff = 0
# These are necessary
ML_BSSN_O2::ML_log_confac_bound = "flat"
ML_BSSN_O2::ML_metric_bound = "flat"
ML_BSSN_O2::ML_Gamma_bound = "flat"
ML_BSSN_O2::ML_trace_curv_bound = "flat"
ML_BSSN_O2::ML_curv_bound = "flat"
ML_BSSN_O2::ML_lapse_bound = "flat"
ML_BSSN_O2::ML_dtlapse_bound = "flat"
ML_BSSN_O2::ML_shift_bound = "flat"
ML_BSSN_O2::ML_dtshift_bound = "flat"
WeylScal4::Psi4r_group_bound = "flat"
WeylScal4::Psi4i_group_bound = "flat"
# Output
IO::out_dir = $parfile
IO::out_fileinfo = "none"
IO::parfile_write = "no"
IOASCII::out1D_every = 1
IOASCII::out1d_zline_x = -0.0625
IOASCII::out1d_zline_y = -0.0625
IOASCII::out1D_vars = "weylscal4::psi4r weylscal4::psi4i"
|