# NOTE: all the values here must be set unless specifically
#       marked as optional


#**********************************
#ALGORITHM CONTROLS (defaults shown; these are both optional)
#**********************************

proj.proj_2 = 1
proj.filter_factor = .1
ns.init_iter=1

#**********************************
#RUN CONTROLS  
#**********************************

# MUST SET ONE OR BOTH OF THESE.  
max_step =         1  # maximum (integer) number of time steps
#stop_time =       2.0 # maximum (real) time (in same units as dt)

amr.n_cell =      64 64 # indices of the domain at level 0

# (this one is optional)
#amr.restart =     # set if starting from a restart file, else comment out

# This will default to file "probin" if not set
amr.probin_file = probin_4

#**********************************
#ADAPTIVITY CONTROLS
#**********************************

amr.max_level      = 2 # maximum number of levels of refinement
amr.regrid_int     = 2 # how often to regrid if max_level > 0

# (this one is optional)
amr.regrid_file    = fixed_grids_6  # file name for specifying fixed grids

# MUST SET ONLY ONE OF THE TWO BELOW (COMMENT OUT THE OTHER)
amr.ref_ratio        = 2 2 # integer refinement ratio 
#amr.ref_ratio_vect         # IntVect refinement ratio 

#**********************************
#VERBOSITY CONTROLS (these are all optional with defaults shown)
#**********************************

ns.v                    = 1
godunov.v               = 0
mac.v                   = 0
proj.v                  = 1
diffuse.v               = 0
amr.v                   = 1

mg.v                    = 0
cg.v                    = 0
proj.Pcode              = 1

#**********************************
#OUTPUT CONTROLS
#**********************************

#THESE FILES WILL BE CREATED ONLY IF YOU SPECIFY THE NAME
#amr.run_log          = runlog  # name of the run log file 
#amr.grid_log          = grdlog  # name of the grid log file
#amr.data_log         = datalog # name of the data log file

#THESE FILES WILL BE CREATED WITH DEFAULT PREFIX IF YOU DONT SPECIFY OTHERWISE
#amr.check_file       # prefix of the checkpoint files (default = "chk")
#amr.plot_file        # prefix of the plot files       (default = "plt")

# (this one is optional - sums will only be done if this is set to positive)
#ns.sum_interval      = - 1# interval at which to sum the specified quantities

# MUST SET ONLY ONE OF THE TWO BELOW (COMMENT OUT THE OTHER)
#amr.check_per        # (real) period between checkpoint files (in same units as dt)
amr.check_int        = 20 # (integer) number of timesteps between checkpoint files
amr.check_int        = 2  # (integer) number of timesteps between checkpoint files

# MUST SET ONLY ONE OF THE TWO BELOW (COMMENT OUT THE OTHER)
#amr.plot_per         # (real) period between plot files (in same units as dt)
amr.plot_int         = 2 # (integer) number of timesteps between plot files

#amr.plot_vars        = x_velocity y_velocity density tracer
amr.derive_plot_vars = avg_pressure

#**********************************
#TIME STEP CONTROLS
#**********************************

ns.cfl                  = 0.9  # CFL number used to set dt
#ns.fixed_dt             = 0.9  # Time step

# (this one is optional)
#ns.init_shrink          = 1.0  # factor which multiplies the very first time step

#**********************************
#PHYSICS
#**********************************

ns.vel_visc_coef        = 0.01 
ns.temp_cond_coef       = 0.0 
ns.scal_diff_coefs      = 0.01
ns.variable_vel_visc    = 0
ns.variable_scal_diff   = 0

# (this one is optional)
#ns.gravity              = 0.0  # forcing term is rho * abs(gravity) "down"

#**********************************
#GEOMETRY
#**********************************

#
# NOTE:  The domain size must be a multiple of 2 for the BCs to work
#        properly for the viscous benchmark.
#
geometry.coord_sys   =  0 # 0 for x-y (x-y-z), 1 for r-z
geometry.prob_lo     =  0. 0. # physical dimensions of the low end of the domain
geometry.prob_hi     =  2. 2. # physical dimensions of the high end of the domain
geometry.is_periodic = 1 0

# >>>>>>>>>>>>>  BC FLAGS <<<<<<<<<<<<<<<<
# 0 = Interior/Periodic  3 = Symmetry
# 1 = Inflow             4 = SlipWall
# 2 = Outflow            5 = NoSlipWall

ns.lo_bc             = 0 4  # boundary conditions on the low end of the domain
ns.hi_bc             = 0 4  # boundary conditions on the high end of the domain

#**********************************
#MORE OUTPUT CONTROLS
#**********************************
#
# Select form of FAB output: default is NATIVE
#
#   ASCII  (this is very slow)
#   NATIVE (native binary form on machine -- the default)
#   IEEE32 (useful if you want 32bit files when running in double precision)
#   8BIT   (eight-bit run-length-encoded)
#
fab.format = NATIVE

#**********************************
#MORE PARALLEL CONTROLS
#**********************************
#
# Initializes DistributionMapping strategy from ParmParse.
#
# ParmParse options are:
#
#   DistributionMapping.strategy = ROUNDROBIN
#   DistributionMapping.strategy = KNAPSACK
#
# The default strategy is ROUNDROBIN.
#
DistributionMapping.strategy = ROUNDROBIN
DistributionMapping.strategy = KNAPSACK

#**********************************
#GRIDDING & SOLVER CONTROLS (these are optional)
#**********************************

amr.blocking_factor  = 4   # factor by which every patch must be coarsenable
amr.n_error_buf      = 2 2 # number of buffer cells around each tagged cell

# necessary for convergence with this problem
mg.usecg = 0              
mg.nu_f = 100

# USE THIS ONLY FOR TESTING
#amr.max_grid_size = 64     # maximum linear dimension of any one grid
#
# StationData.vars     -- Names of StateData components to output
# StationData.coord    -- BL_SPACEDIM array of Reals
# StationData.coord    -- the next one
# StationData.coord    -- ditto ...
#
# e.g.
#
#StationData.vars  = density xmom ymom zmom
#StationData.coord = 0.1 0.1
#StationData.coord = 0.11 0.21
#StationData.coord = 0.5 0.5
#StationData.coord = 1.23 1.23
#StationData.coord = 0.23 1.53
#StationData.coord = 2.34 2.34
