#	Runtime parameters for the Bondi-Hoyle accretion problem
#       with thermal conduction (adapted from windtunnel and 
#       grav_coll problems).
#       FLY units adopted. We choose a box of 2 Mpc, putting at 
#       (x_cb, y_cb, z_cb) a TIS halo, described by a mass (mass_cb) in
#       solar mass units (note: given the mass, a TIS halo is fully determined).
#       Ambient gas temperature is in Kelvin.
#       Ambient pressure, density and inflow velocity are in FLY units.
#       A constant inflow of gas at the left y boundary (as in the wind 
#       tunnel problem) allows one to put oneself in the reference
#       frame in which the halo is at rest.
#
#	Parameters for initial model
#

r_ind_c         = 1.05
h_small         = 0.65
rho_ambient	= 4.7236e-1
T_ambient	= 2.e7
# Note: p_ambient = gascon * rho_ambient * T_ambient
# gascon in FLY units: gascon = 1.7278252e-3
p_ambient	= 3.4658e0
x_cb		= 1.2
y_cb		= 1.
z_cb		= 0.
i_fly_units     = 1

vt_x            = 7.0123
vt_y            = 0.
vt_z            = 0.
#		Gas ratio of specific heats

gamma           = 1.6666667

#	Computational volume parameters

#		Grid geometry

geometry	= "cartesian"
igeomx		= 0
igeomy		= 0
igeomz		= 0

#		Size of computational volume (in Mpc*h^-1)

Nblockx		= 12
Nblocky		= 8
xmin		= 0.
xmax		= 3.
ymin		= 0.
ymax		= 2.

#		Boundary conditions

xl_boundary_type = "user"
xr_boundary_type = "outflow"

yl_boundary_type = "outflow"
yr_boundary_type = "outflow"


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

cfl		= 0.8
lrefine_min	= 2
lrefine_max     = 4
basenm          = "/gpfs/temp/asct001/FLASH/bh_cnd_rcl/m4e11/v2.0/t2e7/bh_"
restart         = .true.
cpnumber	= 0029
ptnumber	= 28
trstrt          = 1.e-2
tplot		= 1.e-2
tmax            = 0.8
dtmax           = 1.e-3
nend            = 30000
eint_switch     = 1.e-8
igrav           = 1
refine_filter	= 1.e-6
# This is the correction factor of the Spitzer thermal conductivity -
# accounts for the redduction/enhancement due to magn. turbulence,
# anomalous diffusion, etc.
c_spitzer	= 1.e-1
diffuse_therm	= .true.
icool		= 1

run_comment     = "2-d Bondi-Hoyle accretion with therm. cond."
log_file        = "/gpfs/temp/asct001/FLASH/bh_cnd_rcl/m4e11/v2.0/t2e7/bh.log"

# Specific to TIS profile

mass_cb		= 4.e11