# ------------------  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.use_fft = true
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         = 256   1      64   # dx=dy=dz=100 m, Straka et al 1993 / Xue et al 2000
erf.fixed_dt       = 1.0       # fixed time step [s] -- Straka et al 1993
erf.fixed_fast_dt  = 0.25      # 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 = plth
erf.plot_int_1  = 50
erf.plot_vars_1 = density x_velocity y_velocity z_velocity pressure theta pert_pres pert_dens dpdx

# 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
prob.U_0 = 0.0

erf.terrain_z_levels = 0. 100. 200. 300. 400. 500. 600. 700. 800. 900. 1000. 1100. 1200. 1300. 1400. 1500. 1600. 1700. 1800. 1900. 2000. 2100. 2200. 2300. 2400. 2500. 2600. 2700. 2800. 2900. 3000. 3100. 3200. 3300. 3400. 3500. 3600. 3700. 3800. 3900. 4000. 4100. 4200. 4300. 4400. 4500. 4600. 4700. 4800. 4900. 5000. 5100. 5200. 5300. 5400. 5500. 5600. 5700. 5800. 5900. 6000. 6100. 6200. 6300. 6400.

################################ 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 ################################