
#define tetradtype 2

#define Cell

#define ghost_width 3



#define GAUGE 0

#define CPBC_ghost_width  (ghost_width)

#define ABV 0

#define EScalar_CC 2

#if 0

define tetradtype
    v:r; u: phi; w: theta
    tetradtype 0
    v^a = (x,y,z)
    orthonormal order: v,u,w
    m = (phi - i theta)/sqrt(2) following Frans, Eq.(8) of  PRD 75, 124018(2007)
    tetradtype 1
    orthonormal order: w,u,v
    m = (theta + i phi)/sqrt(2) following Sperhake, Eq.(3.2) of  PRD 85, 124062(2012)
    tetradtype 2
    v_a = (x,y,z)
    orthonormal order: v,u,w
    m = (phi - i theta)/sqrt(2) following Frans, Eq.(8) of  PRD 75, 124018(2007)

define Cell or Vertex
    Cell center or Vertex center

define ghost_width
    2nd order: 2
    4th order: 3
    6th order: 4
    8th order: 5

define WithShell
    use shell or not

define CPBC
    use constraint preserving boundary condition or not
    only affect Z4c
    CPBC only supports WithShell

define GAUGE
    0: B^i gauge
    1: David puncture gauge
    2: MB B^i gauge
    3: RIT B^i gauge
    4: MB beta gauge (beta gauge not means Eq.(3) of PRD 84, 124006)
    5: RIT beta gauge (beta gauge not means Eq.(3) of PRD 84, 124006)
    6: MGB1 B^i gauge
    7: MGB2 B^i gauge

define CPBC_ghost_width  (ghost_width)
    buffer points for CPBC boundary

define ABV
    0: using BSSN variable for constraint violation and psi4 calculation
    1: using ADM variable for constraint violation and psi4 calculation

define EScalar_CC
Type of Potential and Scalar Distribution in F(R) Scalar-Tensor Theory
    1: Case C of 1112.3928, V=0
    2: shell with   phi(r) = phi0 * a2^2/(1+a2^2), f(R) = R+a2*R^2 induced V
    3: ground state of Schrodinger-Newton system,  f(R) = R+a2*R^2 induced V
    4: a2 = +oo and phi(r) = phi0 * 0.5 * ( tanh((r+r0)/sigma) - tanh((r-r0)/sigma) )
    5: shell with   phi(r) = phi0 * Exp(-(r-r0)**2/sigma), V = 0

#endif

