# ------------------  INPUTS TO MAIN PROGRAM  -------------------
erf.prob_name = "Particles Over Witch of Agnesi Hill"

erf.init_type = Isentropic

max_step =  80

amrex.fpe_trap_invalid = 1

fabarray.mfiter_tile_size = 1024 1024 1024

# PROBLEM SIZE & GEOMETRY
geometry.prob_lo     = 0.   0.  0. 
geometry.prob_hi     = 10.  1.  2.

amr.n_cell           = 256  8   64

geometry.is_periodic = 0 1 0

xlo.type = "Inflow"
xhi.type = "Outflow"
xlo.velocity = 10. 0. 0.
xlo.density  = 1.16
xlo.theta    = 300.
xlo.scalar   = 0.
    
zlo.type = "SlipWall"
zhi.type = "SlipWall"

# PARTICLES
erf.use_tracer_particles = 1
tracer_particles.initial_distribution_type = box
tracer_particles.particle_box_lo = 3.95 -1.0 -1.0
tracer_particles.particle_box_hi = 4.00  2.0  1.0
tracer_particles.place_randomly_in_cells = true

# TIME STEP CONTROL
erf.fixed_dt           = 1E-3
erf.fixed_mri_dt_ratio = 12

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

# REFINEMENT / REGRIDDING
amr.max_level       = 0       # maximum level number allowed

# CHECKPOINT FILES
erf.check_file      = chk     # root name of checkpoint file
erf.check_int       = 100 # number of timesteps between checkpoints

# PLOTFILES
erf.plot_file_1     = plt     # prefix of plotfile name
erf.plot_int_1      = 50      # number of timesteps between plotfiles
erf.plot_vars_1     = density x_velocity y_velocity z_velocity pressure theta pres_hse dens_hse pert_pres pert_dens z_phys detJ dpdx dpdy pres_hse_x pres_hse_y tracer_particles_count tracer_particles_mass_density

# SOLVER CHOICE
erf.use_gravity = true
erf.use_coriolis = false
erf.les_type = "None"

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

erf.refinement_indicators = box1
erf.box1.max_level = 1
erf.box1.in_box_lo =  2. 0.25
erf.box1.in_box_hi =  8. 0.75

# TERRRAIN GRID TYPE
erf.terrain_type = StaticFittedMesh
erf.terrain_smoothing = 0

erf.dycore_horiz_adv_type  = Centered_2nd
erf.dycore_vert_adv_type   = Centered_2nd
erf.dryscal_horiz_adv_type = Centered_2nd
erf.dryscal_vert_adv_type  = Centered_2nd

# Constant diffusion coefficient
erf.molec_diff_type   = "ConstantAlpha"
erf.dynamic_viscosity = 0.0 # [kg/(m-s)]
erf.alpha_T           = 0.0 # [m^2/s]

#erf.abl_driver_type = "PressureGradient"
#erf.abl_pressure_grad = -0.2 0. 0.

# PROBLEM PARAMETERS (optional)
prob.T_0   = 300.0
prob.U_0   = 10.0
prob.rho_0 = 1.16

prob.custom_terrain_type = "WoA"