#==============================================================================
# Cartoon2D, using a bitant grid
#==============================================================================

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

#------------------------------------------------------------------------------
# Run parameters for Cartoon
#------------------------------------------------------------------------------

# Cartoon
cartoon2d::cartoon_active		= "yes"
cartoon2d::verbose              	= "yes"
cartoon2d::stencil              	= "yes"
cartoon2d::order                	= 3
cartoon2d::allow_grid_resize		= "yes"

# This is for 2 processors, you have to
# change manually by processor numbers.

#driver::processor_topology       	= "manual"
#driver::processor_topology_3d_x  	= 1
#driver::processor_topology_3d_y  	= 1
#driver::processor_topology_3d_z  	= 2

# Grid

driver::global_nx                       = 16
driver::global_ny                       = 3
driver::global_nz                       = 32

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

grid::avoid_originy                     = "no"

grid::type                              = "byspacing"
grid::domain				= "bitant"
grid::bitant_plane			= "xy"
grid::dxyz                              = 0.2

time::dtfac                             = 0.25

cactus::cctk_itlast                     = 1

# MODEL 

einstein::initial_data                  = "schwarzschild"

#--- Evolution
einstein::evolution_system 		= "adm_bssn"
adm_bssn::method 			= "leapfrog"
adm_bssn::bound  			= "flat"
adm_bssn::cartoon                	= "yes"
admconstraints::cartoon			= "yes"

#--- gauge
einstein::slicing                       = "1+log"

#--- AHFinder.

ahfinder::ahf_active   			= "yes"
ahfinder::ahf_cartoon			= "yes"
ahfinder::ahf_gaussout 			= "yes"

ahfinder::ahf_phi    			= "no"
ahfinder::ahf_octant 			= "high"

ahfinder::ahf_lmax   			= 0

# Interpolation.

ahfinder::interpolation_order 		= 1

#------------------------------------------------------------------------------
# Output
#------------------------------------------------------------------------------
IO::outdir                              = "test_cartoon_1"
IO::out_fileinfo   			= "none"

IOBasic::outScalar_every                = 1
IOBasic::outScalar_vars                 = "admconstraints::hamiltonian admconstraints::momentum einstein::alp einstein::grr einstein::metric einstein::trace_of_K einstein::betax einstein::betay einstein::betaz adm_bssn::adm_bs_gxx adm_bssn::adm_bs_phi adm_bssn::ADM_BSSN_K adm_bssn::ADM_BSSN_gamma"

IOASCII::out1D_every                    = 1
IOASCII::out1D_vars                     = "admconstraints::hamiltonian admconstraints::momentum einstein::alp einstein::grr einstein::metric einstein::trace_of_K einstein::betax einstein::betay einstein::betaz adm_bssn::adm_bs_gxx adm_bssn::adm_bs_phi adm_bssn::ADM_BSSN_K adm_bssn::ADM_BSSN_gamma"

IOBasic::outInfo_every                  = 1
IOBasic::outInfo_vars                   = "einstein::alp adm_bssn::adm_bs_dtalp"
#==============================================================================