diff options
Diffstat (limited to 'test/bl0.75-025.par')
-rw-r--r-- | test/bl0.75-025.par | 142 |
1 files changed, 142 insertions, 0 deletions
diff --git a/test/bl0.75-025.par b/test/bl0.75-025.par new file mode 100644 index 0000000..52cedeb --- /dev/null +++ b/test/bl0.75-025.par @@ -0,0 +1,142 @@ +# This parameter file sets up Brill-Lindquist 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. + +######################################## + +# +# ***** grid ***** +# + +ActiveThorns = "PUGH" +Driver::global_nx = 43 +Driver::global_ny = 43 +Driver::global_nz = 83 +Driver::ghost_size = 2 + +ActiveThorns = "CoordBase CartGrid3D" +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 = "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 = 0.75 + +IDAnalyticBH::bl_m_2 = 1.0 +IDAnalyticBH::bl_x0_2 = 0.0 +IDAnalyticBH::bl_y0_2 = 0.0 +IDAnalyticBH::bl_z0_2 = -0.75 + +######################################## + +# +# numerical methods +# + +Cactus::cctk_itlast = 0 +ADMBase::metric_type = "static conformal" +StaticConformal::conformal_storage = "factor+derivs+2nd derivs" +SpaceMask::use_mask = "true" + +ADMBase::evolution_method = "none" +ADMBase::initial_lapse = "Cadez" +ADMBase::lapse_evolution_method = "static" +ADMBase::initial_shift = "zero" +ADMBase::shift_evolution_method = "static" + +######################################## + +# +# automagically create output directory +# +ActiveThorns = "IOUtil PUGHSlab IOBasic IOASCII" +IO::out_dir = "bl0.75-025" +IO::parfile_write = "no" +IO::out_fileinfo = "none" + +IOBasic::outscalar_every = 1 +IOBasic::outscalar_vars = "SphericalSurface::sf_valid SphericalSurface::sf_origin SphericalSurface::sf_info" + +IOASCII::out2d_every = 1 +IOASCII::out2d_vars = "SpaceMask::emask SphericalSurface::sf_radius" + +######################################## + +# +# SphericalSurface +# +ActiveThorns = "SphericalSurface" +SphericalSurface::nsurfaces = 3 +SphericalSurface::maxntheta = 62 +SphericalSurface::maxnphi = 62 + +SphericalSurface::nghoststheta[1] = 2 +SphericalSurface::ntheta[1] = 32 +SphericalSurface::nphi [1] = 32 + +SphericalSurface::nghostsphi [2] = 2 +SphericalSurface::ntheta[1] = 32 +SphericalSurface::nphi [1] = 32 + +######################################## + +# +# AHFinderDirect +# + +ActiveThorns = "PUGHInterp PUGHReduce AEILocalInterp AHFinderDirect" + +AHFinderDirect::Theta_norm_for_convergence = 1.0e-12 + +AHFinderDirect::N_horizons = 2 + +AHFinderDirect::set_mask_for_all_horizons = "true" +# only test old-style mask for now, +# since we don't have a thorn to set up the new-style bit fields yet +AHFinderDirect::set_old_style_mask = "true" +AHFinderDirect::set_new_style_mask = "false" +AHFinderDirect::mask_radius_multiplier = 1.0 +AHFinderDirect::mask_buffer_thickness = 4.0 +AHFInderDirect::mask_radius_offset = 0.0 + +AHFinderDirect::which_surface_to_store_info[1] = 1 +AHFinderDirect::N_zones_per_right_angle[1] = 30 +AHFinderDirect::origin_x[1] = 0.0 +AHFinderDirect::origin_y[1] = 0.0 +AHFinderDirect::origin_z[1] = 0.75 +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] = 0.75 +AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.50 + +AHFinderDirect::which_surface_to_store_info[2] = 2 +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 |