amrex.fpe_trap_invalid = 1

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

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#         TIME STEP COMPUTATION         #
#.......................................#
incflo.fixed_dt             =   1.e-5       # Use this constant dt if > 0
incflo.do_initial_proj      =   1
incflo.initial_iterations   =   0

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

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

incflo.fluid_model      =   "newtonian" # Fluid model (rheology)
incflo.mu               =   0.001       # Dynamic viscosity coefficient
incflo.constant_density =   true        #

incflo.advect_tracer = true
incflo.mu_s = 0.1

incflo.advection_type        =  "Godunov" 
incflo.redistribution_type   =  "StateRedist" 


#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#        ADAPTIVE MESH REFINEMENT       #
#.......................................#
amr.n_cell              =   16 34       # Grid cells at coarsest AMRlevel
amr.max_level           =   0           # Max AMR level in hierarchy 
amr.max_grid_size_x     =   1024
amr.max_grid_size_y     =   1024
amr.max_grid_size_z     =   1024
amr.blocking_factor     =   2

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

# Boundary conditions
ylo.type                =   "nsw"
ylo.velocity            =   0.0  0.0  0.0
yhi.type                =   "po"
yhi.pressure            =   0.0

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

box.Lo =  -0.20  0.21  
box.Hi =   1.80  3.60 

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#           INITIAL CONDITIONS          #
#.......................................#
incflo.ic_u             =   0.0         #
incflo.ic_v             =   0.0         #
incflo.ic_p             =   0.0         #

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


#.......................................#
#               EB FLOW                 #
#.......................................#
eb_flow.vel_mag        =  0.1
eb_flow.tracer         =  1.