# Parameters for Rosetta MHD reconnection problem.

# Size of computational volume

xmin		= 0.
xmax		= 12.5
ymin		= 0.
ymax		= 12.5

# Boundary conditions 

#xl_boundary_type      = "outflow"
xl_boundary_type      = "user"
xr_boundary_type      = "outflow"
#yl_boundary_type      = "outflow"
yl_boundary_type      = "user"
yr_boundary_type      = "outflow"

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

run_comment     = "Rosetta MHD reconnection problem"
log_file        = "rosetta.log"
basenm          = "rosetta_"

cpnumber        = 0
ptnumber        = 0
restart         = .FALSE.

nrstrt		= -1
trstrt          = 0.5
tplot           = 0.1
nend            = 1000000000
tmax            = 50.0
cfl		= 0.4

lrefine_max     = 9
lrefine_min     = 3

refine_var_1    = "magx" 
refine_var_1    = "magy"
refine_var_2    = "magz"
refine_var_3	= "divb"

refine_cutoff   = 0.5
derefine_cutoff = 0.1

plot_var_1      = "dens"
plot_var_2      = "velx"
plot_var_3      = "vely"
plot_var_4      = "velz"
plot_var_5      = "pres"
plot_var_6      = "magx"
plot_var_7      = "magy"
plot_var_8      = "magz"
plot_var_9      = "divb"
plot_var_10      = "temp"

gamma           = 1.6666666667

mvisc           = 5.960464478e-4
diff_visc_nu    = 5.960464478e-4

resistive_mhd   = .TRUE.
killdivb        = .FALSE.
eint_switch     = 1.e-4