diff --git a/python/tests/reference/Result/6grains6x7x8_single_phase_tensionY.xdmf b/python/tests/reference/Result/6grains6x7x8_single_phase_tensionY.xdmf
index cc92cd571..80b96fecb 100644
--- a/python/tests/reference/Result/6grains6x7x8_single_phase_tensionY.xdmf
+++ b/python/tests/reference/Result/6grains6x7x8_single_phase_tensionY.xdmf
@@ -26,7 +26,7 @@
 				<Attribute Name="phase/mechanical/L_p / 1/s" Center="Cell" AttributeType="Tensor">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 9">6grains6x7x8_single_phase_tensionY.hdf5:increment_0/phase/pheno_fcc/mechanical/L_p</DataItem>
 				</Attribute>
-				<Attribute Name="phase/mechanical/O / q_0 (q_1 q_2 q_3)" Center="Cell" AttributeType="Matrix">
+				<Attribute Name="phase/mechanical/O / 1" Center="Cell" AttributeType="Matrix">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 4">6grains6x7x8_single_phase_tensionY.hdf5:increment_0/phase/pheno_fcc/mechanical/O</DataItem>
 				</Attribute>
 				<Attribute Name="phase/mechanical/P / Pa" Center="Cell" AttributeType="Tensor">
@@ -177,7 +177,7 @@
 				<Attribute Name="phase/mechanical/L_p / 1/s" Center="Cell" AttributeType="Tensor">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 9">6grains6x7x8_single_phase_tensionY.hdf5:increment_4/phase/pheno_fcc/mechanical/L_p</DataItem>
 				</Attribute>
-				<Attribute Name="phase/mechanical/O / q_0 (q_1 q_2 q_3)" Center="Cell" AttributeType="Matrix">
+				<Attribute Name="phase/mechanical/O / 1" Center="Cell" AttributeType="Matrix">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 4">6grains6x7x8_single_phase_tensionY.hdf5:increment_4/phase/pheno_fcc/mechanical/O</DataItem>
 				</Attribute>
 				<Attribute Name="phase/mechanical/P / Pa" Center="Cell" AttributeType="Tensor">
@@ -328,7 +328,7 @@
 				<Attribute Name="phase/mechanical/L_p / 1/s" Center="Cell" AttributeType="Tensor">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 9">6grains6x7x8_single_phase_tensionY.hdf5:increment_8/phase/pheno_fcc/mechanical/L_p</DataItem>
 				</Attribute>
-				<Attribute Name="phase/mechanical/O / q_0 (q_1 q_2 q_3)" Center="Cell" AttributeType="Matrix">
+				<Attribute Name="phase/mechanical/O / 1" Center="Cell" AttributeType="Matrix">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 4">6grains6x7x8_single_phase_tensionY.hdf5:increment_8/phase/pheno_fcc/mechanical/O</DataItem>
 				</Attribute>
 				<Attribute Name="phase/mechanical/P / Pa" Center="Cell" AttributeType="Tensor">
@@ -479,7 +479,7 @@
 				<Attribute Name="phase/mechanical/L_p / 1/s" Center="Cell" AttributeType="Tensor">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 9">6grains6x7x8_single_phase_tensionY.hdf5:increment_12/phase/pheno_fcc/mechanical/L_p</DataItem>
 				</Attribute>
-				<Attribute Name="phase/mechanical/O / q_0 (q_1 q_2 q_3)" Center="Cell" AttributeType="Matrix">
+				<Attribute Name="phase/mechanical/O / 1" Center="Cell" AttributeType="Matrix">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 4">6grains6x7x8_single_phase_tensionY.hdf5:increment_12/phase/pheno_fcc/mechanical/O</DataItem>
 				</Attribute>
 				<Attribute Name="phase/mechanical/P / Pa" Center="Cell" AttributeType="Tensor">
@@ -630,7 +630,7 @@
 				<Attribute Name="phase/mechanical/L_p / 1/s" Center="Cell" AttributeType="Tensor">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 9">6grains6x7x8_single_phase_tensionY.hdf5:increment_16/phase/pheno_fcc/mechanical/L_p</DataItem>
 				</Attribute>
-				<Attribute Name="phase/mechanical/O / q_0 (q_1 q_2 q_3)" Center="Cell" AttributeType="Matrix">
+				<Attribute Name="phase/mechanical/O / 1" Center="Cell" AttributeType="Matrix">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 4">6grains6x7x8_single_phase_tensionY.hdf5:increment_16/phase/pheno_fcc/mechanical/O</DataItem>
 				</Attribute>
 				<Attribute Name="phase/mechanical/P / Pa" Center="Cell" AttributeType="Tensor">
@@ -781,7 +781,7 @@
 				<Attribute Name="phase/mechanical/L_p / 1/s" Center="Cell" AttributeType="Tensor">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 9">6grains6x7x8_single_phase_tensionY.hdf5:increment_20/phase/pheno_fcc/mechanical/L_p</DataItem>
 				</Attribute>
