blob: d2a20b6116bb0822a59768c7bd4c3c605378dbb8 (
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
|
# 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 CarpetRegrid2 CarpetReduce InitBase CarpetIOAscii IOUtil IOBasic ADMBase ADMMacros ADMConstraints GenericFD TmunuBase ADMCoupling StaticConformal CoordGauge SpaceMask IDLinearWaves LoopControl GSL WeylScal4 ReflectionSymmetry MoL InterpToArray CarpetInterp AEILocalInterp"
# Grid
driver::ghost_size = 2
CartGrid3D::avoid_origin = no
Carpet::max_refinement_levels = 3
Carpet::prolongation_order_space = 3 # 1 seems to be buggy in Carpet/git
Carpet::prolongation_order_time = 0
Carpet::domain_from_coordbase = yes
Carpet::poison_new_timelevels = "yes"
Carpet::check_for_poison = "no"
CarpetLib::poison_new_memory = "yes"
CarpetLib::poison_value = 114
InitBase::initial_data_setup_method = init_some_levels
Carpet::init_fill_timelevels = yes
CarpetRegrid2::min_distance = 0
CarpetRegrid2::ensure_proper_nesting = "yes"
CarpetRegrid2::regrid_every = 0
CarpetRegrid2::num_centres = 1
CarpetRegrid2::num_levels_1 = 3
CarpetRegrid2::position_x_1 = -0.0625
CarpetRegrid2::position_y_1 = -0.0625
CarpetRegrid2::position_z_1 = 0
CarpetRegrid2::radius_1[1] = 2.0
CarpetRegrid2::radius_1[2] = 1.0
Grid::type = "coordbase"
CoordBase::xmin = -4.0625
CoordBase::ymin = -4.0625
CoordBase::zmin = 0.0
CoordBase::xmax = 3.9375
CoordBase::ymax = 3.9375
CoordBase::zmax = 4.0
CoordBase::dx = 0.5
CoordBase::dy = 0.5
CoordBase::dz = 0.5
CoordBase::boundary_size_x_lower = 2
CoordBase::boundary_size_y_lower = 2
CoordBase::boundary_size_z_lower = 2
CoordBase::boundary_size_x_upper = 2
CoordBase::boundary_size_y_upper = 2
CoordBase::boundary_size_z_upper = 2
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 = 0
# Initial data
ADMBase::initial_data = "teukwaves"
ADMBase::initial_shift = "zero"
ADMBase::initial_lapse = "one"
IDLinearWaves::mvalue = 2
# Gauge
ADMBase::lapse_evolution_method = "static"
ADMBase::shift_evolution_method = "static"
ADMBase::dtshift_evolution_method = "static"
# Evolution
ADMBase::evolution_method = "static"
ADMBase::lapse_timelevels = 1
ADMBase::shift_timelevels = 1
ADMBase::metric_timelevels = 1
# Output
IO::out_dir = $parfile
IO::out_fileinfo = "none"
IO::parfile_write = "no"
IOASCII::out1D_every = 1
IOASCII::out1d_yline_x = -0.0625
IOASCII::out1d_xline_y = -0.0625
IOASCII::out1d_zline_x = -0.0625
IOASCII::out1d_zline_y = -0.0625
IOASCII::out1D_vars = "weylscal4::psi4r weylscal4::psi4i"
# Interp pseudoscalar to grid to test pseudoscalar
# interpolation treatment
InterpToArray::narrays2d=1
InterpToArray::array2d_vars[0] = "weylscal4::psi4i"
InterpToArray::array2d_x0 = -2.0
InterpToArray::array2d_y0 = -2.0
InterpToArray::array2d_z0 = -2.0
InterpToArray::array2d_dx_i = 0.3
InterpToArray::array2d_npoints_i = 14
InterpToArray::array2d_dz_j = 0.3
InterpToArray::array2d_npoints_j = 14
IOASCII::out2d_every = 1
IOASCII::out2d_vars = "interptoarray::arrays2d[0]"
WeylScal4::timelevels = 1
|
