Add GABLS1 PBL column inputs

Exec/ABL/inputs_gabls1_mynn25_smag

@@ -0,0 +1,74 @@
+# ------------------  INPUTS TO MAIN PROGRAM  -------------------
+stop_time = 32400.0  # 540 min = 9 h (Cuxart et al. 2006)
+
+amrex.fpe_trap_invalid = 1
+
+fabarray.mfiter_tile_size = 1024 1024 1024
+
+# PROBLEM SIZE & GEOMETRY (Cuxart et al. 2006)
+geometry.prob_extent = 25  25  400
+amr.n_cell           =  4   4   64
+
+geometry.is_periodic = 1 1 0
+
+# MOST BOUNDARY (DEFAULT IS ADIABATIC FOR THETA)
+zlo.type                    = "Most"
+erf.most.z0                 = 0.1  # from Cuxart et al. 2006
+erf.most.zref               = 3.125  # == dz/2
+erf.most.surf_temp_flux     = 0.0
+#erf.most.surf_temp          = 265.0 # initial value, should match input_sounding
+#erf.most.surf_heating_rate  = -0.25 # [K/h] from Cuxart et al. 2006
+
+zhi.type        = "SlipWall"
+zhi.theta_grad  = 0.01  # [K/m] to match the input sounding
+
+# INITIALIZATION (Cuxart et al. 2006)
+erf.init_type           = "input_sounding"
+erf.init_sounding_ideal = 1
+
+# TIME STEP CONTROL
+erf.fixed_dt        = 10.0  # fixed time step (Cuxart et al. 2006)
+erf.no_substepping  = 1
+
+# 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      = 1        # number of timesteps between plotfiles
+erf.plot_vars_1     = density rhoadv_0 x_velocity y_velocity z_velocity pressure temp theta
+
+
+# SOLVER CHOICE
+erf.alpha_T = 0.0
+erf.alpha_C = 1.0
+erf.use_gravity = true
+
+# Coriolis (Cuxart et al. 2006)
+# f = 1.39e-4 s^-1
+erf.use_coriolis = true
+erf.latitude = 73.0
+erf.rotational_time_period = 86455.2516813368
+
+# Geostrophic wind (Cuxart et al. 2006)
+erf.abl_driver_type = "GeostrophicWind"
+erf.abl_geo_wind = 8.0 0.0 0.0
+
+# Turbulence closure
+erf.molec_diff_type = "None"
+erf.les_type        = "Smagorinsky"
+erf.Cs              = 0.1
+erf.rho0_trans      = 1.3223 # from Cuxart et al. 2006
+erf.theta_ref       = 263.5 # from Cuxart et al. 2006
+
+erf.pbl_type    = "MYNN2.5"
+#erf.QKE_0       = 0.8  # == 2*KE_0; note: Cuxart et al. 2006 has an initial TKE profile