# ------------------  INPUTS TO MAIN PROGRAM  -------------------
remora.prob_name = Seamount

remora.max_step = 10
remora.stop_time = 30000.0

# PROBLEM SIZE & GEOMETRY
# dims come from ROMS ana_grid, and must match vals in prob.cpp right now
remora.prob_lo     =   0.       0.  -5000.
remora.prob_hi     =   320000.  320000.       0.

remora.n_cell           =  49   48   13

remora.is_periodic = 1 1 0

# TIME STEP CONTROL
remora.fixed_dt       = 60.0 # Timestep size (seconds)
# NDTFAST  = 30.0 # Number of baratropic steps => 300.0/30.0 = 10.0
#remora.fixed_fast_dt  = 10.0 # Baratropic timestep size (seconds)
# remora.fixed_fast_dt  = 300.0 # Baratropic timestep size (seconds) testing value
remora.fixed_ndtfast_ratio  = 20 # Baratropic timestep size (seconds)

remora.flat_bathymetry=0

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

# 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      = 1          # number of timesteps between plotfiles
remora.plot_vars_3d  = salt temp x_velocity y_velocity z_velocity
remora.plotfile_type = amrex

# SOLVER CHOICE
remora.use_coriolis = true
remora.tracer_horizontal_advection_scheme = "upstream3" # upstream3 or centered4

remora.coriolis_f0 = 1.0e-4 #-8.26e-5

remora.theta_s = 6.5
remora.theta_b = 2.0
remora.tcline = 100.

# PROBLEM PARAMETERS (optional)
remora.R0    = 1027.0  # background density value (Kg/m3) used in Linear Equation of State
remora.S0    = 32.0    # background salinity (nondimensional) constant
remora.T0    = 10.0    # background potential temperature (Celsius) constant
remora.Tcoef = 1.7e-4  # linear equation of state parameter (1/Celsius)
remora.Scoef = 0.0     # linear equation of state parameter (nondimensional)
remora.rho0  = 1025.0  # Mean density (Kg/m3) used when Boussinesq approx is inferred