-				<Attribute Name="phase/mechanical/O / q_0 (q_1 q_2 q_3)" Center="Cell" AttributeType="Matrix">
+				<Attribute Name="phase/mechanical/O / 1" Center="Cell" AttributeType="Matrix">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 4">6grains6x7x8_single_phase_tensionY.hdf5:increment_20/phase/pheno_fcc/mechanical/O</DataItem>
 				</Attribute>
 				<Attribute Name="phase/mechanical/P / Pa" Center="Cell" AttributeType="Tensor">
@@ -932,7 +932,7 @@
 				<Attribute Name="phase/mechanical/L_p / 1/s" Center="Cell" AttributeType="Tensor">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 9">6grains6x7x8_single_phase_tensionY.hdf5:increment_24/phase/pheno_fcc/mechanical/L_p</DataItem>
 				</Attribute>
-				<Attribute Name="phase/mechanical/O / q_0 (q_1 q_2 q_3)" Center="Cell" AttributeType="Matrix">
+				<Attribute Name="phase/mechanical/O / 1" Center="Cell" AttributeType="Matrix">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 4">6grains6x7x8_single_phase_tensionY.hdf5:increment_24/phase/pheno_fcc/mechanical/O</DataItem>
 				</Attribute>
 				<Attribute Name="phase/mechanical/P / Pa" Center="Cell" AttributeType="Tensor">
@@ -1083,7 +1083,7 @@
 				<Attribute Name="phase/mechanical/L_p / 1/s" Center="Cell" AttributeType="Tensor">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 9">6grains6x7x8_single_phase_tensionY.hdf5:increment_28/phase/pheno_fcc/mechanical/L_p</DataItem>
 				</Attribute>
-				<Attribute Name="phase/mechanical/O / q_0 (q_1 q_2 q_3)" Center="Cell" AttributeType="Matrix">
+				<Attribute Name="phase/mechanical/O / 1" Center="Cell" AttributeType="Matrix">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 4">6grains6x7x8_single_phase_tensionY.hdf5:increment_28/phase/pheno_fcc/mechanical/O</DataItem>
 				</Attribute>
 				<Attribute Name="phase/mechanical/P / Pa" Center="Cell" AttributeType="Tensor">
@@ -1234,7 +1234,7 @@
 				<Attribute Name="phase/mechanical/L_p / 1/s" Center="Cell" AttributeType="Tensor">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 9">6grains6x7x8_single_phase_tensionY.hdf5:increment_32/phase/pheno_fcc/mechanical/L_p</DataItem>
 				</Attribute>
-				<Attribute Name="phase/mechanical/O / q_0 (q_1 q_2 q_3)" Center="Cell" AttributeType="Matrix">
+				<Attribute Name="phase/mechanical/O / 1" Center="Cell" AttributeType="Matrix">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 4">6grains6x7x8_single_phase_tensionY.hdf5:increment_32/phase/pheno_fcc/mechanical/O</DataItem>
 				</Attribute>
 				<Attribute Name="phase/mechanical/P / Pa" Center="Cell" AttributeType="Tensor">
@@ -1385,7 +1385,7 @@
 				<Attribute Name="phase/mechanical/L_p / 1/s" Center="Cell" AttributeType="Tensor">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 9">6grains6x7x8_single_phase_tensionY.hdf5:increment_36/phase/pheno_fcc/mechanical/L_p</DataItem>
 				</Attribute>
-				<Attribute Name="phase/mechanical/O / q_0 (q_1 q_2 q_3)" Center="Cell" AttributeType="Matrix">
+				<Attribute Name="phase/mechanical/O / 1" Center="Cell" AttributeType="Matrix">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 4">6grains6x7x8_single_phase_tensionY.hdf5:increment_36/phase/pheno_fcc/mechanical/O</DataItem>
 				</Attribute>
 				<Attribute Name="phase/mechanical/P / Pa" Center="Cell" AttributeType="Tensor">
@@ -1536,7 +1536,7 @@
 				<Attribute Name="phase/mechanical/L_p / 1/s" Center="Cell" AttributeType="Tensor">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 9">6grains6x7x8_single_phase_tensionY.hdf5:increment_40/phase/pheno_fcc/mechanical/L_p</DataItem>
 				</Attribute>
-				<Attribute Name="phase/mechanical/O / q_0 (q_1 q_2 q_3)" Center="Cell" AttributeType="Matrix">
+				<Attribute Name="phase/mechanical/O / 1" Center="Cell" AttributeType="Matrix">
 					<DataItem Format="HDF" NumberType="Float" Precision="8" Dimensions="8 7 6 4">6grains6x7x8_single_phase_tensionY.hdf5:increment_40/phase/pheno_fcc/mechanical/O</DataItem>
 				</Attribute>
 				<Attribute Name="phase/mechanical/P / Pa" Center="Cell" AttributeType="Tensor">