#	Runtime parameters for the Sod shock-tube problem.

#	Parameters for initial model

#		Density, pressure, and velocity on either side of interface

sim_rhoLeft	= 0.00001
sim_rhoRight	= 0.00001
sim_pionLeft    = 0.1
sim_pionRight   = 0.1
sim_peleLeft	= 1000000.
sim_peleRight	= 0.1
#sim_pLeft		= 1.
#sim_pRight		= 0.1
sim_pradLeft	= 0.0
sim_pradRight	= 0.0
sim_uLeft		= 0.
sim_uRight		= 0.
sim_sigma		= 0.00005

#		Angle and position of interface relative to x and y axes

# These settings are better for testing 3-d problem, they are little
# less symmetric.

sim_xangle		= 0.
sim_yangle		= 90.
sim_posn		= 0.0002


# To get a simpler test problem which shows the 1-D Profile
# use the following values.
#sim_xangle		= 0
#sim_yangle		= 90.
#sim_posn		= 0.5

#		Gas ratio of specific heats

#gamma           = 1.4

#	Computational volume parameters

#		Grid dimensionality and geometry

geometry = "cylindrical"

#		Size of computational volume

xmin		= 0.
xmax		= 0.0005
ymin		= 0.
ymax		= 0.0005

#		Boundary conditions

xl_boundary_type = "axisymmetric"
xr_boundary_type = "outflow"

yl_boundary_type = "outflow"
yr_boundary_type = "outflow"


#	Simulation (grid, time, I/O) parameters

cfl		= 0.008
basenm          = "HofiCh_"
restart         = .false.

#	checkpoint file output parameters
checkpointFileIntervalTime = 1.e-11
checkpointFileIntervalStep = 0
checkpointFileNumber = 0

#	plotfile output parameters
plotfileIntervalTime = 0.
plotfileIntervalStep = 0
plotfileNumber = 0
dtmax = 1.e-12
nend            = 100
tmax            = 0.01

run_comment     = "HOFI problem, radial direction"
log_file        = "hofi.log"
eintSwitch     = 1.e-4

eos_singleSpeciesA     = 1.
eos_singleSpeciesZ     = 1.

#gamma                  = 1.66666667
gammaEle               = 1.266666667
gammaIon               = 1.66666667

plot_var_1 = "dens"
plot_var_2 = "pres"
plot_var_3 = "temp"

#    AMR refinement parameters  
lrefine_max     = 6
refine_var_1    = "dens"

# These parameters are used only for the uniform grid

meshCopyCount=1
#iGridSize = 64   #defined as nxb * iprocs
#jGridSize = 64
#kGridSize = 1
iProcs = 1	#number or procs in the i direction	
jProcs = 1
kProcs = 1

# When using UG, iProcs, jProcs and kProcs must be specified.
# These are the processors along each of the dimensions
#FIXEDBLOCKSIZE mode ::
# When using fixed blocksize, iGridSize etc are redundant in
# runtime parameters. These quantities are calculated as 
# iGridSize = NXB*iprocs
# jGridSize = NYB*jprocs
# kGridSize = NZB*kprocs
#NONFIXEDBLOCKSIZE mode ::
# iGridSize etc must be specified. They constitute the global
# number of grid points in the physical domain without taking 
# the guard cell into account. The local blocksize is calculated
# as iGridSize/iprocs  etc.
 
#useDiffuse=.TRUE.
useHeatexchange = .TRUE.
#useConductivity = .TRUE.
#useViscosity = .TRUE.
diff_useRadDiffusion   = .FALSE.
#diff_useEleCond        = .TRUE.
#useHydro               = .TRUE.

## -------------------------------------------------------------##
##  SWITCHES SPECIFIC TO THE UNSPLIT HYDRO SOLVER               ##
#	I. INTERPOLATION SCHEME:
order		= 2      # Interpolation order (first/second/third/fifth order)
slopeLimiter    = "mc"   # Slope limiters (minmod, mc, vanLeer, hybrid, limited)
LimitedSlopeBeta= 1.     # Slope parameter for the "limited" slope by Toro
charLimiting	= .true. # Characteristic limiting vs. Primitive limiting

use_avisc	= .false. # use artificial viscosity (originally for PPM)
#cvisc		= 0.1     # coefficient for artificial viscosity
use_flattening	= .false. # use flattening (dissipative) (originally for PPM)
use_steepening	= .false. # use contact steepening (originally for PPM)
use_upwindTVD	= .false. # use upwind biased TVD slope for PPM (need nguard=6)

#	II. RIEMANN SOLVERS:
RiemannSolver	= "Roe"       # Roe, HLL, HLLC, LLF, Marquina
entropy         = .false.     # Entropy fix for the Roe solver

hy_fPresInMomFlux = 1.0       # Part of pressure gradient term in momentum equation to include in momentum fluxes

#	III. STRONG SHOCK HANDELING SCHEME:
shockDetect	= .false.     # Shock Detect for numerical stability
## -------------------------------------------------------------##

## ---------------------------------------------------------------##
##  SWITCHES SPECIFIC TO THE SUPER-TIME-STEPPING (STS) ALGORITHM  ##
##  NOTE: For details on using STS runtime parameters, please     ##
##        refer to user's guide (Driver chapter).                 ##
useSTS                  = .false.
nstepTotalSTS           = 5
nuSTS                   = 0.2
## ---------------------------------------------------------------##