From 34439126cceb680bbae74cf5972a9e15f7e845a1 Mon Sep 17 00:00:00 2001 From: jthorn Date: Thu, 26 May 2005 10:03:35 +0000 Subject: basic Carpet test of AHFinderDirect git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1409 f88db872-0e4f-0410-b76b-b9085cfa78c5 --- par/Kerr-carpet.par | 135 ++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 135 insertions(+) create mode 100644 par/Kerr-carpet.par (limited to 'par') diff --git a/par/Kerr-carpet.par b/par/Kerr-carpet.par new file mode 100644 index 0000000..cfb3d54 --- /dev/null +++ b/par/Kerr-carpet.par @@ -0,0 +1,135 @@ +# This parameter file sets up Kerr/Kerr-Schild initial data, then +# finds the apparent horizon in it. The local coordinate system origin +# and the initial guess are both deliberately de-centered with respect +# to the black hole, to make this a non-trivial test for the apparent +# horizon finder. + +# flesh +cactus::cctk_itlast = 0 + +######################################## + +# +# ***** 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 = "spacing" +CoordBase::zero_origin_x = true +CoordBase::zero_origin_y = true +CoordBase::zero_origin_z = true +CoordBase::dx = 0.04 +CoordBase::dy = 0.04 +CoordBase::dz = 0.04 +CoordBase::ncells_x = 80 +CoordBase::ncells_y = 80 +CoordBase::ncells_z = 80 + +# +# 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 + +# octant symmetry +ActiveThorns = "ReflectionSymmetry" +ReflectionSymmetry::avoid_origin_x = true +ReflectionSymmetry::avoid_origin_y = true +ReflectionSymmetry::avoid_origin_z = true +ReflectionSymmetry::reflection_x = true +ReflectionSymmetry::reflection_y = true +ReflectionSymmetry::reflection_z = true + +ActiveThorns = "CartGrid3D" +CartGrid3D::type = "CoordBase" + +######################################## + +# +# Kerr initial data +# +ActiveThorns = "ADMBase ADMCoupling StaticConformal CoordGauge SpaceMask Exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_data = "exact" +ADMBase::lapse_evolution_method = "static" +ADMBase::shift_evolution_method = "static" +ADMBase::metric_type = "physical" +Exact::exact_model = "Kerr/Kerr-Schild" +Exact::Kerr_KerrSchild__mass = 1.0 +Exact::Kerr_KerrSchild__spin = 0.6 + +######################################## + +ActiveThorns = "IOUtil" +IOUtil::parfile_write = "no" + +######################################## + +ActiveThorns = "SphericalSurface" +ActiveThorns = "AEILocalInterp AHFinderDirect" + +AHFinderDirect::print_timing_stats = "true" + +AHFinderDirect::output_Theta_every = 1 +AHFinderDirect::h_base_file_name = "Kerr.h" +AHFinderDirect::Theta_base_file_name = "Kerr.Theta" + +AHFinderDirect::N_horizons = 1 +AHFinderDirect::patch_system_type[1] = "+z hemisphere" +AHFinderDirect::origin_x[1] = 0.5 +AHFinderDirect::origin_y[1] = 0.7 +AHFinderDirect::origin_z[1] = 0.0 + +AHFinderDirect::initial_guess_method[1] = "coordinate sphere" +AHFinderDirect::initial_guess__coord_sphere__x_center[1] = -0.2 +AHFinderDirect::initial_guess__coord_sphere__y_center[1] = 0.3 +AHFinderDirect::initial_guess__coord_sphere__z_center[1] = 0.0 +AHFinderDirect::initial_guess__coord_sphere__radius[1] = 2.0 -- cgit v1.2.3