For subch_planar: increase MAX_NPTS, make SDC time integration default (

#2873)
AMReX-Astro · Jun 18, 2024 · 627cc1e · 627cc1e
1 parent 7cf5db7
commit 627cc1e
Show file tree

Hide file tree

Showing 4 changed files with 6 additions and 3 deletions.
diff --git 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
@@ -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
@@ -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
@@ -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;