From 13be62275e9ffad530174b7725ae1dc4454bb4c4 Mon Sep 17 00:00:00 2001 From: Ann Almgren Date: Wed, 16 Oct 2024 21:18:48 -0700 Subject: [PATCH] add ML inputs for DC --- Exec/RegTests/DensityCurrent/inputs_amr | 87 +++++++++++++++++++ Exec/RegTests/DensityCurrent/inputs_anelastic | 86 ------------------ .../DensityCurrent/inputs_crse_halfdomain | 68 +++++++-------- 3 files changed, 121 insertions(+), 120 deletions(-) create mode 100644 Exec/RegTests/DensityCurrent/inputs_amr delete mode 100644 Exec/RegTests/DensityCurrent/inputs_anelastic diff --git a/Exec/RegTests/DensityCurrent/inputs_amr b/Exec/RegTests/DensityCurrent/inputs_amr new file mode 100644 index 000000000..11a5373d4 --- /dev/null +++ b/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 ################################ diff --git a/Exec/RegTests/DensityCurrent/inputs_anelastic b/Exec/RegTests/DensityCurrent/inputs_anelastic deleted file mode 100644 index b555097ca..000000000 --- a/Exec/RegTests/DensityCurrent/inputs_anelastic +++ /dev/null @@ -1,86 +0,0 @@ -# ------------------ INPUTS TO MAIN PROGRAM ------------------- -stop_time = 900.0 - -erf.anelastic = 1 - -erf.use_terrain = false - -amrex.fpe_trap_invalid = 1 - -fabarray.mfiter_tile_size = 1024 1024 1024 - -#SYMMETRY / OUTFLOW VERSION -geometry.prob_lo = 0. 0. 0. -geometry.prob_hi = 25600. 100. 6400. -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 -xlo.type = "Symmetry" -xhi.type = "HO_Outflow" - -#DOUBLY PERIODIC VERSION -#geometry.prob_lo = -25600. 0. 0. -#geometry.prob_hi = 25600. 100. 6400. -#geometry.is_periodic = 1 1 0 -#amr.n_cell = 512 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 - -# DIAGNOSTICS & VERBOSITY -erf.sum_interval = 1 # timesteps between computing mass -erf.v = 1 # verbosity in ERF.cpp -amr.v = 1 # verbosity in Amr.cpp -erf.mg_v = 1 # verbosity in ERF.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 = netcdf # prefix of plotfile name -erf.plotfile_type = amrex # prefix of plotfile name -erf.plot_file_1 = plta # prefix of plotfile name -erf.plot_int_1 = 100 # 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_coriolis = 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.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 ################################ diff --git a/Exec/RegTests/DensityCurrent/inputs_crse_halfdomain b/Exec/RegTests/DensityCurrent/inputs_crse_halfdomain index a1ad57079..ef296b8fd 100644 --- a/Exec/RegTests/DensityCurrent/inputs_crse_halfdomain +++ b/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 ################################