add ML inputs for DC (#1892)

Exec/RegTests/DensityCurrent/inputs_amr

@@ -0,0 +1,87 @@
+# ------------------ INPUTS TO MAIN PROGRAM -------------------
+stop_time = 900.0
+
+amr.blocking_factor_x = 16
+amr.blocking_factor_y = 1
+amr.blocking_factor_z = 16
+
+#erf.anelastic = 1
+#erf.check_file = chka
+#erf.plot_file_1 = plta
+
+erf.anelastic = 0
+erf.check_file = chk
+erf.plot_file_1 = plt
+
+erf.buoyancy_type = 2
+
+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 = 128 1 32 # dx=dy=dz=100 m, Straka et al 1993 / Xue et al 2000
+erf.fixed_dt = 2.0 # fixed time step [s] -- Straka et al 1993
+erf.fixed_fast_dt = 0.5 # fixed time step [s] -- Straka et al 1993
+erf.plot_int_1 = 100 # number of timesteps between plotfiles
+erf.check_int =-1000 # number of timesteps between checkpoints
+
+# DIAGNOSTICS & VERBOSITY
+erf.sum_interval =-1 # timesteps between computing mass
+erf.v = 0 # verbosity in ERF.cpp
+amr.v = 1 # verbosity in Amr.cpp
+
+# CHECKPOINT FILES
+
+# PLOTFILES
+erf.plotfile_type = amrex
+erf.plot_vars_1 = density x_velocity y_velocity z_velocity pressure theta pres_hse dens_hse pert_pres pert_dens
+
+# SOLVER CHOICE
+erf.use_gravity = true
+erf.use_coriolis = false
+erf.use_terrain = 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.rho0_trans = 1.0 # [kg/m^3], used to convert input diffusivities
+erf.dynamicViscosity = 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
+
+################################ MULTILEVEL ################################
+amr.ref_ratio_vect = 2 1 2
+
+erf.coupling_type = "OneWay"
+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 ################################

Exec/RegTests/DensityCurrent/inputs_crse_halfdomain

@@ -1,7 +1,13 @@
 # ------------------ INPUTS TO MAIN PROGRAM -------------------
 stop_time = 900.0
 
-erf.use_terrain = false
+erf.anelastic = 1
+erf.check_file = chka
+erf.plot_file_1 = plta
+
+#erf.anelastic = 0
+#erf.check_file = chk
+#erf.plot_file_1 = plt
 
 erf.buoyancy_type = 2
 
@@ -16,38 +22,49 @@ geometry.prob_hi = 25600. 100. 6400.
 xlo.type = "Symmetry"
 xhi.type = "Outflow"
 
-geometry.is_periodic = 0 1 0
-
-amr.n_cell = 256 1 64 # dx=dy=dz=100 m, Straka et al 1993 / Xue et al 2000
-
 zlo.type = "SlipWall"
 zhi.type = "SlipWall"
 
-# TIME STEP CONTROL
-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
+geometry.is_periodic = 0 1 0
+
+amr.max_level = 0 # maximum level number allowed
+
+#fine
+#amr.n_cell = 512 1 128 # dx=dy=dz=50 m, Straka et al 1993 / Xue et al 2000
+#erf.fixed_dt = 0.5 # fixed time step [s] -- Straka et al 1993
+#erf.fixed_fast_dt = 0.125 # fixed time step [s] -- Straka et al 1993
+#erf.plot_int_1 = 200 # number of timesteps between plotfiles
+#erf.check_int = 200 # number of timesteps between checkpoints
+
+#crse
+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.plot_int_1 = 300 # number of timesteps between plotfiles
+erf.check_int =-1000 # number of timesteps between checkpoints
+
+#crser
+#amr.n_cell = 128 1 32 # dx=dy=dz=100 m, Straka et al 1993 / Xue et al 2000
+#erf.fixed_dt = 2.0 # fixed time step [s] -- Straka et al 1993
+#erf.fixed_fast_dt = 0.5 # fixed time step [s] -- Straka et al 1993
+#erf.plot_int_1 = 150 # number of timesteps between plotfiles
+#erf.check_int = 1 # number of timesteps between checkpoints
 
 # DIAGNOSTICS & VERBOSITY
-erf.sum_interval = 1 # timesteps between computing mass
+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 = -1000 # number of timesteps between checkpoints
 
 # PLOTFILES
 erf.plotfile_type = amrex # prefix of plotfile name
-erf.plot_file_1 = pltc # prefix of plotfile name
-erf.plot_int_1 = 1000 # 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
 
 # SOLVER CHOICE
-erf.use_gravity = true
+erf.use_gravity  = true
 erf.use_coriolis = false
+erf.use_terrain = false
 
 erf.les_type = "None"
 #
@@ -63,20 +80,3 @@ erf.c_p = 1004.0
 # PROBLEM PARAMETERS (optional)
 prob.T_0 = 300.0
 prob.U_0 = 0.0
-
-################################ MULTILEVEL ################################
-amr.max_level = 0
-amr.ref_ratio_vect = 2 1 2
-
-erf.coupling_type = "TwoWay"
-erf.regrid_int = 2
-
-erf.refinement_indicators = lo_theta
-
-erf.lo_theta.max_level = 1
-erf.lo_theta.field_name = theta
-erf.lo_theta.value_less = 299.9
-
-amr.n_error_buf = 5 5
-amr.grid_eff = 0.8
-################################ MULTILEVEL ################################