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

# REQUIRED THORNS

ActiveThorns = "time localreduce pugh pughreduce pughslab CartGrid3D CoordBase Boundary SymBase ADMBase ADMMacros StaticConformal CoordGauge  adm pughinterp timegeodesic idsimple ioascii ioutil iobasic localinterp"

# 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

# INITIAL DATA

admbase::initial_data = "Cartesian Minkowski"

# EVOLUTION

admbase::evolution_method = "adm"

adm::method = "stagleap"
adm::bound  = "flat"

# GAUGE

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

idsimple::gaussian_amplitude = 0.0
idsimple::gaussian_sigma2    = 0.1

# INTERP

#interp::order = 2

# GEODESICS

timegeodesic::num_active = "yes"  #  Integrate numerical geodesics?
timegeodesic::num_test   = "no"   #  Cosistency tests.
timegeodesic::findnum_IJ = "no"   #  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

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

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