diff --git a/PRIVATE b/PRIVATE
index 21105a488..bc6de828c 160000
--- a/PRIVATE
+++ b/PRIVATE
@@ -1 +1 @@
-Subproject commit 21105a488f1acbe119a76d9113a185245b45fd75
+Subproject commit bc6de828cc4ee9c941b37113ca49fcf51abd3512
diff --git a/src/phase_mechanical_elastic.f90 b/src/phase_mechanical_elastic.f90
index 9e5c2653c..45abf48c1 100644
--- a/src/phase_mechanical_elastic.f90
+++ b/src/phase_mechanical_elastic.f90
@@ -91,41 +91,43 @@ module function elastic_C66(ph,en) result(C66)
   integer, intent(in) :: &
     ph, &
     en
+
   real(pReal), dimension(6,6) :: C66
   real(pReal) :: T
 
+
   associate(prm => param(ph))
     C66 = 0.0_pReal
     T = thermal_T(ph,en)
 
-    C66(1,1) =   prm%C_11(1)* 1.0_pReal &
-               + prm%C_11(2)*(T - prm%T_ref)**1 &
-               + prm%C_11(3)*(T - prm%T_ref)**2
+    C66(1,1) = prm%C_11(1) &
+             + prm%C_11(2)*(T - prm%T_ref)**1 &
+             + prm%C_11(3)*(T - prm%T_ref)**2
 
-    C66(1,2) =   prm%C_12(1)* 1.0_pReal &
-               + prm%C_12(2)*(T - prm%T_ref)**1 &
-               + prm%C_12(3)*(T - prm%T_ref)**2
+    C66(1,2) = prm%C_12(1) &
+             + prm%C_12(2)*(T - prm%T_ref)**1 &
+             + prm%C_12(3)*(T - prm%T_ref)**2
 
-    C66(4,4) =   prm%C_44(1)*1.0_pReal &
-               + prm%C_44(2)*(T - prm%T_ref)**1 &
-               + prm%C_44(3)*(T - prm%T_ref)**2
+    C66(4,4) = prm%C_44(1) &
+             + prm%C_44(2)*(T - prm%T_ref)**1 &
+             + prm%C_44(3)*(T - prm%T_ref)**2
 
 
     if (any(phase_lattice(ph) == ['hP','tI'])) then
-      C66(1,3) =   prm%C_13(1)*1.0_pReal &
-                 + prm%C_13(2)*(T - prm%T_ref)**1 &
-                 + prm%C_13(3)*(T - prm%T_ref)**2
+      C66(1,3) = prm%C_13(1) &
+               + prm%C_13(2)*(T - prm%T_ref)**1 &
+               + prm%C_13(3)*(T - prm%T_ref)**2
 
-      C66(3,3) =   prm%C_33(1)*1.0_pReal &
-                 + prm%C_33(2)*(T - prm%T_ref)**1 &
-                 + prm%C_33(3)*(T - prm%T_ref)**2
+      C66(3,3) = prm%C_33(1) &
+               + prm%C_33(2)*(T - prm%T_ref)**1 &
+               + prm%C_33(3)*(T - prm%T_ref)**2
 
     end if
 
     if (phase_lattice(ph) == 'tI') then
-      C66(6,6) =   prm%C_66(1)*1.0_pReal &
-                 + prm%C_66(2)*(T - prm%T_ref)**1 &
-                 + prm%C_66(3)*(T - prm%T_ref)**2
+      C66(6,6) = prm%C_66(1) &
+               + prm%C_66(2)*(T - prm%T_ref)**1 &
+               + prm%C_66(3)*(T - prm%T_ref)**2
     end if
 
     C66 = lattice_symmetrize_C66(C66,phase_lattice(ph))