diff --git a/Exec/SquallLine_2D/inputs_moisture b/Exec/SquallLine_2D/inputs_moisture index 8364107c7..f82dd6612 100644 --- a/Exec/SquallLine_2D/inputs_moisture +++ b/Exec/SquallLine_2D/inputs_moisture @@ -1,5 +1,5 @@ # ------------------ INPUTS TO MAIN PROGRAM ------------------- -max_step = 14400 +max_step = 9000 stop_time = 90000.0 amrex.fpe_trap_invalid = 1 @@ -54,7 +54,7 @@ erf.plot_vars_1 = density rhotheta rhoQt rhoQp x_velocity y_velocity z_v # SOLVER CHOICE erf.use_gravity = true -erf.buoyancy_type = 1 +erf.buoyancy_type = 4 erf.use_coriolis = false erf.use_rayleigh_damping = false @@ -70,7 +70,7 @@ erf.molec_diff_type = "ConstantAlpha" erf.rho0_trans = 1.0 # [kg/m^3], used to convert input diffusivities erf.dynamicViscosity = 200.0 # [kg/(m-s)] ==> nu = 75.0 m^2/s erf.alpha_T = 00.0 # [m^2/s] -erf.alpha_C = 00.0 +erf.alpha_C = 100.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 index 83eb01f20..ad5b642b4 100644 --- a/Exec/SquallLine_2D/inputs_moisture_Gabersek +++ b/Exec/SquallLine_2D/inputs_moisture_Gabersek @@ -1,5 +1,5 @@ # ------------------ INPUTS TO MAIN PROGRAM ------------------- -max_step = 9000 +max_step = 18000 stop_time = 90000.0 amrex.fpe_trap_invalid = 1 @@ -9,7 +9,7 @@ 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 = 200 4 100 # dx=dy=dz=100 m +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 @@ -30,8 +30,8 @@ erf.sponge_z_velocity = 0.0 # TIME STEP CONTROL erf.use_native_mri = 1 -erf.fixed_dt = 1.0 # fixed time step [s] -- Straka et al 1993 -erf.fixed_fast_dt = 0.5 # fixed time step [s] -- Straka et al 1993 +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 @@ -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 = 25.0 # [kg/(m-s)] ==> nu = 75.0 m^2/s -erf.alpha_T = 2.5 # [m^2/s] -erf.alpha_C = 25.0 +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