diff --git a/thermosteam/_graphics.py b/thermosteam/_graphics.py
index 67ba76fd..89a202c9 100644
--- a/thermosteam/_graphics.py
+++ b/thermosteam/_graphics.py
@@ -109,7 +109,7 @@ def get_node_tailored_to_unit(self, unit): # pragma: no coverage
             else:
                 N = None
             label = '\n'.join([tmo.utils.format_title(i) for i in ID.split('.')])
-            label = f"{owner.ID}\n{label}\nAuxiliary"
+            label = f"{owner.ID}\n{label}"
             if N is not None and N > 1: label = f"{label}\n1 of {N}"
             node['label'] = label
         else:
diff --git a/thermosteam/reaction/_reaction.py b/thermosteam/reaction/_reaction.py
index 5ee608a4..94bf3279 100644
--- a/thermosteam/reaction/_reaction.py
+++ b/thermosteam/reaction/_reaction.py
@@ -490,7 +490,10 @@ def __call__(self, material, T=None, P=None, phase=None, time=None):
         values, config, original = as_material_array(
             material, self._basis, self._phases, self.chemicals
         )
-        self._reaction(values)
+        if values.ndim == 2:
+            for i in values: self._reaction(i)
+        else:
+            self._reaction(values)
         if tmo.reaction.CHECK_FEASIBILITY:
             has_negatives = values.has_negatives()
             if has_negatives: