From f1d58e6cd9eab74f429110d6a9fe9d2d961ffb05 Mon Sep 17 00:00:00 2001 From: schnetter Date: Sat, 27 May 2006 00:39:49 +0000 Subject: Example for how to use IDAxiBrillBH without staggering and with mesh refinement. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@97 0a4070d5-58f5-498f-b6c0-2693e757fa0f --- par/octant-refined.par | 203 +++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 203 insertions(+) create mode 100644 par/octant-refined.par diff --git a/par/octant-refined.par b/par/octant-refined.par new file mode 100644 index 0000000..a33b8c8 --- /dev/null +++ b/par/octant-refined.par @@ -0,0 +1,203 @@ + +# Set up the IO Directory +ActiveThorns = "IOUtil" +IO::out_dir = $parfile + +# Checkpoint, if you like +#IO::checkpoint_dir = $parfile + +# TOTAL NUMBER OF TIMESTEPS TO DO +cactus::cctk_itlast = 100 + +# Set up the coordinates and symmetry +ActiveThorns = "CoordBase SymBase" +CoordBase::domainsize = spacing # set domain by grid spacing & number of cells +CoordBase::zero_origin_x = yes # lower bounds are set at 0.0 +CoordBase::zero_origin_y = yes +CoordBase::zero_origin_z = yes +CoordBase::dx = 0.05 # grid spacing +CoordBase::dy = 0.05 +CoordBase::dz = 0.05 +CoordBase::ncells_x = 32 # number of cells (not points) in grid. +CoordBase::ncells_y = 32 +CoordBase::ncells_z = 32 +# These are the offsets between bdry loc and 1st bdry pt at lower x/y/z faces +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 +CoordBase::boundary_size_x_lower = 2 # same as num of ghostzones +CoordBase::boundary_size_y_lower = 2 +CoordBase::boundary_size_z_lower = 2 + +# Set up Carpet stuff. +ActiveThorns = "Carpet CarpetLib CarpetReduce CarpetSlab" +Carpet::domain_from_coordbase = yes # using CoordBase--works better +Carpet::max_refinement_levels = 2 # only coarse and fine. +driver::ghost_size = 2 + +# Here's where the specific coordinates for the fine mesh are set. +ActiveThorns = "CarpetRegrid" +CarpetRegrid::refinement_levels = 2 +CarpetRegrid::refined-regions = manual-coordinate-list +CarpetRegrid::smart_outer_boundaries = yes +CarpetRegrid::coordinates = " [ [ ([0.0, 0.0, 0.0]:[2.0, 2.0, 2.0]:[0.05, 0.05, 0.05]) ] ] " + +# Coordinates? Boundary? +ActiveThorns = "Boundary CartGrid3D" +grid::type = coordbase +grid::avoid_origin = no +grid::domain = octant + +# Set up time stepping. +ActiveThorns = "Time" +Time::dtfac = 0.25 +Cactus::cctk_initial_time = 0.0 + +# Physics, now. +ActiveThorns = "ADMBase StaticConformal PsiKadelia" +ActiveThorns = "SpaceMask ADMCoupling ADMMacros ADMConstraints" +# just to run ADMConstraints! + +# Gauge and Evolution +ActiveThorns = "CoordGauge BSSN_MoL MoL LocalReduce" +ADMBase::lapse_evolution_method = "1+log" +ADMBase::evolution_method = "adm_bssn" # From BSSN_MoL +MoL::ODE_Method = "ICN" +MoL::MoL_Intermediate_Steps = 3 +# Change the above and Carpet::num_integrator_substeps (below) together + +ADM_BSSN::bound = "newrad" + +ADMBase::metric_type = "static conformal" # From StaticConformal + +# The interpolator +ActiveThorns = "LocalInterp" # Needed for IDAxiBrillBH +ActiveThorns = "CarpetInterp" # Needed for IDAxiBrillBH + +Carpet::prolongation_order_space = 3 +Carpet::prolongation_order_time = 2 +Carpet::num_integrator_substeps = 3 # because of ICN +Carpet::use_outer_buffer_zones = yes + + +# Initial Data parameters +ActiveThorns = "IDAxiBrillBH IDAnalyticBH" +admbase::initial_data = "axibrillbh" #From IDAxiBrillBH +idaxibrillbh::amp = 0.5 +idaxibrillbh::eta0 = 0.0 +idaxibrillbh::sigma = 1.0 +idaxibrillbh::etamax = 6.0 +idaxibrillbh::n = 2 +idaxibrillbh::ne = 302 +idaxibrillbh::nq = 54 + + + +# And Excision, because there is a nasty singularity. +ActiveThorns = "LegoExcision" +LegoExcision::fixed_excision = "sphere" +LegoExcision::fixed_size = 0.4 # half of the apparant horizon radius +# the origin is the center, so that's good. +SpaceMask::use_mask = yes # Need these for lego excision +BSSN_MoL::excise = yes +BSSN_MoL::excisionType = lego + + +### Extraction parameters + +ActiveThorns = "Extract" +extract::interpolation_order = 2 + +extract::num_detectors = 2 +extract::itout = 1 + +extract::Nt = 150 +extract::Np = 150 +extract::origin_x = 0. +extract::origin_y = 0. +extract::origin_z = 0. + +extract::all_modes = "no" +extract::l_mode = 2 +extract::m_mode = 0 +extract::detector1 = 5.0 #2 +extract::detector2 = 6.0 #3 + + + +# AH finder. +#ActiveThorns = "AHFinder" +#ahfinder::ahf_active = "yes" +#ahfinder::ahf_flow = "yes" +#ahfinder::ahf_gaussout = "yes" +#ahfinder::ahf_findevery= 10 +#ahfinder::interpolation_order = 2 + + + +# AH finder direct +ActiveThorns = "AHFinderDirect SphericalSurface AEILocalInterp" +#AHFinderDirect::find_every = 10 +AHFinderDirect::find_every = 1 +AHFinderDirect::print_timing_stats = "true" +AHFinderDirect::verbose_level = "algorithm highlights" + +AHFinderDirect::output_Theta_every = 1 +AHFinderDirect::h_base_file_name = "h" +AHFinderDirect::Theta_base_file_name = "Theta" + +AHFinderDirect::N_horizons = 1 +AHFinderDirect::origin_x[1] = 0.0 +AHFinderDirect::origin_y[1] = 0.0 +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] = 1.0 + +AHFinderDirect::geometry_interpolator_name = "Lagrange polynomial interpolation" +AHFinderDirect::geometry_interpolator_pars = "order=4" +AHFinderDirect::surface_interpolator_name = "Lagrange polynomial interpolation" +AHFinderDirect::surface_interpolator_pars = "order=4" + + +# Analysis +ActiveThorns = "ADMAnalysis" +# This += thing doesn't work. +# OUTPUT 1D ASCII +#ioascii::out1D_vars += "admanalysis::spherical_metric" + +# OUTPUT 2D ASCII +#ioascii::out2D_vars += "admanalysis::spherical_metric" + + +# Parameters which will affect performance +# ---------------------------------------- +# Each time admconstraints:: is output thorn ADMConstraints is used +# Each time psikadelia:: is output thorn PsiKadelia is used +# Output can be done in as groups of variables or single variables +# Useful groups are +# einstein::metric - 6 variables +# einstein::conv - 4 variables +# admconstraints::admconstraints - 4 variables + +ActiveThorns = "CarpetIOBasic CarpetIOASCII CarpetIOScalar" + +# OUTPUT TO SCREEN +iobasic::outInfo_every = 2 +#iobasic::outInfo_vars = "ADMBase::alp PsiKadelia::psi4re admconstraints::ham" +IOBasic::outInfo_vars = "ADMBase::alp ADMConstraints::ham" + +# OUTPUT Scalars reductions +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = "ADMBase::alp ADMBase::metric admconstraints::ham PsiKadelia::psi4re" + +# OUTPUT 1D ASCII +ioascii::out1D_every = 1 +ioascii::out1D_vars = "admbase::alp admbase::metric admconstraints::ham PsiKadelia::psi4re" # ADMAnalysis::spherical_metric" + +# OUTPUT 2D ASCII +ioascii::out2D_every = 10 +ioascii::out2D_vars = "admbase::alp admbase::metric admconstraints::ham PsiKadelia::psi4re" # ADMAnalysis::spherical_metric" -- cgit v1.2.3