blob: 6c82c3e0ba0a64f0f1946ec78a0ff24a0dbe9a89 (
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
|
# 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. This run takes about 185 megabytes of memory, and
# 36 cpu seconds on a 1.07GHz PIII laptop.
########################################
#
# ***** grid *****
#
ActiveThorns = "CoordBase SymBase CartGrid3D PUGH"
Driver::global_nx = 83
Driver::global_ny = 83
Driver::global_nz = 115
PUGH::ghost_size = 2
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"
ADMBase::evolution_method = "none"
ADMBase::initial_lapse = "Cadez"
ADMBase::lapse_evolution_method = "static"
ADMBase::initial_shift = "zero"
ADMBase::shift_evolution_method = "static"
########################################
ActiveThorns = "IOUtil"
IOUtil::parfile_write = "no"
########################################
#
# AHFinderDirect
#
ActiveThorns = "AEILocalInterp PUGHInterp PUGHReduce SphericalSurface AHFinderDirect"
AHFinderDirect::print_timing_stats = "true"
##AHFinderDirect::verbose_level = "algorithm details"
AHFinderDirect::h_base_file_name = "misner.h"
AHFinderDirect::BH_diagnostics_base_file_name = "misner.BH_diagnostics"
AHFinderDirect::N_horizons = 2
AHFinderDirect::max_allowable_horizon_radius[1] = 1.0 # no problem
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::max_allowable_horizon_radius[2] = 0.8 # should still allow
# finding this AH
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
|
