From d119c6f8a127ee737fafd9d0f753b7aa5f4ff953 Mon Sep 17 00:00:00 2001 From: knarf Date: Fri, 2 Dec 2005 10:30:12 +0000 Subject: add missing files git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/Dissipation/trunk@18 850bcc8b-0e4f-0410-8c26-8d28fbf1eda9 --- test/test_ah/test_ah.par | 144 ++++++++++++++++++++++++++++++++++++++++++++++ test/test_ob/test_ob.par | 147 +++++++++++++++++++++++++++++++++++++++++++++++ 2 files changed, 291 insertions(+) create mode 100644 test/test_ah/test_ah.par create mode 100644 test/test_ob/test_ob.par diff --git a/test/test_ah/test_ah.par b/test/test_ah/test_ah.par new file mode 100644 index 0000000..15c9835 --- /dev/null +++ b/test/test_ah/test_ah.par @@ -0,0 +1,144 @@ +!DESC "distorted BH, no rotation, fixed mesh refinement evolution with carpet" +# +# $Id$ +# +# Distorted BH +# +# Frank Herrmann 19.Jun.2003 +# +# Edited IH 31.Jul.2003 for Carpet. +# + +ActiveThorns = "ahfinderdirect sphericalsurface symbase CoordBase carpetioascii idaxibrillbh spacemask admcoupling admbase staticconformal admmacros coordgauge admanalysis mol bssn_mol admconstraints carpetinterp localinterp aeilocalinterp boundary time carpet carpetlib carpetregrid carpetreduce carpetslab cartgrid3d ioutil iobasic ellbase nanchecker dissipation" + +#--- grid +time::dtfac = 0.25 +grid::type = "coordbase" +grid::avoid_origin = "no" +grid::domain = "bitant" +carpet::domain_from_coordbase = "yes" +coordbase::xmin = -12.45 +coordbase::xmax = +12.35 +coordbase::ymin = -12.4 +coordbase::ymax = +12.4 +coordbase::zmin = -0.8 +coordbase::zmax = +12.0 +coordbase::dx=0.4 +coordbase::dy=0.4 +coordbase::dz=0.4 +coordbase::boundary_size_x_lower = 2 #ghost size +coordbase::boundary_shiftout_x_lower = 1 #due to avoid_origin="no" +coordbase::boundary_size_z_lower = 2 +coordbase::boundary_shiftout_z_lower = 1 +carpetregrid::smart_outer_boundaries = "yes" +carpet::grid_structure_filename = "grids" +cactus::cctk_itlast = 4 +carpet::ghost_size = 2 +carpet::max_refinement_levels = 2 +carpetregrid::refinement_levels = 2 +carpet::prolongation_order_space = 3 +carpet::prolongation_order_time = 2 +carpet::buffer_width = 4 +carpet::init_3_timelevels = "yes" +carpetregrid::refined_regions = "manual-coordinate-list" +carpetregrid::coordinates = "[ [ ([-0.2,-0.2,-0.2]:[4.0,4.0,4.0]:[0.2,0.2,0.2]) ] ]" + + +#--- axisymmetric brill wave+bh +admbase::initial_data = "axibrillbh" +idaxibrillbh::amp = 0.1 +idaxibrillbh::eta0 = 0.0 +idaxibrillbh::sigma = 1.0 +idaxibrillbh::etamax = 6.0 +idaxibrillbh::n = 4 +idaxibrillbh::ne = 302 +idaxibrillbh::nq = 54 +admbase::initial_shift = "zero" +admbase::evolution_method = "adm_bssn" +admbase::metric_type = "static conformal" +staticconformal::conformal_storage = "factor+derivs+2nd derivs" + +#--- numerics +mol::ode_method = "icn" +mol::initial_data_is_crap = "yes" +adm_bssn::bound = "newrad" +adm_bssn::stencil_size = 2 +adm_bssn::advection = "upwind2" +boundary::radpower = 3 + +#--- gauge +admbase::initial_lapse = "one" +admbase::lapse_evolution_method = "1+log" +adm_bssn::lapsesource = "modified" +adm_bssn::harmonic_f = 2.0 +adm_bssn::force_lapse_positive = "yes" +adm_bssn::LapsePsiPower = 4 +adm_bssn::AlphaDissip = 2.0 +admbase::shift_evolution_method = "gamma0" +adm_bssn::ShiftGammaCoeff = 0.75 +adm_bssn::BetaDriver = 1.1 +adm_bssn::BetaDriverAlpPower = 1 +#adm_bssn::BetaDriverAlpPower = 0 +adm_bssn::ShiftAlpPower = 1 +adm_bssn::ShiftPsiPower = 2 + +#--- output +io::out_dir = $parfile +io::out_fileinfo = "none" +iobasic::outScalar_every = 2 +iobasic::outScalar_vars = "admbase::alp admbase::metric admbase::shift admconstraints::momentum admconstraints::ham dissipation::epsdisA" +iobasic::outInfo_every = 2 +iobasic::outInfo_vars = "dissipation::epsdisA" +ioascii::out1D_every = 1 +ioascii::out1D_vars = "admbase::alp admbase::metric admbase::curv admbase::shift admconstraints::momentum admconstraints::ham dissipation::epsdisA" +## ascii 2D output gets big +#ioascii::out2D_every = 1 +#ioascii::out2D_vars = "admbase::alp admbase::metric admbase::curv admbase::shift admconstraints::momentum admconstraints::ham admconstraints::normalized_hamiltonian adm_bssn::adm_bs_phi adm_bssn::adm_bssn_metric adm_bssn::adm_bs_k adm_bssn::adm_bssn_gamma adm_bssn::adm_bssn_curv dissipation::epsdisA" + +#--- nanchecker +nanchecker::check_every = 2 +nanchecker::check_vars ="admbase::alp admbase::metric admbase::curv admbase::shift" +nanchecker::action_if_found ="terminate" + + +#--- sphericalsurface +sphericalsurface::nsurfaces = 1 +sphericalsurface::maxntheta = 100 +sphericalsurface::maxnphi = 100 +sphericalsurface::ntheta[0] = 100 +sphericalsurface::ntheta[1] = 100 + + +#--- ahfinderdirect +AHFinderDirect::N_horizons = 1 +AHFinderDirect::find_every = 1 +AHFinderDirect::max_Newton_iterations__initial = 50 +AHFinderDirect::max_Newton_iterations__subsequent = 50 +AHFinderDirect::initial_guess__coord_sphere__radius [1] = 2.0 +ahfinderdirect::which_surface_to_store_info [1] = 0 + + + + +#--- dissipation +dissipation::ah_slope = 2.0 +dissipation::epsdis = 0.2 +dissipation::verbose = "no" +dissipation::extra_dissipation_in_horizons = "yes" +dissipation::update_ah_every = 2 +dissipation::surface_number[0] = 0 +dissipation::horizon_number[0] = 1 +Dissipation::vars = " + ADM_BSSN::ADM_BSSN_phi + ADM_BSSN::ADM_BSSN_metric + ADM_BSSN::ADM_BSSN_curv + ADM_BSSN::ADM_BSSN_K + ADM_BSSN::ADM_BSSN_gamma + ADMBase::lapse + ADMBase::shift +" + +#--- timer info +Cactus::cctk_timer_output = "FULL" +IO::print_timing_info = "yes" + diff --git a/test/test_ob/test_ob.par b/test/test_ob/test_ob.par new file mode 100644 index 0000000..a2872c3 --- /dev/null +++ b/test/test_ob/test_ob.par @@ -0,0 +1,147 @@ +!DESC "distorted BH, no rotation, fixed mesh refinement evolution with carpet" +# +# $Id$ +# +# Distorted BH +# +# Frank Herrmann 19.Jun.2003 +# +# Edited IH 31.Jul.2003 for Carpet. +# + +ActiveThorns = "ahfinderdirect sphericalsurface symbase CoordBase carpetioascii idaxibrillbh spacemask admcoupling admbase staticconformal admmacros coordgauge admanalysis mol bssn_mol admconstraints carpetinterp localinterp aeilocalinterp boundary time carpet carpetlib carpetregrid carpetreduce carpetslab cartgrid3d ioutil iobasic ellbase nanchecker dissipation" + +#--- grid +time::dtfac = 0.25 +grid::type = "coordbase" +grid::avoid_origin = "no" +grid::domain = "bitant" +carpet::domain_from_coordbase = "yes" +coordbase::xmin = -12.45 +coordbase::xmax = +12.35 +coordbase::ymin = -12.4 +coordbase::ymax = +12.4 +coordbase::zmin = -0.8 +coordbase::zmax = +12.0 +coordbase::dx=0.4 +coordbase::dy=0.4 +coordbase::dz=0.4 +coordbase::boundary_size_x_lower = 2 #ghost size +coordbase::boundary_shiftout_x_lower = 1 #due to avoid_origin="no" +coordbase::boundary_size_z_lower = 2 +coordbase::boundary_shiftout_z_lower = 1 +carpetregrid::smart_outer_boundaries = "yes" +carpet::grid_structure_filename = "grids" +cactus::cctk_itlast = 4 +carpet::ghost_size = 2 +carpet::max_refinement_levels = 2 +carpetregrid::refinement_levels = 2 +carpet::prolongation_order_space = 3 +carpet::prolongation_order_time = 2 +carpet::buffer_width = 4 +carpet::init_3_timelevels = "yes" +carpetregrid::refined_regions = "manual-coordinate-list" +carpetregrid::coordinates = "[ [ ([-0.2,-0.2,-0.2]:[4.0,4.0,4.0]:[0.2,0.2,0.2]) ] ]" + + +#--- axisymmetric brill wave+bh +admbase::initial_data = "axibrillbh" +idaxibrillbh::amp = 0.1 +idaxibrillbh::eta0 = 0.0 +idaxibrillbh::sigma = 1.0 +idaxibrillbh::etamax = 6.0 +idaxibrillbh::n = 4 +idaxibrillbh::ne = 302 +idaxibrillbh::nq = 54 +admbase::initial_shift = "zero" +admbase::evolution_method = "adm_bssn" +admbase::metric_type = "static conformal" +staticconformal::conformal_storage = "factor+derivs+2nd derivs" + +#--- numerics +mol::ode_method = "icn" +mol::initial_data_is_crap = "yes" +adm_bssn::bound = "newrad" +adm_bssn::stencil_size = 2 +adm_bssn::advection = "upwind2" +boundary::radpower = 3 + +#--- gauge +admbase::initial_lapse = "one" +admbase::lapse_evolution_method = "1+log" +adm_bssn::lapsesource = "modified" +adm_bssn::harmonic_f = 2.0 +adm_bssn::force_lapse_positive = "yes" +adm_bssn::LapsePsiPower = 4 +adm_bssn::AlphaDissip = 2.0 +admbase::shift_evolution_method = "gamma0" +adm_bssn::ShiftGammaCoeff = 0.75 +adm_bssn::BetaDriver = 1.1 +adm_bssn::BetaDriverAlpPower = 1 +#adm_bssn::BetaDriverAlpPower = 0 +adm_bssn::ShiftAlpPower = 1 +adm_bssn::ShiftPsiPower = 2 + +#--- output +io::out_dir = $parfile +io::out_fileinfo = "none" +iobasic::outScalar_every = 2 +iobasic::outScalar_vars = "admbase::alp admbase::metric admbase::shift admconstraints::momentum admconstraints::ham dissipation::epsdisA" +iobasic::outInfo_every = 2 +iobasic::outInfo_vars = "dissipation::epsdisA" +ioascii::out1D_every = 1 +ioascii::out1D_vars = "admbase::alp admbase::metric admbase::curv admbase::shift admconstraints::momentum admconstraints::ham dissipation::epsdisA" +## ascii 2D output gets big +#ioascii::out2D_every = 1 +#ioascii::out2D_vars = "admbase::alp admbase::metric admbase::curv admbase::shift admconstraints::momentum admconstraints::ham admconstraints::normalized_hamiltonian adm_bssn::adm_bs_phi adm_bssn::adm_bssn_metric adm_bssn::adm_bs_k adm_bssn::adm_bssn_gamma adm_bssn::adm_bssn_curv dissipation::epsdisA" + +#--- nanchecker +nanchecker::check_every = 2 +nanchecker::check_vars ="admbase::alp admbase::metric admbase::curv admbase::shift" +nanchecker::action_if_found ="terminate" + + +#--- sphericalsurface +sphericalsurface::nsurfaces = 1 +sphericalsurface::maxntheta = 100 +sphericalsurface::maxnphi = 100 +sphericalsurface::ntheta[0] = 100 +sphericalsurface::ntheta[1] = 100 + + +#--- ahfinderdirect +AHFinderDirect::N_horizons = 1 +AHFinderDirect::find_every = 1 +AHFinderDirect::max_Newton_iterations__initial = 50 +AHFinderDirect::max_Newton_iterations__subsequent = 50 +AHFinderDirect::initial_guess__coord_sphere__radius [1] = 2.0 +ahfinderdirect::which_surface_to_store_info [1] = 0 + + + + +#--- dissipation +dissipation::ah_slope = 2.0 +dissipation::epsdis = 0.2 +dissipation::verbose = "no" +dissipation::extra_dissipation_in_horizons = "yes" +dissipation::update_ah_every = 2 +dissipation::surface_number[0] = 0 +dissipation::horizon_number[0] = 1 +dissipation::extra_dissipation_at_outerbound = "yes" +dissipation::outer_bound_npoints = 5 +dissipation::ob_slope = 1 +Dissipation::vars = " + ADM_BSSN::ADM_BSSN_phi + ADM_BSSN::ADM_BSSN_metric + ADM_BSSN::ADM_BSSN_curv + ADM_BSSN::ADM_BSSN_K + ADM_BSSN::ADM_BSSN_gamma + ADMBase::lapse + ADMBase::shift +" + +#--- timer info +Cactus::cctk_timer_output = "FULL" +IO::print_timing_info = "yes" + -- cgit v1.2.3