blob: d9b55b54d97cb3371d593ca9908696be3c710490 (
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
124
125
126
127
128
|
# This parameter file sets up Misner initial data and finds the
# apparent horizons in it with AHFinderDirect.
#
# Note that the grid is very high-resolution, but small (the outer
# boundaries are very close in). This isn't a problem for the purposes
# of this test.
########################################
#
# ***** grid *****
#
ActiveThorns = "PUGH"
Driver::global_nx = 43
Driver::global_ny = 43
Driver::global_nz = 83
driver::ghost_size = 2
ActiveThorns = "CoordBase CartGrid3D"
Grid::domain = "octant"
Grid::type = "byspacing"
Grid::dxyz = 0.025
Grid::avoid_origin = "true"
########################################
#
# physics
#
ActiveThorns = "ADMBase ADMCoupling CoordGauge StaticConformal SpaceMask"
ActiveThorns = "IDAnalyticBH"
ADMBase::initial_data = "misner_bh"
IDAnalyticBH::mu = 1.2
########################################
#
# numerical methods
#
Cactus::cctk_itlast = 0
ADMBase::metric_type = "static conformal"
StaticConformal::conformal_storage = "factor+derivs+2nd derivs"
SpaceMask::use_mask = "true"
ADMBase::evolution_method = "none"
ADMBase::initial_lapse = "Cadez"
ADMBase::lapse_evolution_method = "static"
ADMBase::initial_shift = "zero"
ADMBase::shift_evolution_method = "static"
########################################
#
# automagically create output directory
#
ActiveThorns = "IOUtil PUGHSlab IOASCII"
IO::out_dir = "misner1.2-025"
IO::parfile_write = "no"
IO::out_fileinfo = "none"
IOASCII::out2d_every = 1
IOASCII::out2d_vars = "SpaceMask::emask SphericalSurface::sf_radius[1] SphericalSUrface::sf_radius[2]"
########################################
#
# SphericalSurface
#
ActiveThorns = "SphericalSurface"
SphericalSurface::nsurfaces = 3
SphericalSurface::maxntheta = 62
SphericalSurface::maxnphi = 62
SphericalSurface::nghoststheta[1] = 2
SphericalSurface::ntheta[1] = 32
SphericalSurface::nphi [1] = 32
SphericalSurface::nghostsphi [2] = 2
SphericalSurface::ntheta[1] = 62
SphericalSurface::nphi [1] = 62
########################################
#
# AHFinderDirect
#
ActiveThorns = "PUGHInterp PUGHReduce AEILocalInterp AHFinderDirect"
AHFinderDirect::Theta_norm_for_convergence = 1.0e-12
AHFinderDirect::N_horizons = 2
AHFinderDirect::set_mask_for_all_horizons = "true"
# only test old-style mask for now,
# since we don't have a thorn to set up the new-style bit fields yet
AHFinderDirect::set_old_style_mask = "true"
AHFinderDirect::set_new_style_mask = "false"
AHFinderDirect::mask_radius_multiplier = 1.0
AHFinderDirect::mask_buffer_thickness = 4.0
AHFInderDirect::mask_radius_offset = 0.0
AHFinderDirect::which_surface_to_store_info[1] = 1
AHFinderDirect::origin_x[1] = 0.0
AHFinderDirect::origin_y[1] = 0.0
AHFinderDirect::origin_z[1] = 1.0
AHFinderDirect::initial_guess_method[1] = "coordinate sphere"
AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 0.0
AHFinderDirect::initial_guess__coord_sphere__y_center[1] = 0.0
AHFinderDirect::initial_guess__coord_sphere__z_center[1] = 1.1
AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.60
AHFinderDirect::which_surface_to_store_info[2] = 2
AHFinderDirect::N_zones_per_right_angle[2] = 30
AHFinderDirect::origin_x[2] = 0.0
AHFinderDirect::origin_y[2] = 0.0
AHFinderDirect::origin_z[2] = 0.0
AHFinderDirect::initial_guess_method[2] = "coordinate ellipsoid"
AHFinderDirect::initial_guess__coord_ellipsoid__x_center[2] = 0.0
AHFinderDirect::initial_guess__coord_ellipsoid__y_center[2] = 0.0
AHFinderDirect::initial_guess__coord_ellipsoid__z_center[2] = 0.0
AHFinderDirect::initial_guess__coord_ellipsoid__x_radius[2] = 0.6
AHFinderDirect::initial_guess__coord_ellipsoid__y_radius[2] = 0.6
AHFinderDirect::initial_guess__coord_ellipsoid__z_radius[2] = 1.8
|
