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

-