diff options
author | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2003-05-27 16:01:36 +0000 |
---|---|---|
committer | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2003-05-27 16:01:36 +0000 |
commit | 1c324095eeb108a7f979081a3a6b397b43946ca3 (patch) | |
tree | 2526415d1b1c9e1adda8084b84835beb6a1392f1 /par | |
parent | f4b924c2a9cdaaa7f595ed199d85722855b1d37f (diff) |
add an example par file showing how to set the interpolator parameters
to allow running with driver::ghost_size = 2
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1071 f88db872-0e4f-0410-b76b-b9085cfa78c5
Diffstat (limited to 'par')
-rw-r--r-- | par/misner-ghost2-mask-noshrink.par | 123 |
1 files changed, 123 insertions, 0 deletions
diff --git a/par/misner-ghost2-mask-noshrink.par b/par/misner-ghost2-mask-noshrink.par new file mode 100644 index 0000000..9868432 --- /dev/null +++ b/par/misner-ghost2-mask-noshrink.par @@ -0,0 +1,123 @@ +# This parameter file sets up Misner initial data and finds the +# apparent horizons in it with AHFinderDirect. +# +# This run shows the AHFinderDirect interpolator parameters you need +# to set to run with PUGH::ghost_size = 2. +# +# 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. This run takes about 198 megabytes of memory, and +# 36 cpu seconds on a 1.07GHz PIII laptop. + +######################################## + +# +# ***** grid ***** +# + +ActiveThorns = "CartGrid3D PUGH" + +Driver::global_nx = 83 +Driver::global_ny = 83 +Driver::global_nz = 115 +PUGH::ghost_size = 2 + +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 = "misner_bh" +IDAnalyticBH::mu = 1.2 + +######################################## + +# +# 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" + +######################################## + +ActiveThorns = "PUGHSlab IOBasic IOUtil IOASCII" +IO::out_dir = "misner-ghost2-mask-noshrink" +IO::parfile_write = "no" +IOASCII::out2D_every = 1 +IOASCII::out2D_style = "gnuplot f(x,y)" +IOASCII::out2D_vars = "SpaceMask::emask" +IOASCII::out3D_dir = "/tmp/misner-ghost2-mask-noshrink" +IOASCII::out3D_every = 1 +IOASCII::out3D_style = "gnuplot f(x,y,z)" +IOASCII::out3D_vars = "SpaceMask::emask" + +######################################## + +# +# AHFinderDirect +# + +ActiveThorns = "LocalInterp PUGHInterp PUGHReduce AHFinderDirect" + +AHFinderDirect::find_AHs_at_postinitial = "true" +AHFinderDirect::find_AHs_at_poststep = "false" +AHFinderDirect::print_timing_stats = "true" +##AHFinderDirect::verbose_level = "algorithm details" + +# this next parameter lowers the geometry interpolator order (default=3) +# so AHFinderDirect can run with driver::ghost_size = 2 +AHFinderDirect::geometry_interpolator_pars = "order=2 boundary_off_centering_tolerance={1.0e-10 1.0e-10 1.0e-10 1.0e-10 1.0e-10 1.0e-10} boundary_extrapolation_tolerance={0.0 0.0 0.0 0.0 0.0 0.0}" + +AHFinderDirect::set_mask = "true" +AHFinderDirect::set_old_style_mask = "true" +AHFinderDirect::set_new_style_mask = "false" + +# inside = horizon - 5 grid points +# buffer = from there out to horizon +AHFinderDirect::mask_radius_multiplier = 1.0 +AHFInderDirect::mask_radius_offset = -5.0 +AHFInderDirect::mask_buffer_thickness = 5.0 +AHFinderDirect::mask_is_noshrink = "true" + +AHFinderDirect::h_base_file_name = "misner-ghost2-mask-noshrink.h" +AHFinderDirect::BH_diagnostics_base_file_name = "misner-ghost2-mask-noshrink.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.1 +AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.60 + +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 |