diff --git a/Exec/SquallLine_2D/inputs_moisture b/Exec/SquallLine_2D/inputs_moisture index 7ede7c674..2e2757fe5 100644 --- a/Exec/SquallLine_2D/inputs_moisture +++ b/Exec/SquallLine_2D/inputs_moisture @@ -68,9 +68,9 @@ erf.les_type = "None" 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.dynamicViscosity = 200.0 # [kg/(m-s)] ==> nu = 75.0 m^2/s erf.alpha_T = 00.0 # [m^2/s] -erf.alpha_C = 100.0 +erf.alpha_C = 25.0 erf.moisture_model = "Kessler" erf.use_moist_background = true diff --git a/Exec/SquallLine_2D/inputs_moisture_Gabersek b/Exec/SquallLine_2D/inputs_moisture_Gabersek deleted file mode 100644 index ad5b642b4..000000000 --- a/Exec/SquallLine_2D/inputs_moisture_Gabersek +++ /dev/null @@ -1,91 +0,0 @@ -# ------------------ 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