From 470b5c17f5737404c7555280f1817be6c9ae3daa Mon Sep 17 00:00:00 2001 From: jthorn Date: Mon, 14 Aug 2006 15:08:43 +0000 Subject: new par files for Carpet tests git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1462 f88db872-0e4f-0410-b76b-b9085cfa78c5 --- par/Kerr-carpet.par | 17 ++--- par/bl-carpet-mask.par | 170 +++++++++++++++++++++++++++++++++++++++++++++++++ par/bl-carpet.par | 149 +++++++++++++++++++++++++++++++++++++++++++ 3 files changed, 325 insertions(+), 11 deletions(-) create mode 100644 par/bl-carpet-mask.par create mode 100644 par/bl-carpet.par (limited to 'par') diff --git a/par/Kerr-carpet.par b/par/Kerr-carpet.par index 7c13a55..423536a 100644 --- a/par/Kerr-carpet.par +++ b/par/Kerr-carpet.par @@ -27,19 +27,17 @@ Carpet::domain_from_coordbase = true ActiveThorns = "SymBase Boundary" ActiveThorns = "CoordBase" -CoordBase::domainsize = "minmax" # specify continuum domain -CoordBase::xmin = -3.0 +CoordBase::domainsize = "minmax" +CoordBase::xmin = -3.0 # specify full domain for "new" symmetries CoordBase::xmax = +3.0 CoordBase::ymin = -3.0 CoordBase::ymax = +3.0 -CoordBase::zmin = 0.0 +CoordBase::zmin = 0.0 # specify active domain for "old" symmetries CoordBase::zmax = +3.0 CoordBase::dx = 0.2 CoordBase::dy = 0.2 CoordBase::dz = 0.2 -Grid::domain = "bitant" - # # ghost zones # @@ -87,10 +85,12 @@ ReflectionSymmetry::avoid_origin_y = true ReflectionSymmetry::avoid_origin_z = true ReflectionSymmetry::reflection_x = true ReflectionSymmetry::reflection_y = true -ReflectionSymmetry::reflection_z = true +ReflectionSymmetry::reflection_z = false # CartGrid3D will take care of + # this reflection symmetry ActiveThorns = "CartGrid3D" CartGrid3D::type = "CoordBase" +Grid::domain = "bitant" ######################################## @@ -110,11 +110,6 @@ Exact::Kerr_KerrSchild__spin = 0.6 ######################################## -ActiveThorns = "IOUtil" -IOUtil::parfile_write = "no" - -######################################## - ActiveThorns = "SphericalSurface" ActiveThorns = "AEILocalInterp AHFinderDirect" diff --git a/par/bl-carpet-mask.par b/par/bl-carpet-mask.par new file mode 100644 index 0000000..7e482ff --- /dev/null +++ b/par/bl-carpet-mask.par @@ -0,0 +1,170 @@ +# 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" +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" +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" + +######################################## + +# +# apparent horizons +# + +ActiveThorns = "SpaceMask" +SpaceMask::use_mask = true + +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::set_mask_for_all_horizons = true + +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 + +######################################## + +# +# ASCII output +# + +######################################## + +ActiveThorns = "CarpetIOASCII" +CarpetIOASCII::out3D_ghosts = true +CarpetIOASCII::out3D_outer_ghosts = true + +IOASCII::out2d_every = 1 +IOASCII::out2d_vars = "SpaceMask::emask" + +######################################## 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 -- cgit v1.2.3