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;