Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

add ML inputs for DC #1892

Merged
merged 1 commit into from
Oct 17, 2024
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
87 changes: 87 additions & 0 deletions Exec/RegTests/DensityCurrent/inputs_amr
Original file line number Diff line number Diff line change
@@ -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 ################################
86 changes: 0 additions & 86 deletions Exec/RegTests/DensityCurrent/inputs_anelastic

This file was deleted.

68 changes: 34 additions & 34 deletions Exec/RegTests/DensityCurrent/inputs_crse_halfdomain
Original file line number Diff line number Diff line change
@@ -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

Expand All @@ -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"
#
Expand All @@ -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 ################################
Loading