From a6be038b102d951e3a535190b966d214c2405ab6 Mon Sep 17 00:00:00 2001 From: jthorn Date: Mon, 13 Jun 2005 15:28:44 +0000 Subject: some more demo par files for this thorn git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@73 0a4070d5-58f5-498f-b6c0-2693e757fa0f --- par/full-040.par | 246 ++++++++++++++++++++++++++++++++++++++++++++++++++ par/octant-040.par | 258 +++++++++++++++++++++++++++++++++++++++++++++++++++++ 2 files changed, 504 insertions(+) create mode 100644 par/full-040.par create mode 100644 par/octant-040.par diff --git a/par/full-040.par b/par/full-040.par new file mode 100644 index 0000000..b693413 --- /dev/null +++ b/par/full-040.par @@ -0,0 +1,246 @@ +# This parameter file +# sets up BH + Brill-wave initial data +# and evolves it with the standard CactusEinstein and AEIThorns thorns, +# finding the apparent horizon frequently + +######################################## + +# +# ***** output setup (also see below for more output stuff) ***** +# +ActiveThorns = "IOUtil IOBasic" +IOUtil::parfile_write = "no" +IO::out_dir = $parfile + +######################################## + +# +# time integration +# + +##Cactus::cctk_itlast = 0 +Cactus::terminate = "time" +Cactus::cctk_final_time = 10 + +ActiveThorns = "Time" +Time::dtfac = 0.5 # dt = 0.025 + +ActiveThorns = "MoL" +MoL::verbose = "register" + +# RK4 time integration +MoL::ODE_Method = "generic" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 3 + +######################################## + +# +# basic grid setup +# + +ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce" +Carpet::domain_from_coordbase = true + +ActiveThorns = "SymBase" +ActiveThorns = "CoordBase" +CoordBase::domainsize = "spacing" +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 1 point outside nominal grid boundary +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_x_upper = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_y_upper = 1 +CoordBase::boundary_shiftout_z_lower = 1 +CoordBase::boundary_shiftout_z_upper = 1 + +# 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" + +######################################## + +# +# BH + Brill-wave initial data +# + +ActiveThorns = "ADMBase StaticConformal SpaceMask" +ADMBase::metric_type = "static conformal" +ADMBase::evolution_method = "none" +ADMBase::initial_lapse = "one" +ADMBase::initial_shift = "zero" +ADMBase::lapse_evolution_method = "static" +ADMBase::shift_evolution_method = "static" + +ActiveThorns = "AEILocalInterp" +ActiveThorns = "IDAxiBrillBH" + +ADMBase::initial_data = "axibrillbh" +IDAxiBrillBH::amp = 1.0 +IDAxiBrillBH::eta0 = 0.0 +IDAxiBrillBH::sigma = 1.0 +IDAxiBrillBH::etamax = 6.0 +IDAxiBrillBH::n = 2 + +IDAxiBrillBH::ne = 1000 +IDAxiBrillBH::nq = 2000 + +IDAxiBrillBH::interpolator_name = "Lagrange polynomial interpolation" +IDAxiBrillBH::interpolator_pars = "order=4" + +######################################## + +# +# BSSN evolution with 1+log slicing +# +ActiveThorns = "CoordGauge ADMMacros ADMCoupling ADMConstraints ADMAnalysis" +ActiveThorns = "Boundary" +ActiveThorns = "BSSN_MoL" + +ADMBase::initial_lapse = "one" +ADMBase::lapse_evolution_method = "1+log" + +BSSN_MoL::lapsesource = "modified" +BSSN_MoL::harmonic_f = 2.0 +BSSN_MoL::LapsePsiPower = 4 +BSSN_MoL::AlphaDissip = 0.002 + +ADMBase::evolution_method = "ADM_BSSN" + +BSSN_MoL::stencil_size = 2 +BSSN_MoL::advection = "upwind2" + +BSSN_MoL::bound = "radiative" +Boundary::radpower = 2 + +######################################## + +# +# gamma-driver shift +# + +ADMBase::initial_shift = "zero" +ADMBase::shift_evolution_method = "gamma2" + +BSSN_MoL::ShiftGammaCoeff = 0.75 +BSSN_MoL::ShiftAlpPower = 1 +BSSN_MoL::ShiftPsiPower = 4 +BSSN_MoL::ShiftPhiPower = 0 + +BSSN_MoL::BetaDriver = 3.0 # scales with 1/M_adm +BSSN_MoL::BetaDriverPsi = 0.0 +BSSN_MoL::BetaDriverPhi = 0.0 + +######################################## + +# +# find apparent horizons +# + +ActiveThorns = "AEILocalInterp" +ActiveThorns = "SphericalSurface AHFinderDirect" + +# run at poststep so excision mask is set for BSSN et al +AHFinderDirect::run_at_CCTK_ANALYSIS = false +AHFinderDirect::run_at_CCTK_POSTSTEP = true + +AHFinderDirect::N_horizons = 1 + +AHFinderDirect::initial_guess_method[1] = "coordinate sphere" +AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.50 + +AHFinderDirect::find_every = 10 +AHFinderDirect::h_base_file_name = "h" + +AHFinderDirect::set_mask_for_all_horizons = true +AHFinderDirect::mask_radius_multiplier = 1.0 +AHFinderDirect::mask_radius_offset = -5.0 +AHFinderDirect::mask_buffer_thickness = 0.0 +AHFinderDirect::mask_is_noshrink = true +AHFinderDirect::set_old_style_mask = true +AHFinderDirect::set_new_style_mask = false + +######################################## + +# +# Lego Excision +# +ActiveThorns = "LegoExcision" +SpaceMask::use_mask = true +BSSN_MoL::excise = true +BSSN_MoL::excisiontype = "Lego" +ADMConstraints::excise = true + +######################################## + +# +# ASCII output +# + +ActiveThorns = "CarpetIOASCII" +CarpetIOASCII::out3D_ghosts = true +CarpetIOASCII::out3D_outer_ghosts = true + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ADMConstraints::ham + ADMBase::alp + SpaceMask::emask + " +IOASCII::out1d_every = 10 +IOASCII::out1d_vars = " + ADMConstraints::ham + ADMBase::alp ADMBase::shift + ADMAnalysis::grr + ADM_BSSN::ADM_BS_phi ADM_BSSN::ADM_BS_K + SpaceMask::emask + " +IOASCII::out2d_every = 1 +##IOASCII::out2d_style = "gnuplot f(x,y)" +IOASCII::out2d_every = 10 +IOASCII::out2d_vars = " + ADMConstraints::ham + ADMBase::alp ADMBase::shift + ADMAnalysis::grr + ADM_BSSN::ADM_BS_phi ADM_BSSN::ADM_BS_K + SpaceMask::emask + " + +######################################## diff --git a/par/octant-040.par b/par/octant-040.par new file mode 100644 index 0000000..3655803 --- /dev/null +++ b/par/octant-040.par @@ -0,0 +1,258 @@ +# This parameter file +# sets up BH + Brill-wave initial data +# and evolves it with the standard CactusEinstein and AEIThorns thorns, +# finding the apparent horizon frequently + +######################################## + +# +# ***** output setup (also see below for more output stuff) ***** +# +ActiveThorns = "IOUtil IOBasic" +IOUtil::parfile_write = "no" +IO::out_dir = $parfile + +######################################## + +# +# time integration +# + +##Cactus::cctk_itlast = 0 +Cactus::terminate = "time" +Cactus::cctk_final_time = 10 + +ActiveThorns = "Time" +Time::dtfac = 0.5 # dt = 0.025 + +ActiveThorns = "MoL" +MoL::verbose = "register" + +# RK4 time integration +MoL::ODE_Method = "generic" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 3 + +######################################## + +# +# basic grid setup +# + +ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce" +Carpet::domain_from_coordbase = true + +ActiveThorns = "SymBase" +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" + +######################################## + +# +# BH + Brill-wave initial data +# + +ActiveThorns = "ADMBase StaticConformal SpaceMask" +ADMBase::metric_type = "static conformal" +ADMBase::evolution_method = "none" +ADMBase::initial_lapse = "one" +ADMBase::initial_shift = "zero" +ADMBase::lapse_evolution_method = "static" +ADMBase::shift_evolution_method = "static" + +ActiveThorns = "AEILocalInterp" +ActiveThorns = "IDAxiBrillBH" + +ADMBase::initial_data = "axibrillbh" +IDAxiBrillBH::amp = 1.0 +IDAxiBrillBH::eta0 = 0.0 +IDAxiBrillBH::sigma = 1.0 +IDAxiBrillBH::etamax = 6.0 +IDAxiBrillBH::n = 2 + +IDAxiBrillBH::ne = 250 +IDAxiBrillBH::nq = 500 + +IDAxiBrillBH::interpolator_name = "Lagrange polynomial interpolation" +IDAxiBrillBH::interpolator_pars = "order=3" + +######################################## + +# +# BSSN evolution with 1+log slicing +# +ActiveThorns = "ADMCoupling CoordGauge ADMMacros ADMConstraints ADMAnalysis" +ActiveThorns = "Boundary" +ActiveThorns = "BSSN_MoL" + +ADMBase::initial_lapse = "one" +ADMBase::lapse_evolution_method = "1+log" + +BSSN_MoL::lapsesource = "modified" +BSSN_MoL::harmonic_f = 2.0 +BSSN_MoL::LapsePsiPower = 4 +BSSN_MoL::AlphaDissip = 0.002 + +ADMBase::evolution_method = "ADM_BSSN" + +BSSN_MoL::stencil_size = 2 +BSSN_MoL::advection = "upwind2" + +BSSN_MoL::bound = "radiative" +Boundary::radpower = 2 + +######################################## + +# +# gamma-driver shift +# + +ADMBase::initial_shift = "zero" +ADMBase::shift_evolution_method = "gamma2" + +BSSN_MoL::ShiftGammaCoeff = 0.75 +BSSN_MoL::ShiftAlpPower = 1 +BSSN_MoL::ShiftPsiPower = 4 +BSSN_MoL::ShiftPhiPower = 0 + +BSSN_MoL::BetaDriver = 3.0 # scales with 1/M_adm +BSSN_MoL::BetaDriverPsi = 0.0 +BSSN_MoL::BetaDriverPhi = 0.0 + +######################################## + +# +# find apparent horizons +# + +ActiveThorns = "AEILocalInterp" +ActiveThorns = "SphericalSurface AHFinderDirect" + +# run at poststep so excision mask is set for BSSN et al +AHFinderDirect::run_at_CCTK_ANALYSIS = false +AHFinderDirect::run_at_CCTK_POSTSTEP = true + +AHFinderDirect::N_horizons = 1 + +AHFinderDirect::initial_guess_method[1] = "coordinate sphere" +AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.50 + +AHFinderDirect::find_every = 10 +AHFinderDirect::h_base_file_name = "h" + +AHFinderDirect::set_mask_for_all_horizons = true +AHFinderDirect::mask_radius_multiplier = 1.0 +AHFinderDirect::mask_radius_offset = -5.0 +AHFinderDirect::mask_buffer_thickness = 0.0 +AHFinderDirect::mask_is_noshrink = true +AHFinderDirect::set_old_style_mask = true +AHFinderDirect::set_new_style_mask = false + +######################################## + +# +# Lego Excision +# +ActiveThorns = "LegoExcision" +SpaceMask::use_mask = true +BSSN_MoL::excise = true +BSSN_MoL::excisiontype = "Lego" +ADMConstraints::excise = true + +######################################## + +# +# ASCII output +# + +ActiveThorns = "CarpetIOASCII" +CarpetIOASCII::out3D_ghosts = true +CarpetIOASCII::out3D_outer_ghosts = true + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ADMConstraints::ham + ADMBase::alp + SpaceMask::emask + " +IOASCII::out1d_every = 10 +IOASCII::out1d_vars = " + ADMConstraints::ham + ADMBase::alp ADMBase::shift + ADMAnalysis::grr + ADM_BSSN::ADM_BS_phi ADM_BSSN::ADM_BS_K + SpaceMask::emask + " +IOASCII::out2d_every = 1 +##IOASCII::out2d_style = "gnuplot f(x,y)" +IOASCII::out2d_every = 10 +IOASCII::out2d_vars = " + ADMConstraints::ham + ADMBase::alp ADMBase::shift + ADMAnalysis::grr + ADM_BSSN::ADM_BS_phi ADM_BSSN::ADM_BS_K + SpaceMask::emask + " + +######################################## -- cgit v1.2.3