#	Runtime parameters for a jet injection problem.

#	Parameters for initial model

# Jet parameters (MH)

jet_vel	= 3.0
p_a	= 1.0
p_j	= 1.0
rho_a   = 1.0
rho_j   = 1.0
rj      = 1.0
mjb     = 25

# --------------

#		Gas ratio of specific heats

gamma           = 1.6666667

#	Computational volume parameters

#		Grid geometry

geometry = "cartesian"

#		Size of computational volume

Nblockx         = 16
Nblocky         = 4
Nblockz         = 1
xmin            = 0.
xmax            = 40.
ymin            = -5.0
ymax            = 5.0
nrefs           = 5
xmin            = 0.
#xmax            = 0.25
zmin            = 0.
zmax            = 0.25
#ymin            = 0.0
#ymax            = 5.0
gconst          =-1.
gdirec          = "y"
igrav           = 1

#		Boundary conditions

xl_boundary_type = "user"
xr_boundary_type = "outflow"

yl_boundary_type = "outflow"
yr_boundary_type = "outflow"


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


plot_var_1      = "dens"
plot_var_2      = "pres"
plot_var_3      = "temp"
plot_var_4      = "velx"
plot_var_5      = "vely"
plot_var_6      = "magx"
plot_var_7      = "magy"
plot_var_8      = "heav"

cfl		= 0.8
lrefine_min     = 1
lrefine_max     = 2
basenm          = "jet_inj_mhd"
nend            = 1000000
tmax            = 50.0
tplot		= 0.1
eint_switch     = 1.e-4

#-------------------------------------------------------------------------------

# Restart

restart                       =  .false.
trstrt                         =  5.0
nrstrt                         =  1000000
#cpnumber			=  0
#ptnumber			=  0
 

#-------------------------------------------------------------------------------

# Refinement criteria

refine_cutoff   = 0.75
derefine_cutoff = 0.25

refine_var_1    = "dens"


# Runtime parameters for the MHD R-T instability problem


lrefine_min     = 1
cpnumber        = 0
ptnumber        = 0
restart         = .FALSE.
tplot           = 0.05
trstrt          = 1.
nend            = 1000000
#tmax            = 4.0
gamma           = 1.66666666666

#xl_boundary_type      = "periodic"
#xr_boundary_type      = "periodic"
zl_boundary_type      = "periodic"
zr_boundary_type      = "periodic"
yl_boundary_type      = "outflow"
yr_boundary_type      = "outflow"


refine_cutoff   = 0.7
derefine_cutoff = 0.1

refine_var_1    = "dens"
refine_var_2    = "heavy tracer"


rho_heavy       = 2.00
rho_light       = 1.00
Bx0             = 0.01
By0             = 0.00
Bz0             = 0.00

epsilon         = 0.025
cfl             = 0.8