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        := 4000
output      := 0-*/10

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

#FP -- second piece of initial profile
amp1_b          := 0.0
r01_b           := 0.0
delta1_b        := 3.0
   # +1  ---  in-going
   #  0  ---  time-symmetric
   # -1  ---  out-going
signum1_b       := 0
rho_center1_b   := 8.0
z_center1_b     := 0.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     := 500
Nrho0       := 128
Nz0         := 256
rhomin      := 0
rhomax      := 40
zmin        := -40
zmax        := 40
ser         := 0
fout        := 1
epsiter     := 1.0e-3
ic_epsiter  := 1.0e-6
tag         := "cyl_0_005_a025_ex_more_"
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     := 1
preswp      := 6
pstswp      := 6

#
# Excision parameters
#

do_excise := 1
exp_out_thresh := -0.0
ex_buffer := 2.00
ex_method := 2
ex_method_am := 2
ex_mg_buffer := 3
ex_test := 0
ex_in_buffer := 2.00
ex_start_t := 0
repop_extrap := 2
ex_fs_psi := 3
ex_fs_bar := 0
ex_fs_phi := 1
can_reduce_ex_zone := 2
ex_alpha_dbc := -1
ex_ri_psi := 0
ex_bz_dt := 0
ex_brho_dt := 0
ex_bz_max_s := 1
ex_brho_max_s := 1
ex_bz_s_am := 1
ex_brho_s_am := 1
ex_ri_iters_pcs := 2
ex_ri_csteps := 10
ex_ri_byrelax := 0

sph_omegabar := 0
extrapolate := 2

dbc_alpha := 1
dbc_betarho := 1
dbc_betaz := 1
dbc_psi := 0
ic_dbc_alpha := 1
ic_dbc_betarho := 1
ic_dbc_betaz := 1
ic_dbc_psi := 0
dbc_hh := 0
dbc_on_all4 := 0
dbc_on_n_finest := 0
b_sweep_x := 1
sym_diff := 1

# The following set of alternative constraint solver
# parameters are (currently) only used in constraints.inc,
# i.e. the initial conditions are still MG with relaxation
# on coarsest.

#bicg instead of MG on finest level
use_bicg := 0
bicg_max_iter := 5000
bicg_eps_iter := 1.0e-5
bicg_newton_iters := 1

use_evo_for_psi := 1

compute_smallest_eig := 0

apph_search_dt := 0.5
apph_range := 1
apph_mxstep := 10
apph_find_z := 0
apph_find_a := 1
apph_num_c  := 0
apph_min_l := 0
apph_apph_delta_A := 0.005
apph_apph_delta_f_a := 0.1

MG_w := 1.0
MG_w_scale := 0.5
MG_w_rescale := 1.025
MG_w_min := 0.05