aboutsummaryrefslogtreecommitdiff
path: root/par/BrillLindDemo.par
blob: b54850e65e2a92e0192ceb02ee5cea2d95bdb6f9 (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
#------------------------------------------------------------------------------
#
# Parameter file for evolutions of colliding black holes (Brill Lindquist)
#
# For use with Cactus 4.0 (Beta 4)
#
# For help on running Cactus, either run the Cactus exectuable with the 
#   command line option -help, or consult the documentation at
#   www.cactuscode.org
#
# Max Planck Institute for Gravitational Physics, Potsdam, Germany, 1999
#
#------------------------------------------------------------------------------

!DESC "Evolution of Colliding Black Holes: Demo"


# These thorns must be compiled

ActiveThorns = "Boundary ADM CartGrid3D Einstein IDanalyticBH IOBasic IOUtil IOASCII PsiKadelia Pugh pughreduce pughslab Time"
 
#------------------------------------------------------------------------------

# Number of grid points in each direction (must be around 100 or larger)
 
driver::global_nsize       = 100

# How many iterations 

cactus::cctk_final_time    = 30.0

# What output and how often

IOBasic::outInfo_vars      = "einstein::alp"
IOBasic::outInfo_every     = 10
IOASCII::out1D_vars        = "einstein::alp einstein::grr psikadelia::psi0re"
IOASCII::out1D_every       = 10


#------------------------------------------------------------------------------
# 
# Evolution and initial data
#
#------------------------------------------------------------------------------


IO::out_dir                 = "BrillLindDemo-Iso"

grid::domain               = "full"
grid::type                 = "byrange"
grid::xyzmin               = -15.0
grid::xyzmax               =  15.0

cactus::cctk_initial_time  = 0.

time::dtfac                = 0.25 
 
einstein::initial_data     = "bl_bh"  
einstein::evolution_system = "ADM"     
einstein::slicing          = "1+log"
einstein::initial_lapse	   = "one"

idanalyticbh::bl_nbh       =  2
idanalyticbh::bl_x0_1      =  0.0
idanalyticbh::bl_y0_1      =  0.0
idanalyticbh::bl_z0_1      =  1.2
idanalyticbh::bl_M_1       =  1.0
idanalyticbh::bl_x0_2      =  0.0
idanalyticbh::bl_y0_2      =  0.0
idanalyticbh::bl_z0_2      = -1.2
idanalyticbh::bl_M_2       =  1.5

adm::harmonic_f            = 2.0
adm::bound                 = "radiative"      
adm::method                = "stagleap"