#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#            SIMULATION STOP            #
#.......................................#
max_step                =  50           # Max number of time steps
steady_state            =  0            # Steady-state solver? 

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#         TIME STEP COMPUTATION         #
#.......................................#
incflo.cfl = 0.45

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#            INPUT AND OUTPUT           #
#.......................................#
amr.plot_int            =   100         # Steps between plot files

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#               PHYSICS                 #
#.......................................#
incflo.gravity          =   0.  0.  0.  # Gravitational force (3D)
incflo.ro_0             =   1.          # Reference density 

incflo.fluid_model      =   "newtonian" # Fluid model (rheology)
incflo.mu               =   0.01        # Dynamic viscosity coefficient
incflo.use_tensor_solve =   false
incflo.advect_momentum  =   true
incflo.constant_density =   false

incflo.advect_tracer    =   true
incflo.ntrac            =   1           # Number of tracers
incflo.mu_s             =   0.001       # Scalar diffusion coefficient

incflo.advection_type        =  "Godunov" 
incflo.diffusion_type        =  1       # Crank-Nicolson
incflo.redistribution_type   =  "StateRedist" 

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#        ADAPTIVE MESH REFINEMENT       #
#.......................................#
amr.n_cell              =   32 32 8     # Grid cells at coarsest AMRlevel
amr.max_level           =   1           # Max AMR level in hierarchy 
amr.max_grid_size_x     =   1024
amr.max_grid_size_y     =   1024
amr.max_grid_size_z     =   1024

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#              GEOMETRY                 #
#.......................................#
geometry.prob_lo        =  -1. -1. -.25 # Lo corner coordinates
geometry.prob_hi        =   1.  1.  .25 # Hi corner coordinates
geometry.is_periodic    =   0   0   1   # Periodicity x y z (0/1)

# Boundary conditions
xlo.type                =   "mixed"
xlo.velocity            =   1.0  0.0  0.0
xlo.pressure            =   0.0
xlo.density             =   2.0
xlo.tracer              =   2.0

xhi.type                =   "mixed"
xhi.velocity            =   1.0  0.0  0.0
xhi.pressure            =   0.0
xhi.density             =   2.0
xhi.tracer              =   2.0

ylo.type                =   "nsw"
yhi.type                =   "nsw"

# Add box
incflo.geometry = "box"
box.internal_flow = false

box.Lo = -0.4  -0.2  -1.1
box.Hi =  1.1   0.2   1.1

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#           INITIAL CONDITIONS          #
#.......................................#
incflo.probtype = 1101

incflo.ic_u             =   0.0         #
incflo.ic_v             =   0.0         #
incflo.ic_w             =   0.0         #
incflo.ic_p             =   0.0         #
incflo.ic_t             =   0.0         #

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#              VERBOSITY                #
#.......................................#
incflo.verbose            =   1         # incflo itself
mac_proj.verbose          =   1         # MAC Projector
nodal_proj.verbose        =   1         # Nodal Projector