#################################################################
# BOOSTED NEUTRON STAR
#################################################################

system        =    "mahc adm"    
ADM_first_step  = "predcorr"
boost_x      = 0.3


nx,ny,nz      =    61,61,61
xmin,ymin,zmin = -30.0,-30.0,-30.0
xmax,ymax,zmax =  30.0, 30.0, 30.0

dtfac         =    0.125            #  courant factor

grid          =    "minmax"                       

itfirst       =     0             #  Evolve from iteration          
itlast        =   480             #  to iteration                   
itout1        =   10
itout3        =   10
itinfo        =   10             #  print info (time,maxvals) each

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

model         =   "boost_star"  

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

hydro         =   "mahc"
use_conformal =   "no"

#SLICING
slicing       =   "harmonic"      #  [geodesic,harmonic,1+log]
shift         =   "zero"
# BOUNDARIES
bound = "static"

outdir = "./boost_star_adm"
showmaxvals   = "gxx ham rho"
slicefields   = "gxx hxx alp rho eps velx vely velz press w_lorentz 
                 dens tau sx sy sz betax betay betaz 
                 gxy gxz gyy gyz gzz hxy hxz hyy hyz hzz 
                 ham_adm momx_adm momy_adm momz_adm "
outfields     = "gxx rho eps ham_adm "