# ------------------  INPUTS TO MAIN PROGRAM  -------------------
remora.max_step = 100

# PROBLEM SIZE & GEOMETRY
remora.prob_lo     =      0.     0.    -47.5
remora.prob_hi     =  100000. 5000.       0.

remora.n_cell           =  39     4      30

remora.is_periodic = 0 1 0

remora.boundary_per_variable = 1
##                      WEST              SOUTH             EAST              NORTH
remora.bc.temp.type   =  slipwall         periodic          orlanski_rad      periodic
remora.bc.salt.type   =  slipwall         periodic          orlanski_rad      periodic
remora.bc.scalar.type =  slipwall         periodic          orlanski_rad      periodic
remora.bc.u.type      =  slipwall         periodic          orlanski_rad      periodic
remora.bc.v.type      =  slipwall         periodic          orlanski_rad      periodic
remora.bc.ubar.type   =  slipwall         periodic          outflow           periodic
remora.bc.vbar.type   =  slipwall         periodic          outflow           periodic
remora.bc.zeta.type   =  slipwall         periodic          outflow           periodic
remora.bc.tke.type    =  slipwall         periodic          outflow           periodic


# TIME STEP CONTROL
remora.fixed_dt       = 50.0 # Timestep size (seconds)
remora.fixed_ndtfast_ratio = 30

# DIAGNOSTICS & VERBOSITY
remora.sum_interval  = 1       # timesteps between integrated/max quantities, if remora.v > 0
remora.v             = 0       # verbosity in REMORA.cpp (0: none, 1: integrated quantities, etc, 2: print boxes)

remora.bottom_stress_type = "quadratic"

# CHECKPOINT FILES
remora.check_file      = chk        # root name of checkpoint file
remora.check_int       = -57600      # number of timesteps between checkpoints

# PLOTFILES
remora.plot_file     = plt        # prefix of plotfile name
remora.plot_int      = 100        # number of timesteps between plotfiles
remora.plot_vars_3d  = salt temp x_velocity y_velocity z_velocity
remora.plotfile_type = amrex

# SOLVER CHOICE
remora.flat_bathymetry = false
remora.tracer_horizontal_advection_scheme = "upstream3" # upstream3 or centered4

remora.horizontal_mixing_type = "constant"
remora.visc2 = 0.0

remora.grid_scale_type = "analytic"

remora.eos_type = "nonlinear"

remora.bulk_fluxes = true

remora.air_temperature = 23.567
remora.air_humidity = 0.776
remora.air_pressure = 1013.48
remora.cloud = 0.5
remora.rain = 0.002
remora.eminusp = true
remora.eminusp_correct_ssh = true

# Linear EOS parameters
remora.R0    = 1027.0  # background density value (Kg/m3) used in Linear Equation of State
remora.S0    = 31.0    # background salinity (nondimensional) constant
remora.T0    = 17.25    # background potential temperature (Celsius) constant
remora.rho0  = 1025.0  # Mean density (Kg/m3) used when Boussinesq approx is inferred

remora.theta_s = 5.0
remora.theta_b = 1.5
remora.tcline = 20.0

remora.vertical_mixing_type="gls"
remora.Akt_bak = 1.0e-6
remora.Akv_bak = 1.0e-5

# Coriolis params
remora.use_coriolis = true
remora.coriolis_type = beta_plane
remora.coriolis_f0 = 9.25e-4
remora.coriolis_beta = 0.0