#==============================================================================
# Cartoon2D, using a full grid, and BSSN_MoL
#==============================================================================

ActiveThorns = "cartoon2d idanalyticbh admbase admcoupling admmacros coordgauge spacemask staticconformal admanalysis bssn_mol mol nanchecker admconstraints pughinterp boundary time pugh pughreduce pughslab cartgrid3d ioascii ioutil iobasic localinterp localreduce coordbase initbase symbase"
## ahfinder

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

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

# 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                       = 19
driver::global_ny                       = 5
driver::global_nz                       = 34

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

grid::avoid_originy                     = "no"

grid::type                              = "byrange"
grid::domain				= "full"
grid::xmin                              = -0.3
grid::xmax                              = +3.3
grid::ymin                              = -0.4
grid::ymax                              = +0.4
grid::zmin                              = -3.3
grid::zmax                              = +3.3

time::dtfac                             = 0.25

cactus::cctk_itlast                     = 1

# MODEL 

admbase::initial_data                   = "schwarzschild"
admbase::metric_type			= "static conformal"

#--- Evolution
admbase::evolution_method 		= "adm_bssn"
methodoflines::ode_method 		= "icn"
adm_bssn::bound  			= "flat"
adm_bssn::cartoon                	= "yes"
admconstraints::cartoon			= "yes"

#--- gauge
admbase::lapse_evolution_method         = "1+log"

## #--- AHFinder.
## 
## ahfinder::ahf_active   			= "yes"
## ahfinder::ahf_cartoon			= "yes"
## 
## ahfinder::ahf_phi    			= "no"
## ahfinder::ahf_octant 			= "high"
## 
## ahfinder::ahf_lmax   			= 0
## 
## # Interpolation.
## 
## ahfinder::interpolation_order 		= 1

#------------------------------------------------------------------------------
# Output
#------------------------------------------------------------------------------
IO::out_dir                             = "test_cartoon_2"
IO::out_fileinfo   			= "none"

IOBasic::outScalar_every                = 1
IOBasic::outScalar_vars                 = "admconstraints::hamiltonian admconstraints::momentum admbase::alp admanalysis::grr admbase::metric admanalysis::trK 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 admbase::alp admanalysis::grr admbase::metric admanalysis::trK admbase::betax admbase::betay admbase::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                   = "admbase::alp adm_bssn::adm_bs_dtalp"
#==============================================================================