diff --git a/t3/simulate/rmg_constantTP.py b/t3/simulate/rmg_constantTP.py
index fdb0b074..69151e28 100755
--- a/t3/simulate/rmg_constantTP.py
+++ b/t3/simulate/rmg_constantTP.py
@@ -359,7 +359,7 @@ def get_species_concentration_lists_from_ranged_params(self) -> List[List[dict]]
             spc_1_vals = get_values_within_range(value_range=self.rmg['species'][spc_indices_w_ranges[1]]['concentration'],
                                                  num=self.t3['options']['num_sa_per_concentration_range'])
             for val_0, val_1 in zip(spc_0_vals, spc_1_vals[::-1]):
-                new_species_list = species_list
+                new_species_list = species_list_wo_ranges  # Not final - needs modification
                 new_species_list.append({'label': self.rmg['species'][spc_indices_w_ranges[0]]['label'],
                                          'concentration': val_0})
                 new_species_list.append({'label': self.rmg['species'][spc_indices_w_ranges[1]]['label'],
@@ -378,7 +378,7 @@ def get_species_concentration_lists_from_ranged_params(self) -> List[List[dict]]
         # 4. Combinations (products)
         else:
             for vals in itertools.product(*species_vals):
-                new_species_list = species_list
+                new_species_list = species_list_wo_ranges  # Not final - needs modification
                 for i, val in enumerate(vals):
                     new_species_list.append({'label': self.rmg['species'][spc_indices_w_ranges[i]]['label'],
                                              'concentration': val})