From 627cc1e0948b9295a0dc408c1e51affc3209b03f Mon Sep 17 00:00:00 2001 From: Ryan Brady Date: Tue, 18 Jun 2024 10:21:20 -0400 Subject: [PATCH] For subch_planar: increase MAX_NPTS, make SDC time integration default (#2873) --- Exec/science/subch_planar/GNUmakefile | 4 +++- Exec/science/subch_planar/inputs_2d | 1 + Exec/science/subch_planar/problem_initialize.H | 2 +- Exec/science/subch_planar/problem_initialize_state_data.H | 2 +- 4 files changed, 6 insertions(+), 3 deletions(-) diff --git a/Exec/science/subch_planar/GNUmakefile b/Exec/science/subch_planar/GNUmakefile index e7eebb65ad..a9c4d6cae1 100644 --- a/Exec/science/subch_planar/GNUmakefile +++ b/Exec/science/subch_planar/GNUmakefile @@ -11,7 +11,9 @@ USE_OMP = FALSE USE_GRAV = TRUE USE_REACT = TRUE -USE_SHOCK_VAR = TRUE +USE_SHOCK_VAR = TRUE + +MAX_NPTS_MODEL = 32678 USE_MODEL_PARSER = TRUE diff --git a/Exec/science/subch_planar/inputs_2d b/Exec/science/subch_planar/inputs_2d index 9340972844..56d6d9854f 100644 --- a/Exec/science/subch_planar/inputs_2d +++ b/Exec/science/subch_planar/inputs_2d @@ -50,6 +50,7 @@ castro.sponge_timescale = 1.0e-3 castro.sum_interval = 1 # timesteps between computing mass castro.v = 1 # verbosity in Castro.cpp amr.v = 1 # verbosity in Amr.cpp +castro.time_integration_method = 3 # uses SDC by default # REFINEMENT / REGRIDDING amr.max_level = 0 # maximum level number allowed diff --git a/Exec/science/subch_planar/problem_initialize.H b/Exec/science/subch_planar/problem_initialize.H index 007cdf644c..01bced54c2 100644 --- a/Exec/science/subch_planar/problem_initialize.H +++ b/Exec/science/subch_planar/problem_initialize.H @@ -14,7 +14,7 @@ void problem_initialize () const Real* problo = dgeom.ProbLo(); const Real* probhi = dgeom.ProbHi(); - + problem::ihe4 = network_spec_index("helium-4"); if (problem::ihe4 < 0) { diff --git a/Exec/science/subch_planar/problem_initialize_state_data.H b/Exec/science/subch_planar/problem_initialize_state_data.H index ea5d3b0b0a..0f29bd07e1 100644 --- a/Exec/science/subch_planar/problem_initialize_state_data.H +++ b/Exec/science/subch_planar/problem_initialize_state_data.H @@ -78,7 +78,7 @@ void problem_initialize_state_data (int i, int j, int k, for (int n = 0; n < NumSpec; n++) { state(i,j,k,UFS+n) = state(i,j,k,URHO) * state(i,j,k,UFS+n); } - + // Initial velocities = 0 state(i,j,k,UMX) = 0.0_rt; state(i,j,k,UMY) = 0.0_rt;