parameters for hlc_graxi
########################################################################
#                                                                      #
#  id       Standard Initial Data File                                 #
#                                                                      #
#           To construct otherwise identical ID files, run             #
#           "mi id" which will create files idd0, idd1, etc            #
#           which then differ only by the resolution.                  #
#                                                                      #
########################################################################

iter        := 1000
output      := 0-*/20

#
#     Matter Parameters
#
amp1          := 0.35
r01           := 0.0
delta1        := 3.0
   # +1  ---  in-going
   #  0  ---  time-symmetric
   # -1  ---  out-going
signum1       := 0
rho_center1   := 0
z_center1     := 0.0
cyl1          := 1

#FP -- second piece of initial profile
amp1_b          := 0.075
r01_b           := 0.0
delta1_b        := 2.0
   # +1  ---  in-going
   #  0  ---  time-symmetric
   # -1  ---  out-going
signum1_b       := 0
rho_center1_b   := 10.0
z_center1_b     := 10.0
cyl1_b          := 1 

#
#     Geometry Parameters
#

# sigma
ampsi         := 0.00000
r0si          := 8.0
deltasi       := 2.0
signumsi      := 0
rho_centersi  := 0.0
z_centersi    := 0.0
cylsi         := 1

# omega
ampo          := 0.00000
r0o           := 7.0
deltao        := 3.0
signumo       := 0
rho_centero   := 0.0
z_centero     := 0.0
cylo          := 1

#
# Tracing level:
#
# > 0  ----  Canonical CN loop output
# > 1  ----  Canonical MG tracing
# > 2  ----  Output relativized errors in CN loop
# > 3  ----  parameter dumping
# > 4  ----  level output in MG
#
trace_level := 1

#
#     General Parameters
#
MAXITER     := 50
maxstep     := 50
Nrho0       := 64
Nz0         := 128
rhomin      := 0
rhomax      := 20
zmin        := -20
zmax        := 20
ser         := 0
fout        := 1
epsiter     := 1.0e-5
tag         := "ex2_"
lambda      := .3
epsdis      := 0.5

#
# Specify form of Initial Data:
#    IDATA:
#      0     ---- general Gaussians with parameters set below
#      1     ---- specified ansatz, and computed source terms
#      2     ---- load in scalar field profiles
#
IDATA := 0

#
# Restrict model to Flat Space?
#     (if one, then metric restricted to flat metric)
#
flat_space := 0

#
#  force_cont_off == 1   do *not* use continuation methods in ID
#                 != 1   do       use them
#
force_cont_off := 0

#
# Output initial time sdf's? (if non-zero)
#
output_init := 0

#
# Solve by relaxation ? (if equal to one)
# nsweeps -- max number of relaxation sweeps
# nout    -- sweep interval for output (according to trace_level)
#
solve_by_relax := 0
nsweeps        := 5000
nout           := 1000

#
# Multigrid parameters
#
nvcycle     := 3
preswp      := 3
pstswp      := 3

#
# Excision parameters
#

do_excise := 1
exp_out_thresh := -0.02

sph_omegabar := 0
in_file :=   "in2.sdf"
out_file := "out2.sdf"
level:=           2