Skip to content

Commit

Permalink
add inputs
Browse files Browse the repository at this point in the history
  • Loading branch information
zingale committed Sep 12, 2023
1 parent 3c9869e commit 26d62a7
Show file tree
Hide file tree
Showing 2 changed files with 86 additions and 5 deletions.
10 changes: 5 additions & 5 deletions Exec/reacting_tests/bubble_convergence/converge_test.sh
Original file line number Diff line number Diff line change
Expand Up @@ -4,16 +4,16 @@ EXEC=./Castro2d.gnu.MPI.TRUESDC.ex

CONV_TOOL=RichardsonConvergenceTest2d.gnu.ex

mpiexec -n 8 ${EXEC} inputs_2d.32 amr.plot_file=bubble_32_sdc4_plt >& 32.out
mpiexec -n 8 ${EXEC} inputs_2d.64 amr.plot_file=bubble_64_sdc4_plt >& 64.out
mpiexec -n 8 ${EXEC} inputs_2d.128 amr.plot_file=bubble_128_sdc4_plt >& 128.out
mpiexec -n 16 ${EXEC} inputs_2d.32 amr.plot_file=bubble_32_sdc4_plt >& 32.out
mpiexec -n 16 ${EXEC} inputs_2d.64 amr.plot_file=bubble_64_sdc4_plt >& 64.out
mpiexec -n 16 ${EXEC} inputs_2d.128 amr.plot_file=bubble_128_sdc4_plt >& 128.out

${CONV_TOOL} coarFile=bubble_32_sdc4_plt00334 mediFile=bubble_64_sdc4_plt00667 fineFile=bubble_128_sdc4_plt01334 > sdc_converge.lo.out

mpiexec -n 8 ${EXEC} inputs_2d.256 amr.plot_file=bubble_256_sdc4_plt >& 256.out
mpiexec -n 16 ${EXEC} inputs_2d.256 amr.plot_file=bubble_256_sdc4_plt >& 256.out

${CONV_TOOL} coarFile=bubble_64_sdc4_plt00667 mediFile=bubble_128_sdc4_plt01334 fineFile=bubble_256_sdc4_plt02667 > sdc_converge.mid.out

mpiexec -n 8 ${EXEC} inputs_2d.512 amr.plot_file=bubble_512_sdc4_plt >& 512.out
mpiexec -n 16 ${EXEC} inputs_2d.512 amr.plot_file=bubble_512_sdc4_plt >& 512.out

${CONV_TOOL} coarFile=bubble_128_sdc4_plt01334 mediFile=bubble_256_sdc4_plt02667 fineFile=bubble_512_sdc4_plt05334 > sdc_converge.hi.out
81 changes: 81 additions & 0 deletions Exec/science/xrb_layered/inputs_2d
Original file line number Diff line number Diff line change
@@ -0,0 +1,81 @@
# ------------------ INPUTS TO MAIN PROGRAM -------------------
max_step = 1000000
stop_time = 1000.0

# PROBLEM SIZE & GEOMETRY
geometry.is_periodic = 1 0
geometry.coord_sys = 0 # 0 => cart, 1 => RZ 2=>spherical
geometry.prob_lo = 0.0 0.0
geometry.prob_hi = 3072 3072
amr.n_cell = 512 512

# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<<
# 0 = Interior 3 = Symmetry
# 1 = Inflow 4 = SlipWall
# 2 = Outflow 5 = NoSlipWall
# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<<
castro.lo_bc = 0 3
castro.hi_bc = 0 2

castro.fill_ambient_bc = 1
castro.ambient_fill_dir = 1
castro.ambient_outflow_vel = 1


# WHICH PHYSICS
castro.do_hydro = 1
castro.do_react = 1
castro.add_ext_src = 0
castro.do_grav = 1
castro.do_sponge = 1

castro.ppm_type = 1
castro.grav_source_type = 2
castro.use_pslope = 1
castro.pslope_cutoff_density = 1.e4

gravity.gravity_type = ConstantGrav
gravity.const_grav = -2.45e14

# TIME STEP CONTROL
castro.cfl = 0.7 # cfl number for hyperbolic system
castro.init_shrink = 0.1 # scale back initial timestep
castro.change_max = 1.1 # max time step growth

# SPONGE
castro.sponge_upper_density = 1.e3
castro.sponge_lower_density = 1.e2
castro.sponge_timescale = 1.0e-4

# DIAGNOSTICS & VERBOSITY
castro.sum_interval = 1 # timesteps between computing mass
castro.v = 1 # verbosity in Castro.cpp
amr.v = 1 # verbosity in Amr.cpp

# REFINEMENT / REGRIDDING
amr.max_level = 0 # maximum level number allowed
amr.ref_ratio = 2 2 2 2 # refinement ratio
amr.regrid_int = 2 2 2 2 # how often to regrid
amr.blocking_factor = 8 # block factor in grid generation
amr.max_grid_size = 64
amr.n_error_buf = 2 2 2 2 # number of buffer cells in error est

# CHECKPOINT FILES
amr.check_file = xrb_chk # root name of checkpoint file
amr.check_int = 1000 # number of timesteps between checkpoints

# PLOTFILES
amr.plot_file = xrb_plt # root name of plotfile
amr.plot_int = 1000 # number of timesteps between plotfiles
amr.derive_plot_vars = ALL

# PROBLEM PARAMETERS
problem.model_name = "toy_atm_3cm.hse"

problem.apply_perturbation = 1

# MICROPHYSICS
integrator.jacobian = 1

integrator.atol_spec = 1.e-6
integrator.rtol_spec = 1.e-6

0 comments on commit 26d62a7

Please sign in to comment.