#------------------------------------------------------------------------------
# cartoon
#------------------------------------------------------------------------------
# This file was automatically generated by from:
#   "../par-files/ivp/schw.par"
#   "../par-files/machines/general.par"
#   "../par-files/grid/cartoon-schw.par"
#   "../par-files/evo/bssn-newrad.par"
#   "../par-files/gauge/1+log.par"
#   "../par-files/gauge/gamma2_b4a2p2.par"
#   "../par-files/output/out_1d.par"
#   "../par-files/misc/nanchecker.par"
#   "../par-files/misc/profile.par"
#------------------------------------------------------------------------------

ActiveThorns = " einstein idanalyticbh time pugh pughreduce pughslab cartgrid3d ioascii ioutil iobasic cartoon2d adm_bssn admconstraints boundary nanchecker "

IOBasic::outScalar_vars = " einstein::metric einstein::curv einstein::lapse einstein::shift einstein::psi adm_bssn::adm_bssn_metric adm_bssn::adm_bssn_curv adm_bssn::adm_bssn_k adm_bssn::adm_bssn_gamma adm_bssn::adm_bssn_phi admconstraints::momy admconstraints::ham "

IOBasic::outInfo_vars = " admconstraints::ham einstein::alp "

IOASCII::out1D_vars = " einstein::metric einstein::curv einstein::lapse einstein::shift einstein::psi adm_bssn::adm_bssn_metric adm_bssn::adm_bssn_curv adm_bssn::adm_bssn_k adm_bssn::adm_bssn_gamma adm_bssn::adm_bssn_phi admconstraints::momy admconstraints::ham "

IOASCII::out2D_vars = ""

#==============================================================================
# schwarzschild initial data
#==============================================================================

einstein::initial_data                  = "bl_bh"
idanalyticbh::bl_M_1                    = 1.0

#==============================================================================
#generic machine outputdir
#==============================================================================

IO::outdir="schw-0050"

#==============================================================================
# Cartoon grid
#==============================================================================

driver::global_nx                       = 51
driver::global_ny                       = 5
driver::global_nz                       = 51

grid::type                              = "byrange"

grid::xmin				= -0.15
grid::xmax				= 4.85

grid::ymin				= -0.2
grid::ymax				= 0.2

grid::zmin				= -4.85
grid::zmax				= 4.85

time::dtfac                             = 0.25

grid::domain                            = "bitant"
grid::bitant_plane			= "xy"

driver::ghost_size_x                    = 2
driver::ghost_size_y                    = 2
driver::ghost_size_z                    = 2

cactus::cctk_itlast                     = 10

grid::avoid_originy 			= "no"

cartoon2d::stencil			= "yes"
cartoon2d::order			= 3
#cartoon2d::verbose			= "yes"

adm_bssn::cartoon			= "yes"
cartoon2d::cartoon_active		= "yes"
admconstraints::cartoon			= "yes"

#==============================================================================
# Standard ADM_BSSN evolution
#==============================================================================

einstein::evolution_system              = "adm_bssn"

adm_bssn::method                        = "icn"
adm_bssn::ICN_itnum                     = 3
adm_bssn::stencil_size                  = 2
einstein::advection                     = "upwind2"
adm_bssn::bound                         = "newrad"
boundary::radpower                      = 3

#==============================================================================
# 1+log
#==============================================================================

einstein::initial_lapse                 = "one"
einstein::slicing                       = "1+log"
adm_bssn::lapsesource                   = "modified"
adm_bssn::harmonic_f                    = 2.0

#==============================================================================
# Gamma2 shift
#==============================================================================

einstein::initial_shift                 = "zero"
einstein::shift                         = "gamma2"

adm_bssn::ShiftGammaCoeff               = 0.75
adm_bssn::BetaDriver                    = 4.0
adm_bssn::ShiftAlpPower                 = 1.0
adm_bssn::ShiftPsiPower                 = 4.0

#==============================================================================
# 1D output
#==============================================================================

IOBasic::outScalar_every                = 1
IOASCII::out1D_every                    = 1
IOBasic::outInfo_every                  = 1

NaNChecker::check_every 	=	1
NaNChecker::check_vars	 	=	"all"
NaNChecker::action_if_found	=	"abort"