#_______________________________________________________________________
# Solver settings

#amrex.fpe_trap_invalid = 1

#! Fluid solver
#!-----------------------------------------------------------//

mfix.cfl = 0.45
mfix.overstep_end_time = 1

mfix.initial_iterations = 3
mfix.max_step = 1

mfix.gravity = 0.0  0.0  0.0

mfix.verbose = 1
mac_proj.verbose   = 3
nodal_proj.verbose = 3

amr.blocking_factor = 1

#_______________________________________________________________________
# Geometry / grids / tiles

# Maximum level in hierarchy (for now must be 0, i.e., one level in total)
amr.max_level         =  1
amr.grid_eff          =  0.75
amr.n_error_buf       =  2

geometry.coord_sys    =   0                  # 0: Cartesian
geometry.is_periodic  =   1    1    1        # Periodicity x y z (0/1)
geometry.prob_lo      =   0.   0.   0.       # Lo corner coordinates
geometry.prob_hi      =   1.   1.   0.25     # Hi corner coordinates

# Number of grid cells in each direction at the coarsest level
amr.n_cell =  32  32  8

#! Grids
#!-----------------------------------------------------------//
# Maximum allowable size of each fluid subdomain in the problem domain;

# Fluid
amr.max_grid_size_x = 64
amr.max_grid_size_y = 64
amr.max_grid_size_z = 64

#! Tiles
#!-----------------------------------------------------------//

# Fluid: Maximum tile size within each grid
fabarray.mfiter_tile_size = 1024 1024 1024

#_______________________________________________________________________
# IO / Checkpointing

mfix.par_ascii_int  = -1
mfix.par_ascii_file ="vis"

mfix.plot_int       = 1

mfix.check_int      = -1
mfix.check_file     ="chk"

mfix.plt_regtest  = 1

#_______________________________________________________________________
# Fluid model settings
#
fluid.solve = fluid

fluid.viscosity.molecular.model          = constant
fluid.viscosity.molecular.constant = 0.0002

#_______________________________________________________________________
# DEM model settings
#
dem.solve = None

#_______________________________________________________________________
# EB geometry

mfix.geometry = "cylinder"

cylinder.internal_flow = false

cylinder.internal_flow  =   false
cylinder.radius         =   0.2
cylinder.direction      =   2
cylinder.center         =   .5  .5  .5

#_______________________________________________________________________
# Regions for defining ICs and BCs
#
regions = full-domain

regions.full-domain.shape = box
regions.full-domain.box.lo = 0.   0.   0.
regions.full-domain.box.hi = 1.   1.   0.25

#_______________________________________________________________________
# Initial Conditions
#
ic.regions = full-domain

ic.full-domain.fluid.volfrac   = 1.0
ic.full-domain.fluid.velocity  = -0.70710678  0.70710678  0.0
ic.full-domain.fluid.density   = 1.0

mfix.stop_for_unused_inputs = true

mfix.redistribute_before_nodal_proj = false
mfix.correction_small_volfrac = 0.
mfix.redistribute_after_initial_nodal_proj = false

# Keys added by fixinputs
mfix.enable.fluid = true