diff options
Diffstat (limited to 'par/bl-carpet.par')
-rw-r--r-- | par/bl-carpet.par | 149 |
1 files changed, 149 insertions, 0 deletions
diff --git a/par/bl-carpet.par b/par/bl-carpet.par new file mode 100644 index 0000000..e2ef125 --- /dev/null +++ b/par/bl-carpet.par @@ -0,0 +1,149 @@ +# This parameter file sets up Brill-Lindquist initial data and finds the +# apparent horizons in it with AHFinderDirect. + +######################################## + +# +# ***** output setup (also see below for more output stuff) ***** +# +ActiveThorns = "IOUtil IOBasic" +IOUtil::parfile_write = "no" +IO::out_dir = $parfile + +######################################## + +# +# basic grid setup +# + +ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce" +Carpet::domain_from_coordbase = true + +ActiveThorns = "SymBase Boundary" +ActiveThorns = "CoordBase" +CoordBase::domainsize = "minmax" +CoordBase::xmin = -3.0 +CoordBase::xmax = +3.0 +CoordBase::ymin = -3.0 +CoordBase::ymax = +3.0 +CoordBase::zmin = -3.0 +CoordBase::zmax = +3.0 +CoordBase::dx = 0.15 +CoordBase::dy = 0.15 +CoordBase::dz = 0.15 + +# +# ghost zones +# + +driver::ghost_size_x = 2 +driver::ghost_size_y = 2 +driver::ghost_size_z = 2 + +# ghost zones are outside the nominal grid +CoordBase::boundary_internal_x_lower = false +CoordBase::boundary_internal_x_upper = false +CoordBase::boundary_internal_y_lower = false +CoordBase::boundary_internal_y_upper = false +CoordBase::boundary_internal_z_lower = false +CoordBase::boundary_internal_z_upper = false + +# ghost zone widths +CoordBase::boundary_size_x_lower = 2 +CoordBase::boundary_size_x_upper = 2 +CoordBase::boundary_size_y_lower = 2 +CoordBase::boundary_size_y_upper = 2 +CoordBase::boundary_size_z_lower = 2 +CoordBase::boundary_size_z_upper = 2 + +# ghost zones start 0.5 point outside nominal grid boundary +CoordBase::boundary_shiftout_x_lower = 0 +CoordBase::boundary_shiftout_x_upper = 0 +CoordBase::boundary_shiftout_y_lower = 0 +CoordBase::boundary_shiftout_y_upper = 0 +CoordBase::boundary_shiftout_z_lower = 0 +CoordBase::boundary_shiftout_z_upper = 0 + +# staggered grids everywhere +CoordBase::boundary_staggered_x_lower = true +CoordBase::boundary_staggered_x_upper = true +CoordBase::boundary_staggered_y_lower = true +CoordBase::boundary_staggered_y_upper = true +CoordBase::boundary_staggered_z_lower = true +CoordBase::boundary_staggered_z_upper = true + +ActiveThorns = "CartGrid3D" +CartGrid3D::type = "CoordBase" + +######################################## +######################################## + +# +# physics +# + +ActiveThorns = "ADMBase ADMCoupling CoordGauge StaticConformal SpaceMask" +ActiveThorns = "IDAnalyticBH" + +ADMBase::initial_data = "bl_bh" +IDAnalyticBH::bl_nbh = 2 + +IDAnalyticBH::bl_m_1 = 1.0 +IDAnalyticBH::bl_x0_1 = 0.0 +IDAnalyticBH::bl_y0_1 = 0.0 +IDAnalyticBH::bl_z0_1 = 1.0 + +IDAnalyticBH::bl_m_2 = 1.0 +IDAnalyticBH::bl_x0_2 = 0.0 +IDAnalyticBH::bl_y0_2 = 0.0 +IDAnalyticBH::bl_z0_2 = -1.0 + +######################################## + +# +# 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" + +######################################## + +# +# AHFinderDirect +# + +ActiveThorns = "SphericalSurface" +ActiveThorns = "AEILocalInterp AHFinderDirect" + +AHFinderDirect::print_timing_stats = "true" +AHFinderDirect::verbose_level = "algorithm highlights" + +AHFinderDirect::h_base_file_name = "bl.h" +AHFinderDirect::BH_diagnostics_base_file_name = "bl.BH_diagnostics" + +AHFinderDirect::N_horizons = 2 + +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.0 +AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.5 + +AHFinderDirect::origin_x[2] = 0.0 +AHFinderDirect::origin_y[2] = 0.0 +AHFinderDirect::origin_z[2] = -1.0 +AHFinderDirect::initial_guess_method[2] = "coordinate sphere" +AHFinderDirect::initial_guess__coord_sphere__x_center[2] = 0.0 +AHFinderDirect::initial_guess__coord_sphere__y_center[2] = 0.0 +AHFinderDirect::initial_guess__coord_sphere__z_center[2] = -1.0 +AHFinderDirect::initial_guess__coord_sphere__radius[2] = 0.5 |