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

!DESC "Binary black hole inspiral, uses TATJacobi for elliptic solve, bitant symmetry"

ActiveThorns = "Boundary HyperCartGrid3D Time PUGH PUGHSlab"
ActiveThorns = "PUGHReduce IOASCII IOBasic IOUtil ADMBase StaticConformal TATelliptic TATJacobi IDPuncture" 
# RotatingSymmetry

Cactus::cctk_itlast		= 0
 
Time::dtfac			= 0.25

Grid::type			= "byrange"
Grid::domain			= "bitant"
Grid::rotation_axis		= "y"
Grid::quadrant_direction	= "z"   # default
Grid::xmin                      =  0
Grid::xmax                      =  4
Grid::ymin                      =  0
Grid::ymax			=  4
Grid::zmin                      = -4
Grid::zmax			=  4
Driver::global_nx		= 32  # the "rotation plane"
Driver::global_ny		= 17  # the bitant plane
Driver::global_nz		= 32

ADMBase::initial_data		= "puncture"

IDPuncture::verbose		= "yes"

IDPuncture::bhm1		= 1.0
IDPuncture::bhz1		= 1.5
IDPuncture::bhpx1		= 2.0

IDPuncture::bhm2		= 1.0
IDPuncture::bhz2		= -1.5
IDPuncture::bhpx2		= -2.0
				   
IDPuncture::bound		= "robin"

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

ADMBase::lapse_evolution_method	= "static"
ADMBase::initial_lapse		= "one"

TATJacobi::verbose		= yes

# TATPETsc::options		= "-snes_type tr -snes_atol 1e-6 -ksp_type cg -ksp_atol 1e-6 -get_total_flops"
  
IO::out_dir			= "puncture-quadrant"
IO::out_fileinfo		= "axis labels"

IOBasic::outScalar_every	= 1
IOBasic::outScalar_vars		= "staticconformal::psi admbase::gxx admbase::gyy admbase::gxy admbase::kxx admbase::kyy admbase::kxy idpuncture::u idpuncture::M idpuncture::N"

IOASCII::out1D_every		= 1
IOASCII::out1D_vars		= "staticconformal::psi admbase::gxx admbase::gyy admbase::gxy admbase::kxx admbase::kyy admbase::kxy idpuncture::u idpuncture::M idpuncture::N"