[Commits] [svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/bbh/ (Rev. 86)
bcmsma at astro.rit.edu
bcmsma at astro.rit.edu
Tue May 17 12:27:29 CDT 2011
User: bmundim
Date: 2011/05/17 12:27 PM
Modified:
/examples/bbh/
BBHLowRes.par
Log:
Parameter file used for this first bbh run.
Comments:
1) According to Peter and Ian comments, I will have to change the
order in WeylScal4::fd_order and maybe increase the angular
resolution in Multipole.
2) This parameter file runs smoothly with the ET curie release
version. When the ET development version is used, it develops
Nans in the AMR boundary closest to Pi symmetry boundary
around 10M of evolution for the second finest grid level.
File Changes:
Directory: /examples/bbh/
=========================
File [modified]: BBHLowRes.par
Delta lines: +74 -69
===================================================================
--- examples/bbh/BBHLowRes.par 2011-05-11 15:34:34 UTC (rev 85)
+++ examples/bbh/BBHLowRes.par 2011-05-17 17:27:28 UTC (rev 86)
@@ -11,12 +11,12 @@
#Cactus::terminate = "never"
#Cactus::terminate = "time"
#Cactus::cctk_final_time = 1500.0
-Cactus::terminate = "iteration"
+#Cactus::terminate = "iteration"
#Cactus::cctk_itlast = 0
-Cactus::cctk_itlast = 256
-#Cactus::terminate="runtime"
+#Cactus::cctk_itlast = 256
+Cactus::terminate="runtime"
#Cactus::max_runtime=28800 # 46 days
-#Cactus::max_runtime=1410 # 23.5 hours
+Cactus::max_runtime=1410 # 23.5 hours
#Cactus::max_runtime=1380 # 23 hours
#Cactus::max_runtime=690 # 11.5 hours
#Cactus::max_runtime=450 # 7.5 hours
@@ -54,6 +54,8 @@
#------------------------------------------------------------------------------
# Diagnostic parameters:
#------------------------------------------------------------------------------
+AHFinderDirect::verbose_level = "physics details"
+
Carpet::output_timers_every = 1024
Carpet::storage_verbose = "no"
#Carpet::timers_verbose = "no"
@@ -68,11 +70,16 @@
#CarpetLib::print_timestats_every = 1024
CarpetMask::verbose = "yes"
-
+CarpetReduce::verbose = "no"
CarpetRegrid2::verbose = "yes"
+#CarpetRegrid2::veryverbose = "yes"
-#NaNChecker::verbose = "all"
+CarpetTracker::verbose = "yes"
+Dissipation::verbose = "no"
+
+NaNChecker::verbose = "all"
+
TimerReport::out_every = 1024
TimerReport::out_filename = "TimerReport"
#TimerReport::output_all_timers = "yes"
@@ -83,18 +90,16 @@
TwoPunctures::verbose = "yes"
-dissipation::verbose = "no"
-#CarpetTracker::verbose = "yes"
-
-AHFinderDirect::verbose_level = "physics details"
+PunctureTracker::verbose = "yes"
QuasiLocalMeasures::verbose = "yes"
+SphericalSurface::verbose = "yes"
#------------------------------------------------------------------------------
# Utility parameters:
#------------------------------------------------------------------------------
-NaNChecker::check_every = 16
-#NaNChecker::check_every = 1
+#NaNChecker::check_every = 16
+NaNChecker::check_every = 1
NaNChecker::check_vars = "
ML_BSSN::ML_log_confac
ML_BSSN::ML_metric
@@ -215,7 +220,7 @@
CarpetInterp::check_tree_search = "no"
CarpetRegrid2::freeze_unaligned_levels = "yes"
-CarpetRegrid2::snap_to_coarse = "yes"
+#CarpetRegrid2::snap_to_coarse = "yes"
CarpetRegrid2::symmetry_rotating180 = "yes"
# System specific Carpet parameters:
@@ -298,8 +303,8 @@
TwoPunctures::npoints_B = 16
TwoPunctures::npoints_phi = 14
TwoPunctures::Newton_tol = 1.e-10
-#TwoPunctures::grid_setup_method = "evaluation"
-TwoPunctures::grid_setup_method = "Taylor expansion"
+TwoPunctures::grid_setup_method = "evaluation"
+#TwoPunctures::grid_setup_method = "Taylor expansion"
TwoPunctures::par_b = 3.0
TwoPunctures::par_m_plus = 0.47656
@@ -316,6 +321,7 @@
InitBase::initial_data_setup_method = "init_all_levels"
ADMMacros::spatial_order = 4
+SummationByParts::order = 4
#SpaceMask::use_mask = "yes"
@@ -374,7 +380,7 @@
ML_BSSN::harmonicN = 1 # 1+log
ML_BSSN::harmonicF = 2.0 # 1+log
ML_BSSN::ShiftGammaCoeff = 0.75
-ML_BSSN::AlphaDriver = 1.0
+#ML_BSSN::AlphaDriver = 1.0
ML_BSSN::BetaDriver = 1.0
ML_BSSN::LapseAdvectionCoeff = 1.0
ML_BSSN::ShiftAdvectionCoeff = 1.0
@@ -432,8 +438,8 @@
IOASCII::one_file_per_group = "yes"
IOASCII::output_symmetry_points = "no"
-IOASCII::output_ghost_points = "no"
-#IOASCII::out3D_ghosts = "no"
+#IOASCII::output_ghost_points = "no" #hg version only!
+IOASCII::out3D_ghosts = "no"
#IOASCII::out3D_outer_ghosts = "no"
@@ -451,32 +457,34 @@
"
#Set these IOASCII options for initial data only:
-#IOASCII::out1D_every = 512000000
+IOASCII::out1D_every = 512000000
+IOASCII::out1D_d = "no"
+IOASCII::out1D_vars = "
+ CarpetReduce::weight
+ ADMBase::lapse
+ ADMBase::metric
+ ADMBase::curv
+ ADMBase::shift
+ ML_BSSN::ML_Ham
+ ML_BSSN::ML_mom
+ ML_ADMConstraints::ML_Ham
+ ML_ADMConstraints::ML_mom
+"
+
+# To draw the AH circunference in the z=0 plane
+# set these IOASCII options (for evolution only):
+# NOTE: it is not possible yet to output array lines
+# as set below.
+#IOASCII::out1D_every = 32
+#IOASCII::out1D_x = "no"
+#IOASCII::out1D_z = "no"
#IOASCII::out1D_d = "no"
+#IOASCII::out1D_yline_xi = 19
+#IOASCII::out1D_yline_zi = 0
#IOASCII::out1D_vars = "
-# CarpetReduce::weight
-# ADMBase::lapse
-# ADMBase::metric
-# ADMBase::curv
-# ADMBase::shift
-# ML_BSSN::ML_Ham
-# ML_BSSN::ML_mom
-# ML_ADMConstraints::ML_Ham
-# ML_ADMConstraints::ML_mom
+# SphericalSurface::sf_radius
#"
-# To draw the AH circunference in the z=0 plane
-# set these IOASCII options (for evolution only):
-IOASCII::out1D_every = 32
-IOASCII::out1D_x = "no"
-IOASCII::out1D_z = "no"
-IOASCII::out1D_d = "no"
-IOASCII::out1D_yline_xi = 19
-IOASCII::out1D_yline_zi = 0
-IOASCII::out1D_vars = "
- SphericalSurface::sf_radius
-"
-
IOASCII::out2D_every = 512000000
IOASCII::out2D_vars = "
ML_ADMConstraints::ML_Ham
@@ -500,6 +508,7 @@
Grid::coordinates
ADMBase::lapse
ADMBase::metric
+ ADMBase::curv
ML_ADMConstraints::ML_Ham
ML_ADMConstraints::ML_Mom
WeylScal4::Psi4r
@@ -512,13 +521,14 @@
#IOHDF5::out3D_outer_ghosts = "no"
IOHDF5::out_every = 512000000
+#IOHDF5::out_every = 256
IOHDF5::out_vars = "
+ CarpetReduce::weight
ADMBase::lapse
ML_ADMConstraints::ML_Ham
ML_ADMConstraints::ML_Mom
"
# Grid::coordinates
-# CarpetReduce::weight
# ADMBase::metric
# ML_BSSN::ML_Ham
# ML_BSSN::ML_mom
@@ -528,22 +538,25 @@
#------------------------------------------------------------------------------
AHFinderDirect::find_every = 32
-AHFinderDirect::run_at_CCTK_ANALYSIS = "yes"
-AHFinderDirect::run_at_CCTK_POSTSTEP = "no"
+#AHFinderDirect::run_at_CCTK_ANALYSIS = "yes"
+#AHFinderDirect::run_at_CCTK_POSTSTEP = "no"
+AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = "no"
-AHFinderDirect::move_origins = yes
+AHFinderDirect::move_origins = "yes"
+AHFinderDirect::reshape_while_moving = "yes"
+AHFinderDirect::predict_origin_movement = "yes"
# Hermite to order 3 to avoid discontinuities in the metric spatial derivatives:
-AHFinderDirect::geometry_interpolator_name = "Hermite polynomial interpolation"
-AHFinderDirect::geometry_interpolator_pars = "order=3"
-AHFinderDirect::surface_interpolator_name = "Hermite polynomial interpolation"
-AHFinderDirect::surface_interpolator_pars = "order=3"
-#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"
+#AHFinderDirect::geometry_interpolator_name = "Hermite polynomial interpolation"
+#AHFinderDirect::geometry_interpolator_pars = "order=3"
+#AHFinderDirect::surface_interpolator_name = "Hermite polynomial interpolation"
+#AHFinderDirect::surface_interpolator_pars = "order=3"
+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"
-AHFinderDirect::output_h_every = 512
+AHFinderDirect::output_h_every = 0
AHFinderDirect::N_horizons = 3
@@ -610,23 +623,15 @@
#------------------------------------------------------------------------------
# Checkpoint/Recovery:
#------------------------------------------------------------------------------
-CarpetIOHDF5::checkpoint = "yes"
-#CarpetIOHDF5::checkpoint = "no"
-
-IO::checkpoint_ID = "yes" #may save you from calculating the initial data again
-IO::checkpoint_every = 6144
-#IO::checkpoint_every = 6912
-##IO::checkpoint_every = 9750
+IOHDF5::checkpoint = "yes"
+IO::checkpoint_dir = $parfile
+IO::checkpoint_ID = "yes"
+IO::checkpoint_every = 6144
+#IO::checkpoint_every_walltime_hours = 6.0
IO::checkpoint_keep=3
-IO::checkpoint_on_terminate = "yes"
-#IO::checkpoint_dir = $parfile
-#IO::checkpoint_file = "checkpoint"
-#IO::verbose = "full"
+IO::checkpoint_on_terminate = "yes"
-#IO::recover_dir = $parfile
-#IO::recover_file = "checkpoint"
-IO::recover = "autoprobe"
-#IO::recover = "auto"
+IO::recover = "autoprobe"
+IO::recover_dir = $parfile
-
More information about the Commits
mailing list