Refactor radiation variable names

Author: James471

integration/VODE/vode_dvode.H

@@ -82,7 +82,7 @@ int dvode (BurnT& state, DvodeT& vstate)
     vstate.ewt(NumSpec+1) = vstate.rtol_enuc * std::abs(vstate.yh(NumSpec+1,1)) + vstate.atol_enuc;
     vstate.ewt(NumSpec+1) = 1.0_rt / vstate.ewt(NumSpec+1);
 #ifdef PHOTOCHEMISTRY
-    for (int i = 1; i <= NumChemActiveRadGroups; i+=MicrophysicsNumRadVarsPerGroup) {
+    for (int i = 1; i <= NumChemBands; i+=MicrophysicsNumRadVarsPerGroup) {
         int rad_n_index = NumSpec + 2 + (i-1) * MicrophysicsNumRadVarsPerGroup;
         vstate.ewt(rad_n_index) = vstate.rtol_rad_num * std::abs(vstate.yh(rad_n_index,1)) + vstate.atol_rad_num;
         for (int j = 2; j <= MicrophysicsNumRadVarsPerGroup; ++j){
@@ -165,7 +165,7 @@ int dvode (BurnT& state, DvodeT& vstate)
            vstate.ewt(NumSpec+1) = vstate.rtol_enuc * std::abs(vstate.yh(NumSpec+1,1)) + vstate.atol_enuc;
            vstate.ewt(NumSpec+1) = 1.0_rt / vstate.ewt(NumSpec+1);
 #ifdef PHOTOCHEMISTRY
-            for (int i = 1; i <= NumChemActiveRadGroups; i+=MicrophysicsNumRadVarsPerGroup) {
+            for (int i = 1; i <= NumChemBands; i+=MicrophysicsNumRadVarsPerGroup) {
                 int rad_n_index = NumSpec + 2 + (i-1) * MicrophysicsNumRadVarsPerGroup;
                 vstate.ewt(rad_n_index) = vstate.rtol_rad_num * std::abs(vstate.yh(rad_n_index,1)) + vstate.atol_rad_num;
                 for (int j = 2; j <= MicrophysicsNumRadVarsPerGroup; ++j){