aboutsummaryrefslogtreecommitdiff
path: root/test/bhns_interp.par
blob: db42acb5102180f85530fceefdbe687ec4238583 (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
ActiveThorns = "Boundary CartGrid3D Time CoordBase SymBase LocalReduce"
ActiveThorns = "Loopcontrol InitBase Carpet CarpetLib CarpetInterp CarpetReduce CarpetSlab CarpetIOASCII CarpetIOBasic IOUtil CarpetIOScalar CarpetIOHDF5"

ActiveThorns = "ADMBase"
ActiveThorns = "StaticConformal"
ActiveThorns = "ADMCoupling ADMMacros"
ActiveThorns = "TmunuBase Constants SpaceMask GRHydro HydroBase"
ActiveThorns = "TOVSolver eos_omni"
ActiveThorns = "TwoPunctures"
ActiveThorns = "coordgauge boundary"

Carpet::init_each_timelevel 		= "no"
InitBase::initial_data_setup_method	= "init_all_levels"

Cactus::cctk_itlast = 0

Time::dtfac			= 0.25

CoordBase::domainsize = minmax

driver::ghost_size                = 3

Carpet::domain_from_coordbase = "yes"
CoordBase::xmin = -60.0
CoordBase::ymin =   0.0
CoordBase::zmin =   0.0
CoordBase::xmax =  60.0
CoordBase::ymax =  60.0
CoordBase::zmax =  60.0
CoordBase::dx   =   3.0
CoordBase::dy   =   3.0
CoordBase::dz   =   3.0

CartGrid3D::type         = "coordbase"
CartGrid3D::domain       = "full"
CartGrid3D::avoid_origin = "no"

ADMBase::initial_data		= "twopunctures"
ADMBase::initial_lapse  = "tov"
ADMBase::initial_shift  = "tov"

TwoPunctures::keep_u_around     = yes
TwoPunctures::schedule_in_ADMBase_InitialData = no

TwoPunctures::npoints_A         = 20
TwoPunctures::npoints_B         = 20
TwoPunctures::npoints_phi       = 4
TwoPunctures::Newton_tol        = 1.e-10
#TwoPunctures::grid_setup_method = "evaluation"
TwoPunctures::TP_Extend_Radius = 1.14

TwoPunctures::par_b             =  30

TwoPunctures::par_m_plus        =  5.0
TwoPunctures::par_m_minus       =  0.0

TwoPunctures::Newton_maxit       =  100
TwoPunctures::use_sources        = 1
TwoPunctures::do_residuum_debug_output = "no"
TwoPunctures::do_initial_debug_output = "no"
TwoPunctures::rescale_sources = "yes"
TwoPunctures::TP_Tiny = 1.e-10

SpaceMask::use_mask = yes

GRHydro::GRHydro_stencil            = 3
GRHydro::GRHydro_eos_type           = "Polytype"
GRHydro::GRHydro_eos_table          = "2D_Polytrope"
GRHydro::recon_method              = "ppm"
GRHydro::rho_abs_min               = 1.e-10
GRHydro::initial_atmosphere_factor = 0.1
#GRHydro::GRHydro_MaxNumSandRVars    = 0
hydrobase::timelevels = 3

TOVSolver::TOV_Enforce_Interpolation = "yes"
TOVSolver::TOV_Num_Radial     = 1000000
TOVSolver::TOV_dr[0]          = 5.e-5
TOVSolver::TOV_Position_x[0]  = -30
TOVSolver::TOV_Rho_Central[0] = 5.e-4
TOVSolver::TOV_Gamma          = 2.0
TOVSolver::TOV_K              = 100.0

IO::out_dir			= $parfile
IO::out_fileinfo= "axis labels"

Carpet::max_refinement_levels = 10
carpet::enable_all_storage        = no
Carpet::prolongation_order_space  = 3
Carpet::prolongation_order_time   = 2
Carpet::convergence_level         = 0
Carpet::init_3_timelevels         = no

CarpetIOScalar::outScalar_reductions	= "minimum maximum norm1 norm2"

IOScalar::outScalar_every	= 512
IOScalar::outScalar_vars		= "admbase::metric admbase::curv TwoPunctures::puncture_u HydroBase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::shift"
IOASCII::out1D_every		= 512
IOASCII::out1D_vars		= "admbase::metric admbase::curv TwoPunctures::puncture_u HydroBase::w_lorentz hydrobase::rho hydrobase::press hydrobase::eps HydroBase::vel GRHydro::dens GRHydro::tau GRHydro::scon admbase::lapse admbase::shift"

admbase::metric_type		= "physical"

boundary::radpower      = 2

admmacros::spatial_order    = 4