From 14073876bc2c1f404de1303a68ed1cd111d9631d Mon Sep 17 00:00:00 2001 From: Mahesh Natarajan Date: Sat, 2 Dec 2023 16:38:09 -0800 Subject: [PATCH] Adding inputs for Gabersek test case: --- Exec/SquallLine_2D/inputs_moisture_Gabersek | 91 +++++++++++++++++++++ 1 file changed, 91 insertions(+) create mode 100644 Exec/SquallLine_2D/inputs_moisture_Gabersek diff --git a/Exec/SquallLine_2D/inputs_moisture_Gabersek b/Exec/SquallLine_2D/inputs_moisture_Gabersek new file mode 100644 index 000000000..ad5b642b4 --- /dev/null +++ b/Exec/SquallLine_2D/inputs_moisture_Gabersek @@ -0,0 +1,91 @@ +# ------------------ INPUTS TO MAIN PROGRAM ------------------- +max_step = 18000 +stop_time = 90000.0 + +amrex.fpe_trap_invalid = 1 + +fabarray.mfiter_tile_size = 2048 1024 2048 + +# PROBLEM SIZE & GEOMETRY +geometry.prob_lo = -60000. 0. 0. +geometry.prob_hi = 60000. 400. 20000. +amr.n_cell = 384 4 128 # dx=dy=dz=100 m + +# periodic in x to match WRF setup +# - as an alternative, could use symmetry at x=0 and outflow at x=25600 +geometry.is_periodic = 1 1 0 +#xlo.type = "Outflow" +#xhi.type = "Outflow" +zlo.type = "SlipWall" +zhi.type = "Outflow" + +erf.sponge_strength = 2.0 +#erf.use_zhi_sponge_damping = true +erf.zhi_sponge_start = 12000.0 + +erf.sponge_density = 1.2 +erf.sponge_x_velocity = 0.0 +erf.sponge_y_velocity = 0.0 +erf.sponge_z_velocity = 0.0 + +# TIME STEP CONTROL +erf.use_native_mri = 1 +erf.fixed_dt = 0.5 # fixed time step [s] -- Straka et al 1993 +erf.fixed_fast_dt = 0.25 # fixed time step [s] -- Straka et al 1993 +#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 +amr.check_file = chk # root name of checkpoint file +amr.check_int = 1000 # number of timesteps between checkpoints +#amr.restart = chk09000 + +# PLOTFILES +erf.plot_file_1 = plt # root name of plotfile +erf.plot_int_1 = 100 # number of timesteps between plotfiles +erf.plot_vars_1 = density rhotheta rhoQt rhoQp x_velocity y_velocity z_velocity pressure theta temp qt qp qv qc qi scalar pert_dens + +# SOLVER CHOICE +erf.use_gravity = true +erf.buoyancy_type = 1 +erf.use_coriolis = false +erf.use_rayleigh_damping = false + +#erf.les_type = "Smagorinsky" +erf.Cs = 0.25 +erf.les_type = "None" + +# +# diffusion coefficient from Straka, K = 75 m^2/s +# +erf.molec_diff_type = "ConstantAlpha" +#erf.molec_diff_type = "Constant" +erf.rho0_trans = 1.0 # [kg/m^3], used to convert input diffusivities +erf.dynamicViscosity = 100.0 # [kg/(m-s)] ==> nu = 75.0 m^2/s +erf.alpha_T = 00.0 # [m^2/s] +erf.alpha_C = 50.0 + +erf.moisture_model = "Kessler" +erf.use_moist_background = true + +erf.moistscal_horiz_adv_string = "Centered_2nd" +erf.moistscal_vert_adv_string = "Centered_2nd" + +# PROBLEM PARAMETERS (optional) +prob.z_tr = 12000.0 +prob.height = 1200.0 +prob.theta_0 = 300.0 +prob.theta_tr = 343.0 +prob.T_tr = 213.0 +prob.x_c = 0.0 +prob.z_c = 2000.0 +prob.x_r = 10000.0 +prob.z_r = 1500.0 +prob.theta_c = 3.0