aboutsummaryrefslogtreecommitdiff
path: root/par/BrillLindDemo.par
blob: 799eb7e186458208df9423670d0464d9c3f275bd (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
#------------------------------------------------------------------------------
#
# 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_BSSN CartGrid3D Einstein IDanalyticBH IOBasic IOUtil IOASCII PsiKadelia Pugh Time"
 

#------------------------------------------------------------------------------


# Number of grid points in each direction
 
driver::global_nsize       = 93

# 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::outdir                 = "BrillLindDemo"

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_BSSN"     
einstein::slicing          = "1+log"
einstein::initial_lapse	   = "one"

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

adm_bssn::harmonic_f       = 2.0
adm_bssn::bound            = "radiative"      
adm_bssn::method           = "icn"
adm_bssn::ICN_itnum        = 3
adm_bssn::evolveG          = "yes"
adm_bssn::RfromH	   = "yes"
adm_bssn::forceKzero	   = "no"
adm_bssn::subtracttrA	   = "no"