#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#            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.  1.  # Gravitational force (3D)
incflo.ro_0             =   1.          # Reference density 

incflo.fluid_model      =   "newtonian" # Fluid model (rheology)
incflo.mu               =   0.001       # Dynamic viscosity coefficient
incflo.use_tensor_solve =   false
incflo.advect_momentum  =   false
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              =   8 32 32     # 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        = -.25  -1. -1. # Lo corner coordinates
geometry.prob_hi        =  .25   1.  1. # Hi corner coordinates
geometry.is_periodic    =    1   0   0  # Periodicity x y z (0/1)

# Boundary conditions
ylo.type                =   "nsw"

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

zlo.type                =   "nsw"
zhi.type                =   "nsw"

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

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

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

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