#################################################################
# twostar.par
#################################################################

system        =    "mahc adm"     
ADM_first_step = "predcorr"

nx,ny,nz      =    131,131,131       #  grid size                      

dx,dy,dz      =    0.6,0.6,0.6    #  resolution
dtfac         =    0.125            #  courant factor

grid          =    "full"                       

itfirst       =     0             #  Evolve from iteration          
itlast        =    1800              #  to iteration                   
itout1         =    5
itout3         = 30
itinfo        =    1             #  print info (time,maxvals) each

method        =   "leapfrog"          
hydromethod   =   "fluxsplit_fluxlim"    
hydro_type    =   "hyperbolic"         
stencil_size  =    2

model         =   "twostar"  

# STAR LOCATIONS
star1x	=      10
star1y	=	0.
star1z	=       0. 
star2x	=     -10.
star2y	=	0.
star2z	=       0.

# MATTER PARAMETERS
eos_gamma     = 1.6666666666666
tov_eos_gamma = 1.6666666666666
eos_k         = 4.349
centden       = 8.1e-4
densmin       = 0.0

hydro         =   "mahc"
use_conformal =   "no"

do_tsivp      =   "yes"
tsivp_lin_tolerance     =       1.e-10
tsivp_nlin_tolerance    =       1.e-10
tsivp_solver            =       "bam"
elliptic_verbose	=       "yes"
dens_conf_pow 	= 8

#SLICING
slicing       =   "maximal full"  
admcons_persist = "yes"
maximal_solver  = "bam"
maximal_every = 1
maximal_thresh = 1.0e-10

# BOUNDARIES
bound = "static bamrobin"


outdir         = "./twostar_bam"
onefileperproc = "no"
showmaxvals   = "rho ham_adm"
outfields     = "rho trK"
slicefields   = " alp rho eps velx vely velz press w_lorentz 
                  momx_adm momy_adm momz_adm vx vy vz trK 
                  hxx hxy hxz hyy hyz hzz gxx gxy gxz gyy gyz gzz 
                  dens tau sx sy sz ham_adm tsivp_psi"