#	Runtime parameters for the gas collapse problem.

print_tstep_loc = .true.

reference_density = 1e6
#jeans_ref_sqr   = 160000
#jeans_dref_sqr  = 64e30
 
eint_switch = 1e-4
#	Parameters for initial model

rho_0		= 9947
R_init		= 6E8
T_ambient	= 1000.
xctr		= 1e5
yctr		= 1e5
zctr		= 1e5

#		Gas ratio of specific heats

gamma           = 1.6666667

#	Computational volume parameters

#		Grid dimensionality and geometry

igeomx		= 0
igeomy		= 0
igeomz		= 0

#		Size of computational volume

xmin		= 0.
xmax		= 2E5
ymin		= 0.
ymax		= 2E5
zmin		= 0.
zmax		= 2E5
Nblockx		= 1
Nblocky		= 1
Nblockz		= 1

#		Boundary conditions

xl_boundary_type = "reflecting"
xr_boundary_type = "reflecting"
yl_boundary_type = "reflecting"
yr_boundary_type = "reflecting"
zl_boundary_type = "reflecting"
zr_boundary_type = "reflecting"
grav_boundary_type = "periodic"

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

cfl		= 0.8
lrefine_min     = 2
lrefine_max     = 5
refine_cutoff   = .8
derefine_cutoff = 0.2
basenm          = "evrard"
restart         = .false.
trstrt          = 2621
nrstrt		= 2500
tplot           = .436
nend            = 9999
tmax            = 7
refine_var_1	= "dens"
#refine_var_2	= "pres"
mpole_lmax	= 5
igrav           = 1

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

run_comment     = "evrard probled"
log_file        = "evrard.log"

smlrho = 1.E-10
smallp = 1.E-10
smalle = 1.E-10
smallt = 1.E-12
dtini  = 1.E-4
dtmin  = 1.E-40
#dtmax  = .002621
nriem  = 100
cvisc  = 0.

#ppm_modifystates = .true.

conserved_var = .true.
monotone      = .true.


mgrid_max_residual_norm = 1e-5
mgrid_solve_max_iter = 10000
#mgrid_nsmooth = 20
mgrid_max_vcycles = 1000