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

some updates for the He merger inputs #3008

Merged
merged 1 commit into from
Dec 6, 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
33 changes: 19 additions & 14 deletions Exec/science/wdmerger/tests/he_double_det/inputs_pakmor_simp_sdc
Original file line number Diff line number Diff line change
Expand Up @@ -61,7 +61,7 @@ castro.hi_bc = 2 2 2
max_step = 10000000

# Simulation end time
stop_time = 200.0
stop_time = 1000.0

# CFL number for hyperbolic system
castro.cfl = 0.5
Expand All @@ -70,10 +70,10 @@ castro.cfl = 0.5
castro.fixed_dt = -1.0

# Scale back initial timestep by this factor
castro.init_shrink = 0.01
castro.init_shrink = 0.1

# Factor by which dt is allowed to change each timestep
castro.change_max = 1.10
castro.change_max = 1.20

# If we regrid on Level 0, compute a new timestep afterward
amr.compute_new_dt_on_regrid = 1
Expand All @@ -91,6 +91,8 @@ castro.abundance_failure_rho_cutoff = 1.0e0
# Default is 10, 16 is recommended value
castro.max_subcycles = 16

# Switching to simplified SDC
castro.time_integration_method = 3

############################################################################################
# Resolution, gridding and AMR
Expand Down Expand Up @@ -292,20 +294,20 @@ integrator.atol_spec = 1.0e-6
integrator.rtol_enuc = 1.0e-6
integrator.atol_enuc = 1.0e-6

# Do not abort or retry on a failed burn (Castro will handle this)
integrator.abort_on_failure = 0

# Renormalize abundances during the burn
integrator.renormalize_abundances = 1

# Maximum temperature allowed in the burn
integrator.MAX_TEMP = 1.0e10

# Use tabular rate evaluation when available
network.use_tables = 1
network.use_tables = 0

# Switching to simplified SDC
castro.time_integration_method = 3
# use retry in the burner and cap the number of steps to make the
# reactions go faster
integrator.use_burn_retry = 1
integrator.retry_swap_jacobian = 1
integrator.ode_max_steps = 10000

############################################################################################
# Gravity
Expand Down Expand Up @@ -380,16 +382,16 @@ amr.check_file = chk
castro.dump_old = 1

# Simulation time between checkpoints
amr.check_per = 1.0
amr.check_per = -1

# Number of timesteps between checkpoints
amr.check_int = -1
amr.check_int = 50

# Root name of plot files
amr.plot_file = plt

# Simulation time between plotfiles
amr.plot_per = 1.0
amr.plot_per = 2.0

# Number of timesteps between plotfiles
amr.plot_int = -1
Expand Down Expand Up @@ -426,17 +428,20 @@ gravity.v = 1
amr.plot_vars = ALL

# Derived variables to add to plot files
amr.derive_plot_vars = pressure
amr.derive_plot_vars = ALL

# State variables to add to small plot files
amr.small_plot_vars = density Temp

# Derived variables to add to small plot files
amr.derive_small_plot_vars = enuc X(He4) X(C12) X(O16)
amr.derive_small_plot_vars = abar enuc MachNumber magvel magvort X(He4) X(C12) X(O16)

# Name of the diagnostic sum output files
amr.data_log = star_diag.out primary_diag.out secondary_diag.out rotation_diag.out

# write plotfiles as single precision
fab.format = NATIVE_32

############################################################################################
# Problem parameters
############################################################################################
Expand Down
24 changes: 15 additions & 9 deletions Exec/science/wdmerger/tests/he_double_det/inputs_pakmor_strang
Original file line number Diff line number Diff line change
Expand Up @@ -292,17 +292,20 @@ integrator.atol_spec = 1.0e-6
integrator.rtol_enuc = 1.0e-6
integrator.atol_enuc = 1.0e-6

# Do not abort or retry on a failed burn (Castro will handle this)
integrator.abort_on_failure = 0

# Renormalize abundances during the burn
integrator.renormalize_abundances = 1

# Maximum temperature allowed in the burn
integrator.MAX_TEMP = 1.0e10

# Use tabular rate evaluation when available
network.use_tables = 1
network.use_tables = 0

# use retry in the burner and cap the number of steps to make the
# reactions go faster
integrator.use_burn_retry = 1
integrator.retry_swap_jacobian = 1
integrator.ode_max_steps = 10000

############################################################################################
# Gravity
Expand Down Expand Up @@ -377,16 +380,16 @@ amr.check_file = chk
castro.dump_old = 1

# Simulation time between checkpoints
amr.check_per = 1.0
amr.check_per = -1

# Number of timesteps between checkpoints
amr.check_int = -1
amr.check_int = 50

# Root name of plot files
amr.plot_file = plt

# Simulation time between plotfiles
amr.plot_per = 1.0
amr.plot_per = 2.0

# Number of timesteps between plotfiles
amr.plot_int = -1
Expand Down Expand Up @@ -423,17 +426,20 @@ gravity.v = 1
amr.plot_vars = ALL

# Derived variables to add to plot files
amr.derive_plot_vars = pressure
amr.derive_plot_vars = ALL

# State variables to add to small plot files
amr.small_plot_vars = density Temp

# Derived variables to add to small plot files
amr.derive_small_plot_vars = enuc X(He4) X(C12) X(O16)
amr.derive_small_plot_vars = abar enuc MachNumber magvel magvort X(He4) X(C12) X(O16)

# Name of the diagnostic sum output files
amr.data_log = star_diag.out primary_diag.out secondary_diag.out rotation_diag.out

# write plotfiles as single precision
fab.format = NATIVE_32

############################################################################################
# Problem parameters
############################################################################################
Expand Down
Loading