######################################################################
#!DESC "Initial data for axisymmetric distored black hole"
######################################################################

ActiveThorns = "coordgauge admanalysis admbase admmacros staticconformal time boundary iobasic ADMconstraints pugh pughreduce pughslab pughinterp cartgrid3d ADM DistortedBHIVP ioascii ioutil"

# GENERAL 

admbase::evolution_method = "ADM"     

driver::global_nx =  32
driver::global_ny =  32
driver::global_nz =  32

grid::type = "byspacing"
grid::dxyz = 0.4
grid::domain = "full"

time::dtfac = 0.25   

cactus::cctk_initial_time = 0.   
cactus::cctk_itlast       = 1

adm::method           = "stagleap" 
adm::bound            = "static"
admconstraints::bound = "static"

# MODEL 

admbase::initial_data    = "distortedbh"
distortedbhivp::amp      = -0.1
distortedbhivp::eta0     = 0.0
distortedbhivp::c        = 0.5
distortedbhivp::sigma    = 1.0
distortedbhivp::etamax   = 6.0
distortedbhivp::n        = 4
distortedbhivp::neta     = 102
distortedbhivp::ntheta   = 27
distortedbhivp::nphi     = 5
distortedbhivp::interpolation_order = 1


# GAUGE

admbase::lapse_evolution_method    = "geodesic"
admbase::metric_type               = "static conformal"
ataticconformal::conformal_storage = "factor+derivs+2nd derivs"

# OUTPUT ########################################################
IO::outdir               = "2test_dbh"
IOASCII::out1d_every     =  1
IOBasic::outScalar_vars  = "admbase::metric admbase::curv admconstraints::ham"
IOASCII::out1D_vars      = "admbase::metric admbase::curv admbase::lapse admconstraints::ham admanalysis::grr staticconformal::psi staticconformal::psix staticconformal::psiy staticconformal::psiz staticconformal::psixx staticconformal::psixy staticconformal::psixz staticconformal::psiyy staticconformal::psiyz staticconformal::psizz"
IOBasic::outInfo_every   = 1
IOBasic::outInfo_vars    = "admconstraints::ham admbase::gxx"

##################################################################