diff --git a/Exec/IdealMiniGrid/inputs_chapman_flather b/Exec/IdealMiniGrid/inputs_chapman_flather
new file mode 100644
index 0000000..9794520
--- /dev/null
+++ b/Exec/IdealMiniGrid/inputs_chapman_flather
@@ -0,0 +1,80 @@
+# ------------------  INPUTS TO MAIN PROGRAM  -------------------
+max_step = 100
+
+amrex.fpe_trap_invalid = 1
+
+fabarray.mfiter_tile_size = 1024 1024 1024
+
+# PROBLEM SIZE & GEOMETRY
+# dims come from ROMS ana_grid, and must match vals in prob.cpp right now
+geometry.prob_lo     =   0.       0.  -5000.
+geometry.prob_hi     =   320000.  320000.       0.
+
+amr.n_cell           =  10 16 20
+
+geometry.is_periodic = 0 0 0
+zlo.type = "SlipWall"
+zhi.type = "SlipWall"
+
+remora.boundary_per_variable = 1
+##             WEST     SOUTH    EAST     NORTH
+temp.type   =  clamped  clamped  clamped  clamped
+salt.type   =  clamped  clamped  clamped  clamped
+scalar.type =  clamped  clamped  clamped  clamped
+u.type      =  clamped  clamped  clamped  clamped
+v.type      =  clamped  clamped  clamped  clamped
+w.type      =  clamped  clamped  clamped  clamped
+ubar.type   =  flather  flather  flather  flather
+vbar.type   =  flather  flather  flather  flather
+zeta.type   =  chapman  chapman  chapman  chapman
+tke.type    =  clamped  clamped  clamped  clamped
+
+# TIME STEP CONTROL
+remora.fixed_dt            = 200.0 # Timestep size (seconds)
+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)
+amr.v                = 1       # verbosity in Amr.cpp
+
+# REFINEMENT / REGRIDDING
+amr.max_level        = 0       # maximum level number allowed
+#amr.max_grid_size    = 600     # maximum level number allowed
+
+# 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     = salt temp x_velocity y_velocity z_velocity
+remora.plotfile_type = amrex
+
+# SOLVER CHOICE
+remora.tracer_horizontal_advection_scheme = "upstream3" # upstream3 or centered4
+remora.spatial_order = 2
+
+remora.use_coriolis  = true
+remora.coriolis_type = real
+remora.coriolis_f0   = 1.0e-4
+
+remora.theta_s = 5.0
+remora.theta_b = 0.5
+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 = 0.0  # 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
+
+remora.ic_bc_type       = "real"
+remora.nc_init_file_0   = "idmini_ini_v1.nc"
+remora.nc_grid_file_0   = "idmini_grd_v1.nc"
+remora.nc_bdry_file     = "idmini_bry_Uvel.nc"
diff --git a/Source/Initialization/REMORA_init_bcs.cpp b/Source/Initialization/REMORA_init_bcs.cpp
index d377d2c..1e990f7 100644
--- a/Source/Initialization/REMORA_init_bcs.cpp
+++ b/Source/Initialization/REMORA_init_bcs.cpp
@@ -156,7 +156,7 @@ void REMORA::init_bcs ()
     set_bcs_by_var = false;
 
     ParmParse pp("remora");
-    pp.query("boundary_per_side", set_bcs_by_var);
+    pp.query("boundary_per_variable", set_bcs_by_var);
     if (!set_bcs_by_var) {
         f_by_side("xlo", Orientation(Direction::x,Orientation::low));
         f_by_side("xhi", Orientation(Direction::x,Orientation::high));