# $Header$
# Denis Pollney <pollney@aei.mpg.de>

ActiveThorns = "Boundary HyperCartGrid3D Time PUGH PUGHReduce PUGHSlab"
ActiveThorns = "IOASCII IOBasic IOUtil"
ActiveThorns = "ADMBase StaticConformal"
ActiveThorns = "ADMCoupling ADMMacros CoordGauge SpaceMask BSSN_MoL MoL NaNChecker"
ActiveThorns = "TATelliptic TATJacobi IDPuncture"

Cactus::cctk_itlast		= 3
 
Time::dtfac			= 0.25

Grid::type			= "byrange"
Grid::domain			= "full"
Grid::xmin                      = -4
Grid::xmax                      =  4
Grid::ymin                      = -4
Grid::ymax			=  4
Grid::zmin                      = -4
Grid::zmax			=  4
Driver::global_nx		= 35
Driver::global_ny		= 66
Driver::global_nz		= 17

PUGH::ghost_size		= 2

ADMBase::initial_data		= "puncture"

IDPuncture::verbose		= "yes"

IDPuncture::bhm1		= 1.0
IDPuncture::bhy1		= 1.0
IDPuncture::bhpx1		= 1.0
IDPuncture::bhpz1		= 1.0

IDPuncture::bhm2		= 1.0
IDPuncture::bhy2		=-1.0
IDPuncture::bhpx2		=-1.0
IDPuncture::bhpz2		= 1.0

IDPuncture::solve_elliptic	= "no"
	   
IDPuncture::bound		= "robin"

ADMBASE::evolution_method	= "adm_bssn"
mol::ode_method			= "icn"
ADM_BSSN::advection		= "upwind2"
ADM_BSSN::stencil_size		= 2
ADM_BSSN::bound			= "newrad"
Boundary::radpower		= 3

ADMBase::initial_lapse		= "one"
ADMBase::lapse_evolution_method = "1+log"
ADMBase::initial_shift		= "none"

ADMBase::metric_type		= "static conformal"
StaticConformal::conformal_storage = "factor+derivs+2nd derivs"

IO::out_dir			= "puncture-nosolve-full"
IO::out_fileinfo		= "axis labels"

IOBasic::outScalar_every	= 1
IOBasic::outScalar_vars		= "staticconformal::psi admbase::metric admbase::curv admbase::shift"
# idpuncture::u idpuncture::M idpuncture::N"
#IOBasic::outScalar_style	= "gnuplot"

IOASCII::out1D_every		= 1
IOASCII::out1D_vars		= "staticconformal::psi adm_bssn::adm_bs_k admbase::metric admbase::curv idpuncture::u"
#IOASCII::out1D_vars		= "staticconformal::psi admbase::gxx admbase::gyy admbase::gxy admbase::kxx admbase::kyy admbase::kxy idpuncture::u idpuncture::M idpuncture::N"
#IOASCII::out1D_style		= "gnuplot f(x)"

IOBasic::outInfo_every		= 1
IOBasic::outInfo_vars		= "admbase::alp"