Add a parameter file for a real puncture evolution. This works on a

coarser resolution, but is untested at this high resolutions.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@52 b2a53a04-0f4f-0410-87ed-f9f25ced00cf

1 files changed, 518 insertions, 0 deletions

diff --git a/par/pi0.par b/par/pi0.par
new file mode 100644
index 0000000..4d4b732
--- /dev/null
+++ b/par/pi0.par
@@ -0,0 +1,518 @@
+# $Header$
+
+Cactus::cctk_run_title     = "Puncture data"
+Cactus::cctk_full_warnings = yes
+Cactus::cctk_timer_output  = full
+
+#Cactus::terminate       = immediately
+Cactus::terminate       = never
+#Cactus::terminate       = time
+#Cactus::cctk_final_time = 1.0
+#Cactus::terminate       = iteration
+#Cactus::cctk_itlast     = 524288
+
+
+
+ActiveThorns = "Formaline"
+
+ActiveThorns = "LocalInterp AEILocalInterp"
+
+## ActiveThorns = "NaNCatcher"
+
+ActiveThorns = "Slab"
+
+ActiveThorns = "TGRtensor"
+
+
+
+ActiveThorns = "SphericalSurface"
+
+SphericalSurface::nsurfaces = 3
+SphericalSurface::maxntheta =  75
+SphericalSurface::maxnphi   = 148
+
+SphericalSurface::ntheta      [0] =  75
+SphericalSurface::nphi        [0] = 148
+SphericalSurface::nghoststheta[0] =   2
+SphericalSurface::nghostsphi  [0] =   2
+
+SphericalSurface::ntheta      [1] =  75
+SphericalSurface::nphi        [1] = 148
+SphericalSurface::nghoststheta[1] =   2
+SphericalSurface::nghostsphi  [1] =   2
+
+SphericalSurface::ntheta      [2] =  75
+SphericalSurface::nphi        [2] = 148
+SphericalSurface::nghoststheta[2] =   2
+SphericalSurface::nghostsphi  [2] =   2
+
+
+
+ActiveThorns = "IOUtil"
+
+IO::out_dir = $parfile
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce CarpetSlab"
+
+## Carpet::verbose = yes
+Carpet::domain_from_coordbase = yes
+Carpet::max_refinement_levels = 20
+
+driver::ghost_size              = 2
+Carpet::use_outer_buffer_zones  = yes
+Carpet::num_integrator_substeps = 3
+
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time  = 2
+
+Carpet::convergence_level = -1
+
+## Carpet::init_3_timelevels = yes
+
+
+
+ActiveThorns = "Boundary CartGrid3D CoordBase ReflectionSymmetry RotatingSymmetry180 SymBase"
+
+CoordBase::domainsize = minmax
+
+CoordBase::xmin =  0.0
+CoordBase::ymin = -8.0
+CoordBase::zmin =  0.0
+CoordBase::xmax =  8.0
+CoordBase::ymax =  8.0
+CoordBase::zmax =  8.0
+CoordBase::dx   =  0.2
+CoordBase::dy   =  0.2
+CoordBase::dz   =  0.2
+
+CoordBase::boundary_size_x_lower     = 2
+CoordBase::boundary_size_z_lower     = 2
+CoordBase::boundary_shiftout_x_lower = 1
+CoordBase::boundary_shiftout_z_lower = 1
+
+grid::type = coordbase
+
+ReflectionSymmetry::reflection_z   = yes
+ReflectionSymmetry::avoid_origin_z = no
+
+
+
+ActiveThorns = "CarpetRegrid"
+
+CarpetRegrid::refinement_levels        = 4
+CarpetRegrid::smart_outer_boundaries   = yes
+
+CarpetRegrid::keep_same_grid_structure = yes
+CarpetRegrid::refined_regions          = manual-coordinate-list
+CarpetRegrid::coordinates = "
+        [ [ ([-00.0,-04.0,-0.0]:[+05.0,+04.0,+04.0]:[0.100,0.100,0.100]) ],
+          [ ([-00.0,-02.0,-0.0]:[+03.0,+02.0,+02.0]:[0.050,0.050,0.050]) ],
+          [ ([-00.0,-01.0,-0.0]:[+02.0,+01.0,+01.0]:[0.025,0.025,0.025]) ] ]
+"
+
+
+
+ActiveThorns = "NaNChecker"
+
+NaNChecker::check_every     = 524288
+NaNChecker::action_if_found = terminate
+NaNChecker::check_vars      = "
+        ADMBase::gxx
+"
+
+
+
+ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase"
+
+ADMBase::metric_type = "physical"
+
+ADMMacros::spatial_order = 4
+
+SpaceMask::use_mask = yes
+
+
+
+## ActiveThorns = "LegoExcision"
+## 
+## LegoExcision::fixed_excision    = sphere
+## LegoExcision::num_fixed_regions = 2
+## 
+## LegoExcision::fixed_size      =  0.076 # 0.101
+## LegoExcision::fixed_origin_x  =  1.168642873
+## LegoExcision::fixed_origin_y  =  0.0
+## LegoExcision::fixed_origin_z  =  0.0
+## 
+## LegoExcision::fixed2_size     =  0.076 # 0.101
+## LegoExcision::fixed2_origin_x = -1.168642873
+## LegoExcision::fixed2_origin_y =  0.0
+## LegoExcision::fixed2_origin_z =  0.0
+
+
+
+## ActiveThorns = "NoExcision"
+## 
+## NoExcision::num_regions = 2
+## 
+## NoExcision::centre_x            [0] =  1.168642873
+## NoExcision::centre_y            [0] =  0.0
+## NoExcision::centre_z            [0] =  0.0
+## NoExcision::radius              [0] =  0.101
+## NoExcision::overwrite_lapse     [0] = no
+## NoExcision::smoothing_zone_width[0] =  0.05
+## 
+## NoExcision::centre_x            [1] = -1.168642873
+## NoExcision::centre_y            [1] =  0.0
+## NoExcision::centre_z            [1] =  0.0
+## NoExcision::radius              [1] =  0.101
+## NoExcision::overwrite_lapse     [1] = no
+## NoExcision::smoothing_zone_width[1] =  0.05
+
+
+
+ActiveThorns = "TwoPunctures"
+
+ADMBase::initial_data  = "twopunctures"
+ADMBase::initial_lapse = "twopunctures-averaged"
+
+TwoPunctures::verbose        = yes
+TwoPunctures::par_b          =  1.168642873
+TwoPunctures::par_m_plus     =  0.453
+TwoPunctures::par_m_minus    =  0.453
+TwoPunctures::par_P_plus [1] =  0.3331917498
+TwoPunctures::par_P_minus[1] = -0.3331917498
+
+TwoPunctures::TP_epsilon     =  0.101
+
+
+
+ActiveThorns = "RigidRotation"
+
+ADMBase::initial_shift = "rigid rotation"
+
+RigidRotation::omega = 0.06
+
+
+
+ActiveThorns = "Time"
+
+Time::dtfac = 0.25
+
+
+
+ActiveThorns = "MoL"
+
+MoL::initial_data_is_crap = yes
+
+MoL::ODE_Method = RK3
+
+
+
+ActiveThorns = "BSSN_MoL"
+
+## BSSN_MoL::excise       = yes
+## BSSN_MoL::excisionType = lego
+
+ADMBase::evolution_method = ADM_BSSN
+
+ADM_BSSN::stencil_size = 2
+ADM_BSSN::advection    = center4
+ADM_BSSN::bound        = newrad
+
+## ADMBase::lapse_evolution_method = 1+log
+## ADM_BSSN::lapsesource           = modified
+## ADM_BSSN::harmonic_f            = 2.0
+ADMBase::lapse_evolution_method = harmonic
+ADM_BSSN::lapsesource           = modified
+ADM_BSSN::harmonic_f            = 1.0
+ADM_BSSN::force_lapse_positive  = yes
+
+ADMBase::shift_evolution_method = gamma0
+ADM_BSSN::ShiftGammaCoeff       = 0.75
+ADM_BSSN::BetaDriver            = 4.0
+
+
+
+## ActiveThorns = "DriftCorrect4"
+## 
+## DriftCorrect4::nsources = 1
+## 
+## DriftCorrect4::position_x[0]                   = 1.168642873
+## DriftCorrect4::position_y[0]                   = 0.0
+## DriftCorrect4::surface_index[0]                = 0
+## DriftCorrect4::first_driftcorrect_iteration[0] = 524288
+## DriftCorrect4::driftcorrect_every[0]           = 524288
+## DriftCorrect4::do_azimuthal_correction[0]      = yes
+## DriftCorrect4::azimuthal_timescale[0]          = 1.0
+## DriftCorrect4::do_radial_correction[0]         = yes
+## DriftCorrect4::radial_timescale[0]             = 1.0
+
+
+
+ActiveThorns = "Dissipation"
+
+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
+        ADM_BSSN::ADM_BSSN_dtlapse
+        ADMBase::shift
+        ADM_BSSN::ADM_BSSN_B
+"
+
+
+
+ActiveThorns = "ADMAnalysis"
+
+
+
+ActiveThorns = "ADMConstraints"
+
+## ADMConstraints::excise = yes
+
+
+
+ActiveThorns = "AHFinderDirect"
+
+AHFinderDirect::N_horizons                        = 3
+AHFinderDirect::find_every                        = 65536
+AHFinderDirect::max_Newton_iterations__initial    = 20
+AHFinderDirect::max_Newton_iterations__subsequent = 20
+
+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"
+
+# individual AH
+AHFinderDirect::origin_x                                 [1] =  1.168642873
+AHFinderDirect::initial_guess__coord_sphere__x_center    [1] =  1.168642873
+AHFinderDirect::initial_guess__coord_sphere__radius      [1] =  0.2
+AHFinderDirect::reset_horizon_after_not_finding          [1] = no
+AHFinderDirect::which_surface_to_store_info              [1] = 0
+
+# outer common AH
+AHFinderDirect::find_after_individual                    [2] = 178257920
+AHFinderDirect::initial_guess__coord_sphere__radius      [2] =  2.2
+AHFinderDirect::max_allowable_horizon_radius             [2] =  4.0
+AHFinderDirect::smoothing_factor                         [2] =  0.25
+AHFinderDirect::reset_horizon_after_not_finding          [2] = no
+AHFinderDirect::which_surface_to_store_info              [2] = 1
+
+# inner common AH
+AHFinderDirect::find_after_individual                    [3] = 178257920
+AHFinderDirect::initial_guess__coord_sphere__radius      [3] =  2.2
+AHFinderDirect::max_allowable_horizon_radius             [3] =  4.0
+AHFinderDirect::smoothing_factor                         [3] = -0.25
+AHFinderDirect::reset_horizon_after_not_finding          [3] = no
+AHFinderDirect::which_surface_to_store_info              [3] = 2
+
+
+
+ActiveThorns = "IsolatedHorizon"
+
+IsolatedHorizon::verbose          = yes
+IsolatedHorizon::num_horizons     = 3
+IsolatedHorizon::surface_index[0] = 0
+IsolatedHorizon::surface_index[1] = 1
+IsolatedHorizon::surface_index[2] = 2
+
+
+
+ActiveThorns = "CarpetIOBasic"
+
+IOBasic::outInfo_every      = 1
+IOBasic::outInfo_reductions = "norm2"
+IOBasic::outInfo_vars       = "
+        ADMConstraints::ham
+        ADMConstraints::momx
+        IsolatedHorizon::ih_spin[0]
+        IsolatedHorizon::ih_radius[0]
+        IsolatedHorizon::ih_mass[0]
+"
+
+
+
+ActiveThorns = "CarpetIOScalar"
+
+IOScalar::outScalar_every = 1
+IOScalar::outScalar_vars  = "
+        grid::coordinates
+        SpaceMask::emask
+        ADM_BSSN::ADM_BSSN_phi
+        ADM_BSSN::ADM_BSSN_metric
+        ADM_BSSN::ADM_BSSN_K
+        ADM_BSSN::ADM_BSSN_curv
+        ADM_BSSN::ADM_BSSN_gamma
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMAnalysis::trK
+        ADMAnalysis::Ricci
+        ADMConstraints::hamiltonian
+        ADMConstraints::normalized_hamiltonian
+        ADMConstraints::momentum
+        SphericalSurface::sf_radius
+        IsolatedHorizon::ih_shapes
+        IsolatedHorizon::ih_coordinates
+        IsolatedHorizon::ih_tetrad_l
+        IsolatedHorizon::ih_tetrad_n
+        IsolatedHorizon::ih_tetrad_m
+        IsolatedHorizon::ih_newman_penrose
+        IsolatedHorizon::ih_weyl_scalars
+        IsolatedHorizon::ih_ricci_scalars
+        IsolatedHorizon::ih_twometric
+        IsolatedHorizon::ih_killing_vector
+        IsolatedHorizon::ih_killed_twometric
+        IsolatedHorizon::ih_invariant_coordinates
+        IsolatedHorizon::ih_3determinant
+        IsolatedHorizon::ih_fluxes
+"
+
+
+
+ActiveThorns = "CarpetIOASCII"
+
+IOASCII::out0D_every = 1
+IOASCII::out0D_vars  = "
+        grid::coordinates
+        SpaceMask::emask
+        ADM_BSSN::ADM_BSSN_phi
+        ADM_BSSN::ADM_BSSN_metric
+        ADM_BSSN::ADM_BSSN_K
+        ADM_BSSN::ADM_BSSN_curv
+        ADM_BSSN::ADM_BSSN_gamma
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMAnalysis::trK
+        ADMAnalysis::Ricci
+        ADMConstraints::hamiltonian
+        ADMConstraints::normalized_hamiltonian
+        ADMConstraints::momentum
+        SphericalSurface::sf_valid
+        SphericalSurface::sf_info
+        SphericalSurface::sf_radius
+        SphericalSurface::sf_origin
+        SphericalSurface::sf_coordinate_descriptors
+        DriftCorrect4::dc4_state
+        DriftCorrect4::dc4_current_time
+        DriftCorrect4::dc4_omega
+        DriftCorrect4::dc4_radius
+        IsolatedHorizon::ih_state
+        IsolatedHorizon::ih_grid_int
+        IsolatedHorizon::ih_grid_real
+        IsolatedHorizon::ih_shapes
+        IsolatedHorizon::ih_coordinates
+        IsolatedHorizon::ih_tetrad_l
+        IsolatedHorizon::ih_tetrad_n
+        IsolatedHorizon::ih_tetrad_m
+        IsolatedHorizon::ih_newman_penrose
+        IsolatedHorizon::ih_weyl_scalars
+        IsolatedHorizon::ih_ricci_scalars
+        IsolatedHorizon::ih_twometric
+        IsolatedHorizon::ih_killing_vector
+        IsolatedHorizon::ih_killed_twometric
+        IsolatedHorizon::ih_scalars
+        IsolatedHorizon::ih_invariant_coordinates
+        IsolatedHorizon::ih_multipole_moments
+        IsolatedHorizon::ih_3determinant
+        IsolatedHorizon::ih_fluxes
+"
+
+IO::out_yline_x = 6.0
+IO::out_zline_x = 6.0
+IOASCII::out1D_every = 1
+IOASCII::out1D_vars  = "
+        grid::coordinates
+        SpaceMask::emask
+        ADM_BSSN::ADM_BSSN_phi
+        ADM_BSSN::ADM_BSSN_metric
+        ADM_BSSN::ADM_BSSN_K
+        ADM_BSSN::ADM_BSSN_curv
+        ADM_BSSN::ADM_BSSN_gamma
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMAnalysis::trK
+        ADMAnalysis::Ricci
+        ADMConstraints::hamiltonian
+        ADMConstraints::normalized_hamiltonian
+        ADMConstraints::momentum
+        SphericalSurface::sf_radius
+        IsolatedHorizon::ih_shapes
+        IsolatedHorizon::ih_coordinates
+        IsolatedHorizon::ih_tetrad_l
+        IsolatedHorizon::ih_tetrad_n
+        IsolatedHorizon::ih_tetrad_m
+        IsolatedHorizon::ih_newman_penrose
+        IsolatedHorizon::ih_weyl_scalars
+        IsolatedHorizon::ih_ricci_scalars
+        IsolatedHorizon::ih_twometric
+        IsolatedHorizon::ih_killing_vector
+        IsolatedHorizon::ih_killed_twometric
+        IsolatedHorizon::ih_invariant_coordinates
+        IsolatedHorizon::ih_3determinant
+        IsolatedHorizon::ih_fluxes
+"
+
+IO::out_yzplane_x        = 6.0
+IOASCII::out2D_criterion = time
+IOASCII::out2D_dt        = 0.8
+IOASCII::out2D_vars      = "
+        SphericalSurface::sf_radius
+        IsolatedHorizon::ih_shapes
+        IsolatedHorizon::ih_coordinates
+        IsolatedHorizon::ih_tetrad_l
+        IsolatedHorizon::ih_tetrad_n
+        IsolatedHorizon::ih_tetrad_m
+        IsolatedHorizon::ih_newman_penrose
+        IsolatedHorizon::ih_weyl_scalars
+        IsolatedHorizon::ih_ricci_scalars
+        IsolatedHorizon::ih_twometric
+        IsolatedHorizon::ih_killing_vector
+        IsolatedHorizon::ih_killed_twometric
+        IsolatedHorizon::ih_invariant_coordinates
+        IsolatedHorizon::ih_3determinant
+        IsolatedHorizon::ih_fluxes
+"
+
+
+
+ActiveThorns = "CarpetIOHDF5"
+
+IOHDF5::out_criterion = time
+IOHDF5::out_dt        = 0.8
+IOHDF5::out_vars      = "
+        SpaceMask::emask
+        ADM_BSSN::ADM_BSSN_phi
+        ADM_BSSN::ADM_BSSN_metric
+        ADM_BSSN::ADM_BSSN_K
+        ADM_BSSN::ADM_BSSN_curv
+        ADM_BSSN::ADM_BSSN_gamma
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMAnalysis::Ricci
+        ADMConstraints::hamiltonian
+        ADMConstraints::normalized_hamiltonian
+        ADMConstraints::momentum
+"
+
+IOHDF5::checkpoint   = yes
+IO::checkpoint_dir   = $parfile
+IO::checkpoint_every = 5242880
+
+
+
+ActiveThorns = "TimerReport"
+
+TimerReport::out_every = 5242880