#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#            SIMULATION STOP            #
#.......................................#
max_step                =   10          # Max number of time steps

incflo.initial_iterations = 2

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#         TIME STEP COMPUTATION         #
#.......................................#
incflo.fixed_dt         =   -1.0        # Use this constant dt if > 0
incflo.cfl              =   0.45        # CFL factor

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#            INPUT AND OUTPUT           #
#.......................................#
amr.plot_int            =   10          # Steps between plot files
amr.check_int           =   1000        # Steps between checkpoint files
amr.restart             =   ""          # Checkpoint to restart from 

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#               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 =  false        #
incflo.advect_tracer    =  true         #
incflo.ntrac            =  2
incflo.trac_is_conservative =  1  0

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

amr.blocking_factor     = 1

incflo.tag_region = true
incflo.tag_region_lo = 0.  1.  
incflo.tag_region_hi = 0.5 2.


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

# Boundary conditions
xlo.type                =   "mi"
xlo.velocity            =   2. 3. 4.
xlo.tracer              =   1. 0.5
xhi.type                =   "po"
xhi.pressure            =   0.0

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

incflo.advection_type = "BDS"

incflo.geometry = "all_regular"

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#           INITIAL CONDITIONS          #
#.......................................#
incflo.probtype                 = 122
incflo.test_tracer_conservation = true

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
#              VERBOSITY                #
#.......................................#
incflo.verbose          =   2           # incflo_level

amr.plt_ccse_regtest    =  1