more updates

integration/VODE/actual_integrator_simplified_sdc.H

@@ -82,6 +82,12 @@ void actual_integrator (BurnT& state, Real dt)
     vode_state.atol_spec = sdc_min_density * atol_spec * sdc_tol_fac;
     vode_state.rtol_spec = rtol_spec * sdc_tol_fac;
 
+    // if we are scaling the system, then update the atols accordingly
+    if (scale_system) {
+        vode_state.atol_enuc /= (state.rho_scale * state.e_scale);
+        vode_state.atol_spec /= state.rho_scale;
+    }
+
     // Call the integration routine.
 
     int istate = dvode(state, vode_state);

integration/integrator_type_simplified_sdc.H

@@ -81,9 +81,21 @@ void burn_to_int(BurnT& state, T& int_state)
 
     // store the original rho and rho e
     state.rho_orig = state.y[SRHO];
-
     state.rhoe_orig = state.y[SEINT];
 
+    // normalize
+    state.e_scale = state.rhoe_orig / state.rho_orig;
+    state.rho_scale = state.rho_orig;
+
+    if (scale_system) {
+        // we are integrating (rho X, rho e), so normalize by the
+        // initial scale
+        int_state.y(net_ienuc) /= (state.rho_scale * state.e_scale);
+        for (int n = 1; n <= NumSpec; ++n) {
+            int_state.y(n) /= state.rho;
+        }
+    }
+
 }