# ------------------  INPUTS TO MAIN PROGRAM  -------------------
erf.prob_name = "Density Current"

erf.init_type = Isentropic

stop_time = 900.0

amr.blocking_factor_x = 16
amr.blocking_factor_y = 1
amr.blocking_factor_z = 16

# Compressible at level 0, anelastic at level 1
erf.anelastic   = 0 1
erf.mg_v = 1

erf.dt_ref_ratio = 1

erf.buoyancy_type = 2
erf.coupling_type = OneWay

amrex.fpe_trap_invalid = 1

fabarray.mfiter_tile_size = 1024 1024 1024

# PROBLEM SIZE & GEOMETRY
geometry.prob_lo     =      0.   0.    0.
geometry.prob_hi     =  25600. 100.  6400.

xlo.type = "Symmetry"
xhi.type = "Outflow"

zlo.type = "SlipWall"
zhi.type = "SlipWall"

geometry.is_periodic = 0 1 0

amr.max_level       = 1       # maximum level number allowed

amr.n_cell         = 128   1      32   # dx=dy=dz=100 m, Straka et al 1993 / Xue et al 2000
erf.fixed_dt       = 2.0       # fixed time step [s] -- Straka et al 1993
erf.fixed_fast_dt  = 0.5       # fixed time step [s] -- Straka et al 1993
erf.check_int      =-100       # number of timesteps between checkpoints

# DIAGNOSTICS & VERBOSITY
erf.sum_interval   =-1       # timesteps between computing mass
erf.v              = 1       # verbosity in ERF.cpp
amr.v              = 1       # verbosity in Amr.cpp

# CHECKPOINT FILES

# PLOTFILES
erf.check_file  = chk
erf.plot_file_1 = plt
erf.plot_int_1  = 50
erf.plot_vars_1 = density x_velocity y_velocity z_velocity pressure theta pert_pres pert_dens

# SOLVER CHOICE
erf.use_gravity  = true
erf.use_coriolis = false

erf.les_type         = "None"
#
# Diffusion coefficient from Straka, K = 75 m^2/s
#
erf.molec_diff_type   = "ConstantAlpha" # where alpha == "K" in Straka et al 1993
erf.dynamic_viscosity = 75.0 # [kg/(m-s)] ==> alpha = 75.0 m^2/s
erf.alpha_T           = 75.0 # [m^2/s]

erf.c_p = 1004.0

# PROBLEM PARAMETERS (optional)
prob.T_0 = 300.0

################################ MULTILEVEL ################################
amr.ref_ratio_vect = 2 1 2

erf.regrid_int = 10

erf.refinement_indicators = lo_theta

erf.lo_theta.max_level     = 1
erf.lo_theta.field_name    = theta
erf.lo_theta.value_less    = 299.0

amr.n_error_buf  = 6 6
amr.grid_eff     = 0.8
################################ MULTILEVEL ################################