#################################################################
# timegeodesic.par
#################################################################

# REQUIRED THORNS

ActiveThorns = "admcoupling LocalReduce admmacros spacemask staticconformal coordgauge time pugh pughreduce pughslab CartGrid3D CoordBase SymBase admbase adm adm_bssn pughinterp timegeodesic ioascii ioutil iobasic boundary"

# GRID

driver::global_nx = 20
driver::global_ny = 20
driver::global_nz = 20

grid::type   = "byspacing"
grid::dxyz   = 0.05
grid::domain = "octant"

time::dtfac = 0.25           

cactus::cctk_initial_time = 0.0
cactus::cctk_itlast       = 50

# METRIC TYPE

#admbase::metric_type      = "Static Conformal"
admbase::metric_type      = "physical"

# INITIAL DATA

#admbase::initial_data = "flat"
admbase::initial_data = "Cartesian Minkowski"

# EVOLUTION

admbase::evolution_method = "adm_bssn"

adm_bssn::method = "stagleap"
adm_bssn::bound  = "flat"

# GAUGE

admbase::lapse_evolution_method       = "geodesic"
#admbase::initial_lapse                = "gaussian"

#admbase::gaussian_amplitude = 0.0
#admbase::gaussian_sigma2    = 0.1

# INTERP

#pughinterp::order = 2

# GEODESICS
timegeodesic::out_xgraph = "yes"

timegeodesic::num_active = "yes"   #  Integrate numerical geodesics?
timegeodesic::num_test   = "yes"   #  Cosistency tests.
timegeodesic::findnum_IJ = "yes"   #  Calculate invariants I and J on geodesics

timegeodesic::exa_active = "no"   #  Integrate analytical geodesics?
timegeodesic::findexa_IJ = "no"   #  Calculate invariants I and J on geodesics

timegeodesic::geopos     = "x00"  #  Initial position of geodesics
timegeodesic::geodesics  = 50     #  Total number of geodesics

# OUTPUT
timegeodesic::out_xgraph  = "yes"
IO::out_dir                = "2timegeodesic"

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

IOBasic::outScalar_every = 1
IOBasic::outScalar_vars  = "admbase::alp admbase::gxx admbase::gyy admbase::gzz admbase::gxy admbase::gxz admbase::gyz"